CN116029518A - Carbon emission planning method for park comprehensive energy system considering generalized energy storage - Google Patents

Carbon emission planning method for park comprehensive energy system considering generalized energy storage Download PDF

Info

Publication number
CN116029518A
CN116029518A CN202310052701.9A CN202310052701A CN116029518A CN 116029518 A CN116029518 A CN 116029518A CN 202310052701 A CN202310052701 A CN 202310052701A CN 116029518 A CN116029518 A CN 116029518A
Authority
CN
China
Prior art keywords
period
unit
energy storage
data center
park
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
CN202310052701.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.)
Hohai University HHU
Original Assignee
Hohai University HHU
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 Hohai University HHU filed Critical Hohai University HHU
Priority to CN202310052701.9A priority Critical patent/CN116029518A/en
Publication of CN116029518A publication Critical patent/CN116029518A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a carbon emission planning method of a park comprehensive energy system considering generalized energy storage, which comprises the following steps: collecting historical photovoltaic output data of a park comprehensive energy system and unit investment cost mainly comprising photovoltaic and storage battery equipment; constructing a multi-energy load increase long-term uncertainty and a photovoltaic output short-term uncertainty scene; establishing an energy hub model taking park hub as a core and park generalized energy storage resource equipment as a main body; establishing a generalized energy storage model; and calculating the operational cost and the extension cost of 365 days of planning years of the park, and outputting the optimal planning and operational strategy. According to the invention, the long-term uncertainty of load increase and the short-term uncertainty of photovoltaic output are fully considered, the time transfer characteristic is utilized to realize the multi-energy cross-period high-efficiency utilization, the generalized energy storage resources taking the data center and the building as main bodies in the park are excavated, the new energy consumption is effectively promoted, and the economical efficiency and the environmental protection performance of the comprehensive energy system planning and scheduling scheme are ensured.

Description

