CN114037292A - Low-carbon optimized scheduling method of electricity-gas integrated energy system considering carbon capture - Google Patents

Low-carbon optimized scheduling method of electricity-gas integrated energy system considering carbon capture Download PDF

Info

Publication number
CN114037292A
CN114037292A CN202111335666.9A CN202111335666A CN114037292A CN 114037292 A CN114037292 A CN 114037292A CN 202111335666 A CN202111335666 A CN 202111335666A CN 114037292 A CN114037292 A CN 114037292A
Authority
CN
China
Prior art keywords
carbon
gas
unit
kcf
power
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202111335666.9A
Other languages
Chinese (zh)
Inventor
余莎
何光层
李新洪
彭文英
聂剑锋
叶志祥
李树东
杨竣皓
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Baoshan Power Supply Bureau of Yunnan Power Grid Co Ltd
Original Assignee
Baoshan Power Supply Bureau of Yunnan Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Baoshan Power Supply Bureau of Yunnan Power Grid Co Ltd filed Critical Baoshan Power Supply Bureau of Yunnan Power Grid Co Ltd
Priority to CN202111335666.9A priority Critical patent/CN114037292A/en
Publication of CN114037292A publication Critical patent/CN114037292A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06313Resource planning in a project environment
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06Electricity, gas or water supply
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Abstract

The invention relates to a low-carbon optimal scheduling method for an electricity-gas integrated energy system considering carbon capture, and belongs to the technical field of optimal scheduling of integrated energy systems. Firstly, constructing a carbon-containing trapped electricity-gas comprehensive energy system; then constructing a low-carbon economic dispatching model of the carbon-containing trapped electricity-gas integrated energy system based on the system; the carbon-containing trapped electricity-gas integrated energy system low-carbon economic dispatching model is based on system operation cost and CO2The minimum comprehensive cost of the sum of the related cost and the wind abandon punishment cost is an objective function, and the constraint conditions are electric power network constraint, natural gas network constraint, electric-to-gas equipment constraint, carbon capture system constraint, carbon storage and hydrogen storage equipment constraint. The invention not only can realize the bidirectional electrical coupling of the power network and the natural gas networkThe wind power is consumed to the maximum extent, carbon circulation can be formed through the carbon collecting device, the carbon emission reduction capability of the system is improved, and the system is easy to popularize and apply.

Description

