CN110474335B - Comprehensive energy system operation method based on evaluation criterion - Google Patents

Comprehensive energy system operation method based on evaluation criterion Download PDF

Info

Publication number
CN110474335B
CN110474335B CN201910883976.0A CN201910883976A CN110474335B CN 110474335 B CN110474335 B CN 110474335B CN 201910883976 A CN201910883976 A CN 201910883976A CN 110474335 B CN110474335 B CN 110474335B
Authority
CN
China
Prior art keywords
heat
electricity
formula
grid
representing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910883976.0A
Other languages
Chinese (zh)
Other versions
CN110474335A (en
Inventor
王志贺
刘元园
沈润
渠卫东
王晓辉
耿伟
张思义
王菲
肖伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Xuzhou Power Supply Co
State Grid Corp of China SGCC
Original Assignee
State Grid Xuzhou Power Supply Co
State Grid Corp of China SGCC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Xuzhou Power Supply Co, State Grid Corp of China SGCC filed Critical State Grid Xuzhou Power Supply Co
Priority to CN201910883976.0A priority Critical patent/CN110474335B/en
Publication of CN110474335A publication Critical patent/CN110474335A/en
Application granted granted Critical
Publication of CN110474335B publication Critical patent/CN110474335B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/008Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Engine Equipment That Uses Special Cycles (AREA)
  • Control Of Steam Boilers And Waste-Gas Boilers (AREA)

Abstract

The invention discloses an evaluation criterion-based operation method of a comprehensive energy system, which comprises the steps of firstly establishing a sub-supply system model, and calculating the electricity purchasing quantity and the fuel consumption quantity of the sub-supply system; then establishing a comprehensive energy system model with two modes of 'fixing power with heat' and 'fixing heat with electricity', and respectively calculating the electricity purchasing amount and the fuel consumption amount of the comprehensive energy system; three evaluation indexes of the system are established: primary energy consumption, operating costs, CO2The method comprises the steps of reducing the discharge capacity, integrating three indexes to establish a comprehensive evaluation index, and determining the operation mode of a comprehensive energy system by calculating the application conditions of two main operation modes of a cooling, heating and power system, so that the utilization rate of system energy is optimal, and the method is more economical and energy-saving than the operation modes of electricity for heat determination and electricity for heat determination. The invention well solves the problem of resource waste caused by inconsistent load randomness and partial load characteristics of the unit when a certain operation mode is adopted.

Description

Comprehensive energy system operation method based on evaluation criterion
Technical Field
The invention belongs to the field of combined cooling heating and power systems, and particularly relates to an evaluation criterion-based operation method of a comprehensive energy system.
Background
An Integrated Energy System (IES) has the characteristics of modularization and decentralization, is a high-efficiency and reliable energy conversion unit arranged around a user side, and a power supply of the IES mainly comprises an internal combustion engine, a micro gas turbine, photovoltaic power generation, wind power generation, biomass power generation and the like, so that the IES can realize near power generation, near grid connection, near conversion and near use of scattered energy, and effectively improve the comprehensive utilization of energy. The development of the comprehensive energy system breaks through the traditional energy system architecture, becomes an important link for linking the energy production links of various forms to the energy consumption links of terminal users, and creates a good technical support platform and a resource allocation center for the coordinated operation and advantage complementation among various energy sources.
The combined cooling heating and power system is a comprehensive energy system with the most potential and prospect, is a comprehensive energy production and energy utilization distributed system established on the basis of energy gradient utilization, is usually installed near an end user end, firstly utilizes primary energy to drive an engine to generate electricity, and then recycles waste heat through various waste heat utilization devices, thereby simultaneously providing power, refrigeration, heating, domestic hot water and the like for users. The operation mode of the comprehensive energy system containing combined supply of cold, heat and electricity mainly comprises a power supply with fixed heat (FEL) and a power supply with fixed heat (FTL). The method comprises the following steps of firstly meeting the heat requirement (cold and hot load) of a system in a heat power ordering mode, determining the power output of the system through a thermoelectric ratio, and meeting the shortage of power by purchasing power from a large power grid; the electric heating mode is used for firstly meeting the power requirement of the system, so that the heat output of the system is determined, and the insufficient heat is supplied by the afterburning boiler. For a combined cooling heating and power system, when the load demand thermoelectric ratio is the same as that of a generator set, the utilization rate of energy is optimal, but the load randomness and the partial load characteristics of the generator set cause inconsistency, and resource waste is caused when a certain operation mode is adopted.
Disclosure of Invention
The invention provides an evaluation criterion-based operation method of a comprehensive energy system, which is characterized in that when the evaluation criterion is fixed, the economical efficiency and the environmental protection property of the operation of the comprehensive energy system are improved by calculating the applicable conditions of two main operation modes of a cooling, heating and power system.
In order to solve the technical problem, the comprehensive energy system operation method based on the evaluation criterion comprises the following steps:
step 10) establishing a model of a sub-supply System (SP), firstly determining the power consumption of the electric refrigerator as follows:
Figure BDA0002205214100000021
in the formula, COPecRepresenting the coefficient of performance, P, of the electric refrigeratorecRepresents the input power of the electric refrigerator, unit: kW; qCLRepresents the system cooling load, unit: kW.
The fuel consumed by the boiler is then determined as:
Figure BDA0002205214100000022
in the formula, FboilerNatural gas, which represents the consumption of a gas boiler, unit: cubic meter; qrhRepresents heat exchanger input heat, in units: kW; η SP builder represents the efficiency of the gas boiler.
Then, determining that the power purchased by the power grid is as follows:
Pgrid=PEL+Pec (3)
in the formula, PgridRepresenting the amount of electricity purchased from the grid, in units: kW; pELIs the electrical load of the system, unit: kW.
The fuel consumed by the system is then determined, only the gas boiler consuming:
Fm=Fboiler (4)
in the formula, FmFuel, expressed as system consumption, unit: cubic meter.
Converting the purchase electricity into corresponding fuel consumption:
Figure BDA0002205214100000023
in the formula, EmRepresenting the conversion of the purchase electricity into the corresponding fuel consumption, unit: standSquare rice.
Step 20) establishing a model of operation with electric fixed heat (FEL), first determining that the power load required by the system is satisfied by the gas turbine:
PMT=PEL (6)
in the formula, PMTRepresenting the turbine output, in units: kW; pELIs the power load unit required by the whole system: kW. In the electric constant heat mode of operation, the electricity generated by the gas turbine is equal to the electrical load required by the building.
The fuel F consumed by the gas turbine is then determinedMTComprises the following steps:
Figure BDA0002205214100000024
in the formula etaMTIs the efficiency of the gas turbine, FMTFuel, expressed as gas turbine consumed, unit: cubic meter.
Then determining the heat quantity Q recovered by the heat recovery systemrhsComprises the following steps:
Qrhs=(FMT-PMTrhs=FMT(1-ηMTrhs (8)
in the formula etarhsIs the efficiency of the heat recovery system, QrhsRepresents the amount of heat recovered by the heat recovery system, in units: kW.
The heat Q required for the adsorption refrigerator is then determinedracComprises the following steps:
Figure BDA0002205214100000031
in the formula, QracRepresents the heat required by the adsorption refrigerator, unit: kW; qCLIs the cooling load of the system, in units: kW; COPacIs the coefficient of performance of an absorption chiller.
The heat Q required by the heat exchanger is then determinedrhComprises the following steps:
Figure BDA0002205214100000032
in the formula, QrhRepresents the heat required by the heat exchanger, in units: kW; qHLIs the heat load, unit: kW; etaheIs the efficiency of the heat exchanger.
The integrated energy system needs to meet the cold and heat demands of the user side while supplying power. If the heat recovered by the heat recovery system is not sufficient, the system will often meet the remaining portion of the heat demand by an augmented gas boiler.
The equivalent thermal load is then determined:
Qreq=Qrac+Qrh (11)
wherein Q isreqIs the equivalent thermal load, in units: kW.
There are two possibilities:
if Qrhs≥Qreq then Qboiler=0 (12)
if Qrhs<Qreq then Qboiler=Qrac+Qrh-Qrhs (13)
determining the fuel consumed by the gas boiler as follows:
Figure BDA0002205214100000033
the fuel consumed by the entire system is then determined to be:
Fm=FMT+Fboiler (15)
and finally, determining that the electricity purchasing quantity from the system to the large power grid is as follows:
Pgrid=0 (16)
in the formula, PgridThe method is characterized in that the power purchasing amount of the system to a large power grid is expressed in a unit: kW.
Step 30) establishing a heat-to-power (FTL) operation model, and firstly determining the heat quantity Q recovered by the heat recovery systemrhsComprises the following steps:
Qrhs=Qreq=Qrac+Qrh (17)
since the waste heat recovered by the system is known, the fuel F consumed by the gas turbineMTComprises the following steps:
Figure BDA0002205214100000041
in the operation mode of the fixed-power by heat, the heat required by the system is satisfied by the waste heat recovered from the gas turbine in each operation period.
Determining the output P of a gas turbineMTComprises the following steps:
PMT=FMTηMT (19)
under this operating strategy, the electricity generated by the gas turbine is unlikely to meet the system requirements, and the shortfall is partially met by buying electricity from the large grid. There are two cases:
if PMT<PEL then Pgrid=PEL-PMT (20)
if PMT≥PEL then Pexcess=PMT-PEL (21)
in the formula, PexcessIs the excess electricity produced by the gas turbine, in units: kW; may be sold to the grid or stored for future use.
Finally, determining the consumption of the system fuel as follows:
Fm=FMT (22)
step 40) from the analysis of step 20), step 30) it can be seen that when the power demand of the system and the heat demand of the system are different, the system may need to purchase electricity or supplement the heat demand by adding a gas boiler, so when electricity needs to be purchased or the boiler needs to be supplemented is a critical issue.
Rated heat-electricity ratio of gas turbine as heat supply Q of unitrhsGenerating capacity P of the unitMTThe ratio is calculated as follows:
Figure BDA0002205214100000051
in the formula, epsilonNIndicating the rated heat to power ratio of the gas turbine.
Determining the load demand thermoelectric ratio:
KL=Qreq/PEL (24)
in the formula, KLRepresenting the load demand thermoelectric ratio.
For the operational mode of the FTL, the electricity purchase situation is: pEL-PMT>0 (i.e. P)EL-Qreq·ηMT/(1-ηMTrhs>0) After finishing, Q is obtainedreq/PELN
For the FEL mode of operation, the case where the boiler requires afterburning is: qreq-Qhrs>0 (i.e. Q)req-PEL·(1-ηMT)·ηrhsMT>0) After finishing, Q is obtainedreq/PELN
Step 50) first establishing a Primary Energy Consumption (PEC) evaluation index:
PEC=kfFm+kePgrid (25)
in the formula, PEC represents a primary energy consumption index, kfRepresents the primary energy consumption coefficient, k, of the natural gas of the systemeAnd representing the primary energy consumption coefficient of the power grid.
Build operating COST (COST):
COST=CfFm+CePgrid (26)
wherein COST represents a system running COST index, CfRepresenting the system operating coefficient, CeRepresenting the grid operating cost factor.
Establishing CO2 emission (CDE) indexes:
CDE=μfFmePgrid (27)
in the formula, CDE represents a system CO2An emission index; mu.sfRepresents the integrated system CO2 emission coefficient, μeRepresenting the grid CO2 emission factor.
Establishing a comprehensive evaluation Index (IPC), synthesizing three indexes of a formula (25), a formula (26) and a formula (27), and determining the comprehensive evaluation Index (IPC):
Figure BDA0002205214100000052
in the formula, ω1、ω2、ω3Representing the weight coefficients. Omega123=1,ω1∈(0,1),ω2∈(0,1),ω3E (0, 1); IPC represents a comprehensive evaluation index.
(1)KLNIn this case, in the FTL case, power is not required to be purchased from the grid, but in the FEL case, a gas boiler needs to be additionally generated.
Figure BDA0002205214100000061
From KLNWhen, KLN>1 is always true, if the FEL mode of operation is better than the FTL, i.e. Δ E>0, need ηboilerrAnd the heating efficiency of the boiler is generally higher than that of the integrated energy system power generation equipment in actual operation, so the operation mode of the FEL is better than that of the FTL in the condition.
(2)KLNAt this time, power is required to be purchased from the power grid in the FTL operation mode, but the power-increasing gas boiler is not required in the FEL mode.
Figure BDA0002205214100000062
When K isLNWhen, KLNIs always less than 1, so in this case if the FEL mode of operation is better than the FTL, i.e. Δ E>0, need (k)m-kgridηMT) I.e. when km/kgridMTIn the meantime, the FEL operation mode is better than the FTL operation mode, whereas the FTL operation mode is better than the FEL.
From the above analysis it was concluded that:
when KLNWhen or when KLNAnd k ism/kgridMTThe FEL operation mode is superior to the FTL operation mode;
when KLNAnd k ism/kgridMTIn time, the FTL mode of operation is superior to the FEL mode of operation.
When the evaluation index is determined in actual operation, a better operation mode of the FEL and the FTL is selected after the judgment, and the system operates in the operation mode, so that the performance of the system is optimal.
Unlike the above single evaluation index, when a plurality of indexes are combined, the combined evaluation index thereof can be written as:
Figure BDA0002205214100000071
derived from equation (31):
Figure BDA0002205214100000072
wherein the content of the first and second substances,
Figure BDA0002205214100000073
let KF=k11+k12+k13,KE=k12+k22+k32If the above formula is converted to IPC ═ KF·Fm+KE·Egrid. It can be seen that similar to the form derived above, similar conclusions can then be drawn from the above analysis:
when KLNWhen or when KLNAnd K isF/KEMTThe FEL operation mode is superior to the FTL operation mode;
when KLNAnd K isF/KEMTIn time, the FTL mode of operation is superior to the FEL mode of operation.
Further, in step 40), the following are assumed:
1) the utilization rate of the waste heat generated by the comprehensive energy system is one hundred percent;
2) the efficiency of each device in actual operation is kept unchanged, and the device is operated without failure in an optimization period;
3) the power output by each device is continuous;
4) in each operation mode, the system can only purchase electricity from the power grid and cannot sell electricity.
Further, in step 20), step 30) and step 40), the equipment parameter is selected to be etaMT=0.25,COPac=0.7,ηrhs=0.8,ηboiler=0.8,COPec=03.0,ηhe=0.8,μe=968,μf=220,ke=3.336,kf=1.047,Cf=0.19,Ce=0.93。
Further, in step 50), it is default that the smaller the value of the evaluation index IPC, the better, but the larger the value of the evaluation index, the better, the opposite is concluded.
Compared with the prior art, the invention has the following advantages:
the invention provides an evaluation criterion-based operation method of a comprehensive energy system, which is based on two modes of 'electricity by heat' and 'heat by electricity', selects a better operation mode of 'electricity by heat' and 'heat by electricity' after judgment, and enables the system to operate in the operation mode so as to optimize the performance of the system. The 'electricity ordering by heat' is based on heat demand, the power output of the system is given, and the insufficient power part is satisfied by purchasing power from a large power grid; "electric constant heat" is the heat output given to the system based on the electricity demand, and the insufficient heat is supplied by the after-burning boiler. Firstly, a separate supply system model is establishedCalculating the electricity purchasing quantity and the fuel consumption quantity of the distribution system; then establishing a comprehensive energy system model with two modes of 'fixing power with heat' and 'fixing heat with electricity', and respectively calculating the electricity purchasing amount and the fuel consumption amount of the comprehensive energy system; three evaluation indexes of the system are established: primary energy consumption, operating costs, CO2The discharge capacity is reduced, the comprehensive evaluation index is established by integrating the three indexes, and the utilization rate of system energy is optimized by calculating the application conditions of two main operation modes of the cooling, heating and power system, so that the system is more economical and energy-saving than the operation modes of electricity for heat determination and electricity for heat determination. The problem that due to the randomness of the load and the partial load characteristics of the unit, the load and the unit are often inconsistent, and when a certain operation mode is adopted, resource waste is often caused is well solved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a working principle diagram of the comprehensive energy system containing combined cooling, heating and power supply of the invention.
Fig. 2 is a schematic diagram of the operation of the dispensing system of the present invention.
FIG. 3 is a thermoelectric ratio curve of the gas turbine of the present invention.
FIG. 4 is a flow chart of the evaluation criteria based operation strategy of the present invention
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in detail with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
The invention relates to an evaluation criterion-based operation method of a comprehensive energy system, belonging to the field of comprehensive energy systems. The working principle is that after fuel is combusted in a combustion chamber, gas passes through a gas turbine to drive the gas turbine to work, high-temperature flue gas after acting is recovered by a heat recovery system, enters an adsorption refrigerator for refrigeration in summer and enters a heat exchanger for heat supply in winter, insufficient heat (the cold and hot load is converted into the heat demand) is supplemented by a gas boiler, and insufficient power load is met by purchasing electricity from a large power grid. The working principle diagram is shown in figure 1. For a conventional split-supply system, the cold load demand is satisfied by an electric chiller, the heat generated by a gas boiler is passed through a heat exchanger to satisfy the system heat demand, and the system's electrical load and the electrical power required by the electric chiller are purchased from a large power grid, the schematic of which is shown in fig. 2.
The graph of the heat-to-power ratio of a gas turbine is shown in fig. 3, ideally with the heat-to-power demand of the system just on the line, where the energy of the system is fully utilized, but it is clear that this is not practical for the system to be at a slope k ═ epsilonNAt the above point, if the system adopts the FTL operation mode, the system will generate redundant electric power; if the FEL operation mode is adopted, the heat recovered by the system cannot meet the heat requirement and needs to be met by an additional gas boiler, such as a point A in the figure. And for a position on the slope k ═ epsilonNIf the system adopts the FTL operation method, the power generated by the system cannot meet the needs of the user side, and therefore, the power needs to be purchased from the large power grid, whereas if the FEL is adopted, the system generates excessive heat, as shown in point B in the figure.
As shown in fig. 4, the method for operating a combined cooling heating and power micro grid based on model predictive control according to the present invention includes the following steps:
step 10) establishing a model of a sub-supply System (SP), firstly determining the power consumption of the electric refrigerator as follows:
Figure BDA0002205214100000091
in the formula, COPecRepresenting the coefficient of performance, P, of the electric refrigeratorecIndicating electric refrigerator inputPower, unit: kW; qCLRepresents the system cooling load, unit: kW.
The fuel consumed by the boiler is then determined as:
Figure BDA0002205214100000092
in the formula, FboilerNatural gas, which represents the consumption of a gas boiler, unit: cubic meter; qrhRepresents heat exchanger input heat, in units: kW; η SP builder represents the efficiency of the gas boiler.
Then, determining that the power purchased by the power grid is as follows:
Pgrid=PEL+Pec (3)
in the formula, PgridRepresenting the amount of electricity purchased from the grid, in units: kW; pELIs the electrical load of the system, unit: kW.
The fuel consumed by the system is then determined, only the gas boiler consuming:
Fm=Fboiler (4)
in the formula, FmFuel, expressed as system consumption, unit: cubic meter.
Converting the purchase electricity into corresponding fuel consumption:
Figure BDA0002205214100000093
in the formula, EmRepresenting the conversion of the purchase electricity into the corresponding fuel consumption, unit: cubic meter.
Step 20) establishing a model of operation with electric fixed heat (FEL), first determining that the power load required by the system is satisfied by the gas turbine:
PMT=PEL (6)
in the formula, PMTRepresenting the turbine output, in units: kW; pELIs the power load unit required by the whole system: kW. In the electric constant-heat operation mode, the electricity generated by the gas turbine is equal to the electricity generated by the generatorThe electrical load required by the building.
The fuel F consumed by the gas turbine is then determinedMTComprises the following steps:
Figure BDA0002205214100000101
in the formula etaMTIs the efficiency of the gas turbine, FMTFuel, expressed as gas turbine consumed, unit: cubic meter.
Then determining the heat quantity Q recovered by the heat recovery systemrhsComprises the following steps:
Qrhs=(FMT-PMTrhs=FMT(1-ηMTrhs (8)
in the formula etarhsIs the efficiency of the heat recovery system, QrhsRepresents the amount of heat recovered by the heat recovery system, in units: kW.
The heat Q required for the adsorption refrigerator is then determinedracComprises the following steps:
Figure BDA0002205214100000102
in the formula, QracRepresents the heat required by the adsorption refrigerator, unit: kW; qCLIs the cooling load of the system, in units: kW; COPacIs the coefficient of performance of an absorption chiller.
The heat Q required by the heat exchanger is then determinedrhComprises the following steps:
Figure BDA0002205214100000103
in the formula, QrhRepresents the heat required by the heat exchanger, in units: kW; qHLIs the heat load, unit: kW; etaheIs the efficiency of the heat exchanger.
The integrated energy system needs to meet the cold and heat demands of the user side while supplying power. If the heat recovered by the heat recovery system is not sufficient, the system will often meet the remaining portion of the heat demand by an augmented gas boiler.
The equivalent thermal load is then determined:
Qreq=Qrac+Qrh (11)
wherein Q isreqIs the equivalent thermal load, in units: kW.
There are two possibilities:
if Qrhs≥Qreq then Qboiler=0 (12)
if Qrhs<Qreq then Qboiler=Qrac+Qrh-Qrhs(13)
determining the fuel consumed by the gas boiler as follows:
Figure BDA0002205214100000111
the fuel consumed by the entire system is then determined to be:
Fm=FMT+Fboiler (15)
and finally, determining that the electricity purchasing quantity from the system to the large power grid is as follows:
Pgrid=0 (16)
in the formula, PgridThe method is characterized in that the power purchasing amount of the system to a large power grid is expressed in a unit: kW.
Step 30) establishing a heat-to-power (FTL) operation model, and firstly determining the heat quantity Q recovered by the heat recovery systemrhsComprises the following steps:
Qrhs=Qreq=Qrac+Qrh (17)
since the waste heat recovered by the system is known, the fuel F consumed by the gas turbineMTComprises the following steps:
Figure BDA0002205214100000112
in the operation mode of the fixed-power by heat, the heat required by the system is satisfied by the waste heat recovered from the gas turbine in each operation period.
Determining the output P of a gas turbineMTComprises the following steps:
PMT=FMTηMT (19)
under this operating strategy, the electricity generated by the gas turbine is unlikely to meet the system requirements, and the shortfall is partially met by buying electricity from the large grid. There are two cases:
if PMT<PEL then Pgrid=PEL-PMT (20)
if PMT≥PEL then Pexcess=PMT-PEL (21)
in the formula, PexcessIs the excess electricity produced by the gas turbine, in units: kW; may be sold to the grid or stored for future use.
Finally, determining the consumption of the system fuel as follows:
Fm=FMT (22)
step 40) from the analysis of step 20), step 30) it can be seen that when the power demand of the system and the heat demand of the system are different, the system may need to purchase electricity or supplement the heat demand by adding a gas boiler, so when electricity needs to be purchased or the boiler needs to be supplemented is a critical issue.
Rated heat-electricity ratio of gas turbine as heat supply Q of unitrhsGenerating capacity P of the unitMTThe ratio is calculated as follows:
Figure BDA0002205214100000121
in the formula, epsilonNIndicating the rated heat to power ratio of the gas turbine.
Determining the load demand thermoelectric ratio:
KL=Qreq/PEL (24)
in the formula, KLRepresenting the load demand thermoelectric ratio.
For the operational mode of the FTL, the electricity purchase situation is: pEL-PMT>0 (i.e. P)EL-Qreq·ηMT/(1-ηMTrhs>0) After finishing, Q is obtainedreq/PELN
For the FEL mode of operation, the case where the boiler requires afterburning is: qreq-Qhrs>0 (i.e. Q)req-PEL·(1-ηMT)·ηrhsMT>0) After finishing, Q is obtainedreq/PELN
Step 50) first establishing a Primary Energy Consumption (PEC) evaluation index:
PEC=kfFm+kePgrid (25)
in the formula, PEC represents a primary energy consumption index, kfRepresents the primary energy consumption coefficient, k, of the natural gas of the systemeAnd representing the primary energy consumption coefficient of the power grid.
Build operating COST (COST):
COST=CfFm+CePgrid (26)
wherein COST represents a system running COST index, CfRepresenting the system operating coefficient, CeRepresenting the grid operating cost factor.
Establishing CO2 emission (CDE) indexes:
CDE=μfFmePgrid (27)
in the formula, CDE represents a system CO2An emission index; mu.sfRepresents the integrated system CO2 emission coefficient, μeRepresenting the grid CO2 emission factor.
Establishing a comprehensive evaluation Index (IPC), synthesizing three indexes of a formula (25), a formula (26) and a formula (27), and determining the comprehensive evaluation Index (IPC):
Figure BDA0002205214100000131
in the formula (I), the compound is shown in the specification,ω1、ω2、ω3representing the weight coefficients. Omega123=1,ω1∈(0,1),ω2∈(0,1),ω3E (0, 1); IPC represents a comprehensive evaluation index.
(1)KLNIn this case, in the FTL case, power is not required to be purchased from the grid, but in the FEL case, a gas boiler needs to be additionally generated.
Figure BDA0002205214100000132
From KLNWhen, KLN>1 is always true, if the FEL mode of operation is better than the FTL, i.e. Δ E>0, need ηboilerrAnd the heating efficiency of the boiler is generally higher than that of the integrated energy system power generation equipment in actual operation, so the operation mode of the FEL is better than that of the FTL in the condition.
(2)KLNAt this time, power is required to be purchased from the power grid in the FTL operation mode, but the power-increasing gas boiler is not required in the FEL mode.
Figure BDA0002205214100000141
When K isLNWhen, KLNIs always less than 1, so in this case if the FEL mode of operation is better than the FTL, i.e. Δ E>0, need (k)m-kgridηMT) I.e. when km/kgridMTIn the meantime, the FEL operation mode is better than the FTL operation mode, whereas the FTL operation mode is better than the FEL.
From the above analysis it was concluded that:
when KLNWhen or when KLNAnd k ism/kgridMTThe FEL operation mode is superior to the FTL operation mode;
when KLNAnd k ism/kgridMTIn time, the FTL mode of operation is superior to the FEL mode of operation.
When the evaluation index is determined in actual operation, a better operation mode of the FEL and the FTL is selected after the judgment, and the system operates in the operation mode, so that the performance of the system is optimal.
Unlike the above single evaluation index, when a plurality of indexes are combined, the combined evaluation index thereof can be written as:
Figure BDA0002205214100000142
derived from equation (31):
Figure BDA0002205214100000143
wherein the content of the first and second substances,
Figure BDA0002205214100000144
let KF=k11+k12+k13,KE=k12+k22+k32If the above formula is converted to IPC ═ KF·Fm+KE·Egrid. It can be seen that similar to the form derived above, similar conclusions can then be drawn from the above analysis:
when KLNWhen or when KLNAnd K isF/KEMTThe FEL operation mode is superior to the FTL operation mode;
when KLNAnd K isF/KEMTIn time, the FTL mode of operation is superior to the FEL mode of operation.
To sum up, the method firstly establishes a distribution system model and calculates the electricity purchasing quantity and the fuel consumption quantity of the distribution system; then, the comprehensive energy of 'fixing the electricity by heat' and 'fixing the heat by electricity' is establishedThe source system model is used for respectively calculating the electricity purchasing amount and the fuel consumption amount of the comprehensive energy system; three evaluation indexes of the system are established: primary energy consumption, operating costs, CO2The method comprises the steps of reducing the discharge capacity, integrating three indexes to establish a comprehensive evaluation index, and determining the operation mode of a comprehensive energy system by calculating the application conditions of two main operation modes of a cooling, heating and power system, so that the utilization rate of system energy is optimal, and the method is more economical and energy-saving than the operation modes of electricity for heat determination and electricity for heat determination. The problem of resource waste caused by the fact that a certain operation mode is adopted due to the fact that the randomness of the load is inconsistent with the partial load characteristics of the unit is well solved.
The foregoing is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, various modifications and decorations can be made without departing from the principle of the present invention, and these modifications and decorations should also be regarded as the protection scope of the present invention.