Carbon emission planning method for park comprehensive energy system considering generalized energy storage
Technical Field
The invention relates to a carbon emission planning method of a park comprehensive energy system considering generalized energy storage, and belongs to the field of comprehensive energy system planning and extension.
Background
The comprehensive energy system of the park is reasonably planned, and the comprehensive energy system has an important supporting effect on improving the energy utilization structure of the park and promoting the low-carbonization transformation of industry. The comprehensive energy system planning around the park has a plurality of defects, such as the research on the planning level focuses on the energy balance for medium and long periods, ignores the problem of the flexibility balance of short-term operation, and does not fully excavate the flexible scheduling potential of various energy resources of the park. The generalized energy storage resource can fully utilize the time transfer characteristic of energy of each energy utilization main body in the park, can store heat in the valley period, release energy in the peak period, and realize the cross period and high value utilization of various energies.
Disclosure of Invention
The invention provides a park comprehensive energy system planning method utilizing a generalized energy storage effect.
The invention adopts the following technical scheme:
a carbon emission planning method of a park comprehensive energy system considering generalized energy storage comprises the following steps:
step 1: collecting historical photovoltaic output data of a park comprehensive energy system and unit investment cost mainly comprising photovoltaic and storage battery equipment;
step 2: according to known park parameters, constructing a scene of long-term uncertainty of multi-energy load growth and short-term uncertainty of photovoltaic output by taking the sum of extension cost and operation cost of a park comprehensive energy system as an objective function; the known park parameters refer to the operating parameters of each plant of the park, such as the efficiency of the units of CCHP, CHP, etc.
Step 3: establishing an energy hub model taking park hub as a core and park generalized energy storage resource equipment as a main body;
step 4: establishing a generalized energy storage model according to the energy storage of the storage battery, the thermal characteristics of the data center and the thermal characteristics of the building;
step 5: and calculating the operational cost and the extension cost of 365 days of planning years of the park, and outputting the optimal planning and operational strategy.
Specifically, the scenario of the long-term uncertainty of the multi-energy load increase in step 2 specifically includes: generating three kinds of upper layer multipotency load long-term uncertainty scenes according to seasonal differences, namely transitional seasons, summer and winter respectively;
the short-term uncertainty scene of the photovoltaic output specifically comprises: based on the photovoltaic historical output data, a k-means clustering method is utilized to generate four types of short-term uncertainty scenes of the lower photovoltaic output.
Specifically, the objective function in step 2 is:
Figure BDA0004058974060000021
in the formula (A-1), a subscript t represents a scheduling time;
Figure BDA0004058974060000022
capacity expansion for photovoltaic units, +.>
Figure BDA0004058974060000023
Capacity expansion for a battery pack>
Figure BDA0004058974060000024
Capacity expansion for ground source heat pump unit +.>
Figure BDA0004058974060000025
The expansion capacity of the combined cooling heating power unit is provided; />
Figure BDA0004058974060000026
Cost for unit expansion of photovoltaic unit, +.>
Figure BDA0004058974060000027
Cost for unit expansion of storage battery unit, < >>
Figure BDA0004058974060000028
The unit expansion cost of the ground source heat pump unit is +.>
Figure BDA0004058974060000029
The unit extension cost of the combined cooling heating power unit is set; omega shape g Omega as a typical day collection g Weighting for each typical day; omega shape PV Omega for the photovoltaic output scene set PV The method comprises the steps of weighing a photovoltaic output scene; p (P) t E For the upper power grid electricity purchasing quantity of t time period, < >>
Figure BDA00040589740600000218
The upper-level gas online shopping amount is t time period; />
Figure BDA00040589740600000219
For the electricity price of the t period, < > for>
Figure BDA00040589740600000220
Gas price for period t; i E Carbon emission cost per unit electricity purchase amount, I G Carbon emission costs per unit amount of purchased gas;
in the formula (A-1), the first four items are annual extension costs of each extension unit, namely a photovoltaic unit, a storage battery unit, a ground source heat pump unit and a combined cooling, heating and power unit, and the fifth item is annual operation costs, wherein carbon emission costs of electricity purchase and gas purchase are taken into account;
each planning variable constraint is:
Figure BDA00040589740600000210
in the formula (A-2),
Figure BDA00040589740600000211
extend the upper capacity limit for the photovoltaic unit, < >>
Figure BDA00040589740600000212
Extend the upper limit of capacity for the battery unit, +.>
Figure BDA00040589740600000213
Expansion capacity upper limit for ground source heat pump unit, < ->
Figure BDA00040589740600000214
The expansion capacity upper limit of the combined cooling heating and power unit is +.>
Figure BDA00040589740600000215
For the original photovoltaic installation capacity, < > a->
Figure BDA00040589740600000216
For the original battery capacity>
Figure BDA00040589740600000217
The capacity of the original cold-heat-electricity triple supply unit is +.>
Figure BDA0004058974060000031
The capacity of the original ground source heat pump unit; the first four items are investment upper limit constraints of the extension unit, and the last four items are operation upper limit constraints of the extension unit.
Specifically, the energy hub model in step 3 is specifically:
the energy level constraint of the park is as follows:
Figure BDA0004058974060000032
wherein P is t LC For the period t of the cold load demand,
Figure BDA0004058974060000037
industrial cooling capacity of GSHP unit in t period>
Figure BDA0004058974060000038
The cooling capacity of the CCHP unit is supplied for the period t; p (P) t LH For t-period heat load demand, < > and->
Figure BDA0004058974060000039
Heat is supplied to the t-period CCHP unit; />
Figure BDA00040589740600000310
Heat is supplied to the GSHP unit in the t period; />
Figure BDA00040589740600000311
Heat supply of the air boiler in the t period of time, +.>
Figure BDA00040589740600000312
Heat is supplied to the electric boiler in the period t; p (P) t PV For the t period of photovoltaic output, P t dis For t period of battery charge, P t LE For t-period electrical load demand, P t GSHP For t-period industrial GSHP unit power consumption, P t EB For t-period electric boiler power consumption, P t ch For the discharge quantity of the storage battery set in the t period, P t BD For the power consumption of the GSHP unit of the building at the period of t, P t DC Power consumption of the data center is t time periods; />
Figure BDA00040589740600000313
The electric quantity provided for the t-period CCHP unit; p (P) t CCHP For the consumption of the CCHP unit in t period, P t GB And the air consumption of the air boiler is t time periods. Constraint (A-2) is a cold bus, a hot bus, an electric bus and an air bus balance equation from top to bottom in sequence;
the operational constraints of conventional multi-energy conversion equipment in a park are as follows:
Figure BDA0004058974060000033
in the method, in the process of the invention,
Figure BDA0004058974060000034
for CCHP unit cold conversion efficiency, P t CCHP,C Natural gas quantity for cold transfer is input into the CCHP unit;
Figure BDA0004058974060000035
for the thermal conversion efficiency of the CCHP unit, P t CCHP,H Natural gas quantity for transferring heat is input into the CCHP unit; />
Figure BDA0004058974060000036
For the electric conversion efficiency of the CCHP unit, P t CCHP,E Natural gas quantity for inputting CCHP units for power conversion; p (P) t TF To input the electric quantity of the transformer, eta TF The electric conversion efficiency of the transformer; η (eta) EB For the heat conversion efficiency of the electric boiler, eta GB The heat conversion efficiency of the gas boiler is; the first three terms are CCHP unit operation equations, and the last three terms are respectively transformer, electric boiler and gas boiler operation equations;
equipment operating constraints related to campus generalized energy storage resources include storage batteries, ground source heat pumps, and data centers, where:
the operation constraint of the storage battery unit is as follows:
Figure BDA0004058974060000041
wherein P is t ch For t period of energy storage charging power, P t dis Energy storage and discharge power is t time period;
Figure BDA0004058974060000043
for the t period of the energy-storing charge state variable, +.>
Figure BDA0004058974060000044
The energy storage and discharge state variable is t time period; e (E) t Energy storage electric quantity is carried out for a period t; lambda (lambda) E For t period of energy storage self-discharge coefficient, eta ch For the energy storage and charging efficiency of t period, eta dis Energy storage and discharge efficiency is t time period; e (E) min For storing energy, E is the lower limit of the electric quantity max Is the upper limit of the energy storage electric quantity;
the operation constraint of the ground source heat pump unit is as follows:
Figure BDA0004058974060000042
wherein P is t BD,C Electric energy for converting heat for t-period building GSHP unit, P t BD,H The GSHP unit is used for converting cold electric energy for the building at the time of t;
Figure BDA0004058974060000045
and->
Figure BDA0004058974060000046
The heat energy and the cold energy which are respectively transmitted to the building by GSHP in t time period; />
Figure BDA0004058974060000047
And->
Figure BDA0004058974060000048
Heating and cooling operation state variables are respectively;
the operation constraint of the data center is as follows;
Figure BDA0004058974060000051
wherein P is t IT For a period t data center IT device power consumption,
Figure BDA0004058974060000052
for each k-shaped suit in t time periodServer power; />
Figure BDA0004058974060000053
For a fixed power consumption of a t-period k-class server, < >>
Figure BDA0004058974060000054
CPU power consumption is t time period; c (C) 1 D 'is the power consumption coefficient of the server' k,t The data load quantity required to be processed for the k-type server of the data center in the t period; a, a k,s,t The working frequency zone bit is s-gear of a k-type server CPU of the data center in the t period; f (f) k,s And the working frequency is s-gear for the CPU of the k-type server of the data center.
According to the storage battery energy storage, the thermal characteristics of a data center and the thermal characteristics of a building, a generalized energy storage model is established, and the method comprises the following steps:
generalized energy storage capacity of storage battery in t period is E t
The thermal characteristics and generalized energy storage constraints of the data center are:
Figure BDA0004058974060000055
in the method, in the process of the invention,
Figure BDA0004058974060000056
generalized energy storage for data center t period, H DC Specific heat capacity of data center, T t dc,in Temperature in data center T time period room, T t dc,in,max An upper limit of the indoor temperature in the t period; η (eta) c Data center refrigeration unit efficiency, Q t Refrigeration T produced by refrigeration equipment of data center T-section t dc,out For t period of external environment temperature, U C The heat transfer coefficient of the data center to the outside;
Figure BDA0004058974060000057
lower limit of temperature change rate of data center, +.>
Figure BDA0004058974060000058
An upper limit for the temperature change rate of the data center; />
Figure BDA0004058974060000059
Is the indoor temperature lower limit of the data center, +.>
Figure BDA00040589740600000510
An upper temperature limit in the data center room;
the thermal characteristics and generalized energy storage constraints of the park building are as follows:
Figure BDA0004058974060000061
in the method, in the process of the invention,
Figure BDA0004058974060000067
generalized energy storage for data center t period, H BD Specific heat capacity of data center, T t BD Indoor temperature in t time period of data center,/o>
Figure BDA0004058974060000062
For the upper limit of the indoor temperature in the t period, < >>
Figure BDA0004058974060000063
The lower limit of the indoor temperature is t time period; />
Figure BDA0004058974060000064
Phi for the solar radiation quantity passing through the glazing of a building t The solar radiation heat is in a t period, and sigma is the solar radiation loss rate; />
Figure BDA0004058974060000068
For building to radiate heat outwards, T t E Building wall temperature for the park at time t; t (T) t The indoor temperature of the park at the period t; u (U) A The heat transfer coefficient is the direct heat transfer coefficient of the outer vertical surface of the building; u (U) E Heat transfer coefficients for indoor to wall; t (T) t a Is the external environment temperature in the t period. Wherein T is t E Building for t periodOuter wall temperature->
Figure BDA0004058974060000065
For the temperature of the building outer wall in t-1 time period, H PE Is the specific heat capacity of the building wall, T t a Is the external environment temperature in the t period.
The generalized energy storage total amount of the park is as follows:
Figure BDA0004058974060000066
calculating the operational cost and the extension cost of 365 days of the planning year of the park, and outputting the optimal planning and the operational strategy, wherein the method comprises the following steps:
in the actual solving process, the objective function is the sum of annual investment cost and planning annual operation cost, wherein the annual investment cost is the total extension cost divided by the planning age, and when the planning annual operation cost is calculated, all loads and photovoltaic output conditions during the planning period need to be met by each equipment operation constraint.
The invention has the beneficial effects that:
according to the invention, the long-term uncertainty of load increase and the short-term uncertainty of photovoltaic output are fully considered, the time transfer characteristic is utilized to realize the multi-energy cross-period high-efficiency utilization, the generalized energy storage resources taking the data center and the building as main bodies in the park are excavated, the new energy consumption is effectively promoted, and the economical efficiency and the environmental protection performance of the comprehensive energy system planning and scheduling scheme are ensured.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Figure 2 is a block diagram of a campus energy hub.
Figure 3 is a graph of typical daily energy load conditions for each campus.
Fig. 4 is a diagram of a campus planning year summer generalized energy storage situation that accounts for generalized energy storage.
Fig. 5 is a short term uncertainty scenario for photovoltaic output.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are only some embodiments of the present invention and are not intended to limit the scope of the present invention, and based on the embodiments of the present invention, all other embodiments obtained by a person skilled in the art without making any creative effort belong to the protection scope of the present invention.
As shown in fig. 1, the invention provides a carbon emission planning method of a park comprehensive energy system considering generalized energy storage, and the method for planning the park comprehensive energy system in the embodiment of the invention comprises the following steps:
step 1: and collecting park operation parameters such as historical photovoltaic output data of a park comprehensive energy system and unit investment cost mainly comprising photovoltaic equipment, storage batteries and the like.
Step 2: generating a scene of long-term uncertainty of multi-energy load growth and short-term uncertainty of photovoltaic output, and taking the sum of annual extension cost and operation cost of system planning as an objective function.
Generating a multi-energy load increase long-term uncertainty and a photovoltaic output short-term uncertainty scene:
generating three kinds of upper layer multipotency load long-term uncertainty scenes according to seasonal differences, namely transitional seasons, summer and winter respectively, as shown in figure 3; based on the photovoltaic historical output data, a k-means clustering method is utilized to generate four types of short-term uncertainty scenes of the lower photovoltaic output, as shown in figure 5.
Based on the uncertainty scene, constructing an objective function with the sum of the extension cost and the running cost as an objective function, wherein the objective function is as follows:
Figure BDA0004058974060000071
in the formula, a subscript t represents a scheduling time;
Figure BDA0004058974060000072
capacity expansion for photovoltaic units, +.>
Figure BDA0004058974060000073
The capacity of the accumulator unit is extended,
Figure BDA0004058974060000074
capacity expansion for ground source heat pump unit +.>
Figure BDA0004058974060000075
The expansion capacity of the combined cooling heating power unit is provided; />
Figure BDA0004058974060000076
Cost for unit expansion of photovoltaic unit, +.>
Figure BDA0004058974060000077
Cost for unit expansion of storage battery unit, < >>
Figure BDA0004058974060000078
The unit expansion cost of the ground source heat pump unit is +.>
Figure BDA0004058974060000081
The unit extension cost of the combined cooling heating power unit is set; omega shape g Omega as a typical day collection g Weighting for each typical day; omega shape PV Omega for the photovoltaic output scene set PV The method comprises the steps of weighing a photovoltaic output scene; p (P) t E For the upper power grid electricity purchasing quantity of t time period, < >>
Figure BDA00040589740600000812
The upper-level gas online shopping amount is t time period; />
Figure BDA00040589740600000813
For the electricity price of the t period, < > for>
Figure BDA00040589740600000814
Gas price for period t; i E Carbon emission cost per unit electricity purchase amount, I G Carbon emission costs per unit amount of purchased gas.
In the formula (A-1), the first four items are annual extension costs of each extension unit respectively, and the annual extension costs comprise a photovoltaic unit, a storage battery unit, a ground source heat pump unit and a combined cooling, heating and power unit, and the fifth item is annual operation costs, wherein carbon emission costs of electricity purchase and gas purchase are taken into account.
Each planning variable constraint is:
Figure BDA0004058974060000082
in the formula (A-2),
Figure BDA0004058974060000083
extend the upper capacity limit for the photovoltaic unit, < >>
Figure BDA0004058974060000084
Extend the upper limit of capacity for the battery unit, +.>
Figure BDA0004058974060000085
Extend the upper capacity limit for GSHP units, < > for>
Figure BDA0004058974060000086
And (5) expanding the upper limit of capacity for the CCHP unit. />
Figure BDA0004058974060000087
For the original photovoltaic installation capacity, < > a->
Figure BDA0004058974060000088
For the original battery capacity>
Figure BDA0004058974060000089
For the original CCHP unit capacity, < > and->
Figure BDA00040589740600000810
Is the original GSHP unit capacity. The first four items are investment upper limit constraints of the extension unit, and the last four items are operation upper limit constraints of the extension unit.
Step 3: building an energy hub model with a park hub as a core:
the energy level constraint of the park is as follows:
Figure BDA00040589740600000811
wherein P is t LC For the period t of the cold load demand,
Figure BDA0004058974060000095
for industrial cold supply of GSHP unit in t period, GSHP unit is ground source heat pump unit, and ∈>
Figure BDA0004058974060000096
The cooling capacity is supplied to a t-period CCHP unit, namely a cold-heat-power cogeneration unit; p (P) t LH For t-period heat load demand, < > and->
Figure BDA0004058974060000097
Heat is supplied to the t-period CCHP unit; />
Figure BDA0004058974060000098
Heat is supplied to the GSHP unit in the t period; />
Figure BDA0004058974060000099
Heat supply of the air boiler in the t period of time, +.>
Figure BDA00040589740600000910
Heat is supplied to the electric boiler in the period t; p (P) t PV For the t period of photovoltaic output, P t dis For t period of battery charge, P t LE For t-period electrical load demand, P t GSHP For t-period industrial GSHP unit power consumption, P t EB For t-period electric boiler power consumption, P t ch For the discharge quantity of the storage battery set in the t period, P t BD For the power consumption of the GSHP unit of the building at the period of t, P t DC Power consumption of the data center is t time periods; />
Figure BDA00040589740600000911
For period tThe electric quantity provided by the CCHP unit; p (P) t CCHP For the consumption of the CCHP unit in t period, P t GB And the air consumption of the air boiler is t time periods. The constraint (A-2) is a cold bus, a hot bus, an electric bus and an air bus balance equation in sequence from top to bottom.
The operational constraints of conventional multi-energy conversion equipment in a park are as follows:
Figure BDA0004058974060000091
in the method, in the process of the invention,
Figure BDA0004058974060000092
for CCHP unit cold conversion efficiency, P t CCHP,C Natural gas quantity for cold transfer is input into the CCHP unit;
Figure BDA0004058974060000093
for the thermal conversion efficiency of the CCHP unit, P t CCHP,H Natural gas quantity for transferring heat is input into the CCHP unit; />
Figure BDA0004058974060000094
For the electric conversion efficiency of the CCHP unit, P t CCHP,E Natural gas quantity for inputting CCHP units for power conversion; p (P) t TF To input the electric quantity of the transformer, eta TF The electric conversion efficiency of the transformer; η (eta) EB For the heat conversion efficiency of the electric boiler, eta GB The heat conversion efficiency of the gas boiler is achieved. The first three terms are CCHP unit operation equations, and the last three terms are respectively transformer, electric boiler and gas boiler operation equations.
Equipment operating constraints related to campus generalized energy storage resources include storage batteries, ground source heat pumps, and data centers, where:
the operation constraint of the storage battery unit is as follows:
Figure BDA0004058974060000101
wherein P is t ch For t period of energy storage charging power, P t dis Energy storage and discharge power is t time period;
Figure BDA0004058974060000107
for the t period of the energy-storing charge state variable, +.>
Figure BDA0004058974060000108
The energy storage and discharge state variable is t time period; e (E) t Energy storage electric quantity is carried out for a period t; lambda (lambda) E For t period of energy storage self-discharge coefficient, eta ch For the energy storage and charging efficiency of t period, eta dis Energy storage and discharge efficiency is t time period; e (E) min For storing energy, E is the lower limit of the electric quantity max Is the upper limit of the stored energy.
The operation constraint of the ground source heat pump unit is as follows:
Figure BDA0004058974060000102
wherein P is t BD,C Electric energy for converting heat for t-period building GSHP unit, P t BD,H The GSHP unit is used for converting cold electric energy for the building at the time of t;
Figure BDA00040589740600001010
and->
Figure BDA0004058974060000109
The heat energy and the cold energy which are respectively transmitted to the building by GSHP in t time period; />
Figure BDA00040589740600001011
And->
Figure BDA00040589740600001012
And the operation state variables are respectively heating and cooling operation state variables. The operating constraints of the industrial GSHP unit are the same as those of the industrial GSHP unit, and therefore the operation constraints are not repeated.
The operation constraint of the data center is as follows;
Figure BDA0004058974060000103
wherein P is t IT For a period t data center IT device power consumption,
Figure BDA0004058974060000104
power for each k-type server for a period t; />
Figure BDA0004058974060000105
For a fixed power consumption of a t-period k-class server, < >>
Figure BDA0004058974060000106
CPU power consumption is t time period; c (C) 1 D 'is the power consumption coefficient of the server' k,t The data load quantity required to be processed for the k-type server of the data center in the t period; a, a k,s,t The working frequency zone bit is s-gear of a k-type server CPU of the data center in the t period; f (f) k,s And the working frequency is s-gear for the CPU of the k-type server of the data center.
Step 4: and building a generalized energy storage model according to the energy storage of the storage battery, the thermal characteristics of the data center and the thermal characteristics of the buildings in the park.
Generalized energy storage capacity of storage battery in t period is E t
The thermal characteristics and generalized energy storage constraints of the data center are:
Figure BDA0004058974060000111
in the method, in the process of the invention,
Figure BDA0004058974060000112
generalized energy storage for data center t period, H DC Specific heat capacity of data center, T t dc,in Temperature in data center T time period room, T t dc,in,max An upper limit of the indoor temperature in the t period; η (eta) c Data center refrigeration unit efficiency, Q t Refrigeration T produced by refrigeration equipment of data center T-section t dc,out At tOutside environmental temperature of the section, U C The heat transfer coefficient of the data center to the outside; />
Figure BDA0004058974060000113
Lower limit of temperature change rate of data center, +.>
Figure BDA0004058974060000114
An upper limit for the temperature change rate of the data center; />
Figure BDA0004058974060000115
Is the indoor temperature lower limit of the data center, +.>
Figure BDA0004058974060000116
An upper temperature limit within the data center room.
The thermal characteristics and generalized energy storage constraints of the park building are as follows:
Figure BDA0004058974060000117
in the method, in the process of the invention,
Figure BDA00040589740600001110
generalized energy storage for data center t period, H BD Specific heat capacity of data center, T t BD Indoor temperature in t time period of data center,/o>
Figure BDA0004058974060000118
For the upper limit of the indoor temperature in the t period, < >>
Figure BDA0004058974060000119
The lower limit of the indoor temperature is t time period; />
Figure BDA0004058974060000121
Phi for the solar radiation quantity passing through the glazing of a building t The solar radiation heat is in a t period, and sigma is the solar radiation loss rate; />
Figure BDA0004058974060000125
For building to radiate heat outwards, T t E Building wall temperature for the park at time t; t (T) t The indoor temperature of the park at the period t; u (U) A The heat transfer coefficient is the direct heat transfer coefficient of the outer vertical surface of the building; u (U) E Heat transfer coefficients for indoor to wall; t (T) t a Is the external environment temperature in the t period. Wherein T is t E For the temperature of the building outer wall in t time period,/for the building outer wall temperature>
Figure BDA0004058974060000122
For the temperature of the building outer wall in t-1 time period, H PE Is the specific heat capacity of the building wall, T t a Is the external environment temperature in the t period.
The generalized energy storage total amount of the park is as follows:
Figure BDA0004058974060000126
step 5: and calculating the operational cost and the extension cost of 365 days of planning years of the park, and outputting the optimal planning and operational strategy.
In the actual solving process, the objective function is the sum of annual investment cost and planning annual operation cost, wherein the annual investment cost is the total extension cost divided by the planning age, and when the planning annual operation cost is calculated, all loads and photovoltaic output conditions during the planning period need to be met by each equipment operation constraint.
Example analysis
And selecting each typical daily energy load of an industrial park in Jiangsu province as a test load, wherein the photovoltaic input data is derived from a database of solar forecast and intelligent power grid projects in the Netherlands. Typical daily selection is shown in table 1, and time-sharing energy prices are shown in table 2.
Table 1 typical day selection
Figure BDA0004058974060000123
TABLE 2 price of time-sharing electricity and gas purchase [ this/kWh ]
Figure BDA0004058974060000124
In this example we compared the proposed method with respect to investment cost, total cost and carbon emissions with a campus planning method that does not take into account generalized energy storage, where the building and data center temperatures are controlled to a constant value without peak clipping and valley filling capabilities, as shown in table 3. The investment cost of the planning method for the park comprehensive energy system considering the flexible operation mode of the generalized energy storage resources is slightly higher than that of the conventional method, but the total annual planning cost is reduced by 170.7 ten thousand yuan, the annual carbon emission is reduced by 0.15 kiloton, and the temperature of a building and a data center can be changed in a specific comfort zone to realize energy cross-period transfer, so that the planning method considering the flexible operation of the park generalized energy storage main body improves the operation economy and environmental protection of the park comprehensive energy system.
Table 3 investment costs and carbon emissions comparisons with and without park generalized energy storage resources
Figure BDA0004058974060000131
In this example, a campus planning year summer generalized energy storage scenario that accounts for generalized energy storage is shown in fig. 4. As can be seen from fig. 4, in a typical day in summer, electricity prices are lower in 2-7 and 15-18 periods, the generalized energy storage preferentially stores energy, that is, the energy storage histogram has a "double-peak" characteristic, wherein the storage battery purchases electric energy in the flat valley period and discharges electric energy in 12-14 and 19-22 periods, so as to meet the load demand at peak time. The data center and the building perform thermal energy storage in the peak period and the valley period of new energy output, and release part of energy storage in the 18-23 period. From the above results, the generalized energy storage utilizes the flexible charge/discharge characteristics thereof to realize peak clipping and valley filling of energy loads for park diversification.

