CN113393077B - Method for configuring an electric-gas multi-energy storage system taking into account the uncertainty of the energy used by the user - Google Patents
Method for configuring an electric-gas multi-energy storage system taking into account the uncertainty of the energy used by the user Download PDFInfo
- Publication number
- CN113393077B CN113393077B CN202110451877.2A CN202110451877A CN113393077B CN 113393077 B CN113393077 B CN 113393077B CN 202110451877 A CN202110451877 A CN 202110451877A CN 113393077 B CN113393077 B CN 113393077B
- Authority
- CN
- China
- Prior art keywords
- gas
- energy storage
- energy
- power
- storage system
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004146 energy storage Methods 0.000 title claims abstract description 117
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000003860 storage Methods 0.000 claims abstract description 82
- 230000005540 biological transmission Effects 0.000 claims abstract description 16
- 238000005457 optimization Methods 0.000 claims abstract description 5
- 239000007789 gas Substances 0.000 claims description 216
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 34
- 238000006243 chemical reaction Methods 0.000 claims description 22
- 230000005611 electricity Effects 0.000 claims description 17
- 239000003345 natural gas Substances 0.000 claims description 17
- 238000007599 discharging Methods 0.000 claims description 9
- 238000002347 injection Methods 0.000 claims description 9
- 239000007924 injection Substances 0.000 claims description 9
- 238000009826 distribution Methods 0.000 claims description 7
- 150000001875 compounds Chemical class 0.000 claims description 6
- 238000009827 uniform distribution Methods 0.000 claims description 5
- 238000002485 combustion reaction Methods 0.000 claims description 3
- 238000012423 maintenance Methods 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 3
- 239000011248 coating agent Substances 0.000 claims description 2
- 238000000576 coating method Methods 0.000 claims description 2
- 238000002360 preparation method Methods 0.000 claims description 2
- 230000000295 complement effect Effects 0.000 abstract description 3
- 238000012795 verification Methods 0.000 abstract 1
- 238000011160 research Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 238000006467 substitution reaction Methods 0.000 description 3
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 239000002131 composite material Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000011049 filling Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000005338 heat storage Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000011068 loading method Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 210000000352 storage cell Anatomy 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06313—Resource planning in a project environment
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Water Supply & Treatment (AREA)
- Computer Hardware Design (AREA)
- General Health & Medical Sciences (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Public Health (AREA)
- Biodiversity & Conservation Biology (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Primary Health Care (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses an electric-gas multi-energy storage system configuration method considering user energy uncertainty. Firstly, establishing a regional comprehensive energy system model containing an electricity-gas multi-energy storage system; then, analyzing the influence of different load characteristics and uncertainty of user energy types on the configuration of the energy storage system; on the basis, with the economy as a target, an electricity-gas multi-energy storage system optimization configuration model containing an energy storage battery, P2G equipment (power to gas), gas storage equipment and a gas transmission pipeline is established; finally, calling CPLEX in Matlab to solve the verification method is provided, and the complementary characteristics of the multi-energy storage system are proved to be utilized, so that the cost can be further reduced.
Description
Technical Field
The invention belongs to the research field of a multi-energy storage system, and particularly relates to an electric-gas multi-energy storage system configuration method considering user energy uncertainty.
Background
With the exhaustion of fossil energy and the increasing severity of environmental pollution, the energy internet has been widely recognized as an effective way to solve the problem. The energy internet is a novel technical information energy fusion open system based on the internet concept, can greatly improve the energy use efficiency, and promotes the large-scale development of renewable energy.
The regional comprehensive energy system generally comprises energy forms such as cold, heat, electricity and gas as an important component of an energy internet, and all energy supply equipment sources in a region are integrated and scheduled in a unified manner by using an internet of things technology and an information technology so as to achieve the effects of optimizing energy supply on regional cold, heat and electricity loads and improving energy utilization efficiency.
Energy storage technology has received wide attention as a support technology for integrated energy systems. The comprehensive energy system introduces a multi-energy storage system based on the traditional energy storage battery, and comprises various forms of electricity storage, heat storage, cold storage, gas storage, composite energy storage and the like. The energy storage system with a single energy form is installed in the comprehensive energy system to have a certain peak clipping and valley filling effect, but when the multi-energy type load needs to be supplemented by energy storage output, certain energy loss can be caused by certain interconversion among multiple energy sources, for example, the heat load cannot be supplied by the energy storage system in time. And the electricity storage system is relatively high in manufacturing cost, and the loss of large-scale energy storage is large. The heat and gas storage system can store energy in a large scale, but the inertia of the system is larger than that of the electricity storage system, and the energy transfer is slow. Thus, the different types of energy storage systems complement each other in performance. In addition, the multi-energy storage system plays a key role in electric energy replacement, and the multi-energy storage system with certain configuration can ensure complete electric energy replacement. Therefore, in order to enhance the flexibility of the integrated energy network, a multi-energy storage system is introduced to achieve the purpose of complementary advantages.
The research on the multi-energy storage system is concerned by students recently, and if the students establish an optimized operation model of a regional comprehensive energy system taking the minimum comprehensive cost as a target and considering the multi-energy storage; scholars propose a novel real-time energy management method to control an island direct-current microgrid driven by a photovoltaic array, and the microgrid is simultaneously provided with 2 different types of energy storage systems of electricity and hydrogen; researchers study the configuration of a hydrogen storage system in the process of converting electricity into natural gas; researchers have proposed an energy storage optimization configuration method for a regional comprehensive energy system considering electric/thermal flexible load;
however, the research does not research the optimal configuration of the electric-gas multi-energy system, does not research the optimal configuration model of various devices of the electric-gas multi-energy storage system such as an energy storage battery, a P2G device, a gas storage device and a gas transmission pipeline, does not consider the uncertainty of the user energy types in the optimal configuration of the multi-energy storage system, and has little discussion on the feasibility of electric energy substitution after the multi-energy storage system is added.
Disclosure of Invention
The present invention is directed to overcoming the disadvantages of the prior art and providing a method for configuring an electric-gas multi-energy storage system that takes into account user energy uncertainty.
The purpose of the invention is realized by the following technical scheme: an electric-gas multi-energy storage system configuration method considering user energy uncertainty, comprising the steps of:
(1) Establishing a regional comprehensive energy system model containing an electricity-gas multi-energy storage system;
(2) Analyzing the influence of different load characteristics and uncertainty of user energy types on the configuration of the energy storage system and adjusting the load:
and (2.1) determining different supply sources of the loads.
And (2.2) determining the proportion of the selectable load.
(3) With the economy as a target, an electricity-gas multi-energy storage system optimization configuration model containing an energy storage battery, a P2G device, a gas storage device and a gas transmission pipeline is established:
and (3.1) planning the rated power and the capacity of the multi-energy storage system in the established integrated energy system to minimize the total cost in the whole life cycle. Meanwhile, the time value of the capital is considered, and all the cost is reduced to the equal-year value.
And (3.2) adding gas turbine output constraint, gas boiler constraint, electric boiler constraint, energy storage battery constraint and gas storage constraint in the established model.
(4) And (4) optimizing and solving the model established in the step (3).
Further, the step (1) specifically comprises:
(1.1) establishing a physical model of the regional comprehensive energy system based on the electricity-gas multi-energy storage system. The multifunctional storage system comprises an energy storage battery, a P2G device, a gas storage device and input and output pipelines of the gas storage device. The loads are divided into electrical and thermal loads, and other devices are gas turbines, wind generators, gas boilers and electric boilers.
And (1.2) establishing a mathematical model of the regional comprehensive energy system based on the electricity-gas multi-energy storage system. The method specifically comprises the following steps:
(1.2.1) external energy injection: injection power of distribution network
And net injection power
(1.2.2) electric load
Power generated by fan WT Energy storage battery output P Bat Part for power purchasing of distribution network
And the output P of the gas turbine GT Supplying:
wherein, mu and (1-mu) respectively represent the shunt coefficients of the power used for the P2G equipment and the direct power supply; λ and (1- λ) represent the split coefficients of natural gas input to the gas turbine and the gas boiler. Eta P2G 、η GF 、η GT The conversion efficiencies of the P2G plant, gas boiler, gas turbine and electric boiler, respectively.
(1.2.3) Heat load
Supplied by a gas boiler:
wherein, P GS Outputting power to the gas storage device.
(1.2.4) Power P of P2G device P2G Indicate electric power conversion natural gas power, the P2G equipment directly gets into gas storage system through gas transmission pipeline as a part of gas storage device, from gas storage system through gas transmission pipeline supply load again when needing:
(1.2.5) the gas source of the gas turbine is a part of the output of the gas storage device and the input of the gas network:
wherein eta is GS The conversion efficiency of the gas storage equipment.
Further, in step (2.2), the optional load is partially powered by the gas boiler and partially powered by the electric boiler and is treated as an electric load. The load type of each time period is distributed with a proportional load uniform distribution.
Further, the step (3.1) is specifically:
(3.1) the objective function of the plan can be expressed as:
min C＝k pa IC+OC+FC
wherein C is the total cost converted to annuity; IC is the initial investment cost; OC is the annual operating cost; FC is a fixed cost and a fixed value, and mainly refers to annual maintenance cost; k is a radical of pa When the annual interest rate is r and the energy storage life cycle is y years, the expression is as follows:
further, in step (3.1):
(3.1.1) the investment cost mainly comprises two parts of an energy storage battery and a gas storage device:
in the formula:
respectively representing the maximum power and the maximum capacity of the storage battery;
represents the maximum power of the P2G device;
indicating the maximum capacity of the gas storage device;
respectively representing the maximum power of the input and output gas storage devices; k is a radical of 1 、k 2 Respectively representing the unit power investment coefficient and the unit capacity investment coefficient of the storage battery; k is a radical of 3 Expressing the unit power investment coefficient of the P2G equipment; k is a radical of 4 Representing the unit capacity investment coefficient of the gas storage device; k is a radical of 5 The unit power investment coefficient of the gas transmission pipeline is shown.
(3.1.2) operating costs refer to electricity purchase from the main grid and gas purchase from the gas grid, and are expressed as follows:
in the formula, S and T respectively represent typical days and scheduling hours of a day; d s Represents the total number of days for a typical day s;
the electricity purchase price for the t-th scheduled hour;
representing the purchased electric power of the tth scheduling hour of the s typical day; lambda g Representing a uniform price of the natural gas network; r represents a calorific value coefficient of combustion of natural gas;
and scheduling the purchased gas power of the hour for the tth typical day.
Further, since the natural gas produced by the P2G equipment is directly transferred into the gas storage device, the method is approximately considered as
Further, in step (3.2):
(3.2.1) device constraints
(3.2.1.1) gas turbine output constraint
0≤P GT ≤P GT，max
In the formula: p is GT，max Is the maximum output of the gas turbine; eta GT The gas-electricity conversion efficiency of the gas turbine is obtained;
for input of gas turbine gas power.
(3.2.1.2) gas boiler constraint
0≤P GF ≤P GF，max
In the formula: p GF，max The maximum thermal output of the gas boiler; eta GF The gas heat conversion efficiency of the gas boiler is obtained;
for inputting the gas power of the gas boiler.
(3.2.1.3) electric boiler restraint
0≤P EB ≤P EB，max
In the formula: p EB，max The maximum thermal output of the electric boiler; eta EB The electric heat conversion efficiency of the electric boiler is obtained;
is the input of electric boiler power.
(3.2.2) energy storage System restraint
(3.2.2.1) energy storage battery restraint:
in the formula (I), the compound is shown in the specification,
representing the state of charge of the energy storage battery at the tth scheduling hour of the sth typical day; eta e,C 、η e,F Respectively representing the charge and discharge efficiency of the energy storage battery; delta e Representing the self-discharge efficiency of the energy storage battery;
the charge state upper and lower limits of the energy storage battery are set;
the charging and discharging power in the period of the t scheduled hour of the s typical day cannot be charged and discharged at the same time in the same period, and the product of the charging and discharging power and the discharging power is 0; Δ t represents a time interval.
(3.2.2.2) gas storage system constraint:
in the formula (I), the compound is shown in the specification,
representing the charge energy state of the gas storage system at the tth scheduling hour of the sth typical day; eta g,C 、η g,F Respectively representing the charge and discharge efficiency of the energy storage battery; delta g The gas consumption efficiency of the gas storage system is shown;
the upper limit and the lower limit of the charge energy state of the gas storage system;
representing the charge energy state of the gas storage system at the end of the t-1 time period of the s-th typical day;
the gas storage and gas consumption power of the time interval of the t scheduled hour of the s typical day is represented; eta P2G The electric gas conversion efficiency of the P2G equipment is obtained;
electric power input for the P2G device.
Further, the SOC of the last period of the multi-energy storage system is equal to the SOC of the first period for each typical day:
further, in the step (4), the optimal configuration model established in the step (3) is solved by using CPLEX. Wherein, for the condition that the product of [0,1] continuous variable and [0, M ] continuous variable appears in the model, the product is linearized by adopting a relaxation method to obtain a linear programming model with the following form:
obj.min f(x)
s.t.A eq x＝b eq Ax≤b
wherein x is a control variable and a state variable, and f (x) is a planning total cost objective function. A. The eq For equality-constrained independent variable coefficients, b eq Is a constant in the equality constraint; a is an independent variable coefficient in an inequality constraint, and b is a constant in the inequality constraint.
Furthermore, the control variable and the state variable x comprise rated power and capacity of the energy storage device, output of the device in each scheduling period, system electricity and gas purchasing quantity, a shunt coefficient and the like.
The beneficial effects of the invention are:
(1) The invention plans and constructs the multi-energy storage system in the comprehensive energy system to influence the total cost of the whole system, and the method of the invention selects more reasonable energy storage system configuration to increase the system profit and reduce the total cost of the system;
(2) The planning construction of the multi-energy storage system is beneficial to electric energy substitution, and the multi-energy storage system with certain configuration can realize complete electric energy substitution.
Detailed Description
The invention relates to a configuration method of an electricity-gas multi-energy storage system considering user energy uncertainty, which researches the configuration of the electricity-gas multi-energy storage system in a regional comprehensive energy system and aims to further establish the relation between the electricity-gas systems through the electricity-gas multi-energy storage system. Firstly, establishing a comprehensive energy network model containing an electricity-gas multi-energy storage system; then, analyzing the influence of different load characteristics and uncertainty of user energy types on the configuration of the energy storage system; on the basis, with the economy as a target, an electricity-gas multi-energy storage system optimization configuration model containing energy storage batteries, P2G equipment, gas storage equipment, gas transmission pipelines and other equipment is established; finally, it is solved by CPLEX. The method comprises the following specific steps:
(1) Establishing a regional comprehensive energy system model containing an electric-gas multi-energy storage system; the method specifically comprises the following steps:
(1.1) establishing a physical model of the regional comprehensive energy system based on the electricity-gas multi-energy storage system. The multifunctional storage system comprises an energy storage battery (Bat), a P2G device (power to gas, electricity to gas), a gas storage device (GS) and input and output pipelines of the gas storage device. The loads are divided into electrical and thermal loads, and other devices include Gas Turbines (GT), wind Turbines (WT), gas boilers (GF), and Electric Boilers (EB). The energy conversion device comprises P2G, GF, GT and EB.
(1.2) establishing a mathematical model of a regional comprehensive energy system based on the electricity-gas containing multi-energy storage system; the method comprises the following specific steps:
(1.2.1) external energy injection: injection power of distribution network
And net injection power
(1.2.2) electric load
Power generated by fan P WT Energy storage battery output P Bat Part for power purchase of distribution network
And the output P of the gas turbine GT Supplying:
wherein, mu and (1-mu) respectively represent the shunt coefficients of the power used for the P2G equipment and the direct power supply, and the sum of the shunt coefficients is 1; λ and (1- λ) represent the split coefficients of natural gas input to the gas turbine and the gas boiler, the sum of which is 1; eta P2G 、η GF 、η GT The conversion efficiencies of the P2G plant, gas boiler, gas turbine and electric boiler, respectively.
(1.2.3) Heat load
Supplied by a gas boiler:
wherein, P GS And outputting power to the gas storage device. In particular, the electric boiler can also be supplied with a thermal load, which the model of the invention treats as an electrical load.
(1.2.4) Power P of P2G device P2G Indicate the natural gas power of electric power conversion, the output of P2G equipment directly gets into gas storage system through gas transmission pipeline, and from gas storage system through gas transmission pipeline supply load again when needing:
(1.2.5) gas turbine gas source is gas storage device output P GS Part of the inputs to the air grid:
wherein eta GS The conversion efficiency of the gas storage equipment.
(2) Analyzing the influence of different load characteristics and uncertainty of user energy types on the configuration of the energy storage system and adjusting the load; the method specifically comprises the following steps:
(2.1) determining the supply sources of different loads, for example, the supply of heat load can be selected from an electric boiler or a gas boiler.
And (2.2) determining the proportion of the selectable load.
The invention mainly aims at the uncertainty of selecting different heating equipment according to different preferences of users, and a certain proportion of heat loads are used as selectable supply loads. In the invention, part of the selective load is powered by a gas boiler, and part of the selective load is powered by an electric boiler and is treated as an electric load. The load type of each time period is distributed with a proportional load uniform distribution as follows:
in the formula, P ran.e For selectable electrical load portions of the load, P [0,1] Is in [0,1]]Probability of medium uniform distribution, η h To select the ratio of load to thermal load, P ran.h The portion of the load that remains as a thermal load is selected.
(3) Establishing an electricity-gas multi-energy storage system optimal configuration model containing an energy storage battery, a P2G device, a gas storage device and a gas transmission pipeline by taking economy as a target; the method comprises the following steps:
(3.1) establishing an objective function
In the integrated energy system that has been built, the power rating and capacity of the multi-energy storage system are planned to minimize its total cost over the life cycle. Meanwhile, the time value of the capital is considered, and all the cost is reduced to the equal-year value. The planned objective function can be expressed as:
min C＝k pa IC+OC+FC
wherein C is the total cost converted to annuity; IC is the initial investment cost; OC is the annual operating cost; FC is a fixed cost and a fixed value, and mainly refers to annual maintenance cost; k is a radical of pa When the annual interest rate is r and the energy storage life cycle is y years, the expression is as follows:
(3.1.1) investment cost
As other equipment is built, the investment plan of multi-energy storage in the comprehensive energy network mainly comprises two parts of an energy storage battery and a gas storage device:
in the formula:
respectively representing the maximum power and the maximum capacity of the storage battery;
represents the maximum power of the P2G device;
indicating the maximum capacity of the gas storage device;
respectively representing the maximum power of the input and output gas storage devices; wherein, because the natural gas produced by the P2G equipment is directly transmitted into the gas storage device, the natural gas is approximately considered
k 1 、k 2 Respectively representing the unit power investment coefficient and the unit capacity investment coefficient of the storage battery; k is a radical of formula 3 Expressing the unit power investment coefficient of the P2G equipment; k is a radical of 4 Representing the unit capacity investment coefficient of the gas storage device; k is a radical of 5 The unit power investment coefficient of the gas transmission pipeline is shown.
(3.1.2) running cost
In the comprehensive energy network, the operation cost refers to the electricity purchasing cost from the main network and the gas purchasing cost from the gas network, and the expression of the operation cost is as follows:
in the formula, S and T respectively represent typical days and scheduling hours of one day; d is a radical of s Represents the total number of days for a typical day s;
the electricity purchase price for the t-th scheduled hour;
representing the purchased electric power of the t scheduled hour on the s typical day; lambda [ alpha ] g The unit of the unit is Yuan/m and represents the unified price of the natural gas network 3 (ii) a r represents the calorific value coefficient of combustion of natural gas in kW/m 3 ；
And scheduling the gas purchasing power of the hour for the tth typical day.
(3.2) constraint conditions: the built model is additionally provided with gas turbine output constraint, gas boiler constraint, electric boiler constraint, energy storage battery constraint, gas storage constraint and the like; the method specifically comprises the following steps:
(3.2.1) device constraints:
(3.2.1.1) gas turbine output constraint
0≤P GT ≤P GT，max
In the formula: p GT，max Is the maximum output of the gas turbine; eta GT The gas-electricity conversion efficiency of the gas turbine is obtained;
for inputting gas turbine gas power.
(3.2.1.2) gas boiler constraint
0≤P GF ≤P GF，max
In the formula: p GF，max The maximum thermal output of the gas boiler; eta GF The gas heat conversion efficiency of the gas boiler is obtained;
for inputting the gas power of the gas boiler.
(3.2.1.3) electric boiler restraint
0≤P EB ≤P EB，max
In the formula: p is EB，max The maximum thermal output of the electric boiler; eta EB The electric heat conversion efficiency of the electric boiler is obtained;
for inputting electric power to the electric boiler.
(3.2.2) energy storage System constraints
(3.2.2.1) energy storage cell restraint:
in the formula (I), the compound is shown in the specification,
representing the state of charge of the energy storage battery at the tth scheduling hour of the sth typical day; eta e,C 、η e,F Respectively representing the charge and discharge efficiency of the energy storage battery; delta. For the preparation of a coating e Representing the self-discharge efficiency of the energy storage battery;
the charge state upper and lower limits of the energy storage battery are set;
the charging and discharging power in the period of the t scheduled hour of the s typical day cannot be charged and discharged at the same time in the same period, and the product of the charging and discharging power and the discharging power is 0; at represents the time interval, and the time-of-use electricity price is adopted in the invention, namely, the at is taken as 1 hour.
(3.2.2.2) gas storage system constraint:
in the formula (I), the compound is shown in the specification,
representing the charge energy state of the gas storage system at the tth scheduling hour of the sth typical day; eta g,C 、η g,F Respectively representing the charge and discharge efficiency of the energy storage battery; delta g The gas consumption efficiency of the gas storage system is shown;
the upper and lower limits of the charge energy state of the gas storage system;
the energy loading state of the gas storage system at the end time of the t-1 time period of the s-th typical day is shown;
the gas storage and gas consumption power of the time interval representing the t scheduled hour of the s typical day is worth explaining thatThe gas and the gas consumption respectively pass through different gas transmission pipelines, so that the gas storage and the gas consumption can be carried out simultaneously, and the gas storage power is approximately equal to the output power of the P2G equipment; eta P2G The electric gas conversion efficiency of the P2G equipment is obtained;
electric power input for the P2G device; and in order not to generate errors, the SOC of the last period of the multi-energy storage system is specified to be equal to the SOC of the first period for each typical day:
(4) And calling CPLEX in Matlab to solve.
For the occurrence of [0,1 in the model]Continuous variable x 1 And [0,M]Continuous variable y 1 In the case of the product, it is linearized by a relaxation method (second-order cone relaxation method); wherein, [0,1]Continuous variable x 1 Refers to the shunt coefficients μ and λ, [0,M]Continuous variable y 1 Refers to the power of a device multiplied by μ and λ in the partial equality and inequality constraints (equations (6) and (7)), which, since both are variables, results in model non-linearity; m is a positive real number. After a T variable is introduced, a CPLEX solver can be called in Matlab to solve the optimized configuration model established in the step (3) after the nonlinear part is subjected to linearization processing through linearization processing. The method specifically comprises the following steps:
introducing a new variable T, let
T＝x 1 y 1 x 1 ∈[0，1]y 1 ∈[0，M]
The linearizable process is:
T-0·x 1 ≥0
T-M·x 1 ＜0
T-y+0·(1-x 1 )≤0
T-y+M·(1-x 1 )≥0
in order to clearly describe the flow of the algorithm, the model is unified and compact to obtain a linear programming model in the following form:
obj.min f(x) (5)
s.t.A eq x＝b eq (6)
Ax≤b (7)
wherein, the formula (5) is a planning total cost objective function; x is a control variable and a state variable, and comprises rated power and capacity of energy storage equipment, equipment output force in each scheduling period, system electricity and gas purchasing quantity, a shunting coefficient and the like; the equations (6) and (7) are respectively linear equality constraint and inequality constraint, and comprise balance of electric and gas output, constraint of equipment output, conservation of energy storage equipment and the like; a. The eq For equality-constrained independent variable coefficients, b eq Is a constant in the equality constraint; a is an independent variable coefficient in an inequality constraint, and b is a constant in the inequality constraint.
Claims (9)
1. An electric-gas multi-energy storage system configuration method considering uncertainty of user energy, comprising the steps of:
(1) Establishing a regional comprehensive energy system model containing an electricity-gas multi-energy storage system, comprising the following steps of:
(1.1) establishing a physical model of a regional comprehensive energy system based on the electricity-gas-containing multi-energy storage system; the multi-energy storage system comprises an energy storage battery, a P2G device, a gas storage device and input and output pipelines of the gas storage device; the load is divided into an electric load and a heat load, and other equipment comprises a gas turbine, a wind driven generator, a gas boiler and an electric boiler;
(1.2) establishing a mathematical model of a regional comprehensive energy system based on the electricity-gas containing multi-energy storage system; the method specifically comprises the following steps:
(1.2.1) external energy injection: injection power of distribution network
And net injection power
(1.2.2) electric load
Power generated by fan P WT Energy storage battery output P Bat Part for power purchasing of distribution network
And the output P of the gas turbine GT Supplying:
wherein, mu and (1-mu) respectively represent the shunt coefficients of the power used for the P2G equipment and the direct power supply; λ and (1- λ) represent the split coefficients of natural gas input to the gas turbine and the gas boiler; eta P2G 、η GF 、η GT And η EB The conversion efficiencies of the P2G equipment, the gas boiler, the gas turbine and the electric boiler are respectively;
(1.2.3) Heat load
Supplied by a gas boiler:
wherein, P GS Outputting power to the gas storage device;
(1.2.4) Power P of P2G device P2G Indicate electric power conversion natural gas power, the P2G equipment directly gets into gas storage system through gas transmission pipeline as a part of gas storage device, from gas storage system through gas transmission pipeline supply load again when needing:
(1.2.5) the gas source of the gas turbine is a part of the output of the gas storage device and the input of the gas network:
wherein eta GS The conversion efficiency of the gas storage equipment.
(2) Analyzing the influence of different load characteristics and user energy type uncertainty on the configuration of the energy storage system and adjusting the load:
(2.1) determining different supply sources of loads;
(2.2) determining a duty ratio of the selectable load;
(3) With the economy as a target, an electricity-gas multi-energy storage system optimization configuration model containing an energy storage battery, a P2G device, a gas storage device and a gas transmission pipeline is established:
(3.1) planning the rated power and the capacity of the multi-energy storage system in the built integrated energy system to minimize the total cost in the whole life cycle; meanwhile, the time value of capital is considered, and all the cost is reduced to the equal-year value;
(3.2) adding output constraint, gas boiler constraint, electric boiler constraint, energy storage battery constraint and gas storage constraint of the gas turbine in the established model;
(4) And (4) optimizing and solving the model established in the step (3).
2. The method for configuring an electric-gas multi-energy storage system in consideration of uncertainty of user energy according to claim 1, wherein in the step (2.2), the alternative load is partially powered by the gas boiler and partially powered by the electric boiler to be treated as the electric load; the load type of each time period is distributed with a proportional load uniform distribution as follows:
in the formula, P ran.e For selectable electrical load portions of the load, P [0,1] Is in [0,1]]Probability of medium uniform distribution, η h To select the ratio of load to thermal load, P ran.h The portion of the selectable load that remains as the thermal load is selected.
3. The method for configuring an electric-gas multi-energy storage system taking account of uncertainty in user energy according to claim 1, wherein the step (3.1) is specifically:
(3.1) the planned objective function can be expressed as:
min C＝k pa IC+OC+FC；
wherein C is the total cost converted to annuity; IC is the initial investment cost; OC is the annual operating cost; FC is a fixed cost and a fixed value, and mainly refers to annual maintenance cost; k is a radical of pa When the annual interest rate is r and the energy storage life cycle is y years, the expression is as follows:
4. the method for configuring an electro-pneumatic multi energy storage system considering uncertainty of user energy according to claim 3, wherein in the step (3.1):
(3.1.1) the investment cost mainly comprises two parts of an energy storage battery and a gas storage device:
in the formula:
respectively representing the maximum power and the maximum capacity of the storage battery;
represents the maximum power of the P2G device;
indicating the maximum capacity of the gas storage device;
respectively representing the maximum power of the input and output gas storage devices; k is a radical of formula 1 、k 2 Respectively representing the unit power investment coefficient and the unit capacity investment coefficient of the storage battery; k is a radical of formula 3 Expressing the investment coefficient of unit power of the P2G equipment; k is a radical of 4 The unit capacity investment coefficient of the gas storage device is represented; k is a radical of 5 Expressing the unit power investment coefficient of the gas transmission pipeline;
(3.1.2) operating costs refer to electricity purchase from the main grid and gas purchase from the gas grid, and are expressed as follows:
in the formula, S and T respectively represent typical days and scheduling hours of a day; d s Represents the total number of days for a typical day s;
the electricity purchase price for the t-th scheduled hour;
representing the purchased electric power of the t scheduled hour on the s typical day; lambda g Representing a natural gas network uniform price; r represents a calorific value coefficient of combustion of natural gas;
and scheduling the gas purchasing power of the hour for the tth typical day.
5. The method for configuring an electric-gas multi-energy storage system taking account of uncertainty in user energy according to claim 4, wherein natural gas produced by the P2G plant is approximately considered as being directly transferred to the gas storage facility
6. The method for configuring an electric-gas multi-energy storage system considering uncertainty of user energy according to claim 3, wherein in the step (3.2):
(3.2.1) device constraints
(3.2.1.1) gas turbine output constraint
0≤P GT ≤P GT，max ；
In the formula: p GT，max Is the maximum output of the gas turbine; eta GT The gas-electricity conversion efficiency of the gas turbine is obtained;
inputting gas power of a gas turbine;
(3.2.1.2) gas boiler constraint
0≤P GF ≤P GF，max ；
In the formula: p GF，max The maximum thermal output of the gas boiler; eta GF The gas heat conversion efficiency of the gas boiler is obtained;
inputting the gas power of the gas boiler;
(3.2.1.3) electric boiler restraint
0≤P EB ≤P EB，max ；
In the formula: p EB，max The maximum thermal output of the electric boiler; eta EB The electric heat conversion efficiency of the electric boiler is obtained;
inputting electric power of an electric boiler;
(3.2.2) energy storage System restraint
(3.2.2.1) energy storage battery restraint:
in the formula (I), the compound is shown in the specification,
representing the state of charge of the energy storage battery at the tth scheduling hour of the sth typical day; eta e,C 、η e,F Respectively representing stored energy electricityThe charge-discharge efficiency of the cell; delta. For the preparation of a coating e Representing the self-discharge efficiency of the energy storage battery;
the charge state upper and lower limits of the energy storage battery are set;
the charging and discharging power in the period of the tth scheduling hour of the sth typical day cannot be charged and discharged simultaneously in the same period, and the product of the charging and discharging power and the discharging power is 0; Δ t represents a time interval;
(3.2.2.2) gas storage system constraint:
in the formula (I), the compound is shown in the specification,
representing the charge energy state of the gas storage system at the tth scheduling hour of the sth typical day; eta g,C 、η g,F Respectively representing the charge and discharge efficiency of the energy storage battery; delta g The gas consumption efficiency of the gas storage system is shown;
the upper and lower limits of the charge energy state of the gas storage system;
representing the charge energy state of the gas storage system at the end of the t-1 time period of the s-th typical day;
the gas storage and gas consumption power of the time interval representing the t scheduled hour of the s typical day; eta P2G The electric gas conversion efficiency of the P2G equipment is obtained;
electric power input for the P2G device.
7. The method of configuring an electric-to-gas multi-energy storage system taking into account user energy uncertainty as claimed in claim 5, wherein the SOC of the last period of the multi-energy storage system is equal to the SOC of the first period every typical day:
8. the method for configuring an electro-pneumatic multi-energy storage system in consideration of user energy uncertainty as claimed in claim 1, wherein in the step (4), the optimized configuration model established in the step (3) is solved by CPLEX; wherein, for the condition that the product of [0,1] continuous variable and [0, M ] continuous variable appears in the model, the product is linearized by adopting a relaxation method to obtain a linear programming model with the following form:
obj.min f(x)；
s.t.A eq x＝b eq ；Ax≤b；
wherein x is a control variable and a state variable, and f (x) is a planning total cost objective function; a. The eq For equality-constrained independent variable coefficients, b eq Is a constant in the equality constraint; a is the independent variable coefficient in the inequality constraint, and b is the constant in the inequality constraint.
9. The method according to claim 8, wherein the control variables and the state variables x comprise rated power of the energy storage device, capacity of the energy storage device, output of the device in each scheduling period, system power purchase amount, system gas purchase amount and flow distribution coefficient.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110451877.2A CN113393077B (en) | 2021-04-26 | 2021-04-26 | Method for configuring an electric-gas multi-energy storage system taking into account the uncertainty of the energy used by the user |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110451877.2A CN113393077B (en) | 2021-04-26 | 2021-04-26 | Method for configuring an electric-gas multi-energy storage system taking into account the uncertainty of the energy used by the user |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113393077A CN113393077A (en) | 2021-09-14 |
| CN113393077B true CN113393077B (en) | 2023-03-14 |
Family
ID=77617700
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110451877.2A Active CN113393077B (en) | 2021-04-26 | 2021-04-26 | Method for configuring an electric-gas multi-energy storage system taking into account the uncertainty of the energy used by the user |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113393077B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114169137B (en) * | 2021-11-10 | 2024-07-19 | 华中科技大学 | Function space optimization method and system for coordinating multi-energy storage of gas-electricity combined system |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109885009A (en) * | 2019-03-19 | 2019-06-14 | 广东电网有限责任公司电网规划研究中心 | Meter and electricity turn the garden energy source optimization configuration method of providing multiple forms of energy to complement each other of gas planning |
| CN111210054A (en) * | 2019-12-22 | 2020-05-29 | 上海电力大学 | Micro-energy network optimization scheduling method considering direct load control uncertainty |
| CN111445090A (en) * | 2020-04-21 | 2020-07-24 | 清华大学 | Double-layer planning method for off-grid type comprehensive energy system |
| CN111476509A (en) * | 2020-05-28 | 2020-07-31 | 南方电网科学研究院有限责任公司 | User side comprehensive energy system planning method and device based on IGDT model |
| CN111861783A (en) * | 2020-04-26 | 2020-10-30 | 国网江苏省电力有限公司经济技术研究院 | Comprehensive energy system multi-objective optimization configuration method considering load transfer |
| CN112036652A (en) * | 2020-09-06 | 2020-12-04 | 华北电力大学 | Photovoltaic-energy storage integrated energy system planning method based on opportunity constraint planning |
| CN112417651A (en) * | 2020-10-29 | 2021-02-26 | 国网浙江省电力有限公司温州供电公司 | Regret avoidance-based user-level comprehensive energy system optimization method |
| CN112464477A (en) * | 2020-11-27 | 2021-03-09 | 国网山东省电力公司青岛供电公司 | Multi-energy coupling comprehensive energy operation simulation method considering demand response |
| CN112580938A (en) * | 2020-12-03 | 2021-03-30 | 国家电网有限公司 | Multi-uncertainty-oriented optimization scheduling method and device for comprehensive energy system |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100217642A1 (en) * | 2009-02-26 | 2010-08-26 | Jason Crubtree | System and method for single-action energy resource scheduling and participation in energy-related securities |
| US20110246381A1 (en) * | 2010-03-30 | 2011-10-06 | Aide Audra Fitch | Systems and methods of modeling energy consumption of buildings |
| CN112529244B (en) * | 2020-10-23 | 2022-08-23 | 河海大学 | Comprehensive energy system collaborative optimization operation method considering electric load demand response |
-
2021
- 2021-04-26 CN CN202110451877.2A patent/CN113393077B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109885009A (en) * | 2019-03-19 | 2019-06-14 | 广东电网有限责任公司电网规划研究中心 | Meter and electricity turn the garden energy source optimization configuration method of providing multiple forms of energy to complement each other of gas planning |
| CN111210054A (en) * | 2019-12-22 | 2020-05-29 | 上海电力大学 | Micro-energy network optimization scheduling method considering direct load control uncertainty |
| CN111445090A (en) * | 2020-04-21 | 2020-07-24 | 清华大学 | Double-layer planning method for off-grid type comprehensive energy system |
| CN111861783A (en) * | 2020-04-26 | 2020-10-30 | 国网江苏省电力有限公司经济技术研究院 | Comprehensive energy system multi-objective optimization configuration method considering load transfer |
| CN111476509A (en) * | 2020-05-28 | 2020-07-31 | 南方电网科学研究院有限责任公司 | User side comprehensive energy system planning method and device based on IGDT model |
| CN112036652A (en) * | 2020-09-06 | 2020-12-04 | 华北电力大学 | Photovoltaic-energy storage integrated energy system planning method based on opportunity constraint planning |
| CN112417651A (en) * | 2020-10-29 | 2021-02-26 | 国网浙江省电力有限公司温州供电公司 | Regret avoidance-based user-level comprehensive energy system optimization method |
| CN112464477A (en) * | 2020-11-27 | 2021-03-09 | 国网山东省电力公司青岛供电公司 | Multi-energy coupling comprehensive energy operation simulation method considering demand response |
| CN112580938A (en) * | 2020-12-03 | 2021-03-30 | 国家电网有限公司 | Multi-uncertainty-oriented optimization scheduling method and device for comprehensive energy system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113393077A (en) | 2021-09-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108229025B (en) | Economic optimization scheduling method for cooling, heating and power combined supply type multi-microgrid active power distribution system | |
| CN106372742A (en) | Power-to-gas multi-source energy storage type microgrid day-ahead optimal economic dispatching method | |
| CN110689189A (en) | Combined cooling heating and power supply and demand balance optimization scheduling method considering energy supply side and demand side | |
| CN106849190A (en) | A kind of microgrid real-time scheduling method of providing multiple forms of energy to complement each other based on Rollout algorithms | |
| Wu et al. | Optimal generation scheduling of a microgrid | |
| CN113131513B (en) | Method for optimizing operation of electric, thermal and gas conversion system with consideration of carbon emission and storage medium | |
| CN105186583B (en) | Energy router and its energy dispatching method based on multi-agent modeling | |
| CN112990523A (en) | Regional comprehensive energy system layered optimization operation method based on multi-objective model predictive control | |
| CN112036652A (en) | Photovoltaic-energy storage integrated energy system planning method based on opportunity constraint planning | |
| CN117081143A (en) | Method for promoting coordination and optimization operation of park comprehensive energy system for distributed photovoltaic on-site digestion | |
| CN113393077B (en) | Method for configuring an electric-gas multi-energy storage system taking into account the uncertainty of the energy used by the user | |
| Yu et al. | Optimization of urban multi-energy flow systems considering seasonal peak shaving of natural gas | |
| Deckmyn et al. | Multi-objective optimization for environomic scheduling in microgrids | |
| Li et al. | Planning model of integrated energy system considering P2G and energy storage | |
| CN116663818A (en) | Low-carbon economic dispatching method for virtual power plant under ladder carbon transaction mechanism | |
| Xing et al. | Multi-energy simulation and optimal scheduling strategy based on digital twin | |
| Zhou et al. | Dynamic pricing of integrated energy service providers based on master-slave game | |
| CN115906456A (en) | Hydrogen-containing energy IES scheduling optimization model considering response uncertainty of demand side | |
| CN115483718A (en) | Electric-qi deficiency simulated power plant double-layer optimal scheduling method considering energy storage and demand response | |
| Pazouki et al. | Short term scheduling of multi carrier systems through interruptible load and Energy Storage toward future sustainable energy needs | |
| Rahmanzadeh et al. | Optimal energy management of microgrid based on fcchp in the presence of electric and thermal loads considering energy storage systems | |
| Ma et al. | Optimization of Integrated Energy System Scheduling Considering Stepped Carbon Emission Trading Mechanism and Electricity Heat Demand Response | |
| CN111668882A (en) | Method and device for optimizing output of micro power supply in intelligent energy ring network | |
| Shihab et al. | Optimal resources planning algorithm for a coordinated multi-energy system with industrial loads | |
| Song et al. | Optimal Energy for Grid-Connected Microgrid with Battery Swapping Station and Wind Photovoltaic and Energy Storage |
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 |