Low-carbon optimized scheduling method of electricity-gas integrated energy system considering carbon capture
Technical Field
The invention belongs to the technical field of optimization scheduling of an integrated energy system, and particularly relates to a carbon capture considered low-carbon optimization scheduling method for an electricity-gas integrated energy system.
Background
Under the background that the global warming problem is becoming more serious, a series of green economic concepts such as 'carbon neutralization', 'carbon peak-reaching action' and the like are complied with and proposed. The electricity-gas comprehensive energy system plays an important role in green transformation of a future energy system due to the characteristics of supply and demand interaction, high comprehensive energy utilization rate and high operation efficiency. However, the traditional electric comprehensive energy system only couples the power network and the natural gas network through the electric-to-gas equipment singly, and still has the problem of high carbon emission, so that the research on how to operate the comprehensive energy system in a low-carbon and economic manner is of great significance to green transformation of the power industry.
The power generation of a thermal power generating unit at the side of a power grid in the comprehensive energy system is the largest carbon emission source, and the comprehensive energy system can adopt an effective method to increase the grid access rate of green energy such as wind power, photoelectricity and the like, change the production mode and develop a carbon capture technology. For the two methods, the following disadvantages exist: firstly, the wind power output is influenced by natural factors such as seasons, climate and the like, uncertainty exists, and as the scale of wind power grid connection is continuously increased, electricity is generatedThe stability and the safety of the operation of the power system are also influenced, and the popularization and the utilization of wind power in the comprehensive energy system are restricted by the problems of wind power access and consumption. And secondly, a carbon capture system is introduced into the traditional electricity-gas comprehensive energy system, so that the problem of large carbon emission is solved to a certain extent, but at present, the flexible operation between the carbon capture system and the electricity-gas conversion equipment and the mechanism for jointly consuming wind power by the carbon capture system and the electricity-gas conversion equipment are rarely researched after the carbon capture system is introduced. The coupling of a pure carbon capture system with an electric gas conversion device has certain disadvantages. Firstly, the methanation reaction of the electric-to-gas equipment in the electric comprehensive energy system consumes CO by using the residual electric energy only in the period of wind abandonment2Synthesis of methane, but carbon capture system captures CO2Possibly throughout the scheduling period, CO2The imbalance of gas flow supply and demand causes unnecessary cost and waste of gas resources. Secondly, the application of the hydrogen is wide, and compared with the direct participation in methanation reaction, the hydrogen is used for the electric automobile market, a fuel cell and a hydrogen gas turbine to generate electricity or possibly bring better economic benefit.
Therefore, how to overcome the defects of the prior art is a problem to be solved urgently in the technical field of the optimization scheduling of the comprehensive energy system at present.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provides a low-carbon optimal scheduling method of an electricity-gas integrated energy system considering carbon capture.
In order to achieve the purpose, the technical scheme adopted by the invention is as follows:
a low-carbon optimized dispatching method for an electricity-gas integrated energy system considering carbon capture adopts the electricity-gas integrated energy system with carbon capture;
the carbon-containing trapped electrical comprehensive energy system comprises a conventional thermal power generating unit, a carbon trapping unit, carbon storage equipment, methanation equipment, hydrogen storage equipment, a gas turbine and electrolysis equipment; the carbon capture unit is a conventional thermal power unit additionally provided with a carbon capture system;
the conventional thermal power generating unit is connected with a carbon capture system of the carbon capture unit;
the carbon capture unit is respectively connected with the carbon storage equipment and the methanation equipment;
the methanation equipment is also respectively connected with the carbon storage equipment and the hydrogen storage equipment;
the hydrogen storage equipment is also respectively connected with the gas turbine and the electrolysis equipment;
the scheduling method comprises the following steps:
constructing a low-carbon economic dispatching model of the carbon-containing trapped electrical comprehensive energy system; the carbon-containing captured electrical comprehensive energy system low-carbon economic dispatching model is based on system operation cost and CO2The minimum comprehensive cost of the sum of the related cost and the wind abandon punishment cost is a target function, and the constraint conditions are electric power network constraint, natural gas network constraint, electric-to-gas equipment constraint, carbon capture system constraint, carbon storage and hydrogen storage equipment constraint; said CO2Associated costs include CO emissions from fossil fuel units2Carbon tax cost FCT、CO2Transmission storage cost FTsCarbon trapping unit for trapping CO from air2Cost F ofCC(ii) a Then, solving is carried out, and scheduling is carried out according to the solving result.
Further, it is preferable that the objective function is:
min F=FO+FC+FW (1)
in the formula, Fo$ for system operating cost, $; fcIs CO2Associated cost, $; fwPunishing cost for wind abandonment;
(1) the system operation cost comprises the operation cost F of the conventional thermal power generating unitGGas production cost F of natural gas sourceAnd natural gas storage
Gas cost F for gas facilityNS,$;
Figure BDA0003350397790000021
In the formula, NG、Nω、NGsThe number of the conventional coal burners, the number of natural gas sources and the number of natural gas storage devices are respectively; n is a radical ofTIs a scheduling period; beta is ai、βω、βsThe unit is the fuel price of a conventional coal-fired unit i, the unit gas production cost of a natural gas source omega and the unit gas consumption cost of a natural gas storage device s, and the units are $/MBtu, $/kcf and $/kcf; a isGi、bGi、cGiThe cost coefficient of the ith conventional thermal power generating unit; pGi,tThe output power, MW, of the conventional thermal power generating unit i at the moment t; qω,tThe gas production amount of the natural gas source omega at the time t is kcf:
Figure BDA0003350397790000022
dividing the natural gas storage device s into the natural gas injection and extraction amount at the time t, kcf;
(2)CO2associated costs include CO emissions from fossil fuel units2Carbon tax cost FCT,$、CO2Transmission storage cost FTSAnd the carbon capture system captures CO from the air2Cost F ofCC,$;
Figure BDA0003350397790000023
In the formula, NFU、NCC、NCSRespectively representing the number of fossil fuel units, the number of carbon gathering units and the number of carbon storage equipment; beta is aCT、βTS、βCCRespectively carbon tax price and unit CO2Transportation and storage price of, capture of unit CO from air2The required price cost in $/kcf; u. ofiCO for conventional coal-fired unit i2The discharge intensity; pi,tThe output power, MW, of the fossil fuel unit i at the moment t;
Figure BDA0003350397790000024
for CO trapped by carbon trapping machine set at time t2Amount captured, kfc;
Figure BDA0003350397790000025
for supplying carbon from airs carbon dioxide flow, kfc/h;
Figure BDA0003350397790000026
the carbon dioxide flow rate for supplying the carbon capture unit with the carbon storage device s is kfc/h; the fossil fuel unit comprises a conventional thermal power unit and a gas unit;
(3) wind abandon penalty cost:
Figure BDA0003350397790000027
in the formula, NWThe number of wind power plants; cw is a wind abandon punishment cost coefficient, $ MWh; pw,tThe actual output value, MW, of the wind power plant in the time period t; and P' w and t are output predicted values, MW, of the wind power plant in the t period.
Further, preferably, the power network constraints include a power balance constraint, a rotation standby constraint, a unit operation constraint and a power grid operation safety and reliability constraint; the method specifically comprises the following steps:
1) power balance constraint
Figure BDA0003350397790000031
In the formula, Ab(.)Representing an association matrix of a fossil fuel unit, a wind power plant, an electric-to-gas device and a power transmission line associated with a power grid node b; n is a radical ofFU、NW、NPtG、NTLThe total number of fossil fuel units, wind power plants, power-to-gas equipment and power transmission lines is respectively;
Figure BDA0003350397790000032
power, MW, to the input electrical to gas device; pED,b,tRepresents the electrical load at node b at time t, MW; PLl,tThe transmission power of the transmission line l at the time t, MW;
2) rotational back-up restraint
Figure BDA0003350397790000033
In the formula, τi,tThe variable is 0-1 and represents the start-stop state of the fossil fuel unit i at the moment t, the value is 1 when the unit i starts to operate, and the value is 0 when the unit i stops operating; pi,min、Pi,maxRespectively representing the lower limit, the upper limit and the MW of the output of the fossil fuel unit i;
Figure BDA0003350397790000034
respectively representing the up-regulation reserve capacity demand value and the down-regulation reserve capacity demand value MW which are at least required to be met by the power system at the moment t;
3) the unit operation constraints comprise wind power output constraints, fossil fuel unit output constraints, climbing rate constraint minimum startup and shutdown time constraints and startup and shutdown energy consumption constraints, and specifically comprise the following steps:
wind power output constraint:
0≤Pw,t≤P′w,t (7)
output constraint of fossil fuel units:
τi,tPi,min≤Pi,t≤τi,tPi,max (8)
third, climbing speed constraint:
Figure BDA0003350397790000035
in the formula, Sup,i、Sdn,iRespectively representing the upward and downward climbing rates, MW, of the fossil fuel unit i; pi,t+1The output power, MW, of the fossil fuel unit i at the moment t + 1; tau isi,t+1The variable is 0-1 and represents the start-stop state of the fossil fuel unit i at the moment of t + 1;
and fourthly, constraint of minimum start-up and shut-down time:
Figure BDA0003350397790000041
in the formula (I), the compound is shown in the specification,
Figure BDA0003350397790000042
respectively representing the minimum time interval number h of starting and stopping the fossil fuel unit i;
Figure BDA0003350397790000043
respectively indicating the number h of the periods when the fossil fuel unit i is continuously started and stopped in the period t-1; tau isi,t-1The variable is 0-1 and represents the start-stop state of the fossil fuel unit i at the time t-1;
energy consumption constraint of starting and stopping
Figure BDA0003350397790000044
In the formula, sui、sdiRespectively representing the starting and stopping energy consumption, MBtu, of the fossil fuel unit i; SUi,t、SDi,tRespectively representing the starting-up and stopping energy consumption, MBtu, of the fossil fuel unit i at the time t;
4) power grid operation safety and reliability constraint
PLl,t=(πl1,tl2,t)/xl (12)
In the formula, pil1,t、πl2,tThe voltage phase angle values, rad, of the head end and the tail end of the power transmission line l at the time t are respectively; x is the number oflIs the reactance value of line l, p.u.;
the line transmission power constraint and the node voltage angle constraint are as follows:
Figure BDA0003350397790000045
in the formula, PLl,min、PLl,maxRespectively the minimum and maximum active power, MW, transmitted by the line l; piref、πb,min、πb,max、πb,tRespectively representing the voltage phase angle value of the reference node and the minimum and maximum voltage phase angle values of the node b, the actual voltage phase at the time tAngle value, rad.
Further, preferably, the natural gas network constraints include gas source gas production constraints, natural gas pipeline flow constraints, natural gas storage equipment constraints, node pressure constraints and node supply and demand balance constraints;
1) gas source gas production constraint
Qω,min≤Qω,t≤Qω,max (14)
In the formula, Qω,min、Qω,maxThe unit of the lower limit value and the upper limit value of the omega gas production of the natural gas source is kcf;
2) natural gas pipeline flow restriction
For a natural gas pipeline mn, gas flows from m node to n node, and the pipeline gas flow is expressed as:
Figure BDA0003350397790000051
in the formula, Qmn,tThe gas flow rate of the natural gas pipeline mn flowing at the time t is kcf/h; sgn (p)m,t,pn,t) The direction of gas flow is shown, when the direction is 1, the gas flows from the m node to the n node, and when the direction is-1, the gas flows from the n node to the m node; pm,t、pn,tPressure values at nodes m and n at the time t, Psig; cmnIs the pipeline constant;
for a natural gas pipeline mn with a compressor j, the compression ratio between the output node and the input node is constrained by:
Figure BDA0003350397790000052
in the formula, Rj,max、Rj,minMaximum and minimum values of the allowable compressor compression ratio;
3) natural gas storage equipment restraint
The gas storage amount of the gas storage equipment at any time t in the operation process
Figure BDA0003350397790000053
Expressed as:
Figure BDA0003350397790000054
in the formula (I), the compound is shown in the specification,
Figure BDA0003350397790000055
indicating the gas storage amount of the gas storage device s at the time t and t-1, respectively, kcf;
Figure BDA0003350397790000056
dividing the natural gas storage device s into the natural gas injection and extraction amount at t moment, kcf/h; Δ t is the time step;
Figure BDA0003350397790000057
respectively representing the minimum and maximum gas storage capacities of the gas storage device s, kcf;
Figure BDA0003350397790000058
respectively representing the maximum gas injection and extraction flow rates of natural gas storage equipment s, kcf/h;
Figure BDA0003350397790000059
respectively representing the gas capacity of the gas storage device s at the initial and ending moments of a scheduling period, kcf;
4) nodal pressure constraint
For the natural gas node m, the pressure magnitude should satisfy the following constraint:
pm,min≤pm,t≤pm,max (18)
in the formula, pm,min、pm,maxRespectively representing the lower limit and the upper limit of the pressure value of the node m;
5) node supply and demand balance constraints
For natural gas node m, the equilibrium equation is as follows:
Figure BDA00033503977900000510
in the formula, Nω、NGs、NPtG、NNU、NGPRespectively indicating the quantity of the gas source, the gas storage facility, the electric gas conversion equipment, the gas unit and the gas transmission pipeline; b ism(.)Representing a correlation matrix associated with a natural gas node m;
Figure BDA0003350397790000061
the amount of the methane gas synthesized by PtG equipment is kcf/h; qNi,tThe natural gas consumption of the gas unit i at the time t is kcf/h; qGL,m,tThe air load at the time t at the node m is kcf/h; qmn,tThe natural gas flow rate flowing through the pipeline mn in the period t is kcf/h.
Further, preferably, the electric transfer device is constrained by:
the hydrogen quantity generated by the electrolytic hydrogen production reaction of the electric gas conversion equipment
Figure BDA0003350397790000062
Expressed as:
Figure BDA0003350397790000063
in the formula (I), the compound is shown in the specification,
Figure BDA0003350397790000064
the unit of the hydrogen quantity generated by the electrolytic hydrogen production reaction of the electrolysis equipment is kcf/h;
Figure BDA0003350397790000065
power, MW, for input to the electrical to gas equipment;
Figure BDA0003350397790000066
hydrogen production efficiency in the electrolysis reaction; xie,gMBtu/MWh, which is the equivalent thermal energy coefficient;
Figure BDA0003350397790000067
high heating value of hydrogen, MBtu/kcf;
The upper and lower limits of the power input into the electric gas conversion equipment in the electrolytic reaction are respectively 20% and 5% of the designed capacity, so that the input power of the electric gas conversion equipment can meet the following constraint:
Figure BDA0003350397790000068
in the formula, τptgThe variable is 0-1, which represents the running state of the electric gas conversion device, and the value is 1 when the electric gas conversion device runs and 0 when the electric gas conversion device stops running;
Figure BDA0003350397790000069
the lower limit value, the upper limit value and the MW of the allowable input power of the electric gas conversion device are respectively set;
CO consumed by methanation2The amounts and the amount of methane gas produced are expressed as:
Figure BDA00033503977900000610
in the formula (I), the compound is shown in the specification,
Figure BDA00033503977900000611
CO consumed for methanation reactions2Amount, kcf/h;
Figure BDA00033503977900000612
the amount of methane gas synthesized by methanation reaction is kcf/h;
Figure BDA00033503977900000613
the hydrogen flow taken off for the hydrogen storage device s for the synthesis of methane, kcf/h;
Figure BDA00033503977900000614
the reaction coefficients between hydrogen and carbon dioxide, hydrogen and methane;
meanwhile, the heat energy generated by methanation reaction is as follows:
Figure BDA00033503977900000615
in the formula (I), the compound is shown in the specification,
Figure BDA00033503977900000616
the recovered thermal power, MW, for the electric gas-converting equipment; ΨheatMWh/kcf for heat recovery coefficient;
further, it is preferred that the carbon capture system is constrained by:
the carbon capture unit operation model is represented as follows:
Figure BDA0003350397790000071
in the formula, PcEquivalent output power, MW; pc netNet output power, MW; pc conFor capturing CO2Required energy consumption, MW; maintenance energy consumption Pc mIs a fixed value, MW; pc oEnergy consumption for operation, MW;
Figure BDA0003350397790000072
CO that can be processed for carbon capture systems2Amount, kcf/h; gamma raycFor the carbon capture unit c to capture unit CO2The energy consumption generated is kfc/MWh;
the carbon capture unit generates equivalent output power and simultaneously generates carbon dioxide Q according to a certain proportionc
Figure BDA0003350397790000073
In the formula, mucCO for carbon capture train c2Emission intensity, kfc/MWh;
Figure BDA0003350397790000074
CO for carbon capture train c2Trapping capacity, kcf/h; beta is acCO for carbon capture train c2The capture rate;
in summary, the net output power of the carbon capture unit is expressed as:
Pc net=Pc-Pc m-Pc o (26)
and the net output power meets the output limit constraint:
Figure BDA0003350397790000075
in the formula (I), the compound is shown in the specification,
Figure BDA0003350397790000076
respectively the minimum and maximum net power values allowed to be output by the carbon capture unit;
further, preferably, the constraints of the carbon storage and hydrogen storage device are as follows:
1) hydrogen storage apparatus
Figure BDA0003350397790000077
In the formula (I), the compound is shown in the specification,
Figure BDA0003350397790000078
the hydrogen flow rate injected for the hydrogen storage device s, kcf/h;
Figure BDA0003350397790000079
the hydrogen storage amounts of the hydrogen storage equipment s at the times t and t-1, respectively, kcf;
Figure BDA00033503977900000710
dividing into kcf/h hydrogen injection and extraction quantity of hydrogen storage equipment s at t moment;
Figure BDA00033503977900000711
respectively representing the minimum and maximum hydrogen storage capacities of the hydrogen storage device s, kcf;
Figure BDA00033503977900000712
respectively representing the maximum hydrogen injection and extraction flow rates of the hydrogen storage device s, kcf/h; n is a radical ofTIs the end time of a scheduling period, h;
Figure BDA00033503977900000713
respectively representing the hydrogen capacities of the hydrogen storage equipment at the initial and end moments of a scheduling period, wherein the hydrogen capacities of the hydrogen storage equipment and the hydrogen capacity are equal, kcf;
the hydrogen required by the methanation equipment and the gas turbine is from the hydrogen storage equipment, so the hydrogen flow rate should satisfy the following constraints:
Figure BDA0003350397790000081
in the formula
Figure BDA0003350397790000082
Representing the hydrogen flow rate of the hydrogen storage device to the gas turbine, kcf/h;
2) carbon storage equipment
Carbon storage amount of carbon storage equipment at any moment in operation process
Figure BDA0003350397790000083
Expressed as:
Figure BDA0003350397790000084
in the formula (I), the compound is shown in the specification,
Figure BDA0003350397790000085
respectively representing CO of the carbon storage equipment at t and t-12Capacity, kfc;
Figure BDA0003350397790000086
respectively representing CO of carbon storage equipment s at time t2The injection and extraction flow rates are kcf/h;
corresponding carbon storage facility capacity constraint, CO2The injection quantity and the extraction quantity of the flow are restricted in one scheduling periodInitial, end CO2The capacity constraints are as follows:
Figure BDA0003350397790000087
in the formula (I), the compound is shown in the specification,
Figure BDA0003350397790000088
respectively representing minimum and maximum CO of carbon storage facility s2Storage capacity, kfc;
Figure BDA0003350397790000089
CO of carbon storage equipment s at t moment2The injection and extraction amount is kcf/h;
Figure BDA00033503977900000810
respectively representing CO of carbon storage units s2The maximum injection and extraction flow rate is kcf/h;
Figure BDA00033503977900000811
respectively representing CO of carbon storage equipment at the initial and end time of a scheduling period2Capacity, kcf:
CO trapped by carbon trapping system in carbon removal trapping unit2CO of carbon storage equipment2The injected amount also comes from CO captured in the air2Expressed as:
Figure BDA00033503977900000812
in the formula (I), the compound is shown in the specification,
Figure BDA00033503977900000813
the flow rate of carbon dioxide supplied from the carbon capture unit to the carbon storage facility s at time t is shown at kcf.
The invention also provides the carbon-containing trapping electricity-gas comprehensive energy system.
The invention also provides an electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor executes the program to implement the steps of the low-carbon optimal scheduling method for the electricity-gas integrated energy system considering carbon capture as described above.
The invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the steps of the method for low-carbon optimal scheduling of an electric-gas integrated energy system in consideration of carbon capture as described above.
In the present invention, aGi、bGi、cGiThe cost coefficient of the ith conventional thermal power generating unit can be found from an electric power regulation manual, a network or a factory nameplate of an actual generator manufacturer.
In the present invention, the fossil fuel power plant includes a conventional thermal power plant NGAnd gas turbine unit NNUPower calculation system-P involving bothiAnd (4) showing.
In the present invention, CmnIs a pipe constant, the size depending on diameter, length, coefficient of friction, gas composition, etc., kcf/Psig; the specific calculation is only needed by the calculation of the prior art content and the related formula. And delta t is a time step, and is preferably 1 h. N is a radical ofTThe value is 24h for one scheduling period.
Figure BDA0003350397790000091
The preferred value is 0.73 for the hydrogen production efficiency in the electrolysis reaction; zetae,gThe value is 3.41, MBtu/MWh;
Figure BDA0003350397790000092
the high heat value of the hydrogen is preferably 0.342;
Figure BDA0003350397790000093
the reaction coefficients of hydrogen and carbon dioxide and hydrogen and methane are all 0.25; ΨheatThe MWh/kcf value is 0.0126 for the heat recovery coefficient. Maintenance energy consumption Pc mIs a fixed value, is the lowest energy consumption for maintaining the operation of the carbon capture unit, andthe running state of the carbon capture unit is irrelevant, and the value of the actual running condition can be generally 10-20% of the rated power of the carbon capture unit. Pc oThe value of MW, which is the operation energy consumption, changes along with the change of the operation state of the carbon capture unit and is positively correlated with the carbon dioxide capture capacity of the carbon capture unit. Beta is acCO for carbon capture train c2The trapping rate, which is 0.9:
the invention discloses a low-carbon optimal scheduling method of an electricity-gas integrated energy system considering carbon capture, aiming at the problem of wind abandon caused by new energy access in the electricity-gas integrated energy system and the problem of carbon emission of a thermal power generating unit. Firstly, on the basis of a traditional electricity-gas integrated energy system, a carbon capture system is additionally arranged on a conventional thermal power generating unit to form a carbon capture unit, and CO is captured2And the carbon emission of the unit is reduced by trapping.
And carbon dioxide (CO) capture based on carbon capture system2) While methanation reactions in the electric gas-conversion plant (P2G) consume CO2The system forms carbon cycle in the scheduling period by taking the joint operation mode of the carbon storage and hydrogen storage equipment into consideration. Reduction of P2G facility for purchasing CO from methanation reactions2And increases the net output power of the carbon capture plant. Secondly, in this combined mode of operation, the cost of operation of the system, CO2The minimum comprehensive cost of the sum of the related cost and the wind abandon punishment cost is an objective function, an environment economic dispatching model of the carbon-containing capture electric comprehensive energy system is constructed, and the constraints of an electric power network, a natural gas network, an electric gas conversion device, a carbon capture system and carbon and hydrogen storage devices are considered. Finally, the model is compared and verified under different simulation scenes, and the result shows that: compared with the traditional comprehensive energy system architecture, the system of the invention not only can realize the bidirectional electrical coupling of the power network and the natural gas network by introducing the carbon capture unit, the carbon storage equipment and the hydrogen storage equipment 3, and furthest consume the wind power, but also can form carbon circulation through the carbon capture device, improve the carbon emission reduction capability of the system, ensure that the system has both economy and environmental protection, and has important theoretical performance under the background of carbon neutralizationValue and practical significance.
Preferably, the model is solved by using MATLAB to call GUROBI.
According to the invention, the carbon storage equipment is additionally arranged between the carbon capture system and the methanation equipment, so that CO required by the methanation reaction2The gas quantity is provided by the carbon storage equipment, the transmission of the gas flow in the carbon circulation process is balanced, the resource waste is reduced, and the flexible operation capacity of the carbon capture power plant is improved.
According to the invention, the hydrogen storage equipment is arranged between the electrolysis equipment and the methanation equipment, so that the two processes of electricity-to-natural gas and electricity-to-hydrogen are decoupled. The electrolytic hydrogen production by wind power can be absorbed to the maximum extent through the electrolytic reaction, and the hydrogen extraction amount is supplied for methanation reaction, and the hydrogen can also be used for power generation of a gas turbine, so that better economic benefit is created.
Under the combined operation mode, the wind power can be basically fully consumed in the scheduling period, and the optimal economic scheduling in the whole day is realized. Meanwhile, the comprehensive operation cost and the carbon emission of the system are optimal, and the economical efficiency and the environmental protection performance of the system are greatly improved.
Compared with the prior art, the invention has the beneficial effects that:
1. the invention is improved on the basis of the traditional electricity-gas integrated energy system, introduces the carbon capture system, provides a combined operation mode of carbon storage and hydrogen storage equipment between the carbon capture and electricity-gas conversion equipment, realizes the carbon circulation of the electricity-gas integrated energy system in a scheduling period, and reduces the purchase of CO by the electricity-gas conversion equipment due to methanation reaction2The cost of the method increases the net output power of the power plant of the carbon-containing capture system by nearly 10-15%.
2. In a combined operation mode, the flexible operation capacity of the carbon capture power plant and the electricity-to-gas equipment is remarkably improved, the full consumption of wind power can be basically realized in a scheduling period, the penalty cost of wind abandoning is basically 0, and the optimal economic scheduling in the whole day is realized.
3. In combined operation mode, the carbon capture unit captures CO2The carbon emission of the electricity-gas integrated energy system is reduced, and compared with the carbon capture system which is not additionally arranged, the carbon emission reduction of the power plant in a combined operation modeThe height is about 31.2%; compared with the method that the carbon capture system is additionally arranged but the combined operation mode is not considered, the carbon emission reduction of the power plant is improved by nearly 7.5 percent, and the environmental cost of the system is greatly reduced.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
FIG. 1 is a block diagram of a 6-node power system and a 7-node natural gas system;
FIG. 2 is a graph of predicted values of power load, natural gas load, and wind power;
FIG. 3 is a wind power consumption comparison diagram under three scenes;
FIG. 4 is a comparison graph of the output of the carbon capture unit 5 under three scenes;
FIG. 5 is a diagram of the carbon capture and recovery heat energy of the unit 5 in scenario 3;
FIG. 6 shows CO under three scenarios2Emission and CO2A collection amount graph;
FIG. 7 is a schematic diagram of the configuration of an electric-gas integrated energy system for carbon capture according to the present invention;
FIG. 8 is a schematic view of the operation of the integrated electric-gas energy system for carbon capture according to the present invention;
fig. 9 is a schematic structural diagram of an electronic device according to the present invention.
Detailed Description
The present invention will be described in further detail with reference to examples.
It will be appreciated by those skilled in the art that the following examples are illustrative of the invention only and should not be taken as limiting the scope of the invention. The examples are given without reference to the specific techniques or conditions, according to the techniques or conditions described in the literature in the field or according to the product description. The materials or equipment used are not indicated by manufacturers, and all are conventional products available by purchase.
A low-carbon optimal scheduling method of an electricity-gas integrated energy system considering carbon capture comprises the following steps:
one, construct the electricity that contains the carbon entrapment-gas comprehensive energy system
The system comprises an electric power network, a natural gas network, a carbon capture system and an electric gas conversion device. Specifically, as shown in fig. 7 and 8, the carbon-containing capture electricity-gas integrated energy system comprises a conventional thermal power generating unit, a carbon capture unit, a carbon storage device, a methanation device, a hydrogen storage device, a gas turbine and an electrolysis device; the carbon capture unit is a conventional thermal power unit additionally provided with a carbon capture system;
the conventional thermal power generating unit is connected with a carbon capture system of the carbon capture unit;
the carbon capture unit is respectively connected with the carbon storage equipment and the methanation equipment;
the methanation equipment is also respectively connected with the carbon storage equipment and the hydrogen storage equipment;
the hydrogen storage equipment is also respectively connected with the gas turbine and the electrolysis equipment;
secondly, a combined operation mode of the carbon storage and hydrogen storage equipment is provided
A carbon storage device is additionally arranged between the carbon capture system and the methanation device, so that CO required by the methanation reaction2The gas quantity is provided by the carbon storage equipment, the transmission of the gas flow in the carbon circulation process is balanced, the resource waste is reduced, and the flexible operation capacity of the carbon capture power plant is improved. The hydrogen storage equipment is arranged between the electrolysis equipment and the methanation equipment, so that the decoupling of the two processes of electricity-to-natural gas and electricity-to-hydrogen is realized. The electrolytic hydrogen production by wind power can be absorbed to the maximum extent through the electrolytic reaction, and the hydrogen extraction quantity is supplied to methanation reaction, and the hydrogen can also be used for power generation of a gas turbine, so that better economic benefit is created.
Thirdly, constructing a low-carbon economic dispatching model of the carbon-trapping electricity-gas integrated energy system
The model uses the system running cost, CO2The minimum comprehensive cost of the sum of the related cost and the wind curtailment cost is an objective function, and the electric power network constraint, the natural gas network constraint and the coupling are consideredSynthesizing unit constraints and considering carbon storage and carbon capture system constraints of the hydrogen storage device. The method comprises the following specific steps:
1. objective function
The carbon capture-containing electrical integrated energy system eco-economic dispatch is optimized with minimum integrated cost in a combined mode of operation.
min F=FO+FC+FW (1)
In the formula, FO$ for system operating cost, $; fCIs CO2Associated cost, $; fWPunishment cost for wind abandon, $.
(1) The system operation cost comprises the operation cost F of the conventional thermal power generating unitGAnd the sound and gas cost F of the natural gas sourceGas cost F of natural gas storage facilityNS,$:
Figure BDA0003350397790000111
In the formula, NG、Nω、NGsThe number of the conventional coal burners, the number of natural gas sources and the number of natural gas storage devices are respectively; n is a radical ofTTaking the value of 24h for one scheduling period; beta is ai、βω、βsThe unit is the fuel price of a conventional coal-fired unit i, the unit gas production cost of a natural gas source omega and the unit gas consumption cost of a natural gas storage device s, and the units are $/MBtu, $/kcf and $/kcf; a isGi、bGi、cGiThe cost coefficient of the ith conventional thermal power generating unit; pGi,tThe output power, MW, of the conventional thermal power generating unit i at the moment t; qω,tBeta gas quantity at time t for natural gas source omega, kcf;
Figure BDA0003350397790000112
dividing the natural gas storage device s into the natural gas injection and extraction amount at the time t, kcf;
(2)CO2associated costs include CO emissions from fossil fuel units2Carbon tax cost FCT,$、CO2Is transmitted and storedCost FTSAnd the carbon capture system captures CO from the air2Cost F ofCC,$;
Figure BDA0003350397790000113
In the formula, NFU、NCC、NCSRespectively representing the number of fossil fuel units, the number of carbon gathering units and the number of carbon storage equipment; beta is aCT、βTS、βCCRespectively carbon tax price and unit CO2Transportation and storage price of, capture of unit CO from air2The required price cost, in $/kcf: u. ofiCO for conventional coal-fired unit i2The discharge intensity; pi,tThe output power, MW, of the fossil fuel unit i at the moment t;
Figure BDA0003350397790000114
for CO trapped by carbon trapping machine set at time t2Amount captured, kfc;
Figure BDA0003350397790000115
the flow rate of carbon dioxide supplied from the air to the carbon storage facility s, kfc/h;
Figure BDA0003350397790000116
the carbon dioxide flow rate for supplying the carbon capture unit with the carbon storage device s is kfc/h; the fossil fuel unit comprises a conventional thermal power unit and a gas unit;
(3) wind abandon penalty cost:
Figure BDA0003350397790000121
in the formula, NWThe number of wind power plants; cw is a wind abandon punishment cost coefficient, $ MWh; pw,tThe actual output value, MW, of the wind power plant in the time period t; p' w, t is the output predicted value, MW of the wind power plant in the period t;
2. constraint conditions
Including power network constraints, natural gas network constraints, coupling unit constraints, and carbon capture system constraints that take into account carbon storage, hydrogen storage facilities. The method comprises the following specific steps:
(1) power network constraints
The power network constraint comprises a power balance constraint, a rotation standby constraint, a unit operation constraint and a power grid operation safety and reliability constraint; the concrete is as follows:
1) power balance constraint
Figure BDA0003350397790000122
In the formula, Ab(.)Representing an association matrix of a fossil fuel unit, a wind power plant, an electric-to-gas device and a power transmission line associated with a power grid node b; n is a radical ofFU、NW、NPtG、NTLThe total number of fossil fuel units, wind power plants, power-to-gas equipment and power transmission lines is respectively;
Figure BDA0003350397790000123
power, MW, to the input electrical to gas device; pED,b,tRepresents the electrical load at node b at time t, MW; PLl,tThe transmission power of the transmission line l at the time t, MW;
2) rotational back-up restraint
Figure BDA0003350397790000124
In the formula, τi,tThe variable is 0-1 and represents the start-stop state of the fossil fuel unit i at the moment t, the value is 1 when the unit i starts to operate, and the value is 0 when the unit i stops operating; pi,min、Pi,maxRespectively representing the lower limit, the upper limit and the MW of the output of the fossil fuel unit i;
Figure BDA0003350397790000125
respectively representing the up-regulation reserve capacity demand value and the down-regulation reserve capacity demand value MW which are at least required to be met by the power system at the moment t;
3) the unit operation constraints comprise wind power output constraints, fossil fuel unit output constraints, climbing rate constraint minimum startup and shutdown time constraints and startup and shutdown energy consumption constraints, and specifically comprise the following steps:
wind power output constraint:
0≤Pw,t≤P′w,t (7)
output constraint of fossil fuel units:
τi,tPi,min≤Pi,t≤τt,tPt,max (8)
third, climbing speed constraint:
Figure BDA0003350397790000132
in the formula, Sup,i、sdn,iRespectively representing the upward and downward climbing rates, MW, of the fossil fuel unit i; pi,t+1The output power, MW, of the fossil fuel unit i at the moment t + 1; tau isi,t+1The variable is 0-1 and represents the start-stop state of the fossil fuel unit i at the moment of t + 1;
and fourthly, constraint of minimum start-up and shut-down time:
Figure BDA0003350397790000133
in the formula (I), the compound is shown in the specification,
Figure BDA0003350397790000134
respectively representing the minimum time interval number h of starting and stopping the fossil fuel unit i;
Figure BDA0003350397790000135
respectively indicating the number h of the periods when the fossil fuel unit i is continuously started and stopped in the period t-1; tau isi,t-1And the variable is a 0-1 variable and represents the start-stop state of the fossil fuel unit i at the time t-1.
Energy consumption constraint of starting and stopping
Figure BDA0003350397790000136
In the formula, sui、sdiRespectively representing the starting and stopping energy consumption, MBtu, of the fossil fuel unit i; SUi,t、SDi,tRespectively representing the starting-up and stopping energy consumption, MBtu, of the fossil fuel unit i at the time t;
4) power grid operation safety and reliability constraint
PLl,t=(πl1,tl2,t)/xl (12)
In the formula, pil1,t、πl2,tThe voltage phase angle values, rad, of the head end and the tail end of the power transmission line l at the time t are respectively; x is the number oflIs the reactance value of line l, p.u.;
the line transmission power constraint and the node voltage angle constraint are as follows:
Figure BDA0003350397790000137
in the formula, PLl,min、PLl,maxRespectively the minimum and maximum active power, MW, transmitted by the line l; piref、πb,min、πb,max、πb,tRespectively, the voltage phase angle value of the reference node and the minimum and maximum voltage phase angle values of node b, the actual voltage phase angle value at time t, rad.
(2) Natural gas network constraints
The natural gas network constraints comprise gas source gas production constraints, natural gas pipeline flow constraints, natural gas storage equipment constraints, node pressure constraints and node supply and demand balance constraints;
1) beta gas confinement of gas source
Qω,min≤Qω,t≤Qω,max (14)
In the formula, Qω,min、Qω,maxThe unit of the lower limit value and the upper limit value of the gas production rate of the natural gas source omega is kcf:
2) natural gas pipeline flow restriction
For a natural gas pipeline mn, gas flows from m node to n node, and the pipeline gas flow is expressed as:
Figure BDA0003350397790000141
in the formula, Qmn,tThe gas flow rate of the natural gas pipeline mn flowing at the time t is kcf/h; sgn (p)m,t,pn,t) The direction of gas flow is shown, when the direction is 1, the gas flows from the m node to the n node, and when the direction is-1, the gas flows from the n node to the m node; p is a radical ofm,t、pn,tPressure values at nodes m and n at the time t, Psig; cmnThe size is a pipe constant and depends on diameter, length, coefficient of friction, gas composition, etc., kcf/Psig.
For a natural gas pipeline mn with a compressor j, the compression ratio between the output node and the input node is constrained by:
Figure BDA0003350397790000142
in the formula, Rj,max、Rj,minMaximum and minimum values of the allowable compressor compression ratio;
3) natural gas storage equipment restraint
The gas storage capacity of the gas storage equipment at any time t in the operation process
Figure BDA0003350397790000143
Expressed as:
Figure BDA0003350397790000144
in the formula (I), the compound is shown in the specification,
Figure BDA0003350397790000145
indicating the gas storage amount of the gas storage device s at the time t and t-1, respectively, kcf;
Figure BDA0003350397790000146
dividing the natural gas storage device s into the natural gas injection and extraction amount at t moment, kcf/h; Δ t is the time step;
Figure BDA0003350397790000147
respectively representing the minimum and maximum gas storage capacities of the gas storage device s, kcf;
Figure BDA0003350397790000148
respectively representing the maximum gas injection and extraction flow rates of natural gas storage equipment s, kcf/h;
Figure BDA0003350397790000149
respectively representing the gas capacity of the gas storage device s at the initial and ending moments of a scheduling period, kcf;
4) nodal pressure constraint
For the natural gas node m, the pressure magnitude should satisfy the following constraint:
pm,min≤pm,t≤pm,max (18)
in the formula, pm,min、pm,maxRespectively representing the lower limit and the upper limit of the pressure value of the node m;
5) node supply and demand balance constraints
For natural gas node m, the equilibrium equation is as follows:
Figure BDA0003350397790000151
in the formula, Nω、NGs、NPtG、NNU、NGPRespectively indicating the quantity of the gas source, the gas storage facility, the electric gas conversion equipment, the gas unit and the gas transmission pipeline; b ism(.)Representing a correlation matrix associated with a natural gas node m;
Figure BDA0003350397790000152
the amount of the methane gas synthesized by PtG equipment is kcf/h; qNi,tThe natural gas consumption of the gas unit i at the time t is kcf/h; qGL,m,tThe air load at the time t at the node m is kcf/h; qmn,tThe natural gas flow rate flowing through the pipeline mn in the period t is kef/h.
(3) Electric gas conversion equipment restraint
The hydrogen quantity generated by the electrolytic hydrogen production reaction of the electric gas conversion equipment
Figure BDA0003350397790000153
Expressed as:
Figure BDA0003350397790000154
in the formula (I), the compound is shown in the specification,
Figure BDA0003350397790000155
the unit of the hydrogen generated by the electrolysis hydrogen production reaction beta of the electrolysis equipment is kcf/h;
Figure BDA0003350397790000156
to input power to the electrical conversion device, MW:
Figure BDA0003350397790000157
hydrogen production efficiency in the electrolysis reaction; zetae,gMBtu/MWh, which is the equivalent thermal energy coefficient;
Figure BDA0003350397790000158
high heating value for hydrogen, MBtu/kcf;
the upper and lower limits of the power input into the electric gas conversion equipment in the electrolytic reaction are respectively 20% and 5% of the designed capacity, so that the input power of the electric gas conversion equipment can meet the following constraint:
the upper and lower limits of the power input into the electric gas conversion device in the electrolytic reaction are respectively 20% and 5% of the designed capacity, so the input power of the PtG device should satisfy the following constraint:
Figure BDA0003350397790000159
in the formula, τptgThe variable is 0-1, which represents the running state of the electric gas conversion device, and the value is 1 when the electric gas conversion device runs and 0 when the electric gas conversion device stops running;
Figure BDA00033503977900001510
the lower limit value, the upper limit value and the MW of the allowable input power of the electric gas conversion device are respectively set;
CO consumed by methanation2The amounts and the amount of methane gas produced are expressed as:
Figure BDA00033503977900001511
in the formula (I), the compound is shown in the specification,
Figure BDA0003350397790000161
CO consumed for methanation reactions2Amount, kcf/h;
Figure BDA0003350397790000162
the amount of methane gas synthesized by methanation reaction is kcf/h;
Figure BDA0003350397790000163
the hydrogen flow taken off for the hydrogen storage device s for the synthesis of methane, kcf/h;
Figure BDA0003350397790000164
the reaction coefficients between hydrogen and carbon dioxide, hydrogen and methane are shown;
meanwhile, the heat energy generated by methanation reaction is as follows:
Figure BDA0003350397790000165
in the formula (I), the compound is shown in the specification,
Figure BDA0003350397790000166
the recovered thermal power, MW, for the electric gas-converting equipment; ΨheatMWh/kcf for heat recovery coefficient;
(4) carbon capture system constraints
The carbon capture unit operation model is represented as follows:
Figure BDA0003350397790000167
in the formula, PcEquivalent output power, MW; pc netNet output power, MW; pc conEnergy consumption, MW, required to capture CO 2; maintenance energy consumption Pc mIs a fixed value, MW; pc oEnergy consumption for operation, MW;
Figure BDA0003350397790000168
CO that can be processed for carbon capture systems2Amount, kcf/h; gamma raycFor the carbon capture unit c to capture unit CO2The energy consumption generated is kfc/MWh;
the carbon capture unit generates equivalent output power and simultaneously generates carbon dioxide Q according to a certain proportionc
Figure BDA0003350397790000169
In the formula, mucCO for carbon capture train c2Emission intensity, kfc/MWh;
Figure BDA00033503977900001610
CO for carbon capture train c2Trapping capacity, kcf/h; beta is acCO for carbon capture train c2The capture rate;
in summary, the net output power of the carbon capture unit is expressed as:
Pc net=Pc-Pc m-Pc o (26)
and the net output power meets the output limit constraint:
Figure BDA00033503977900001611
in the formula (I), the compound is shown in the specification,
Figure BDA00033503977900001612
respectively the minimum and maximum net power values allowed to be output by the carbon capture unit;
(5) carbon and hydrogen storage equipment restraint
1) Hydrogen storage apparatus
Figure BDA0003350397790000171
In the formula (I), the compound is shown in the specification,
Figure BDA0003350397790000172
the hydrogen flow rate injected for the hydrogen storage device s, kcf/h;
Figure BDA0003350397790000173
the hydrogen storage amounts of the hydrogen storage equipment s at the times t and t-1, respectively, kcf;
Figure BDA0003350397790000174
dividing into kcf/h hydrogen injection and extraction quantity of hydrogen storage equipment s at t moment;
Figure BDA0003350397790000175
respectively representing the minimum and maximum hydrogen storage capacities of the hydrogen storage device s, kcf;
Figure BDA0003350397790000176
respectively representing the maximum hydrogen injection and extraction flow rates of the hydrogen storage device s, kcf/h; n is a radical ofTIs the end time of a scheduling period, h;
Figure BDA0003350397790000177
respectively representing the hydrogen capacities of the hydrogen storage equipment at the initial and end moments of a scheduling period, wherein the hydrogen capacities of the hydrogen storage equipment and the hydrogen capacity are equal, kcf;
the hydrogen required by the methanation equipment and the gas turbine is from the hydrogen storage equipment, so the hydrogen flow rate should satisfy the following constraints:
Figure BDA0003350397790000178
in the formula
Figure BDA0003350397790000179
Representing the hydrogen flow rate of the hydrogen storage device to the gas turbine, kcf/h;
2) carbon storage equipment
Carbon storage amount of carbon storage equipment at any moment in operation process
Figure BDA00033503977900001710
Expressed as:
Figure BDA00033503977900001711
in the formula (I), the compound is shown in the specification,
Figure BDA00033503977900001712
respectively representing CO of the carbon storage equipment at t and t-12Capacity, kfc;
Figure BDA00033503977900001713
respectively representing CO of carbon storage equipment s at time t2The injection and extraction flow rates are kcf/h;
corresponding carbon storage facility capacity constraint, CO2The injection quantity and the extraction quantity of the flow are restricted, and the CO at the initial and the end in a scheduling period2The capacity constraints are as follows:
Figure BDA00033503977900001714
in the formula (I), the compound is shown in the specification,
Figure BDA00033503977900001715
respectively representing minimum and maximum CO of carbon storage facility s2Storage capacity, kfc;
Figure BDA00033503977900001716
CO of carbon storage equipment s at t moment2The injection and extraction amount is kcf/h;
Figure BDA0003350397790000181
respectively representing CO of carbon storage units s2The maximum injection and extraction flow rate is kcf/h;
Figure BDA0003350397790000182
respectively representing CO of carbon storage equipment at the initial and end time of a scheduling period2Capacity, kcf;
CO trapped by carbon trapping system in carbon removal trapping unit2CO of carbon storage equipment2The injected amount also comes from CO captured in the air2Expressed as:
Figure BDA0003350397790000183
in the formula (I), the compound is shown in the specification,
Figure BDA0003350397790000184
the flow rate of carbon dioxide supplied from the carbon capture unit to the carbon storage facility s at time t is shown at kcf.
Fig. 9 is a schematic structural diagram of an electronic device according to an embodiment of the present invention, and referring to fig. 9, the electronic device may include: a processor (processor)201, a communication Interface (communication Interface)202, a memory (memory)203 and a communication bus 204, wherein the processor 201, the communication Interface 202 and the memory 203 complete communication with each other through the communication bus 204. The processor 201 may call logic instructions in the memory 203 to perform the following method: constructing a low-carbon economic dispatching model of the carbon-containing trapped electrical comprehensive energy system; the low-carbon economic dispatching model of the carbon-containing captured electrical comprehensive energy system is based on system operation cost and CO2The minimum comprehensive cost of the sum of the related cost and the wind abandoning penalty cost is a target function, and the power network constraint, the natural gas network constraint and the electricity-to-gas conversion are usedThe equipment constraint, the carbon capture system constraint and the carbon storage and hydrogen storage equipment constraint are constraint conditions; said CO2Associated costs include CO emissions from fossil fuel units2Carbon tax cost FCT、CO2Transmission storage cost FTSCarbon trapping unit for trapping CO from air2Cost F ofCC(ii) a Then, solving is carried out, and scheduling is carried out according to the solving result.
In addition, the logic instructions in the memory 203 may be implemented in the form of software functional units and stored in a computer readable storage medium when the logic instructions are sold or used as independent products. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
In another aspect, an embodiment of the present invention further provides a non-transitory computer-readable storage medium, on which a computer program is stored, which when executed by a processor, is implemented to perform a method for low-carbon optimal scheduling of an electricity-gas integrated energy system considering carbon capture according to the embodiments, for example, the method includes: constructing a low-carbon economic dispatching model of the carbon-containing trapped electrical comprehensive energy system; the low-carbon economic dispatching model of the carbon-containing captured electrical comprehensive energy system is based on system operation cost and CO2The minimum comprehensive cost of the sum of the related cost and the wind abandon punishment cost is an objective function, and the constraint conditions are power network constraint, natural gas network constraint, electric-to-gas equipment constraint, carbon capture system constraint, carbon storage and hydrogen storage equipment constraint; said CO2Associated costs include fossil fuel unit emissionsCO2Carbon tax cost FCT、CO2Transmission storage cost FTSCarbon trapping unit for trapping CO from air2Cost F ofCC(ii) a Then, solving is carried out, and scheduling is carried out according to the solving result.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium, such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute the method according to the embodiments or some parts of the embodiments.
Examples of the applications
In the description, based on a 6-node power system and a 7-node natural gas system network model in fig. 1, three simulation scenes are set for comparative simulation analysis, and the feasibility and the effectiveness of the method are verified.
Scene 1: environmental economic dispatch of typical integrated energy system
Scene 2: integrated energy system environmental economic dispatch of carbonaceous capture system
Scene 3: comprehensive energy system environment economic dispatch in combined operation mode of carbon storage and hydrogen storage equipment
And setting the parameters of the system as follows: the prediction curves of the power load, the natural gas load and the wind power output are shown in fig. 2. In addition, the rotating standby requirement value in the system scheduling period in the power network is 30% of the total load; the lower limit value of the allowable input power of the electric gas conversion equipment is 10MW, and the upper limit value of the allowable input power of the electric gas conversion equipment is 100 MW; maintaining energy consumption of the carbon capture system to be a fixed value, and taking 0.5% of the installed capacity of the unit; carbon tax value of 20$/t, transmission storage unit CO2The price of the carbon dioxide is 5$/t, and the unit CO is captured from the air2The required price cost takes the value of 200 $/t. And (4) carrying out optimization analysis on the electricity-gas comprehensive energy system within 24h in one day by taking 1h as a time interval.
The method comprises the steps of calling Gurobi optimization software through MATLAB to carry out simulation solution on environment economic optimization scheduling of an electricity-gas comprehensive energy network in three scenes within 24h a day, and obtaining wind power consumption, net output power of a carbon capture unit, system carbon emission and target optimization results in the three scenes. The results show that:
1) wind power consumption condition and analysis: as can be seen from the wind power consumption comparison graphs in the three scenes in fig. 3, a large amount of abandoned wind occurs in the scene 1, so that the penalty cost of the system is increased, a small amount of abandoned wind still exists in the scene 2, and the carbon storage and hydrogen storage units are added in the scene 3, so that the flexible operation capacity of the carbon capture power plant and the electricity-to-gas conversion equipment can be improved, and the full consumption of the wind power can be basically realized in the scheduling period.
2) Analyzing the output condition of the carbon capture unit: fig. 4 shows the net output power curve of the unit 5 in three scenarios, in scenario 2, the unit itself provides the capture energy required for the operation of the capture system, so the net output of the unit 5 in scenario 2 is lower than that in scenario 1. In the combined operation mode of scene 3, the carbon capture system recycles the heat released by the methanation reaction (as shown in fig. 5), and reduces the operation energy consumption of the carbon capture system, so that the net output power of the unit 5 is higher than that of scene 2. In addition, the recovered heat energy cannot fully compensate for the energy consumption of the carbon capture operation, and the unit itself still needs to provide a part of the energy consumption, so the net output of the unit 5 in the scene 3 is still lower than that in the scene 1 as a whole.
3) Carbon emission analysis under three scenarios: CO under three scenes in FIG. 62Emission and CO2The capture amount graph shows that the carbon capture unit 5 can capture CO in scene 32The most amount, the lowest carbon emission.
4) Target optimization results and analysis: CO under different scenes2The associated cost pairs are shown in table 1. CO Capture from air in scenario 32Greatly reduced cost, CO2The associated cost is optimal under 3 scenes. The system cost and the comprehensive cost under the three scenes are shown in a table 2, and the table shows that the system comprehensive cost, the system operation cost, the CO2 related cost and the wind abandon penalty cost in the scene 3 are the optimal values under the three scenes.
TABLE 1 CO under different scenarios2Correlation cost comparison table
Scene CO2Associated cost/$ Carbon tax cost/$ CO capture from air2Cost/$ Cost of transfer storage/$
Scene 1 196336.64 190950.00 5386.64 0
Scene 2 134472.14 131405.42 1532.45 1534.27
Scene 3 132969.90 131364.98 0 1601.92
TABLE 2 comparison table of comprehensive cost of system under different scenes
Scene Integrated cost/$ System operating cost/$ CO2Associated cost/$ Wind curtailment penalty cost/$
Scene 1 498563.89 267576.90 196336.64 34650.35
Scene 2 422075.47 282461.58 134472.14 5141.75
Scene 3 387701.92 254732.02 132969.90 0
Compared with the traditional electricity-gas comprehensive energy system, the invention provides a combined operation mode of considering carbon storage and hydrogen storage equipment between the carbon capture equipment and the electricity-gas conversion equipment by introducing the carbon capture system, thereby realizing the carbon circulation of the comprehensive energy system in a scheduling period, reducing the carbon emission and lowering the environmental cost of the system; the wind power is fully consumed in the dispatching period, the net output power of the carbon capture power plant is improved, and the optimal economic dispatching in all days is realized.
The foregoing shows and describes the general principles, essential features, and advantages of the invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are given by way of illustration of the principles of the present invention, and that various changes and modifications may be made without departing from the spirit and scope of the invention as defined by the appended claims. The scope of the invention is defined by the appended claims and equivalents thereof.