Claims (6)

1. A carbon emission planning method of a park comprehensive energy system considering generalized energy storage is characterized by comprising the following steps:
step 1: collecting historical photovoltaic output data of a park comprehensive energy system and unit investment cost mainly comprising photovoltaic and storage battery equipment;
step 2: according to known park parameters, constructing a scene of long-term uncertainty of multi-energy load growth and short-term uncertainty of photovoltaic output by taking the sum of extension cost and operation cost of a park comprehensive energy system as an objective function;
step 3: establishing an energy hub model taking park hub as a core and park generalized energy storage resource equipment as a main body;
step 4: establishing a generalized energy storage model according to the energy storage of the storage battery, the thermal characteristics of the data center and the thermal characteristics of the building;
step 5: and calculating the operational cost and the extension cost of 365 days of planning years of the park, and outputting the optimal planning and operational strategy.
2. The method for planning carbon emissions in a campus integrated energy system involving generalized energy storage according to claim 1, wherein constructing a scenario of long-term uncertainty of multi-energy load increase in step 2 specifically comprises: generating three kinds of upper layer multipotency load long-term uncertainty scenes according to seasonal differences, namely transitional seasons, summer and winter respectively;
the short-term uncertainty scene of the photovoltaic output specifically comprises: based on the photovoltaic historical output data, a k-means clustering method is utilized to generate four types of short-term uncertainty scenes of the lower photovoltaic output.
3. A method for planning carbon emissions in a campus integrated energy system involving generalized energy storage according to claim 2, wherein the objective function in step 2 is:
Figure FDA0004058974050000011
in the formula (A-1), a subscript t represents a scheduling time;
Figure FDA0004058974050000012
capacity expansion for photovoltaic units, +.>
Figure FDA0004058974050000013
The capacity of the accumulator unit is extended,
Figure FDA0004058974050000014
capacity expansion for ground source heat pump unit +.>
Figure FDA0004058974050000015
The expansion capacity of the combined cooling heating power unit is provided; />
Figure FDA0004058974050000016
Cost for unit expansion of photovoltaic unit, +.>
Figure FDA0004058974050000017
Cost for unit expansion of storage battery unit, < >>
Figure FDA0004058974050000018
The unit expansion cost of the ground source heat pump unit is +.>
Figure FDA0004058974050000019
The unit extension cost of the combined cooling heating power unit is set; omega shape g Omega as a typical day collection g Weighting for each typical day; omega shape PV Omega for the photovoltaic output scene set PV The method comprises the steps of weighing a photovoltaic output scene; p (P) t E For the upper power grid electricity purchasing quantity of t time period, < >>
Figure FDA00040589740500000213
The upper-level gas online shopping amount is t time period; />
Figure FDA0004058974050000021
For the electricity price of the t period, < > for>
Figure FDA0004058974050000022
Gas price for period t; i E Carbon emission cost per unit electricity purchase amount, I G Carbon emission costs per unit amount of purchased gas;
in the formula (A-1), the first four items are annual extension costs of each extension unit, namely a photovoltaic unit, a storage battery unit, a ground source heat pump unit and a combined cooling, heating and power unit, and the fifth item is annual operation costs, wherein carbon emission costs of electricity purchase and gas purchase are taken into account;
each planning variable constraint is:
Figure FDA0004058974050000023
in the formula (A-2),
Figure FDA0004058974050000024
extend the upper capacity limit for the photovoltaic unit, < >>
Figure FDA0004058974050000025
Extend the upper limit of capacity for the battery unit, +.>
Figure FDA0004058974050000026
Expansion capacity upper limit for ground source heat pump unit, < ->
Figure FDA0004058974050000027
The expansion capacity upper limit of the combined cooling heating and power unit is +.>
Figure FDA0004058974050000028
For the original photovoltaic installation capacity, < > a->
Figure FDA0004058974050000029
For the original battery capacity>
Figure FDA00040589740500000210
The capacity of the original cold-heat-electricity triple supply unit is +.>
Figure FDA00040589740500000211
The capacity of the original ground source heat pump unit; the first four items are investment upper limit constraints of the extension unit, and the last four items are operation upper limit constraints of the extension unit.
4. The method for planning carbon emissions in a campus integrated energy system involving generalized energy storage according to claim 1, wherein step 3 specifically comprises:
the energy level constraint of the park is as follows:
Figure FDA00040589740500000212
wherein P is t LC For the period t of the cold load demand,
Figure FDA00040589740500000214
industrial cooling capacity of GSHP unit in t period>
Figure FDA00040589740500000215
The cooling capacity of the CCHP unit is supplied for the period t; p (P) t LH For t-period heat load demand, < > and->
Figure FDA00040589740500000216
Heat is supplied to the t-period CCHP unit; />
Figure FDA0004058974050000039
Heat is supplied to the GSHP unit in the t period; />
Figure FDA0004058974050000031
Heat supply of the air boiler in the t period of time, +.>
Figure FDA0004058974050000032
Heat is supplied to the electric boiler in the period t; p (P) t PV For the t period of photovoltaic output, P t dis For t period of battery charge, P t LE For t-period electrical load demand, P t GSHP For t-period industrial GSHP unit power consumption, P t EB For t-period electric boiler power consumption, P t ch For the discharge quantity of the storage battery set in the t period, P t BD For the power consumption of the GSHP unit of the building at the period of t, P t DC Power consumption of the data center is t time periods; />
Figure FDA0004058974050000033
The electric quantity provided for the t-period CCHP unit; p (P) t CCHP For the consumption of the CCHP unit in t period, P t GB And the air consumption of the air boiler is t time periods. Constraint (A-2) is a cold bus, a hot bus, an electric bus and an air bus balance equation from top to bottom in sequence;
the operational constraints of conventional multi-energy conversion equipment in a park are as follows:
Figure FDA0004058974050000034
in the method, in the process of the invention,
Figure FDA0004058974050000035
for CCHP unit cold conversion efficiency, P t CCHP,C Natural gas quantity for cold transfer is input into the CCHP unit;
Figure FDA0004058974050000036
for the thermal conversion efficiency of the CCHP unit, P t CCHP,H Natural gas quantity for transferring heat is input into the CCHP unit; />
Figure FDA0004058974050000037
For the electric conversion efficiency of the CCHP unit, P t CCHP,E Natural gas quantity for inputting CCHP units for power conversion;P t TF to input the electric quantity of the transformer, eta TF The electric conversion efficiency of the transformer; η (eta) EB For the heat conversion efficiency of the electric boiler, eta GB The heat conversion efficiency of the gas boiler is; the first three terms are CCHP unit operation equations, and the last three terms are respectively transformer, electric boiler and gas boiler operation equations;
equipment operating constraints related to campus generalized energy storage resources include storage batteries, ground source heat pumps, and data centers, where:
the operation constraint of the storage battery unit is as follows:
Figure FDA0004058974050000038
wherein P is t ch For t period of energy storage charging power, P t dis Energy storage and discharge power is t time period;
Figure FDA0004058974050000041
for the t period of the energy-storing charge state variable, +.>
Figure FDA0004058974050000042
The energy storage and discharge state variable is t time period; e (E) t Energy storage electric quantity is carried out for a period t; lambda (lambda) E For t period of energy storage self-discharge coefficient, eta ch For the energy storage and charging efficiency of t period, eta dis Energy storage and discharge efficiency is t time period; e (E) min For storing energy, E is the lower limit of the electric quantity max Is the upper limit of the energy storage electric quantity;
the operation constraint of the ground source heat pump unit is as follows:
Figure FDA0004058974050000043
wherein P is t BD,C Electric energy for converting heat for t-period building GSHP unit, P t BD,H The GSHP unit is used for converting cold electric energy for the building at the time of t;
Figure FDA0004058974050000044
and->
Figure FDA0004058974050000045
The heat energy and the cold energy which are respectively transmitted to the building by GSHP in t time period; />
Figure FDA0004058974050000046
And->
Figure FDA0004058974050000047
Heating and cooling operation state variables are respectively;
the operation constraint of the data center is as follows;
Figure FDA0004058974050000048
wherein P is t IT For a period t data center IT device power consumption,
Figure FDA0004058974050000049
power for each k-type server for a period t; />
Figure FDA00040589740500000410
For a fixed power consumption of a t-period k-class server, < >>
Figure FDA00040589740500000411
CPU power consumption is t time period; c (C) 1 D 'is the power consumption coefficient of the server' k,t The data load quantity required to be processed for the k-type server of the data center in the t period; a, a k,s,t The working frequency zone bit is s-gear of a k-type server CPU of the data center in the t period; f (f) k,s And the working frequency is s-gear for the CPU of the k-type server of the data center.
5. The method for planning carbon emissions of a campus integrated energy system with generalized energy storage according to claim 1, wherein the building thermal characteristics and thermal characteristics of the data center are used to build the generalized energy storage model according to the energy storage of the storage battery, comprising:
generalized energy storage capacity of storage battery in t period is E t
The thermal characteristics and generalized energy storage constraints of the data center are:
Figure FDA0004058974050000051
in the method, in the process of the invention,
Figure FDA0004058974050000052
generalized energy storage for data center t period, H DC Specific heat capacity of data center, T t dc,in Temperature in data center T time period room, T t dc,in,max An upper limit of the indoor temperature in the t period; η (eta) c Data center refrigeration unit efficiency, Q t Refrigeration T produced by refrigeration equipment of data center T-section t dc,out For t period of external environment temperature, U C The heat transfer coefficient of the data center to the outside; />
Figure FDA0004058974050000053
Lower limit of temperature change rate of data center, +.>
Figure FDA0004058974050000054
An upper limit for the temperature change rate of the data center; />
Figure FDA0004058974050000055
Is the indoor temperature lower limit of the data center, +.>
Figure FDA0004058974050000056
An upper temperature limit in the data center room;
the thermal characteristics and generalized energy storage constraints of the park building are as follows:
Figure FDA0004058974050000057
in the method, in the process of the invention,
Figure FDA00040589740500000512
generalized energy storage for data center t period, H BD Specific heat capacity of data center, T t BD Indoor temperature in t time period of data center,/o>
Figure FDA0004058974050000058
For the upper limit of the indoor temperature in the t period, < >>
Figure FDA0004058974050000059
The lower limit of the indoor temperature is t time period; />
Figure FDA00040589740500000510
Phi for the solar radiation quantity passing through the glazing of a building t The solar radiation heat is in a t period, and sigma is the solar radiation loss rate; />
Figure FDA00040589740500000511
For building to radiate heat outwards, T t E Building wall temperature for the park at time t; t (T) t The indoor temperature of the park at the period t; u (U) A The heat transfer coefficient is the direct heat transfer coefficient of the outer vertical surface of the building; u (U) E Heat transfer coefficients for indoor to wall; t (T) t a Is the external environment temperature in the t period. Wherein T is t E For the temperature of the building outer wall in t time period,/for the building outer wall temperature>
Figure FDA0004058974050000061
For the temperature of the building outer wall in t-1 time period, H PE Is the specific heat capacity of the building wall, T t a Is the external environment temperature in the t period.
The generalized energy storage total amount of the park is as follows:
Figure FDA0004058974050000062
6. the method for planning carbon emissions in a campus integrated energy system with generalized energy storage according to claim 1, wherein calculating the operational cost and the extension cost for 365 days of the campus planning, outputting the optimal planning and the operational strategy comprises:
in the actual solving process, the objective function is the sum of annual investment cost and planning annual operation cost, wherein the annual investment cost is the total extension cost divided by the planning age, and when the planning annual operation cost is calculated, all loads and photovoltaic output conditions during the planning period need to be met by each equipment operation constraint.
CN202310052701.9A 2023-02-02 2023-02-02 Carbon emission planning method for park comprehensive energy system considering generalized energy storage Pending CN116029518A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310052701.9A CN116029518A (en) 2023-02-02 2023-02-02 Carbon emission planning method for park comprehensive energy system considering generalized energy storage

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310052701.9A CN116029518A (en) 2023-02-02 2023-02-02 Carbon emission planning method for park comprehensive energy system considering generalized energy storage