Claims (6)

1. An integrated energy system operation method based on evaluation criteria is characterized in that: firstly, establishing a distribution system model, and calculating the electricity purchasing quantity and the fuel consumption quantity of a distribution system;
then establishing a comprehensive energy system model with two modes of 'fixing power by heat FTL' and 'fixing power by heat FEL', and respectively calculating the electricity purchasing amount and the fuel consumption of the comprehensive energy system;
three evaluation indexes of the system are established: primary energy consumption, operating costs, CO2Reducing the discharge capacity, integrating the three indexes to establish a comprehensive evaluation index, and determining the operation mode of the comprehensive energy system by calculating the application conditions of the two operation modes of the cold and hot electric system so as to optimize the utilization rate of system energy;
establishing a sub-supply system model, firstly determining the power consumption of the electric refrigerator as follows:
Figure RE-160275DEST_PATH_IMAGE001
(1)
in the formula (I), the compound is shown in the specification,COP ec indicating electric refrigerationThe coefficient of performance of the machine is,P ec represents the input power of the electric refrigerator, unit: kW;Q CL represents the system cooling load, unit: kW;
the fuel consumed by the boiler is then determined as:
Figure RE-890464DEST_PATH_IMAGE002
(2)
in the formula (I), the compound is shown in the specification,F boiler natural gas, which represents the consumption of a gas boiler, unit: cubic meter;Q rh represents heat exchanger input heat, in units: kW;η boiler representing the efficiency of the gas boiler;
then, determining that the power purchased by the power grid is as follows:
Figure RE-670202DEST_PATH_IMAGE003
(3)
in the formula (I), the compound is shown in the specification,P grid representing the amount of electricity purchased from the grid, in units: kW;P EL is the electrical load of the system, unit: kW;
the fuel consumed by the system is then determined, only the gas boiler consuming:
Figure RE-127728DEST_PATH_IMAGE004
(4)
in the formula, FmFuel, expressed as system consumption, unit: cubic meter;
converting the purchase electricity into corresponding fuel consumption:
Figure RE-253684DEST_PATH_IMAGE005
(5)
in the formula (I), the compound is shown in the specification,E m representing the conversion of the purchase electricity into the corresponding fuel consumption, unit: cubic meter;
establishing an electric heating operation model, firstly determining that the power load required by the system is met by the gas turbine:
Figure RE-325415DEST_PATH_IMAGE006
(6)
in the formula (I), the compound is shown in the specification,P MT representing the turbine output, in units: kW;P EL is the power load unit required by the whole system: kW; in the electric constant-heat operation mode, the electricity generated by the gas turbine is equal to the electric load required by the building;
the fuel consumed by the gas turbine is then determinedF MT Comprises the following steps:
Figure RE-592448DEST_PATH_IMAGE007
(7)
in the formula (I), the compound is shown in the specification,η MT is the efficiency of the gas turbine engine and,F MT fuel, expressed as gas turbine consumed, unit: cubic meter;
the heat recovered by the heat recovery system is then determinedQ rhs Comprises the following steps:
Figure RE-666714DEST_PATH_IMAGE008
(8)
in the formula (I), the compound is shown in the specification,η rhs is the efficiency of the heat recovery system and,Q rhs represents the amount of heat recovered by the heat recovery system, in units: kW;
the heat required by the adsorption chiller is then determinedQ rac Comprises the following steps:
Figure RE-397910DEST_PATH_IMAGE009
(9)
in the formula (I), the compound is shown in the specification,Q rac representing the amount of heat required by the adsorption chiller,unit: kW;Q CL is the cooling load of the system, in units: kW;COP ac is the coefficient of performance of the absorption chiller;
the heat required by the heat exchanger is then determinedQ rh Comprises the following steps:
Figure RE-771034DEST_PATH_IMAGE010
(10)
in the formula (I), the compound is shown in the specification,Q rh represents the heat required by the heat exchanger, in units: kW;Q HL is the heat load, unit: kW;η he is the efficiency of the heat exchanger;
the comprehensive energy system needs to meet the cold demand and the heat demand of a user side while supplying power; if the heat recovered by the heat recovery system cannot meet the requirement, the system often meets the heat requirement of the rest part by an additional gas boiler;
the equivalent thermal load is then determined:
Figure RE-525364DEST_PATH_IMAGE011
(11)
wherein the content of the first and second substances,Q req is the equivalent thermal load, in units: kW;
establishing a primary energy consumption evaluation index:
Figure RE-75425DEST_PATH_IMAGE012
(25)
in the formula, PEC represents a primary energy consumption index,k f the coefficient of primary energy consumption of the natural gas of the system is shown,k e representing the primary energy consumption coefficient of the grid, FmWhich is indicative of the fuel consumed by the system,P grid representing the purchase of electricity from the grid;
establishing operation cost:
Figure RE-661127DEST_PATH_IMAGE013
(26)
wherein COST represents a system running COST index,C f the coefficient of operation of the system is represented,C e representing a grid operating cost coefficient;
establishing CO2 emission indexes:
Figure RE-464873DEST_PATH_IMAGE014
(27)
in the formula, CDE represents a system CO2An emission index; u. offRepresents the integrated system CO2 emission coefficient, ueRepresents the grid CO2 emission coefficient;
establishing a comprehensive evaluation index, synthesizing three indexes of a formula (25), a formula (26) and a formula (27), and determining the comprehensive evaluation index IPC:
Figure RE-519548DEST_PATH_IMAGE015
(28)
in the formula, ω1、ω2、ω3Representing the weight coefficient, ω123=1,ω1∈(0,1),ω2∈(0,1),ω3E (0, 1); IPC denotes the evaluation index, PECIESRepresenting a primary energy consumption index, PEC, of an integrated energy systemSpIndicating a primary energy consumption index, COST, of the separate supply systemIESRepresenting the COST of operation, COST, of the integrated energy systemSpIndicating the operating cost index, CDE, of the separate supply systemIESCO representing an integrated energy system2Emission index, CDESPCO indicating separate supply system2An emission index;
(1) when in useK L >ε N At this time, in the case of FTL, the power is not required to be purchased from the power grid, but in the case of FEL, a gas boiler is required to be additionally generated;
Figure RE-794671DEST_PATH_IMAGE016
(29)
in the formula, FmWhich is indicative of the fuel consumed by the system,E m representing the conversion of the purchase electricity into the corresponding fuel consumption,P MT the output of the turbine is indicated,P EL is the unit of electrical load required by the entire system,η MT is the efficiency of the gas turbine engine and,F MT which represents the fuel consumed by the gas turbine,Q req is the equivalent thermal load of the heat exchanger,F boiler which represents the natural gas consumed by the gas boiler,η boiler which represents the efficiency of the gas-fired boiler,η r is the efficiency of the heat recovery system and,K L which is indicative of the heat-to-electricity ratio demanded by the load,ε N representing the rated heat to power ratio of the gas turbine;
byK L >ε N When the temperature of the water is higher than the set temperature,K L /ε N >1 is always true, if the FEL mode of operation is better than the FTL, ΔE>0, need toη boiler >η r In actual operation, the heating efficiency of the boiler is higher than that of the power generation equipment of the comprehensive energy system, so the operation mode of the FEL is better than that of the FTL in the condition;
(2)K L <ε N at the moment, electricity needs to be purchased from a power grid in the FTL operation mode, but an additional gas boiler is not needed in the FEL mode;
Figure RE-500459DEST_PATH_IMAGE017
(30)
in the formula, FmWhich is indicative of the fuel consumed by the system,E m representing the conversion of the purchase electricity into the corresponding fuel consumption,P MT the output of the turbine is indicated,P EL is the unit of electrical load required by the entire system,η MT is the efficiency of the gas turbine engine and,F MT which represents the fuel consumed by the gas turbine,Q req is the equivalent thermal load of the heat exchanger,F boiler which represents the natural gas consumed by the gas boiler,η boile which represents the efficiency of the gas-fired boiler,η r is the efficiency of the heat recovery system and,K L which is indicative of the heat-to-electricity ratio demanded by the load,ε N representing the rated heat to power ratio of the gas turbine;
when in useK L <ε N When, KL/ε N Is always less than 1, so if the FEL mode of operation is better than the FTL in this case, ΔE>0, whenk m /k grid <η MT In the meantime, the FEL operation mode is superior to the FTL operation mode, otherwise, the FTL operation mode is superior to the FEL;
from the above analysis it was concluded that:
firstly, whenK L >ε N When or whenK L <ε N And isk m /k grid <η MT The FEL operation mode is superior to the FTL operation mode;
② whenK L <ε N And isk m /k grid >η MT Meanwhile, the FTL operation mode is superior to the FEL operation mode;
when the evaluation index is determined in actual operation, selecting a more optimal operation mode of the FEL and the FTL after the judgment, and enabling the system to operate in the operation mode so as to enable the performance of the system to be optimal;
when a plurality of indexes are combined, the comprehensive evaluation index IPC is written as:
Figure RE-147210DEST_PATH_IMAGE018
(31)
in the formula, ω1、ω2、ω3Representing the weight coefficient, ω123=1,ω1∈(0,1),ω2∈(0,1),ω3E (0, 1); IPC denotes the evaluation index, PECIESRepresenting a primary energy consumption index, PEC, of an integrated energy systemSpIndicating a primary energy consumption index, COST, of the separate supply systemIESRepresenting the COST of operation, COST, of the integrated energy systemSpIndicating the operating cost index, CDE, of the separate supply systemIESCO representing an integrated energy system2Emission index, CDESPCO indicating separate supply system2An emission index;
derived from equation (31):
Figure RE-DEST_PATH_IMAGE019
(32)
Figure RE-938448DEST_PATH_IMAGE020
order toK F =k 11 +k 12 +k 13 K E =k 12 +k 22 +k 32 Then the above formula can be converted intoIPC=K F ∙F m +K E ∙E grid (ii) a The following conclusions can be drawn from the above analysis:
when in useK L >ε N When or whenK L <ε N And isK F /K E <η MT The FEL operation mode is superior to the FTL operation mode;
when in useK L <ε N And isK F /K E >η MT Meanwhile, the FTL operation mode is superior to the FEL operation mode;
the default is that the smaller the value of the evaluation index IPC, the better, but the larger the value of the evaluation index, the better, the opposite is concluded.
2. The method of claim 1, wherein the method comprises: there are two possibilities for equivalent thermal loading:
Figure RE-DEST_PATH_IMAGE021
determining the fuel consumed by the gas boiler as follows:
Figure RE-266530DEST_PATH_IMAGE022
(14)
the fuel consumed by the entire system is then determined to be:
Figure RE-764508DEST_PATH_IMAGE023
(15)
and finally, determining that the electricity purchasing quantity from the system to the large power grid is as follows:
Figure RE-411521DEST_PATH_IMAGE024
(16)
in the formula (I), the compound is shown in the specification,P grid the method is characterized in that the power purchasing amount of the system to a large power grid is expressed in a unit: kW.
3. The method of claim 1, wherein the method comprises: establishing a heat-to-power operation model, and firstly determining the heat recovered by the heat recovery systemQ rhs Comprises the following steps:
Figure RE-691322DEST_PATH_IMAGE025
(17)
since the waste heat recovered by the system is known, the fuel consumed by the gas turbineF MT Comprises the following steps:
Figure RE-636144DEST_PATH_IMAGE026
(18)
in the operation mode of electricity by heat, the heat required by the system is satisfied by the waste heat recovered from the gas turbine in each operation period;
determining the output of a gas turbineP MT Comprises the following steps:
Figure RE-988628DEST_PATH_IMAGE027
(19)。
4. the method of claim 3, wherein the method comprises: under the operation strategy, the electricity generated by the gas turbine cannot necessarily meet the system requirement, and the shortage part is met by buying electricity to a large power grid;
there are two cases:
Figure RE-744226DEST_PATH_IMAGE028
in the formula (I), the compound is shown in the specification,P excess is the excess electricity produced by the gas turbine, in units: kW; can be sold to the grid or stored for future use;
finally, determining the consumption of the system fuel as follows:
Figure RE-244477DEST_PATH_IMAGE029
(22)。
5. the method of claim 3, wherein the method comprises:
rated heat-electricity ratio of gas turbine as heat supply of unitQ rhs Generating capacity of the unitP MT The ratio is calculated as follows:
Figure RE-665094DEST_PATH_IMAGE030
(23)
wherein, the rated thermoelectric ratio of the gas turbine is expressed;
determining the load demand thermoelectric ratio:
K L =Q req /P EL (24)
in the formula (I), the compound is shown in the specification,K L representing the load demand thermoelectric ratio.
6. The method of claim 5, wherein the method comprises:
for the operational mode of the FTL, the electricity purchase situation is:P EL -P MT >0 (i.e. theP EL -Q req η MT /(1-η MT )η rhs >0) Is finished to obtainQ req /P EL <ε N
For the FEL mode of operation, the case where the boiler requires afterburning is:Q req -Q rhs >0 (i.e. theQ req -P EL ∙(1-η MT )∙η rhs /η MT >0) Is finished to obtainQ req /P EL >ε N
CN201910883976.0A 2019-09-18 2019-09-18 Comprehensive energy system operation method based on evaluation criterion Active CN110474335B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910883976.0A CN110474335B (en) 2019-09-18 2019-09-18 Comprehensive energy system operation method based on evaluation criterion

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910883976.0A CN110474335B (en) 2019-09-18 2019-09-18 Comprehensive energy system operation method based on evaluation criterion