Claims (10)

1. A low-carbon optimized dispatching method of an electricity-gas integrated energy system considering carbon capture is characterized in that the electricity-gas integrated energy system with carbon capture is adopted;
the carbon-containing trapped electricity-gas comprehensive energy system comprises a conventional thermal power generating unit, a carbon trapping unit, carbon storage equipment, methanation equipment, hydrogen storage equipment, a gas turbine and electrolysis equipment; the carbon capture unit is a conventional thermal power unit additionally provided with a carbon capture system;
the conventional thermal power generating unit is connected with a carbon capture system of the carbon capture unit;
the carbon capture unit is respectively connected with the carbon storage equipment and the methanation equipment;
the methanation equipment is also respectively connected with the carbon storage equipment and the hydrogen storage equipment;
the hydrogen storage equipment is also respectively connected with the gas turbine and the electrolysis equipment;
the scheduling method comprises the following steps:
constructing a low-carbon economic dispatching model of the carbon-containing trapped electricity-gas integrated energy system; the carbon-containing trapped electricity-gas integrated energy system low-carbon economic dispatching model is based on system operation cost and CO2The minimum comprehensive cost of the sum of the related cost and the wind abandon punishment cost is an objective function, and the constraint conditions are power network constraint, natural gas network constraint, electric-to-gas equipment constraint, carbon capture system constraint, carbon storage and hydrogen storage equipment constraint; said CO2Associated costs include CO emissions from fossil fuel units2Carbon tax cost FCT、CO2Transmission storage cost FTSCarbon trapping unit for trapping CO from air2Cost F ofCC(ii) a Then, solving is carried out, and scheduling is carried out according to the solving result.
2. The carbon capture-considered low-carbon optimal scheduling method for the electricity-gas integrated energy system according to claim 1, wherein the objective function is as follows:
minF=FO+FC+FW(1)
in the formula, FO$ for system operating cost, $; fCIs CO2Associated cost, $; fWPunishing cost for wind abandonment;
(1) the system operation cost comprises the operation cost F of the conventional thermal power generating unitGGas production cost F of natural gas sourceGas cost F of natural gas storage facilityNS,$;
Figure FDA0003350397780000011
In the formula, NG、Nω、NGsThe number of the conventional coal burners, the number of natural gas sources and the number of natural gas storage devices are respectively; n is a radical ofTIs a scheduling period; beta is ai、βω、βsThe unit is the fuel price of a conventional coal-fired unit i, the unit gas production cost of a natural gas source omega and the unit gas consumption cost of a natural gas storage device s, and the units are $/MBtu, $/kcf and $/kcf; a isGi、bGi、cGiThe cost coefficient of the ith conventional thermal power generating unit; pGi,tThe output power, MW, of the conventional thermal power generating unit i at the moment t; qω,tThe gas production rate of the natural gas source omega at the time t, kcf;
Figure FDA0003350397780000026
dividing the natural gas storage device s into the natural gas injection and extraction amount at the time t, kcf;
(2)CO2associated costs include CO emissions from fossil fuel units2Carbon tax cost FCT,$、CO2Transmission storage cost FTSAnd the carbon capture system captures CO from the air2Cost F ofCC,$;
Figure FDA0003350397780000021
In the formula, NFU、NCC、NCSRespectively representing the number of fossil fuel units, the number of carbon gathering units and the number of carbon storage equipment; beta is aCT、βTS、βCCRespectively carbon tax price and unit CO2Transportation and storage price of, capture of unit CO from air2The required price cost in $/kcf; u. ofiCO for conventional coal-fired unit i2The discharge intensity; pi,tThe output power, MW, of the fossil fuel unit i at the moment t;
Figure FDA0003350397780000022
for CO trapped by carbon trapping unit at time t2Amount captured, kfc;
Figure FDA0003350397780000023
the flow rate of carbon dioxide supplied from the air to the carbon storage facility s, kfc/h;
Figure FDA0003350397780000024
the carbon dioxide flow rate for supplying the carbon capture unit with the carbon storage device s is kfc/h; the fossil fuel unit comprises a conventional thermal power unit and a gas unit;
(3) wind abandon penalty cost:
Figure FDA0003350397780000025
in the formula, NWThe number of wind power plants; cw is a wind abandon punishment cost coefficient, $ MWh; pw,tThe actual output value, MW, of the wind power plant in the time period t is obtained; and P' w and t are output predicted values, MW, of the wind power plant in the t period.
3. The carbon capture-considered low-carbon optimal scheduling method for the electricity-gas integrated energy system according to claim 1, wherein the power network constraints comprise power balance constraints, rotating standby constraints, unit operation constraints and grid operation safety and reliability constraints; the method specifically comprises the following steps:
1) power balance constraint
Figure FDA0003350397780000031
In the formula, Ab(.)Representing an incidence matrix of a fossil fuel unit, a wind power plant, an electric-to-gas device and a power transmission line which are associated with a power grid node b; n is a radical ofFU、NW、NPtG、NTLThe total number of fossil fuel units, wind power plants, power-to-gas equipment and power transmission lines is respectively;
Figure FDA0003350397780000032
power, MW, input to the electric gas-to-gas device; pED,b,tRepresents the electrical load at node b at time t, MW; PLl,tThe transmission power of the transmission line l at the moment t, MW;
2) rotational back-up restraint
Figure FDA0003350397780000033
In the formula, τi,tThe variable is 0-1 and represents the start-stop state of the fossil fuel unit i at the moment t, the value is 1 when the unit i starts to operate, and the value is 0 when the unit i stops operating; pi,min、Pi,maxRespectively representing the lower limit, the upper limit and the MW of the output of the fossil fuel unit i;
Figure FDA0003350397780000034
respectively representing the up-regulation reserve capacity demand value and the down-regulation reserve capacity demand value MW which are at least required to be met by the power system at the moment t;
3) the unit operation constraints comprise wind power output constraints, fossil fuel unit output constraints, climbing rate constraint minimum startup and shutdown time constraints and startup and shutdown energy consumption constraints, and specifically comprise the following steps:
wind power output constraint:
0≤Pw,t≤P′w,t(7)
output constraint of fossil fuel units:
τi,tPi,min≤Pi,t≤τi,tPi,max(8)
third, climbing speed constraint:
Figure FDA0003350397780000035
in the formula, Sup,i、Sdn,iRespectively representing the upward and downward climbing rates, MW, of the fossil fuel unit i;Pi,t+1the output power, MW, of the fossil fuel unit i at the moment t + 1; tau isi,t+1The variable is 0-1 and represents the start-stop state of the fossil fuel unit i at the moment of t + 1;
and fourthly, constraint of minimum start-up and shut-down time:
Figure FDA0003350397780000036
in the formula (I), the compound is shown in the specification,
Figure FDA0003350397780000041
respectively representing the minimum time interval number h of starting and stopping the fossil fuel unit i;
Figure FDA0003350397780000042
Figure FDA0003350397780000043
respectively representing the number h of the continuous startup and shutdown time periods of the fossil fuel unit i in the t-1 period; tau isi,t-1The variable is 0-1 and represents the start-stop state of the fossil fuel unit i at the time t-1;
energy consumption constraint of starting and stopping
Figure FDA0003350397780000044
In the formula, sui、sdiRespectively representing the starting and stopping energy consumption, MBtu, of the fossil fuel unit i; SUi,t、SDi,tRespectively representing the starting-up and stopping energy consumption, MBtu, of the fossil fuel unit i at the time t;
4) power grid operation safety and reliability constraint
PLl,t=(πl1,tl2,t)/xl(12)
In the formula, pil1,t、πl2,tThe voltage phase angle values, rad, of the head end and the tail end of the power transmission line l at the time t are respectively; x is the number oflIs the reactance of the line lValue, p.u.;
the line transmission power constraint and the node voltage angle constraint are as follows:
Figure FDA0003350397780000045
in the formula, PLl,min、PLl,maxRespectively the minimum and maximum active power, MW, transmitted by the line l; piref、πb,min、πb,max、πb,tRespectively, the voltage phase angle value of the reference node and the minimum and maximum voltage phase angle values of node b, the actual voltage phase angle value at time t, rad.
4. The carbon capture-considered low-carbon optimal scheduling method for the electricity-gas integrated energy system is characterized in that the natural gas network constraints comprise gas source gas production constraints, natural gas pipeline flow constraints, natural gas storage equipment constraints, node pressure constraints and node supply and demand balance constraints;
1) gas source gas production constraint
Qω,min≤Qω,t≤Qω,max(14)
In the formula, Qω,min、Qω,maxThe unit of the lower limit value and the upper limit value of the omega gas production of the natural gas source is kcf;
2) natural gas pipeline flow restriction
For a natural gas pipeline mn, gas flows from m node to n node, and the pipeline gas flow is expressed as:
Figure FDA0003350397780000051
in the formula, Qmn,tThe gas flow rate of the natural gas pipeline mn flowing at the time t is kcf/h; sgn (p)m,t,pn,t) The direction of gas flow is shown, when the direction is 1, the gas flows from the m node to the n node, and when the direction is-1, the gas flows from the n node to the m node; p is a radical ofm,t、pn,tPressure values at nodes m and n at the time t, Psig; cmnIs the pipeline constant;
for a natural gas pipeline mn with a compressor j, the compression ratio between the output node and the input node is constrained by:
Figure FDA0003350397780000052
in the formula, Rj,max、Rj,minMaximum and minimum values of the allowable compressor compression ratio;
3) natural gas storage equipment restraint
The gas storage capacity of the gas storage equipment at any time t in the operation process
Figure FDA0003350397780000053
Expressed as:
Figure FDA0003350397780000054
in the formula (I), the compound is shown in the specification,
Figure FDA0003350397780000055
indicating the gas storage amount of the gas storage device s at the time t and t-1, respectively, kcf;
Figure FDA0003350397780000056
dividing the natural gas storage device s into the natural gas injection and extraction amount at t moment, kcf/h; Δ t is the time step;
Figure FDA0003350397780000057
respectively representing the minimum and maximum gas storage capacities of the gas storage device s, kcf;
Figure FDA0003350397780000058
respectively representing the maximum gas injection and extraction flow rates of natural gas storage equipment s, kcf/h;
Figure FDA0003350397780000059
respectively representing the gas capacity of the gas storage device s at the initial and ending moments of a scheduling period, kcf;
4) nodal pressure constraint
For the natural gas node m, the pressure magnitude should satisfy the following constraint:
pm,min≤pm,t≤pm,max(18)
in the formula, pm,min、pm,maxRespectively representing the lower limit and the upper limit of the pressure value of the node m;
5) node supply and demand balance constraints
For natural gas node m, the equilibrium equation is as follows:
Figure FDA0003350397780000061
in the formula, Nω、NGs、NPtG、NNU、NGPRespectively indicating the quantity of the gas source, the gas storage facility, the electric gas conversion equipment, the gas unit and the gas transmission pipeline; b ism(.)Representing a correlation matrix associated with a natural gas node m;
Figure FDA0003350397780000062
the amount of the methane gas synthesized by PtG equipment is kcf/h; qNi,tThe natural gas consumption of the gas unit i at the time t is kcf/h; qGL,m,tThe air load at the time t at the node m is kcf/h; qmn,tThe natural gas flow rate flowing through the pipeline mn in the period t is kcf/h.
5. The carbon capture-considered low-carbon optimal scheduling method for the electricity-gas integrated energy system according to claim 1, wherein the constraints of the electricity-to-gas equipment are as follows:
the hydrogen quantity generated by the electrolytic hydrogen production reaction of the electric gas conversion equipment
Figure FDA0003350397780000063
Expressed as:
Figure FDA0003350397780000064
in the formula (I), the compound is shown in the specification,
Figure FDA0003350397780000065
the unit of the hydrogen quantity generated by the electrolytic hydrogen production reaction of the electrolysis equipment is kcf/h;
Figure FDA0003350397780000066
power, MW, for input to the electrical to gas equipment;
Figure FDA0003350397780000067
hydrogen production efficiency in the electrolysis reaction; zetae,gMBtu/MWh, which is the equivalent thermal energy coefficient;
Figure FDA00033503977800000615
high heating value for hydrogen, MBtu/kcf;
the upper and lower limits of the power input into the electric gas conversion equipment in the electrolytic reaction are respectively 20% and 5% of the designed capacity, so that the input power of the electric gas conversion equipment can meet the following constraint:
Figure FDA0003350397780000068
in the formula, τptgThe variable is 0-1, which represents the running state of the electric gas conversion device, and the value is 1 when the electric gas conversion device runs and 0 when the electric gas conversion device stops running;
Figure FDA0003350397780000069
the lower limit value, the upper limit value and the MW of the allowable input power of the electric gas conversion device are respectively set;
CO consumed by methanation2The amounts and the amount of methane gas produced are expressed as:
Figure FDA00033503977800000610
in the formula (I), the compound is shown in the specification,
Figure FDA00033503977800000611
CO consumed for methanation reactions2Amount, kcf/h;
Figure FDA00033503977800000612
the amount of methane gas synthesized by methanation reaction is kcf/h;
Figure FDA00033503977800000613
the hydrogen flow taken off for the hydrogen storage device s for the synthesis of methane, kcf/h;
Figure FDA00033503977800000614
the reaction coefficients between hydrogen and carbon dioxide, hydrogen and methane;
meanwhile, the heat energy generated by methanation reaction is as follows:
Figure FDA0003350397780000071
in the formula (I), the compound is shown in the specification,
Figure FDA0003350397780000072
the recovered thermal power, MW, for the electric gas-converting equipment; ΨheatMWh/kcf is the heat recovery factor.
6. The carbon capture-considered low-carbon optimal scheduling method for the electric-gas integrated energy system according to claim 1, wherein the carbon capture system constraints are as follows:
the carbon capture unit operation model is represented as follows:
Figure FDA0003350397780000073
in the formula, PcEquivalent output power, MW;
Figure FDA0003350397780000074
net output power, MW;
Figure FDA0003350397780000075
for capturing CO2Required energy consumption, MW; maintenance of energy consumption
Figure FDA0003350397780000076
Is a fixed value, MW;
Figure FDA0003350397780000077
energy consumption for operation, MW;
Figure FDA0003350397780000078
CO that can be processed for carbon capture systems2Amount, kcf/h; gamma raycFor the carbon capture unit c to capture unit CO2The energy consumption generated is kfc/MWh;
the carbon capture unit generates equivalent output power and simultaneously generates carbon dioxide Q according to a certain proportionc
Figure FDA0003350397780000079
In the formula, mucCO for carbon capture train c2Emission intensity, kfc/MWh;
Figure FDA00033503977800000710
CO for carbon capture train c2Trapping capacity, kcf/h; beta is acCO for carbon capture train c2The capture rate;
in summary, the net output power of the carbon capture unit is expressed as:
Figure FDA00033503977800000711
and the net output power meets the output limit constraint:
Figure FDA00033503977800000712
in the formula (I), the compound is shown in the specification,
Figure FDA00033503977800000713
respectively, the minimum and maximum net power values allowed to be output by the carbon capture unit.
7. The carbon capture-considered low-carbon optimal scheduling method for the electricity-gas integrated energy system is characterized in that the carbon storage and hydrogen storage equipment constraints are as follows:
1) hydrogen storage apparatus
Figure FDA0003350397780000081
In the formula (I), the compound is shown in the specification,
Figure FDA0003350397780000082
the hydrogen flow rate injected for the hydrogen storage device s, kcf/h;
Figure FDA0003350397780000083
the hydrogen storage amounts of the hydrogen storage equipment s at the times t and t-1, respectively, kcf;
Figure FDA0003350397780000084
dividing into kcf/h hydrogen injection and extraction quantity of hydrogen storage equipment s at t moment;
Figure FDA0003350397780000085
respectively representing the minimum and maximum of the hydrogen storage apparatus sHydrogen storage capacity, kcf;
Figure FDA0003350397780000086
respectively representing the maximum hydrogen injection and extraction flow rates of the hydrogen storage equipment s, kcf/h; NT is the ending time of a scheduling period, h;
Figure FDA0003350397780000087
respectively representing the hydrogen capacities of the hydrogen storage equipment at the initial and end moments of a scheduling period, wherein the hydrogen capacities of the hydrogen storage equipment and the hydrogen capacity are equal, kcf;
the hydrogen required by the methanation equipment and the gas turbine is from the hydrogen storage equipment, so the hydrogen flow rate should satisfy the following constraints:
Figure FDA0003350397780000088
in the formula
Figure FDA0003350397780000089
Representing the hydrogen flow rate of the hydrogen storage device to the gas turbine, kcf/h;
2) carbon storage equipment
Carbon storage amount of carbon storage equipment at any moment in operation process
Figure FDA00033503977800000810
Expressed as:
Figure FDA00033503977800000811
in the formula (I), the compound is shown in the specification,
Figure FDA00033503977800000812
respectively representing CO of the carbon storage equipment at t and t-12Capacity, kfc;
Figure FDA00033503977800000813
Figure FDA00033503977800000814
respectively representing CO of carbon storage equipment s at time t2The injection and extraction flow rates are kcf/h;
corresponding carbon storage facility capacity constraint, CO2The injection quantity and the extraction quantity of the flow are restricted, and the CO at the initial and the end in a scheduling period2The capacity constraints are as follows:
Figure FDA0003350397780000091
in the formula (I), the compound is shown in the specification,
Figure FDA0003350397780000092
respectively representing minimum and maximum CO of carbon storage facility s2Storage capacity, kfc;
Figure FDA0003350397780000093
CO of carbon storage equipment s at t moment2The injection and extraction amount is kcf/h;
Figure FDA0003350397780000094
respectively representing CO of carbon storage units s2The maximum injection and extraction flow rate is kcf/h;
Figure FDA0003350397780000095
respectively representing CO of carbon storage equipment at the initial and end time of a scheduling period2Capacity, kcf;
CO trapped by carbon trapping system in carbon removal trapping unit2CO of carbon storage equipment2The injected amount also comes from CO captured in the air2Expressed as:
Figure FDA0003350397780000096
in the formula (I), the compound is shown in the specification,
Figure FDA0003350397780000097
the flow rate of carbon dioxide supplied from the carbon capture unit to the carbon storage facility s at time t is shown at kcf.
8. The carbon capture-containing electro-gas integrated energy system of claim 1.
9. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the program performs the steps of the method for carbon capture-aware low carbon optimal scheduling of an electric-gas integrated energy system according to any one of claims 1 to 7.
10. A non-transitory computer readable storage medium having stored thereon a computer program, wherein the computer program when executed by a processor implements the steps of the carbon capture aware low carbon optimal scheduling method for an electricity-gas integrated energy system according to any of claims 1 to 7.
CN202111335666.9A 2021-11-11 2021-11-11 Low-carbon optimized scheduling method of electricity-gas integrated energy system considering carbon capture Pending CN114037292A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111335666.9A CN114037292A (en) 2021-11-11 2021-11-11 Low-carbon optimized scheduling method of electricity-gas integrated energy system considering carbon capture

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111335666.9A CN114037292A (en) 2021-11-11 2021-11-11 Low-carbon optimized scheduling method of electricity-gas integrated energy system considering carbon capture

Publications (1)

Publication Number Publication Date
CN114037292A true CN114037292A (en) 2022-02-11

Family

ID=80144139

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111335666.9A Pending CN114037292A (en) 2021-11-11 2021-11-11 Low-carbon optimized scheduling method of electricity-gas integrated energy system considering carbon capture

Country Status (1)

Country Link
CN (1) CN114037292A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114545878A (en) * 2022-02-22 2022-05-27 山东大学 Optimized scheduling method and system for comprehensive energy system
CN115293518A (en) * 2022-07-13 2022-11-04 河海大学 Low-carbon economic dispatching method for gas-electric coupling comprehensive energy system considering flexible climbing
CN117314708A (en) * 2023-10-26 2023-12-29 南方电网能源发展研究院有限责任公司 Carbon emission reduction system optimization method, device, equipment, medium and program product
CN117639069A (en) * 2023-11-28 2024-03-01 国网青海省电力公司清洁能源发展研究院 Light-hydrogen-storage comprehensive energy system suitable for remote agriculture and animal husbandry area

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114545878A (en) * 2022-02-22 2022-05-27 山东大学 Optimized scheduling method and system for comprehensive energy system
CN114545878B (en) * 2022-02-22 2023-09-01 山东大学 Comprehensive energy system optimal scheduling method and system
CN115293518A (en) * 2022-07-13 2022-11-04 河海大学 Low-carbon economic dispatching method for gas-electric coupling comprehensive energy system considering flexible climbing
CN115293518B (en) * 2022-07-13 2023-05-16 河海大学 Low-carbon economic dispatching method of gas-electricity coupling comprehensive energy system considering flexible climbing
CN117314708A (en) * 2023-10-26 2023-12-29 南方电网能源发展研究院有限责任公司 Carbon emission reduction system optimization method, device, equipment, medium and program product
CN117639069A (en) * 2023-11-28 2024-03-01 国网青海省电力公司清洁能源发展研究院 Light-hydrogen-storage comprehensive energy system suitable for remote agriculture and animal husbandry area

Similar Documents

Publication Publication Date Title
CN114037292A (en) Low-carbon optimized scheduling method of electricity-gas integrated energy system considering carbon capture
CN109919399B (en) Day-ahead economic dispatching method and system for comprehensive energy system
CN115018260A (en) Peak regulation optimal scheduling method of cascade water-light complementary power generation system
CN109741110A (en) A kind of wind hydrogen system combined optimization modeling method based on chance constrained programming
Huang et al. Assessment of an integrated energy system embedded with power-to-gas plant
CN113256045A (en) Park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty
CN113205273A (en) Low-carbon power supply planning method and system considering off-site electric energy transaction
CN114583725A (en) Hydrogen-based near-zero carbon emission comprehensive energy system and operation optimization method thereof
CN113806952B (en) Cold-hot electricity comprehensive energy system considering source-charge-storage and optimal operation method thereof
CN114282828A (en) Carbon circulation system and application thereof
CN117254503A (en) Source-net-storage-load double-layer collaborative low-carbon scheduling method, equipment and medium based on wind power consumption
CN112418537A (en) Optimized scheduling method for multi-energy cloud energy storage system
CN112270096A (en) Flexible resource contribution degree solving method for promoting new energy consumption
CN109687423B (en) Energy storage configuration method of direct current power grid
CN116646984A (en) Gas-water coordination peak shaving method considering system operation risk
CN111210119A (en) Establishment method of VPP electricity, heat and gas optimized scheduling model in various markets
CN114387129B (en) Comprehensive energy scheduling method, device, equipment and medium considering natural gas pipe network
CN115293518A (en) Low-carbon economic dispatching method for gas-electric coupling comprehensive energy system considering flexible climbing
Chen et al. Integrated energy system day-ahead dispatch considering supply-side and demand-side flexibility
CN114881328A (en) Comprehensive energy system economic dispatching method considering gas network hydrogen mixing and low-carbon reward
CN114091917A (en) Dynamic environment-friendly economic dispatching method and system for cold-heat-power combined supply type micro-grid
CN113742944A (en) Virtual power plant modeling method considering electric hydrogen production system
CN117216963B (en) Comprehensive energy system operation method for carbon-hydrogen utilization
CN218673288U (en) Park new forms of energy supply system
CN117171944A (en) Low-carbon economic dispatching method for electric coupling system considering hydrogen loading of natural gas network

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