Publications (1)

Publication Number Publication Date
CN116029518A true CN116029518A (en) 2023-04-28

Family

ID=86075878

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310052701.9A Pending CN116029518A (en) 2023-02-02 2023-02-02 Carbon emission planning method for park comprehensive energy system considering generalized energy storage

Country Status (1)

Country Link
CN (1) CN116029518A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117094507A (en) * 2023-08-21 2023-11-21 四川大学 Method and system for planning agricultural industry building comprehensive energy based on biomass resources

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117094507A (en) * 2023-08-21 2023-11-21 四川大学 Method and system for planning agricultural industry building comprehensive energy based on biomass resources
CN117094507B (en) * 2023-08-21 2024-02-20 四川大学 Method and system for planning agricultural industry building comprehensive energy based on biomass resources

Similar Documents

Publication Publication Date Title
Luo et al. Multi-objective capacity optimization of a distributed energy system considering economy, environment and energy
CN111445090B (en) Double-layer planning method for off-grid type comprehensive energy system
Liu et al. Two-phase collaborative optimization and operation strategy for a new distributed energy system that combines multi-energy storage for a nearly zero energy community
CN111738502B (en) Multi-energy complementary system demand response operation optimization method for promoting surplus wind power consumption
CN108629462A (en) Comprehensive energy microgrid Method for optimized planning containing energy storage and comprehensive energy micro-grid system
CN111144620A (en) Electricity-hydrogen comprehensive energy system considering seasonal hydrogen storage and robust planning method thereof
CN114722591B (en) Planning method for electric heating hydrogen multi-energy flow energy supply equipment of building with zero energy consumption
Liu et al. Co-optimization of a novel distributed energy system integrated with hybrid energy storage in different nearly zero energy community scenarios
CN112907098B (en) Multi-stage capacity configuration method and configuration system of park comprehensive energy system
CN111144707A (en) Multi-energy system collaborative planning modeling method based on energy hub
CN110956314B (en) Improved particle swarm optimization-based capacity planning method for hybrid energy storage system
CN110363353A (en) The optimization design and dispatching method and system of a kind of Distributed Integration energy resource system
CN112446141B (en) Double-layer planning method for electric heating comprehensive energy system
CN109993345B (en) Garden-oriented dynamic economic dispatching method for multi-energy complementary system for island operation
Fan et al. Two-layer collaborative optimization for a renewable energy system combining electricity storage, hydrogen storage, and heat storage
CN111639819B (en) Multi-stage optimization control method for comprehensive energy park
CN115099007B (en) Comprehensive energy system optimized operation method based on comprehensive cost-energy consumption curve
CN113240204A (en) Energy station capacity optimal configuration method and system considering renewable energy consumption area
CN115857348A (en) Distributed energy system capacity optimization method considering comfortable energy supply of two-combined heat pump
Guan et al. Optimal configuration and operation of multi-energy complementary distributed energy systems
CN116029518A (en) Carbon emission planning method for park comprehensive energy system considering generalized energy storage
CN117081143A (en) Method for promoting coordination and optimization operation of park comprehensive energy system for distributed photovoltaic on-site digestion
CN114818078A (en) Planning method for net zero energy consumption building energy system considering long-term and short-term uncertainty
CN114742276A (en) Multi-objective optimization scheduling method of park comprehensive energy system with ORC (organic Rankine cycle) considering exergy efficiency
Yuan et al. An advanced multicarrier residential energy hub system based on mixed integer linear programming

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