Publications (2)

Publication Number Publication Date
CN110474335A CN110474335A (en) 2019-11-19
CN110474335B true CN110474335B (en) 2021-12-07

Family

ID=68516275

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910883976.0A Active CN110474335B (en) 2019-09-18 2019-09-18 Comprehensive energy system operation method based on evaluation criterion

Country Status (1)

Country Link
CN (1) CN110474335B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111520841A (en) * 2020-03-30 2020-08-11 国网天津市电力公司电力科学研究院 Cooling, heating and power combined supply system regulation and control strategy based on efficient low-carbon emission criterion
CN114035434B (en) * 2021-11-22 2023-09-01 西南石油大学 Operation optimization method of gas-steam combined cycle power generation system
CN114779750B (en) * 2022-06-21 2022-11-11 成都安美科燃气技术股份有限公司 Control method for determining power by heat for gas internal combustion engine combined heat and power unit

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102236349A (en) * 2010-04-30 2011-11-09 新奥科技发展有限公司 System energy efficiency controller, energy efficiency grain device and intelligent energy service system for energy utilization
CN106503914A (en) * 2016-11-04 2017-03-15 国电南瑞科技股份有限公司 A kind of polynary economy of energy method for optimizing scheduling towards cooling heating and power generation system
CN108052722B (en) * 2017-12-08 2021-06-08 中国电力工程顾问集团西南电力设计院有限公司 Distributed cooling, heating and power hybrid energy system design method oriented to comprehensive energy efficiency optimization
CN108717594B (en) * 2018-04-16 2020-02-18 东南大学 Economic optimization scheduling method for combined cooling heating and power type multi-microgrid system
CN109255560B (en) * 2018-11-20 2022-10-28 成都大学 CCHP system evaluation optimization method based on cooling, heating and power load proportion
CN109946977A (en) * 2019-04-24 2019-06-28 华北电力大学 Wisdom garden energy source optimization dispatching method containing cold, heat and power triple supply system

