CN113131513A - Method for optimizing operation of electric, thermal and gas conversion system with consideration of carbon emission and storage medium - Google Patents

Method for optimizing operation of electric, thermal and gas conversion system with consideration of carbon emission and storage medium Download PDF

Info

Publication number
CN113131513A
CN113131513A CN202110348701.4A CN202110348701A CN113131513A CN 113131513 A CN113131513 A CN 113131513A CN 202110348701 A CN202110348701 A CN 202110348701A CN 113131513 A CN113131513 A CN 113131513A
Authority
CN
China
Prior art keywords
representing
power
gas
cost
time
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110348701.4A
Other languages
Chinese (zh)
Other versions
CN113131513B (en
Inventor
马万里
王旭东
段敬
陈伟
景峰
张雪芹
汪小闯
郑强
蔡苏华
马骏
王昌平
马少松
赵亮
其他发明人请求不公开姓名
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Nari Technology Co Ltd
NARI Nanjing Control System Co Ltd
Information and Telecommunication Branch of State Grid Shanxi Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Nari Technology Co Ltd
NARI Nanjing Control System Co Ltd
Information and Telecommunication Branch of State Grid Shanxi Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, Nari Technology Co Ltd, NARI Nanjing Control System Co Ltd, Information and Telecommunication Branch of State Grid Shanxi Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202110348701.4A priority Critical patent/CN113131513B/en
Publication of CN113131513A publication Critical patent/CN113131513A/en
Application granted granted Critical
Publication of CN113131513B publication Critical patent/CN113131513B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381—Dispersed generators
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
    • H02J3/14—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by switching loads on to, or off from, network, e.g. progressively balanced loading
    • H02J3/144—Demand-response operation of the power transmission or distribution network
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/28—Arrangements for balancing of the load in a network by storage of energy
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/40—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation wherein a plurality of decentralised, dispersed or local energy generation technologies are operated simultaneously
    • 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
    • Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
    • Y02B70/00—Technologies for an efficient end-user side electric power management and consumption
    • Y02B70/30—Systems integrating technologies related to power network operation and communication or information technologies for improving the carbon footprint of the management of residential or tertiary loads, i.e. smart grids as climate change mitigation technology in the buildings sector, including also the last stages of power distribution and the control, monitoring or operating management systems at local level
    • Y02B70/3225—Demand response systems, e.g. load shedding, peak shaving
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S20/00—Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
    • Y04S20/20—End-user application control systems
    • Y04S20/222—Demand response systems, e.g. load shedding, peak shaving

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses an optimal operation configuration method and a storage medium for electricity, heat and gas conversion considering carbon emission, which comprises the following steps: establishing a mathematical model, an operation cost estimation model and a carbon emission model of the electricity-heat-gas conversion component; constructing an objective function and establishing constraint conditions according to the established mathematical model, the operation cost estimation model and the carbon emission model of the electricity-heat-gas conversion component; and solving and calculating to obtain the optimal configuration scheme of electricity, heat and gas according to the objective function and the constraint condition. By optimizing the conversion configuration of the electricity-heat-gas energy, the aims of maximizing social benefits and minimizing carbon emission are achieved, the utilization efficiency of an energy system is improved, and the practical problems of insufficient current energy, environmental deterioration and the like are solved.

Description

Method for optimizing operation of electric, thermal and gas conversion system with consideration of carbon emission and storage medium
Technical Field
The invention belongs to the technical field of comprehensive energy, and relates to an optimal operation configuration method for electricity, heat and gas conversion considering carbon emission and a storage medium.
Background
The traditional energy system is limited to single energy form systems such as electricity, gas, heat (cold) and the like, and complementary advantages and synergistic benefits between the systems cannot be fully played. The energy internet is a novel energy industry development form combining information network and energy links. The energy internet connects the power network with the natural gas network and the heat supply network on the basis of the information network, and finally realizes a novel energy system with multiple complementary energy sources such as electricity, heat, gas and the like. The development of the electricity-heat-gas comprehensive energy system has important significance for promoting the construction of new energy patterns in China and improving the proportion of renewable energy in power generation, and is beneficial to promoting the clean and efficient utilization of energy such as coal, crude oil and the like. The electricity-heat-gas comprehensive energy system comprises an electric power system network, a heat power network and a natural gas network, and various energy coupling units connect different energy networks at the terminals of the networks to realize the mutual conversion of different energy forms. By researching the electricity-heat-gas conversion optimization operation configuration method, the efficiency of the multi-energy system can be improved, the operation cost is reduced, the stability of the multi-energy system is improved, the carbon emission is reduced, and the method has great significance to energy crisis and environmental crisis.
Disclosure of Invention
The invention aims to provide an optimal operation configuration method and a storage medium for conversion of electricity, heat and gas with consideration of carbon emission, which connect different energy networks, realize interconversion of different energy forms, improve the efficiency of a multi-energy system, and reduce the carbon emission to deal with energy crisis and environmental crisis.
In order to achieve the purpose, the invention adopts the following technical scheme: an optimal operation configuration method for electricity, heat and gas conversion considering carbon emission comprises the following steps:
establishing a mathematical model, an operation cost estimation model and a carbon emission model of the electricity-heat-gas conversion component;
constructing an objective function and establishing constraint conditions according to the established mathematical model, the operation cost estimation model and the carbon emission model of the electricity-heat-gas conversion component;
and solving and calculating to obtain the optimal configuration scheme of electricity, heat and gas according to the objective function and the constraint condition.
Further, a mathematical model of an electro-thermal-gas conversion assembly, comprising:
the mathematical model of the ac power grid includes:
the bus voltage and circuit node admittance matrix is:
Vi＝|Vi|∠θi＝|Vi|(cosθi+jsinθi)
Yij＝|Yij|∠θij＝|Yij|(cosθij+jsinθij)＝Gij+jBij
in the formula: viRepresenting the voltage on bus i, YijDenotes the mutual admittance, θ, between circuit node i and circuit node jij＝θi-θj，θiRepresenting the vector angle of the voltage of the bus i, thetajRepresenting the vector angle of the voltage of the bus j, GijRepresenting the conductance between circuit node i and circuit node j, BijRepresents the susceptance between circuit node i and circuit node j;
the injected active and reactive power on the different buses is:
Figure BDA0003001536080000021
Figure BDA0003001536080000022
in the formula: piRepresenting the active power, Q, of circuit node iiRepresenting reactive power, V, of circuit node ijRepresenting the voltage on the bus j, wherein N is the total number of circuit nodes;
the power balance equation of the circuit node is as follows:
Figure BDA0003001536080000023
Figure BDA0003001536080000024
in the formula: pg,iAnd Pd,iActive power, Δ P, representing the amount of generated power on bus i and the active power of the amount of consumed power, respectivelyiRepresenting the difference between the active power emitted and the active power consumed on the bus i, Qg,iAnd Qd,iReactive power, Δ Q, representing the generated energy on the bus i and the consumed power respectivelyiRepresenting the difference between the reactive power emitted and the reactive power consumed on the bus i;
the mathematical model of the heating system comprises:
the flow continuity is expressed as:
A×mpipe＝mnode
in the formula: m ispipeIs a vector representing the mass flow rate in each pipe, mnodeA vector representing the mass flow rate through each heat network node, a representing a heat network node correlation matrix;
HP＝cpmnode(Ts-To)
in the formula: HP represents thermal power, cpRepresents the specific heat capacity of water, TsAnd ToRepresenting the supply temperature and the outlet temperature, respectively;
the mathematical model of the natural gas system includes:
Figure BDA0003001536080000031
in the formula: qmnRepresenting gas flows at natural gas node m and natural gas node n, KmnRepresenting a natural gas characteristic factor, PmAnd PnThe pressure of the natural gas node m and the pressure of the natural gas node n are shown.
Further, the operation cost estimation model comprises:
the operation cost of the power grid generator set is as follows:
Figure BDA0003001536080000032
in the formula: cgIs the operating cost of the grid generator, aG，bGAnd cGIs the cost coefficient of the grid generator, gGAnd eGA coefficient representing the effect of the grid generator valve point loading,
Figure BDA0003001536080000033
the representation is the amount of power exchanged by the utility grid at time t,
Figure BDA0003001536080000034
representing a minimum value of the electric energy exchange quantity of the public power grid;
the cost function of the gas generator is as follows:
Figure BDA0003001536080000035
in the formula: cgfIs the operating cost of the gas generator, aGf，bGfAnd cGfAre all cost factors, g, of gas generatorsGfAnd eGfIs a coefficient representing the valve point load effect of the gas generator,
Figure BDA0003001536080000041
the representation is the amount of power exchanged by the gas generator at time t,
Figure BDA0003001536080000042
is the minimum value of the electric energy exchange quantity of the gas generator;
the natural gas supply cost function is:
Figure BDA0003001536080000043
in the formula: cgsRepresents the total cost of the natural gas supply,
Figure BDA0003001536080000044
representing the cost factor for the vth natural gas supply facility,
Figure BDA0003001536080000045
representing the supply quantity of a vth natural gas supply device at the time t;
the cost of wind power generation is:
Figure BDA0003001536080000046
in the formula:
Figure BDA0003001536080000047
represents the total cost of the w-th wind power generation,
Figure BDA0003001536080000048
represents a cost coefficient of planned wind power generation of the w-th wind power generator,
Figure BDA0003001536080000049
representing the active power produced by the w-th wind turbine at time t,
Figure BDA00030015360800000410
represents a penalty cost for the w-th wind turbine power,
Figure BDA00030015360800000411
representing the available work produced by the w-th wind farm at time tThe ratio of the total weight of the particles,
Figure BDA00030015360800000412
representing the spare cost of the w wind turbine;
the operating cost of the photovoltaic generator is as follows:
Figure BDA00030015360800000413
in the formula: cpvRepresents the total cost of the photovoltaic power generation,
Figure BDA00030015360800000414
represents the power generation coefficient, pv, of the qth photovoltaic generatorq,tRepresenting the active power of the qth photovoltaic generator at time t,
Figure BDA00030015360800000415
represents the penalty cost of the power of the qth photovoltaic generator,
Figure BDA00030015360800000416
representing the available power produced by the qth photovoltaic generator at time t,
Figure BDA00030015360800000417
representing the standby cost of the qth photovoltaic generator;
the cost of the energy storage device is:
Figure BDA00030015360800000418
in the formula: cSDRepresents the total cost of the energy storage device, CSDD，CSDCRespectively representing the charging cost of the charging device of the second SDD station and the discharging cost of the discharging device of the second SDC station,
Figure BDA0003001536080000051
and
Figure BDA0003001536080000052
respectively showing charging and discharging power of the r-th energy storage equipment at the time t, NSDD，NSDCRespectively representing the number of charging and discharging devices;
the total cost function of the heating unit is as follows:
Figure BDA0003001536080000053
in the formula: cHOURepresents the total operating cost of the heating unit,
Figure BDA0003001536080000054
all represent cost factors for the xth heating unit,
Figure BDA0003001536080000055
representing the thermal power of the x-th heating unit at the time t;
the CHP cost of the cogeneration unit is as follows:
Figure BDA0003001536080000056
in the formula: cCHPRepresents the overall cost of the CHP generator,
Figure BDA0003001536080000057
and
Figure BDA0003001536080000058
are the y-th CHP generator cost coefficients,
Figure BDA0003001536080000059
shows the active power generated by the y-th CHP generator at time t,
Figure BDA00030015360800000510
the heat power generated by the y-th CHP generator at the time t is shown;
the cost of the electric conversion P2G device is as follows:
Figure BDA00030015360800000511
in the formula: cP2GRepresents the total operating cost of the P2G plant, cP2GIs the cost factor of the P2G device,
Figure BDA00030015360800000512
representing the amount of conversion of electrical energy to natural gas at time t in the z-th P2G plant,
Figure BDA00030015360800000513
representing the real power of the conversion of electrical energy to natural gas at time t in the z-th P2G plant,
Figure BDA00030015360800000514
representing the conversion factor of electrical energy to natural gas in the z-th P2G plant.
Further, the carbon emission model is as follows:
Figure BDA00030015360800000515
in the formula: eGRepresents the total pollutant emission, aE，bE，dE，γE，δEIs the emission coefficient of the thermal power generating unit.
Further, the objective function F is:
F＝SW-CE
SW＝RD-cost
RD＝λePD+λheatPheat+λgasPgas
Figure BDA0003001536080000061
cost＝Cg+Cgf+Cgs+CPw+Cpv+CSD+CHOU+CCHP+CP2G
SWmaxis the maximum value of the social benefit,
Figure BDA0003001536080000062
maximum amount of carbon emission, λeIs the cost of the electrical energy consumed, λheatIs the cost of heat energy, λgasIs the cost of consuming gas energy, PD，PheatAnd PgasRespectively representing the demand for electrical energy, the demand for thermal energy and the demand for natural gas.
Further, the constraint conditions include:
active power balance:
Figure BDA0003001536080000063
in the formula:
Figure BDA0003001536080000064
representing the active power at time t of the f-th gas generator, NGIndicates the number of the gas generators,
Figure BDA0003001536080000065
representing the active power at time t of the vth natural gas supply facility, NgfIndicating the number of natural gas supply facilities,
Figure BDA0003001536080000066
representing the active power at time t of the w-th wind turbine, NwIndicating the number of wind generators, pvq,tRepresenting the active power of the qth photovoltaic generator at time t, NpvIndicates the number of the photovoltaic generators,
Figure BDA0003001536080000067
indicating the active power at time t of the y-th CHP device, NCHPIndicates the number of the CHP devices,
Figure BDA0003001536080000068
indicating the r-th energy storage equipment at time tWork power, NBIndicates the number of the energy storage devices,
Figure BDA0003001536080000069
represents the active power, N, of the z-th station P2G device at time tP2GDenotes the number of P2G devices, PDtRepresenting the total power demand at time t, PlossRepresents the power loss;
Figure BDA00030015360800000610
and (3) constraint of bus voltage and branch flow:
Figure BDA0003001536080000071
in the formula: vi max，Vi minRepresents the ith bus voltage maximum and minimum, Sflow,iShowing the tidal flow distribution on the ith branch,
Figure BDA0003001536080000072
represents the maximum value of the tidal flow distribution on the ith branch;
and (3) restricting the climbing rate of the wind driven generator:
Figure BDA0003001536080000073
in the formula: UR and DR respectively represent the rising rate constraint and the falling rate constraint of any wind turbine generator set, Pw,tRepresenting the power of the w-th wind driven generator set at the time t;
output power constraint:
Figure BDA0003001536080000074
in the formula:
Figure BDA0003001536080000075
respectively representing the maximum value and the minimum value of the output power of the f-th gas generator,
Figure BDA0003001536080000076
showing the output power of the f-th gas generator at the moment t,
Figure BDA0003001536080000077
respectively representing the maximum value and the minimum value of the output power of the y-th cogeneration generator,
Figure BDA0003001536080000078
represents the maximum value of active power generated by the wind turbine,
Figure BDA0003001536080000079
denotes the maximum value of the output power of the z-th stage P2G device, PwmaxRepresents the maximum power value, pv, of the w-th wind turbinemaxRepresenting the maximum power of the qth photovoltaic generator;
thermal node balance constraint:
Figure BDA00030015360800000710
in the formula:
Figure BDA00030015360800000711
indicating the thermal power, N, at time t of the y-th CHP deviceCHPIndicates the number of the CHP devices,
Figure BDA00030015360800000712
indicating the thermal power, N, of the xth heating unit at time tHOUIndicates the number of the heat supply units,
Figure BDA00030015360800000713
indicating the thermal power, N, of the s-th other device at time tHSNumber of other devices, HDtIndicating the heat demand at time t, HlossRepresents heat loss;
heat restraint of cogeneration and heating units:
Figure BDA0003001536080000081
in the formula:
Figure BDA0003001536080000082
represents the maximum value and the minimum value of the heat production quantity of the y-th cogeneration generator,
Figure BDA0003001536080000083
the maximum value and the minimum value of the heat production quantity of the x-th heat supply unit set are represented;
gas node balance constraint:
Figure BDA0003001536080000084
in the formula:
Figure BDA0003001536080000085
indicating the gas output of the z-th station P2G device at time t,
Figure BDA0003001536080000086
representing the gas output of the vth natural gas plant at time t,
Figure BDA0003001536080000087
indicating the gas consumption of the y-th CHP device at time t,
Figure BDA0003001536080000088
denotes the gas consumption, Q, of the f-th gas generator unit at time tD.tIndicating the gas demand at time t, QlossRepresents the amount of gas lost;
pressure and airflow constraints:
Figure BDA0003001536080000089
in the formula:
Figure BDA00030015360800000810
representing the maximum value and the minimum value of the pressure intensity between the natural gas node m and the natural gas node n;
Figure BDA00030015360800000811
the pressure intensity between a natural gas node m and a natural gas node n at the moment t of the natural gas equipment is represented;
charge and discharge power constraint:
Figure BDA00030015360800000812
in the formula: emaxThe maximum energy storage device capacity of different storage types is represented, and delta t represents time difference;
energy restraint in the energy storage device:
Figure BDA00030015360800000813
in the formula:
Figure BDA0003001536080000091
representing the maximum and minimum values of energy in the r-th energy storage device,
Figure BDA0003001536080000092
representing the initial energy of the r-th energy storage device.
A computer readable storage medium storing one or more programs, the one or more programs comprising instructions, which when executed by a computing device, cause the computing device to perform any of the above methods of electrical, thermal, gas conversion optimization operating configuration that accounts for carbon emissions.
A computing device comprising one or more processors, memory, and one or more programs stored in the memory and configured to be executed by the one or more processors, the one or more programs including instructions for performing any of the above electrical, thermal, gas shift optimization operational configuration methods that take into account carbon emissions.
This beneficial effect:
by optimizing the conversion configuration of the electricity-heat-gas energy, the aims of maximizing social benefits and minimizing carbon emission are achieved, the utilization efficiency of an energy system is improved, the carbon emission is reduced to deal with energy crisis and environmental crisis, and the practical problems of insufficient energy, environmental deterioration and the like are relieved.
Drawings
FIG. 1 is a flow chart of a method for configuring an electric-to-heat-gas shift optimization operation in consideration of carbon emissions;
FIG. 2 is a schematic diagram of electric-to-heat-to-gas energy conversion;
FIG. 3 is a schematic diagram of an IEEE 69 node standard test system;
FIG. 4 is a schematic diagram of a heating system configuration;
FIG. 5 is a schematic of a natural gas system configuration;
FIG. 6 is a schematic diagram of electrical power output;
FIG. 7 is a schematic of thermal energy output;
fig. 8 is a schematic illustration of natural gas output.
Detailed Description
The technical scheme of the invention is further explained in detail by combining the attached drawings:
example 1:
as shown in fig. 1, a method for optimizing operation and configuration of electric, thermal and gas conversion considering carbon emission comprises the following steps:
step 1, establishing a mathematical model, an operation cost estimation model and a carbon emission model of the electricity-heat-gas conversion component:
the conversion component refers to a component of an electric-heat-gas energy junction of an alternating current power grid and a heat supply system.
(1) The mathematical model of the electro-thermal-gas conversion assembly comprises:
1) AC power network
The invention provides a bus voltage V and a circuit node admittance matrix Y of a system in a polar coordinate form, which are shown as the following formula:
Vi＝|Vi|∠θi＝|Vi|(cosθi+jsinθi) (4)
Yij＝|Yij|∠θij＝|Yij|(cosθij+jsinθij)＝Gij+jBij (5)
in the formula: the bus and the circuit nodes are in one-to-one correspondence, the total number is the same, ViRepresenting the voltage on bus i, YijDenotes the mutual admittance, θ, between circuit node i and circuit node jij＝θi-θj，θiRepresenting the vector angle of the voltage of the bus i, thetajRepresenting the vector angle of the voltage of the bus j, GijRepresenting the conductance between circuit node i and circuit node j, BijRepresenting the susceptance between circuit node i and circuit node j.
The injected active and reactive power on the different buses can be expressed as:
Figure BDA0003001536080000101
Figure BDA0003001536080000102
in the formula: piRepresenting the active power, Q, of circuit node iiRepresenting the reactive power at circuit node i,
Vjrepresenting the voltage on bus j, and N is the total number of circuit nodes.
The circuit node power balance equation can be expressed as follows:
Figure BDA0003001536080000103
Figure BDA0003001536080000104
in the formula: pg,iAnd Pd,iThe active power representing the amount of generated power on the bus i and the active power representing the amount of consumed power are shown. Delta PiRepresenting the difference between the real power drawn and the real power consumed on the bus i. Qg,iAnd Qd,iAnd reactive power representing the generated energy on the bus i and the reactive power representing the consumed power. Delta QiRepresenting the difference between the reactive power drawn and the reactive power consumed on the bus i.
2) Heating system
The mass flow around any heat network node is equal to the sum of the mass flow entering the heat network node, the mass flow leaving the heat network node and the flow consumption at the heat network node. The flow continuity is expressed as:
A×mpipe＝mnode (10)
in the formula: m ispipeIs a vector representing the mass flow rate (kg/s) in each pipe, mnodeA vector representing the mass flow rate (kg/s) of each heat network node passing through, and a represents the heat network node correlation matrix. Head loss in a pipe is the change in pipe pressure (in meters) due to pipe friction.
The thermal power can be calculated by:
HP＝cpmnode(Ts-To) (11)
in the formula: HP represents thermal power, cpRepresents the specific heat capacity (J/Kg ℃) of water, TsAnd ToRepresenting the supply temperature and the outlet temperature (. degree. C.) respectively.
3) Natural gas system
Modeling of natural gas systems is similar to the analysis of heating systems, with the following two assumptions: (1) the gas flow temperature remains constant, assuming no change in the temperature of the gas piping. (2) The difference in the height of the pipes is ignored, i.e. the two pipes are horizontal pipes. The invention gives the flow equation of natural gas through any pipeline as shown in the following formula:
Figure BDA0003001536080000111
in the formula: qmnGas flow (m3/h), K, representing natural gas node m and natural gas node nmnRepresenting a natural gas characteristic factor, PmAnd PnRepresenting the Pressure (PSIA) of natural gas node m and natural gas node n.
(2) The operating cost estimation model includes:
1) cost of the grid
The operating cost of the grid generator set can be expressed as:
Figure BDA0003001536080000121
in the formula: cgIs the operating cost of the grid generator, aG，bGAnd cGIs the cost coefficient of the grid generator, gGAnd eGA coefficient representing the effect of the grid generator valve point loading,
Figure BDA0003001536080000122
the representation is the amount of power exchanged by the utility grid at time t,
Figure BDA0003001536080000123
representing the minimum value of the electric energy exchange of the utility grid.
2) Cost of gas-fired generator
The gas generator cost function can be expressed as:
Figure BDA0003001536080000124
in the formula: cgfIs the operating cost of the gas generator, aGf，bGfAnd cGfAre all cost factors, g, of gas generatorsGfAnd eGfIs a coefficient representing the valve point load effect of the gas generator,
Figure BDA0003001536080000125
the expression is the amount of electric energy exchanged, i.e. the active power,
Figure BDA0003001536080000126
is the minimum value of the electric energy exchange quantity of the gas generator.
3) Cost of natural gas supply
The natural gas supply cost function may be expressed as:
Figure BDA0003001536080000127
in the formula: cgsRepresents the total cost of the natural gas supply,
Figure BDA0003001536080000128
representing the cost factor for the vth natural gas supply facility,
Figure BDA0003001536080000129
indicating the supply of the vth natural gas supply facility at time t.
4) Cost of wind power generation
The cost of wind power generation can be expressed as:
Figure BDA00030015360800001210
in the formula:
Figure BDA00030015360800001211
represents the total cost of the w-th wind power generation,
Figure BDA00030015360800001212
represents a cost coefficient of planned wind power generation of the w-th wind power generator,
Figure BDA00030015360800001213
representing the active power produced by the w-th wind turbine at time t,
Figure BDA00030015360800001214
represents a penalty cost for the w-th wind turbine power,
Figure BDA00030015360800001215
it is usually quantified by the overestimated output, underestimated output and penalty cost coefficients of the wind power.
Figure BDA00030015360800001216
Representing the available power (MW) generated by the w-th wind farm at time t,
Figure BDA00030015360800001217
representing the spare cost for the w wind turbine.
5) Cost of photovoltaic power generation
The operating cost of a photovoltaic generator is expressed as follows:
Figure BDA0003001536080000131
in the formula: cpvRepresents the total cost of the photovoltaic power generation,
Figure BDA0003001536080000132
represents the power generation coefficient, pv, of the qth photovoltaic generatorq,tRepresenting the active power of the qth photovoltaic generator at time t,
Figure BDA0003001536080000133
represents the penalty cost of the power of the qth photovoltaic generator,
Figure BDA0003001536080000134
representing the available power produced by the qth photovoltaic generator at time t,
Figure BDA0003001536080000135
representing the standby cost of the qth photovoltaic generator.
6) Cost of energy storage device
The cost of the energy storage device can be calculated by:
Figure BDA0003001536080000136
in the formula: cSDRepresents the total cost of the energy storage device, CSDD，CSDCRespectively represent the charge cost and the discharge cost,
Figure BDA0003001536080000137
and
Figure BDA0003001536080000138
respectively showing charging and discharging power of the r-th energy storage equipment at the time t, NSDD，NSDCRespectively representing the number of charging and discharging devices.
7) Cost of heat supply unit
The role of heating units (HOU) in energy hubs is mainly to provide large amounts of energy to district heating systems. In the presence of cogeneration units, the heating unit is normally only used during periods of high demand, and the total cost function can be expressed as:
Figure BDA0003001536080000139
in the formula: cHOURepresents the total operating cost of the heating unit,
Figure BDA00030015360800001310
all represent cost factors for the xth heating unit,
Figure BDA00030015360800001311
indicating the thermal power of the x-th heating unit at time t.
8) Combined heat and power generation unit (CHP) cost
The CHP device has a convex cost function with a variable of PCHPAnd QCHPAs shown in the following formula:
Figure BDA00030015360800001312
in the formula: cCHPRepresents the overall cost of the CHP generator,
Figure BDA00030015360800001313
and
Figure BDA00030015360800001314
are the cost coefficients of the y-th CHP generator, and the size of the coefficients is determined according to empirical values.
Figure BDA0003001536080000141
Shows the active power generated by the y-th CHP generator at time t,
Figure BDA0003001536080000142
and the heat power generated by the y-th CHP generator at the time t is shown.
9) Cost of electric gas (P2G) conversion equipment
The total operating cost of the P2G plant includes the cost of electricity obtained minus the profit of selling the natural gas produced by P2G, as shown by the following equation:
Figure BDA0003001536080000143
in the formula: cP2GRepresents the total operating cost of the P2G plant, cP2GIs the cost factor of the P2G device,
Figure BDA0003001536080000144
representing the amount of conversion of electrical energy to natural gas at time t in the z-th P2G plant,
Figure BDA0003001536080000145
representing the activity of the conversion of electric energy to natural gas at time t in the z-th P2G plantThe power of the electric motor is controlled by the power controller,
Figure BDA0003001536080000146
representing the conversion factor of electrical energy to natural gas in the z-th P2G plant.
(3) Carbon emission model
The total pollutant emissions can be expressed as a function of the power generation as shown in the following equation:
Figure BDA0003001536080000147
in the formula: eGRepresents the total pollutant emission, aE，bE，dE，γE，δEIs the emission coefficient of the thermal power generating unit.
Step 2, constructing an objective function and establishing constraint conditions according to the established mathematical model, the operation cost estimation model and the carbon emission model of the electricity-heat-gas conversion component;
(1) objective function
The overall goal of system optimization is to achieve more social benefits by reducing the total cost of system operation. In addition, the total carbon emissions should be minimized. The total operating cost function can be expressed as:
cost＝Cg+Cgf+Cgs+CPw+Cpv+CSD+CHOU+CCHP+CP2G (23)
the revenue obtained by selling energy to the consumer can be expressed as:
RD＝λePD+λheatPheat+λgasPgas (24)
in the formula: lambda [ alpha ]eIs the cost of the electrical energy consumed, λheatIs the cost of heat energy, λgasIs the cost of consuming gas energy, PD，PheatAnd PgasRespectively representing the demand for electrical energy, the demand for thermal energy and the demand for natural gas. Social benefits can be calculated from the following formula:
SW＝RD-cost (25)
the carbon emissions are converted to emission costs by a maximum penalty price factor h. Therefore, on the premise of meeting the load demand and the operation constraint, the social benefit obtained by the power plant is maximized and the carbon emission is minimized through solving the optimization problem. The highest penalty price factor h is the maximum value SW of social benefitmaxAnd maximum amount of carbon emission
Figure BDA0003001536080000151
The ratio therebetween is shown by the following formula:
Figure BDA0003001536080000152
cost C corresponding to carbon emissionsECan be calculated from the following formula:
CE＝hEG (27)
the overall objective function F is as follows;
Max→F＝SW-CE (28)
(2) the constraint conditions include:
1) active power balance:
Figure BDA0003001536080000153
in the formula:
Figure BDA0003001536080000154
representing the active power at time t of the f-th gas generator, NGIndicates the number of the gas generators,
Figure BDA0003001536080000155
representing the active power at time t of the vth natural gas supply facility, NgfIndicating the number of natural gas supply facilities,
Figure BDA0003001536080000156
representing the active power at time t of the w-th wind turbine, NwIndicating the number of wind generators, pvq,tRepresenting the active power of the qth photovoltaic generator at time t, NpvIndicates the number of the photovoltaic generators,
Figure BDA0003001536080000157
indicating the active power at time t of the y-th CHP device, NCHPIndicates the number of the CHP devices,
Figure BDA0003001536080000158
representing the active power at time t of the r-th energy storage device, NBIndicates the number of the energy storage devices,
Figure BDA0003001536080000159
represents the active power, N, of the z-th station P2G device at time tP2GDenotes the number of P2G devices, PDtRepresenting the total power demand at time t, PlossRepresents the power loss, PlossThe calculation formula is as follows:
Figure BDA0003001536080000161
in the formula: gijRepresenting the conductance, V, between circuit node i and circuit node jiAnd VjRepresenting the voltages on bus i and bus j, respectively, thetaijIs the voltage vector angular difference of the generatrices i and j.
2) And (3) constraint of bus voltage and branch flow:
Figure BDA0003001536080000162
in the formula: vi max，Vi minRepresents the ith bus voltage maximum and minimum, Sflow,iShowing the tidal flow distribution on the ith branch,
Figure BDA0003001536080000163
represents the maximum value of the tidal flow distribution on the ith branch. The bus and the branch are in one-to-one relation;
3) and (3) restricting the climbing rate of the wind driven generator:
Figure BDA0003001536080000164
in the formula: UR and DR respectively represent the rising rate constraint and the falling rate constraint of any wind turbine generator set, Pw,tAnd the power of the w-th wind power generator set at the time t is shown.
4) Output power constraint:
Figure BDA0003001536080000165
in the formula:
Figure BDA0003001536080000166
respectively representing the maximum value and the minimum value of the output power of the f-th gas generator,
Figure BDA0003001536080000167
Figure BDA0003001536080000168
respectively representing the maximum value and the minimum value of the output power of the y-th cogeneration generator,
Figure BDA0003001536080000169
denotes the maximum value of the output power of the z-th stage P2G device, PwmaxRepresents the maximum power value, pv, of the w-th wind turbinemaxRepresents the maximum power of the qth photovoltaic generator.
5) Thermal node balance constraint:
Figure BDA0003001536080000171
in the formula:
Figure BDA0003001536080000172
indicating the thermal power, N, at time t of the y-th CHP deviceCHPIndicates the number of the CHP devices,
Figure BDA0003001536080000173
indicating the thermal power, N, of the xth heating unit at time tHOUIndicates the number of the heat supply units,
Figure BDA0003001536080000174
represents the thermal power at time t of the s-th other plant (gas, natural gas, wind, etc.), NHSNumber of other devices, HDtIndicating the heat demand at time t, HlossIndicating heat loss.
6) Heat restraint of cogeneration and heating units:
Figure BDA0003001536080000175
in the formula:
Figure BDA0003001536080000176
represents the maximum value and the minimum value of the heat production quantity of the y-th cogeneration generator,
Figure BDA0003001536080000177
and the maximum value and the minimum value of the heat production quantity of the x-th heat supply unit set are represented.
7) Gas node balance constraint:
Figure BDA0003001536080000178
in the formula:
Figure BDA0003001536080000179
indicating the gas output of the z-th station P2G device at time t,
Figure BDA00030015360800001710
representing the gas output of the vth natural gas plant at time t,
Figure BDA00030015360800001711
indicating the gas consumption of the y-th CHP device at time t,
Figure BDA00030015360800001712
denotes the gas consumption, Q, of the f-th gas generator unit at time tD.tIndicating the gas demand at time t, QlossIndicating the amount of gas lost.
8) Pressure and airflow constraints:
Figure BDA00030015360800001713
in the formula:
Figure BDA00030015360800001714
and the maximum value and the minimum value of the pressure intensity between the natural gas node m and the natural gas node n are shown.
Figure BDA0003001536080000181
And the pressure intensity between the natural gas node m and the natural gas node n at the moment t of the natural gas equipment is shown.
9) Charge and discharge power constraint:
Figure BDA0003001536080000182
in the formula: emaxRepresenting the maximum energy storage device capacity for different storage types and deltat representing the time difference.
10) Energy restraint in the energy storage device:
Figure BDA0003001536080000183
in the formula:
Figure BDA0003001536080000184
representing the maximum and minimum values of energy in the r-th energy storage device,
Figure BDA0003001536080000185
representing the initial energy of the r-th energy storage device.
And 3, solving and calculating to obtain an optimal configuration scheme of electricity, heat and gas according to the objective function and the constraint conditions.
Example 2:
a computer readable storage medium storing one or more programs, the one or more programs comprising instructions, which when executed by a computing device, cause the computing device to perform any of the electrical, thermal, gas shift optimization operating configuration methods of embodiment 1 that take into account carbon emissions.
A computing device comprising one or more processors, memory, and one or more programs stored in the memory and configured to be executed by the one or more processors, the one or more programs comprising instructions for performing any of the electrical, thermal, gas shift optimization operation configuration methods of embodiment 1 that take into account carbon emissions.
The feasibility of the proposed method was verified by calculation:
in order to verify the applicability and validity of the proposed algorithm, it is applied to an IEEE 69 node standard test system. The system consists of 9 conventional generators, 48 loads and 68 branches (transmission line and transformer) as shown in fig. 3. Wherein the equipment parameters are shown in table 1. In addition, an energy storage device having a capacity of 60MWh is connected to node 60. The heating system consists of two CHP units, an HOU and a heat storage facility, wherein the heat storage facility supplies heat to a local heating network, as shown in fig. 4. The natural gas system configuration is shown in fig. 5, GF representing a gas turbine.
TABLE 1 Equipment parameters Table
Figure BDA0003001536080000191
The present invention considers 5 examples to illustrate the role of renewable energy, energy storage devices and P2G devices in energy hub operations.
Example 1: the original energy hub configuration was investigated without including any renewable energy or energy storage devices.
Example 2: an energy storage device is added on the basis of the formula 1 for research.
Example 3: the energy hub includes renewable energy and P2G devices, but does not include energy storage devices.
Example 4: the energy hub contains only renewable energy sources.
Example 5: the energy hub comprises all energy output.
The invention calculates the total operation cost, the carbon emission and the system loss in each calculation example respectively. The electric, thermal, and gas costs of each node are calculated with the electric, thermal, and gas costs of the bus 2 as references and the electric, thermal, and gas costs of each node as indexes, taking the loss into consideration, as shown in table 2.
TABLE 2 operating results under different energy hub configurations
Figure BDA0003001536080000192
Figure BDA0003001536080000201
As can be seen from table 2:
1. in example 1, the CHP unit is the main source of electricity and heat supply during the day due to the low price of natural gas, and in addition, the loss of natural gas is increased due to the addition of the feeding path of the gas system. However, as daytime loads increase, CHP cannot meet these loads and can only be purchased from HOU, thus increasing operating costs and carbon emissions.
2. In example 2, with the addition of the energy storage device, the system performance is greatly improved. The carbon emission and the energy consumption are slightly reduced. In this case, the excess energy in the energy hub is stored for use when the energy source is insufficient.
3. Compared with the embodiments 1 and 2, the embodiment 3 has the advantages that due to the addition of renewable energy, social benefits, carbon emission, energy consumption and the like are obviously improved, but because the energy storage device is not arranged in the energy hub, the loss is increased.
4. The absence of the P2G device of example 4 resulted in a reduction in social benefit, carbon emissions and energy consumption compared to example 3.
5. In example 5, the energy storage device and the P2G device are operated simultaneously due to the renewable energy, thereby reducing the electric energy purchased from the power grid and the thermal power generating unit. Thus, the performance parameters as well as the stability of the energy hub are improved. The simulation results of the optimal output in the embodiment 5 are shown in FIGS. 6 to 8.
As can be seen from the results of comparative analysis of the examples, when the thermal power generating unit becomes the main source of energy in the energy hub, both the operating cost and the carbon emission increase. When energy storage devices are connected into the energy hub, the energy hub is improved in terms of carbon emission, system loss and social benefits. Similarly, when renewable energy is added to the energy hub, various aspects of the system are greatly improved. Without the P2G unit, the total load is reduced, thereby increasing the power supplied to the grid and reducing social benefits, carbon emissions and losses. In example 5, when the energy hub includes renewable energy, energy storage device and P2G device, the parameters and stability of the energy hub are greatly improved. The simulation results of the optimal output in the embodiment 5 are shown in FIGS. 6 to 8.
The above description is only an embodiment of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art can understand that the modifications or substitutions within the technical scope of the present invention are included in the scope of the present invention, and therefore, the scope of the present invention should be subject to the protection scope of the claims.

Claims (8)

1. An optimal operation configuration method for electricity, heat and gas conversion considering carbon emission is characterized by comprising the following steps:
establishing a mathematical model, an operation cost estimation model and a carbon emission model of the electricity-heat-gas conversion component;
constructing an objective function and establishing constraint conditions according to the established mathematical model, the operation cost estimation model and the carbon emission model of the electricity-heat-gas conversion component;
and solving and calculating to obtain the optimal configuration scheme of electricity, heat and gas according to the objective function and the constraint condition.
2. The method of claim 1, wherein the mathematical model of the electric-to-heat-gas conversion module comprises:
the mathematical model of the ac power grid includes:
the bus voltage and circuit node admittance matrix is:
Vi＝|Vi|∠θi＝|Vi|(cosθi+jsinθi)
Yij＝|Yij|∠θij＝|Yij|(cosθij+jsinθij)＝Gij+jBij
in the formula: viRepresenting the voltage on bus i, YijDenotes the mutual admittance, θ, between circuit node i and circuit node jij＝θi-θj，θiRepresenting the vector angle of the voltage of the bus i, thetajRepresenting the vector angle of the voltage of the bus j, GijRepresenting the conductance between circuit node i and circuit node j, BijRepresents the susceptance between circuit node i and circuit node j;
the injected active and reactive power on the different buses is:
Figure FDA0003001536070000011
Figure FDA0003001536070000012
in the formula: piRepresenting the active power, Q, of circuit node iiRepresenting reactive power, V, of circuit node ijRepresenting the voltage on the bus j, wherein N is the total number of circuit nodes;
the power balance equation of the circuit node is as follows:
Figure FDA0003001536070000013
Figure FDA0003001536070000021
in the formula: pg,iAnd Pd,iActive power, Δ P, representing the amount of generated power on bus i and the active power of the amount of consumed power, respectivelyiRepresenting the difference between the active power emitted and the active power consumed on the bus i, Qg,iAnd Qd,iReactive power, Δ Q, representing the generated energy on the bus i and the consumed power respectivelyiRepresenting the difference between the reactive power emitted and the reactive power consumed on the bus i;
the mathematical model of the heating system comprises:
the flow continuity is expressed as:
A×mpipe＝mnode
in the formula: m ispipeIs a vector representing the mass flow rate in each pipe, mnodeA vector representing the mass flow rate through each heat network node, a representing a heat network node correlation matrix;
HP＝cpmnode(Ts-To)
in the formula: HP represents thermal power, cpRepresents the specific heat capacity of water, TsAnd ToRepresenting the supply temperature and the outlet temperature, respectively;
the mathematical model of the natural gas system includes:
Figure FDA0003001536070000022
in the formula: qmnRepresenting gas flows at natural gas node m and natural gas node n, KmnRepresenting a natural gas characteristic factor, PmAnd PnThe pressure of the natural gas node m and the pressure of the natural gas node n are shown.
3. The method of claim 2, wherein the operation cost estimation model comprises:
the operation cost of the power grid generator set is as follows:
Figure FDA0003001536070000023
in the formula: cgIs the operating cost of the grid generator, aG，bGAnd cGIs the cost coefficient of the grid generator, gGAnd eGA coefficient representing the effect of the grid generator valve point loading,
Figure FDA0003001536070000031
the representation is the amount of power exchanged by the utility grid at time t,
Figure FDA0003001536070000032
representing a minimum value of the electric energy exchange quantity of the public power grid;
the cost function of the gas generator is as follows:
Figure FDA0003001536070000033
in the formula: cgfIs the operating cost of the gas generator, aGf，bGfAnd cGfAre all cost factors, g, of gas generatorsGfAnd eGfIs a coefficient representing the valve point load effect of the gas generator,
Figure FDA0003001536070000034
the representation is the amount of power exchanged by the gas generator at time t,
Figure FDA0003001536070000035
is the minimum value of the electric energy exchange quantity of the gas generator;
the natural gas supply cost function is:
Figure FDA0003001536070000036
in the formula: cgsRepresents the total cost of the natural gas supply,
Figure FDA0003001536070000037
representing the cost factor for the vth natural gas supply facility,
Figure FDA0003001536070000038
representing the supply quantity of a vth natural gas supply device at the time t;
the cost of wind power generation is:
Figure FDA0003001536070000039
in the formula:
Figure FDA00030015360700000310
represents the total cost of the w-th wind power generation,
Figure FDA00030015360700000311
indicating the w-th wind turbineThe cost factor for planning the wind power generation,
Figure FDA00030015360700000312
representing the active power produced by the w-th wind turbine at time t,
Figure FDA00030015360700000313
represents a penalty cost for the w-th wind turbine power,
Figure FDA00030015360700000314
representing the available power generated by the w-th wind farm at time t,
Figure FDA00030015360700000315
representing the spare cost of the w wind turbine;
the operating cost of the photovoltaic generator is as follows:
Figure FDA00030015360700000316
in the formula: cpvRepresents the total cost of the photovoltaic power generation,
Figure FDA00030015360700000317
represents the power generation coefficient, pv, of the qth photovoltaic generatorq,tRepresenting the active power of the qth photovoltaic generator at time t,
Figure FDA00030015360700000318
represents the penalty cost of the power of the qth photovoltaic generator,
Figure FDA0003001536070000041
representing the available power produced by the qth photovoltaic generator at time t,
Figure FDA0003001536070000042
represents the standby cost of the qth photovoltaic generator；
The cost of the energy storage device is:
Figure FDA0003001536070000043
in the formula: cSDRepresents the total cost of the energy storage device, CSDD，CSDCRespectively representing the charging cost of the charging device of the second SDD station and the discharging cost of the discharging device of the second SDC station,
Figure FDA0003001536070000044
and
Figure FDA0003001536070000045
respectively showing charging and discharging power of the r-th energy storage equipment at the time t, NSDD，NSDCRespectively representing the number of charging and discharging devices;
the total cost function of the heating unit is as follows:
Figure FDA0003001536070000046
in the formula: cHOURepresents the total operating cost of the heating unit,
Figure FDA0003001536070000047
all represent cost factors for the xth heating unit,
Figure FDA0003001536070000048
representing the thermal power of the x-th heating unit at the time t;
the CHP cost of the cogeneration unit is as follows:
Figure FDA0003001536070000049
in the formula: cCHPRepresents the overall cost of the CHP generator,
Figure FDA00030015360700000410
and
Figure FDA00030015360700000411
are the y-th CHP generator cost coefficients,
Figure FDA00030015360700000412
shows the active power generated by the y-th CHP generator at time t,
Figure FDA00030015360700000413
the heat power generated by the y-th CHP generator at the time t is shown;
the cost of the electric conversion P2G device is as follows:
Figure FDA00030015360700000414
in the formula: cP2GRepresents the total operating cost of the P2G plant, cP2GIs the cost factor of the P2G device,
Figure FDA00030015360700000415
representing the amount of conversion of electrical energy to natural gas at time t in the z-th P2G plant,
Figure FDA00030015360700000416
representing the real power of the conversion of electrical energy to natural gas at time t in the z-th P2G plant,
Figure FDA00030015360700000417
representing the conversion factor of electrical energy to natural gas in the z-th P2G plant.
4. The method for optimizing operation and configuration of power, heat and gas conversion considering carbon emission according to claim 3, wherein the carbon emission model is as follows:
Figure FDA0003001536070000051
in the formula: eGRepresents the total pollutant emission, aE，bE，dE，γE，δEIs the emission coefficient of the thermal power generating unit.
5. The method for optimizing operation of an electric, thermal and gas conversion system according to claim 4, wherein the objective function F is:
F＝SW-CE
SW＝RD-cost
RD＝λePD+λheatPheat+λgasPgas
CE＝hEG,
Figure FDA0003001536070000052
cost＝Cg+Cgf+Cgs+CPw+Cpv+CSD+CHOU+CCHP+CP2G
SWmaxis the maximum value of the social benefit,
Figure FDA0003001536070000053
maximum amount of carbon emission, λeIs the cost of the electrical energy consumed, λheatIs the cost of heat energy, λgasIs the cost of consuming gas energy, PD，PheatAnd PgasRespectively representing the demand for electrical energy, the demand for thermal energy and the demand for natural gas.
6. The method of claim 5, wherein the constraints include:
active power balance:
Figure FDA0003001536070000054
in the formula:
Figure FDA0003001536070000055
representing the active power at time t of the f-th gas generator, NGIndicates the number of the gas generators,
Figure FDA0003001536070000056
representing the active power at time t of the vth natural gas supply facility, NgfIndicating the number of natural gas supply facilities,
Figure FDA0003001536070000057
representing the active power at time t of the w-th wind turbine, NwIndicating the number of wind generators, pvq,tRepresenting the active power of the qth photovoltaic generator at time t, NpvIndicates the number of the photovoltaic generators,
Figure FDA0003001536070000061
indicating the active power at time t of the y-th CHP device, NCHPIndicates the number of the CHP devices,
Figure FDA0003001536070000062
representing the active power at time t of the r-th energy storage device, NBIndicates the number of the energy storage devices,
Figure FDA0003001536070000063
represents the active power, N, of the z-th station P2G device at time tP2GDenotes the number of P2G devices, PDtRepresenting the total power demand at time t, PlossRepresents the power loss;
Figure FDA0003001536070000064
and (3) constraint of bus voltage and branch flow:
Figure FDA0003001536070000065
in the formula: vi max，Vi minRepresents the ith bus voltage maximum and minimum, Sflow,iShowing the tidal flow distribution on the ith branch,
Figure FDA0003001536070000066
represents the maximum value of the tidal flow distribution on the ith branch;
and (3) restricting the climbing rate of the wind driven generator:
Figure FDA0003001536070000067
in the formula: UR and DR respectively represent the rising rate constraint and the falling rate constraint of any wind turbine generator set, Pw,tRepresenting the power of the w-th wind driven generator set at the time t;
output power constraint:
Figure FDA0003001536070000068
in the formula:
Figure FDA0003001536070000069
respectively representing the maximum value and the minimum value of the output power of the f-th gas generator,
Figure FDA00030015360700000610
showing the output power of the f-th gas generator at the moment t,
Figure FDA00030015360700000611
are respectively provided withRepresents the maximum value and the minimum value of the output power of the y-th cogeneration generator,
Figure FDA00030015360700000612
represents the maximum value of active power generated by the wind turbine,
Figure FDA0003001536070000071
denotes the maximum value of the output power of the z-th stage P2G device, PwmaxRepresents the maximum power value, pv, of the w-th wind turbinemaxRepresenting the maximum power of the qth photovoltaic generator;
thermal node balance constraint:
Figure FDA0003001536070000072
in the formula:
Figure FDA0003001536070000073
indicating the thermal power, N, at time t of the y-th CHP deviceCHPIndicates the number of the CHP devices,
Figure FDA0003001536070000074
indicating the thermal power, N, of the xth heating unit at time tHOUIndicates the number of the heat supply units,
Figure FDA0003001536070000075
indicating the thermal power, N, of the s-th other device at time tHSNumber of other devices, HDtIndicating the heat demand at time t, HlossRepresents heat loss;
heat restraint of cogeneration and heating units:
Figure FDA0003001536070000076
in the formula:
Figure FDA0003001536070000077
represents the maximum value and the minimum value of the heat production quantity of the y-th cogeneration generator,
Figure FDA0003001536070000078
the maximum value and the minimum value of the heat production quantity of the x-th heat supply unit set are represented;
gas node balance constraint:
Figure FDA0003001536070000079
in the formula:
Figure FDA00030015360700000710
indicating the gas output of the z-th station P2G device at time t,
Figure FDA00030015360700000711
representing the gas output of the vth natural gas plant at time t,
Figure FDA00030015360700000712
indicating the gas consumption of the y-th CHP device at time t,
Figure FDA00030015360700000713
denotes the gas consumption, Q, of the f-th gas generator unit at time tD.tIndicating the gas demand at time t, QlossRepresents the amount of gas lost;
pressure and airflow constraints:
Figure FDA00030015360700000714
in the formula:
Figure FDA0003001536070000081
representing the maximum value and the minimum value of the pressure intensity between the natural gas node m and the natural gas node n;
Figure FDA0003001536070000082
the pressure intensity between a natural gas node m and a natural gas node n at the moment t of the natural gas equipment is represented;
charge and discharge power constraint:
Figure FDA0003001536070000083
in the formula: emaxThe maximum energy storage device capacity of different storage types is represented, and delta t represents time difference;
energy restraint in the energy storage device:
Figure FDA0003001536070000084
in the formula:
Figure FDA0003001536070000085
representing the maximum and minimum values of energy in the r-th energy storage device,
Figure FDA0003001536070000086
representing the initial energy of the r-th energy storage device.
7. A computer readable storage medium storing one or more programs, characterized in that: the one or more programs include instructions which, when executed by a computing device, cause the computing device to perform any of the carbon emissions aware electricity, heat, gas conversion optimization running configuration methods of claims 1 to 6.
8. A computing device, characterized by: comprises the steps of (a) preparing a mixture of a plurality of raw materials,
one or more processors, memory, and one or more programs stored in the memory and configured to be executed by the one or more processors, the one or more programs including instructions for performing any of the carbon emissions aware electricity, heat, gas shift optimization operational configuration methods of claims 1-6.
CN202110348701.4A 2021-03-31 2021-03-31 Method for optimizing operation of electric, thermal and gas conversion system with consideration of carbon emission and storage medium Expired - Fee Related CN113131513B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110348701.4A CN113131513B (en) 2021-03-31 2021-03-31 Method for optimizing operation of electric, thermal and gas conversion system with consideration of carbon emission and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110348701.4A CN113131513B (en) 2021-03-31 2021-03-31 Method for optimizing operation of electric, thermal and gas conversion system with consideration of carbon emission and storage medium

Publications (2)

Publication Number Publication Date
CN113131513A true CN113131513A (en) 2021-07-16
CN113131513B CN113131513B (en) 2022-08-02

Family

ID=76775055

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110348701.4A Expired - Fee Related CN113131513B (en) 2021-03-31 2021-03-31 Method for optimizing operation of electric, thermal and gas conversion system with consideration of carbon emission and storage medium

Country Status (1)

Country Link
CN (1) CN113131513B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116681545A (en) * 2023-01-30 2023-09-01 兰州理工大学 Facility agriculture park near-zero carbon implementation method considering biomass-P2G coupling

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110322051A (en) * 2019-06-06 2019-10-11 国网浙江省电力有限公司经济技术研究院 Consider the integrated energy system Optimal Configuration Method of N-1 security constraint
CN112332412A (en) * 2020-10-26 2021-02-05 三峡大学 Optimized dispatching method for electricity-gas interconnected comprehensive energy system with P2G wind-abandoning start-stop function
CN112488525A (en) * 2020-12-01 2021-03-12 燕山大学 Electric heating rolling scheduling method and system considering source-charge side response under carbon transaction mechanism

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110322051A (en) * 2019-06-06 2019-10-11 国网浙江省电力有限公司经济技术研究院 Consider the integrated energy system Optimal Configuration Method of N-1 security constraint
CN112332412A (en) * 2020-10-26 2021-02-05 三峡大学 Optimized dispatching method for electricity-gas interconnected comprehensive energy system with P2G wind-abandoning start-stop function
CN112488525A (en) * 2020-12-01 2021-03-12 燕山大学 Electric heating rolling scheduling method and system considering source-charge side response under carbon transaction mechanism

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116681545A (en) * 2023-01-30 2023-09-01 兰州理工大学 Facility agriculture park near-zero carbon implementation method considering biomass-P2G coupling

Also Published As

Publication number Publication date
CN113131513B (en) 2022-08-02

Similar Documents

Publication Publication Date Title
CN109327042B (en) Multi-energy joint optimization scheduling method for micro-grid
CN108229025B (en) Economic optimization scheduling method for cooling, heating and power combined supply type multi-microgrid active power distribution system
CN110866627B (en) Multi-zone electricity-gas coupling comprehensive energy system optimal scheduling method considering step gas price
CN109919399B (en) Day-ahead economic dispatching method and system for comprehensive energy system
CN112035984B (en) Collaborative planning method for comprehensive energy system of electricity-gas-storage area
CN110391655B (en) Multi-energy-coupling micro-energy-network economic optimization scheduling method and device
CN110163415B (en) Multi-objective fuzzy cooperative optimization method for multi-energy flow system under variable working condition characteristic
Zhai et al. Optimization of integrated energy system considering photovoltaic uncertainty and multi-energy network
CN113131513B (en) Method for optimizing operation of electric, thermal and gas conversion system with consideration of carbon emission and storage medium
Li et al. Planning model of integrated energy system considering P2G and energy storage
CN113128868A (en) Regional comprehensive energy system scheduling optimization method and device
CN112836957A (en) Regional comprehensive energy system planning method considering power supply reliability
Liu et al. Low carbon economic dispatch for multi-energy distribution network with compressed air energy storage system as energy hub
CN111724026A (en) Optimization method for coupling operation of multi-energy network and water distribution network
Liu et al. Influence Evaluation of Integrated Energy System on the Unit Commitment in Power System
Meng et al. Economic optimization operation approach of integrated energy system considering wind power consumption and flexible load regulation
CN114139837A (en) Regional multi-system double-layer distributed optimization scheduling method considering double-layer carbon emission optimization distribution model
Chenrui et al. Optimal Configuration of Energy Storage for Integrated Energy Stations Considering Multiple Scenarios
Liu et al. The Research and Method on Planning Model of Integrated Energy Internet
CN112949903A (en) Day-ahead optimization scheduling method for regional comprehensive energy system with power-to-gas and energy storage functions
CN111668882A (en) Method and device for optimizing output of micro power supply in intelligent energy ring network
Sun et al. Research on Equilibrium Carbon Price Method-based Carbon Trade Pricing of Electric-Natural Gas Multi-Energy System
CN113393077B (en) Method for configuring an electric-gas multi-energy storage system taking into account the uncertainty of the energy used by the user
Bayasgalan et al. The Utilization of an Energy Storage Facility in Dispatching Regulation of Power System
Sui et al. Optimization of Monthly Power Generation Plan for Thermal Power Units Considering Access of Large-scale New Energy

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220802

CF01 Termination of patent right due to non-payment of annual fee