Also Published As

Publication number Publication date
CN110474335A (en) 2019-11-19

Similar Documents

Publication Publication Date Title
CN104716644B (en) Renewable energy source cooling, heating and power microgrid system and control method
CN108717594A (en) A kind of more micro-grid system economic optimization dispatching methods of supply of cooling, heating and electrical powers type
CN110474335B (en) Comprehensive energy system operation method based on evaluation criterion
CN112736939B (en) Optimized capacity configuration method for hydrogen production and storage device of hydrogen-doped natural gas comprehensive energy system
CN109696891B (en) Micro energy network system comprising air source heat pump and energy storage and operation control method thereof
CN110503241B (en) Multi-objective optimization method of cold-heat-electricity comprehensive energy system
Feng et al. Analysis of energy-matching performance and suitable users of conventional CCHP systems coupled with different energy storage systems
CN110676847A (en) Optimal scheduling method considering wind power-heat storage unit-electric boiler combined operation
CN110361969B (en) Optimized operation method of cooling, heating and power comprehensive energy system
CN112329988A (en) Demand side response calculation method based on park adjustability analysis
CN106593771A (en) Power-heat decoupling system in heat-power cogeneration system and control method thereof
CN105305472B (en) A kind of substation capacity optimization method based on multiple-energy-source collaboration energy supply
CN110555618A (en) Networked comprehensive energy system optimization scheduling method based on improved goblet sea squirt algorithm
Wang et al. A hybrid operating strategy of combined cooling, heating and power system for multiple demands considering domestic hot water preferentially: A case study
CN113313305A (en) Non-cooperative game-based comprehensive energy system optimization scheduling method
CN116167483A (en) Park comprehensive energy system robust scheduling method considering stepped demand response
CN109255487A (en) A kind of integrated energy system optimization method based on normalized matrix model
CN106096248A (en) The selection method of the optimal heating system under varying environment
CN106440902B (en) A kind of heat accumulation regulating system and its adjusting method based on cogeneration
CN111625961A (en) Comprehensive energy system collaborative optimization operation regulation and control method
CN110516863A (en) A kind of more microgrid active distribution system dual blank-holders of supply of cooling, heating and electrical powers type
CN113239550B (en) Energy storage capacity configuration method for combined cooling heating and power system and combined cooling and power system
Zeng et al. The integrated operation of multiple energy carrier used to enhance the reliability of the integrated energy system
CN112529728A (en) Non-cooperative game-based comprehensive energy system optimization scheduling method
CN114881416B (en) Comprehensive energy inertial support method considering flexibility of cogeneration unit

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant