Journal of Modern Power Systems and Clean Energy

ISSN 2196-5625 CN 32-1884/TK

网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

Optimal Scheduling of Integrated Energy Systems in Low-carbon Communities Considering Flexibility of Resources and Segmental Control of Solid Oxide Fuel Cells  PDF

  • Mao Yang
  • Yuxin Wang
  • Jinxin Wang
  • Dongxu Liu
  • Weihang Xu
Key Laboratory of Modern Power System Simulation and Control & Renewable Energy Technology, Ministry of Education, Northeast Electric Power University, Jilin 132012, China

Updated:2025-05-21

DOI:10.35833/MPCE.2023.000888

  • Full Text
  • Figs & Tabs
  • References
  • Authors
  • About
CITE
OUTLINE

Abstract

To address the strong thermoelectric coupling of the combined heat and power (CHP) units, the low utilization rate of energy storage, and the underexploitation of load-side resource flexibility in integrated energy systems (IESs), this paper proposes an optimal scheduling model of an IES in low-carbon communities considering flexibility of resources and the segmental control of solid oxide fuel cells (SOFCs). Firstly, by replacing the gas turbine (GT) in the CHP unit with an SOFC array to reduce carbon emissions and simultaneously weakening the thermoelectric coupling of the CHP unit, the segmental control method is used to control the SOFC array to improve the overall efficiency of the CHP unit. Secondly, coupled interactions among different types of energy storage equipments are mobilized through the integrated energy storage system to make full use of the remaining space in the heat and natural gas storage tanks. Finally, load-side flexible resources are utilized by considering transferable, substitutable, and heat loads, taking into account the thermal inertia of the building and categorizing rooms based on floors, orientations, and room area. Additionally, different user characteristics are characterized, and the flexible resources of building heating periods in northern cities in China are tapped in depth according to the actual factors. Compared with the traditional model, the optimal scheduling model proposed in this paper can reduce the wind abandonment rate and the carbon emission of community-integrated energy system (CIES) by 4.54% and 70.63%, respectively, and increase the utilization rate of heat and natural gas storage tanks by 12.34% and 30.52%, respectively, and lower the total cost by ¥2183.6 under the premise of ensuring user comfort during energy consumption, which promotes the economic and low-carbon operation of the CIES.

I. Introduction

WITH the introduction of the target “carbon emission peaking and carbon neutrality”, the renewable energy installations in China are increasing [

1]. However, this trend highlights several problems in integrated energy systems (IESs) in terms of the sizeable thermoelectric coupling of the combined heat and power (CHP) units [2], the low utilization of energy storage [3], and the underexploitation of flexibility of load-side resources [4], which directly contributes to the wind abandonment, whereby a large amount of wind resources are not effectively utilized. Therefore, it is necessary to strengthen the technical transformation and optimization of CHP units, improve the utilization rate of energy storage, fully tap into load-side flexibility, realize the sustainable and efficient utilization of new energy sources, and contribute to the transformation and upgrade of the energy infrastructure of China.

In the context of the energy crisis and environmental protection, conventional CHP technology has become a vital energy utilization method [

5]. Compared with conventional gas turbines (GTs), solid oxide fuel cells (SOFCs) have a higher electricthermal conversion efficiency, which can more fully utilize chemical energy of fuel, reduce energy waste, and improve energy efficiency [6]-[8]. Owing to their operation principles and application of unique materials, SOFCs produce less pollutant emissions during the electrochemical process, which are more environmentally friendly [9]. However, the SOFC-CHP plants using a single SOFC are not suitable for scenarios with high electrothermal power demands because of their power rating limitations, and few studies address the combined operation and control of multiple SOFCs.

Existing studies on how IESs consume wind power suggest that IESs should be equipped with energy storage, which can maximize wind power resources by storing wind power during high wind power output periods and releasing it during low wind power output periods. Reference [

10] applies a battery energy storage system for load shaving and peak filling to establish an optimization model to maximize the net benefit. A battery energy storage system configuration method is proposed in [11]. In addition, a thermoelectric hybrid energy storage system is used to increase the wind power accommodation of IES [11]. Part of the IES is equipped with energy storage, but could not entirely consume wind power. There still exists the wind abandonment phenomenon, primarily due to the lack of storage capacity [12]. However, too much energy storage capacity leads to high investment costs and low utilization rates of energy storage, and few studies explore how to balance the relationship between storage capacity and utilization rates of energy storage.

For load-side flexible resources, three main types of demand response exist: price-based demand response, which is based on the impact of time-of-energy prices; incentive-based demand response, which is based on incentive compensation instruments; and substitution-based demand response, which is based on the substitution of energy consumption by customers due to differences in heterogeneous energy price signals [

13]. Existing studies demonstrate the economic and low-carbon benefits of integrated demand response for IES operations [14]. Among them, the flexibility of heat load regulation is also a critical factor in changing the IES ability to consume wind power, and tapping the potential of heat load regulation plays a vital role in realizing the economic and low-carbon operation of IES.

The heat network storage and discharge characteristics and the thermal inertia of the heating area under a CHP system are analyzed in [

15] to exploit the elastic change in heat load resources relative to internal and external temperatures, promoting the wind power accommodation of IES. Reference [16] demonstrates that an integrated demand response for electricity and heat based on time-of-use tariffs and heating comfort ambiguity balances the IES economics and low-carbon nature. Notably, few studies address the mining of flexibility of the community building architecture, which has the potential to be mined in most northern cities in China with a relatively high proportion of residential areas and long heating time [17]. Moreover, existing models of thermal inertia of buildings are more general and cannot accurately portray the different characteristics of buildings for typical urban community users.

In summary, this paper focuses on combining multiple SOFCs to form a SOFC/GT-CHP system. A segmental control method is proposed to improve the overall efficiency of the CHP units and weaken the thermoelectric coupling of the CHP units. At the same time, to improve the utilization rate of energy storage, the remaining space of heat and natural gas storage is used to improve the utilization rate of wind power. The primary focus of this paper is to analyze the load-side flexible resources in community buildings during heating period, accounting for the actual factors involved in buildings and analyzing load-side flexible resources in depth.

This paper proposes an optimal scheduling model of IES in low-carbon communities considering flexibility of resources and the segmental control of SOFC to promote the wind power consumption, realize the economic and low-carbon operation of the community-integrated energy system (CIES), and ensure the stable and efficient consumption of energy.

II. Architecture of Operation of CIES

Figure 1 shows the block diagram of the CIES. As shown in Fig. 1, a CIES consists of an upper energy network containing power grid and natural gas network. Wind power is a renewable energy source, and coupling equipment includes an SOFC/GT-CHP system, a power-to-gas (P2G), a gas boiler (GB), and a ground-source heat pump (GSHP). Among these resources, the electricity grid and wind turbine provide electricity and natural gas supplies, respectively; the natural gas networks and P2G provide natural gas supplies; and the SOFC/GT-CHP system, GB, and GSHP provide heat supplies. The combined energy storage system (CESS) stores excess energy in the CIES, converts, and releases it when needed. The energy demand side contains electric, heat, and natural gas loads.

Fig. 1  Block diagram of CIES.

A. SOFC/GT-CHP System

The operating principles of SOFC/GT-CHP system are shown in Fig. 2. The SOFC/GT-CHP system combines the power of an SOFC and a GT to improve the overall efficiency of CHP unit. The GT, which drives high-speed electric motor, can generate electricity from the exhaust waste heat generated by SOFCs at high temperature. Despite the leading role of the SOFC in power generation, the GT can be used as a backup device when needed.

Fig. 2  Operating principles of SOFC/GT-CHP system.

B. SOFC Array Control Strategy Based on Segmental Control Method

Considering that a single SOFC has a small generating power output, multiple SOFCs are required to produce energy simultaneously to meet high electric power demands. To control SOFCs, the simplest method is to start other SOFCs sequentially when a power shortage occurs. However, since the operating efficiency of an SOFC is related to the output power, the specific relationship of operating efficiency and output power is shown in Fig. 3. Thus, by controlling a single SOFC to maintain the optimal output power, the operating efficiency of each SOFC can be significantly improved, thus improving the overall efficiency.

Fig. 3  Specific relationship of operating efficiency and output power.

As shown in Fig. 4, the segmental controller consists of a segmental control selector, a segmental selection switch, n SOFCs, and an actuator. Figure 5 shows the control diagram of SOFC based on segmental control. The specific process of the segmental control is as follows.

Fig. 4  Block diagram structure of segmental controller.

Fig. 5  Control diagram of SOFC based on segmental control.

Step 1:   determine the output power of the SOFC array for each period according to the CIES scheduling results.

Step 2:   determine the power output interval of the SOFC array.

Step 3:   allocate the specific output power of each SOFC via Step 1, which is determined as follows: when PSOFCall is in power interval N1, all the SOFCs do not operate; when PSOFCall is in power interval N2, there are k SOFCs in the SOFC array with n SOFCs operating at the optimal output power, and the others do not operate; when PSOFCall is in power interval N3, all the SOFCs in the SOFC array can operate at the optimal output power; when PSOFCall is in power interval N4, all the SOFCs are operating, with l SOFCs operating at the rated power and the rest operating at the optimal output power; and when PSOFCall is in power interval N5, all the SOFCs are operating at the rated power.

According to the allocation above, the actuator execution is input to control each SOFC.

PSOFCall=0                                      PSOFCallN1kPSOFCo                            k<n,PSOFCallN2nPSOFCo                            PSOFCallN3lPSOFCrated+(n-l)PSOFCo    l<n,PSOFCallN4nPSOFCrated                            PSOFCallN5 (1)

where PSOFCo is the optimal output power of each SOFC; and PSOFCrated is the rated power of each SOFC.

C. CESS

As a result of the increase in energy storage device for wind power accommodation, traditional energy storage device has limitations due to the limited capacity of the coupling device and the weak coupling link, thus reducing the utilization of the energy storage space.

The CESS includes control systems, batteries, heat storage tanks, and natural gas storage tanks, as well as storage tanks applicable for microconversion device. The microconversion device in the CESS does not operate when each energy storage state does not reach its peak value. When the battery reaches its limit and cannot hold the excess electric energy, the CESS can convert the excess electric energy to a microconversion device, which then convert the excess electric energy into natural gas and heat energy, where the remaining space in the natural gas and heat storage tanks is used to store natural gas and heat energy, respectively. As part of a CESS, the control system monitors and manages the conversion and storage processes among electric, thermal, and natural gas energy. The control system can intelligently optimize the way energy is converted and stored according to supply and demand, resulting in the highest possible efficiency of energy consumption.

The model representation is given as:

PinCESS=PechaCESS+PhchaCESS+PpchaCESSPoutCESS=PedisCESS+PhdisCESS+PpdisCESSPpme2p=Peme2pλme2pPhmeh=Pemehλmeh (2)
PedisCESSPhdisCESSPpdisCESS=PechaCESSPhchaCESSPpchaCESS+-1-100λmeh0λme2p00Peme2pPemeh0 (3)

where PinCESS and PoutCESS are the input and output power of the CESS, respectively; PedisCESS, PhdisCESS, and PpdisCESS are the discharging power, heat releasing power, and natural gas release power of the CESS, respectively; PechaCESS, PhchaCESS, and PpchaCESS are the charging power, heating power, and inflation power of the CESS, respectively; Peme2p and Ppme2p are the electric power consumed to charge the storage tanks with excess battery power and the natural gas power charged into the CESS, respectively; Pemeh and Phmeh are the electric power consumed to charge the heat storage tank with excess battery power and the heating power as well as charged to the heat storage tank in the CESS, respectively; and λme2p and λmeh are the operating efficiencies of the microconversion devices in the CESS.

To evaluate the efficiency of utilization of energy storage, this paper takes the capacity utilization rate of energy storage as an evaluation index.

θc,𝓁=Cac,𝓁Cra,𝓁×100% (4)

where subscript 𝓁 is the index of energy storage type; θc,𝓁 is the capacity utilization rate; Cac,𝓁 is the actual capacity; and Cra,𝓁 is the total capacity.

III. Flexible Resources

A. Flexible Heat Loads Considering Thermal Inertia of Building

Traditionally, centralized heating systems cause significant fluctuations in indoor temperature, which makes users uncomfortable. Within a building, distributed heating strategies are thermally inert, resulting in small fluctuations in indoor temperature. In this way, the stability and comfort of the indoor temperature can be improved, allowing for the use of operating strategies based on the temperature comfort zones of the users. Furthermore, it provides more room for the consumption of renewable energy.

1) Level of Temperature Comfort Zone

Individuals have different indoor temperature comfort requirements during different time periods. Hence, their work and rest patterns are considered via a time-sharing model. On the basis of their indoor temperature comfort requirements, the time in a day can be classified into two time periods: daytime and nighttime [

18]. Table I shows the indoor temperature comfort range of different time periods.

TABLE I  Indoor Temperature Comfort Range of Different Time Periods
Time periodLower comfort limit (℃)Upper comfort limit (℃)
07:00-22:00 20 24
22:00-07:00 17 21

2) Equivalent Model of Room Temperature Variability Process

The temperature variability process in the building is determined by both the heating power of the electric heating device and the outdoor temperature. The schematic of the equivalent model of room temperature variability process shown in Fig. 6 can approximate the room temperature variability process [

19]. The set of total parameter ordinary differential equations describing the room temperature change process is given as:

Cindθin(t)dt=Pheat(t)-θin(t)-θwall(t)R1Cwalldθwall(t)dt=θin(t)-θwall(t)R1-θwall(t)-θout(t)R2 (5)

Fig. 6  Schematic of equivalent model of room temperature variability process.

where Cin and Cwall are the equivalent heat capacity of the indoor air and the equivalent heat capacity of the wall, respectively; R1 and R2 are the equivalent thermal resistance of the indoor air and the inside of the wall and the equivalent thermal resistance of the outdoor air and the outside of the wall, respectively; θin(t), θwall(t), and θout(t) are the indoor, wall, and outdoor temperatures during period t, respectively; and Pheat is the heating power of the heat source in the room.

Figure 7 shows the heating power-temperature variation characteristics. Heating device has two operating modes: on (i.e., maintenance, full power, etc.) and off. At the lower comfort limit, the device is turned on to maintain the temperature, and the power is less than the rated power. The fully rated state increases the room temperature, and when the temperature increases to the upper comfort limit, the device is turned off or maintained at the upper comfort limit.

Fig. 7  Heating power-temperature variation characteristics.

B. Substitutable Loads

To meet the energy demands and improve economics, users with multiple energy loads intelligently select the most cost-effective energy form on the basis of differences in the temporal distribution of energy price. When energy prices are low, users can adjust their energy sources to choose more appropriate ones. They can also convert energy by switching to more economical options during peak periods. This paper investigates the conversion of electric heat into electric energy while considering heterogeneous energy price factors and deviations in the conversion efficiency of coupling CHP. By comparing electricity and natural gas prices, users can choose the energy supply mode directly, and improve their energy economy by comparing the integrated electricity and natural gas prices. Compared with a lower natural gas price, a higher composite electricity price encourages users with multiple energy loads to use natural gas as a means to reduce the electricity consumption, whereas a relatively low composite electricity price encourages users with multiple energy loads to increase their electricity demand and reduce natural gas consumption to maximize their benefits, as well as switch to electric heat. The substitutable load modeling method is as follows.

Pi,loadal(t)=Pi,loadals0(t)+ΔPi,loadal(t)ΔPi,loadal(t)=Sialin(t)Pi,loadalin(t)-Sialout(t)Pi,loadalout(t)Sialin(t)+Sialout(t)1t=1T(ΔPe,loadal(t)+ΔPh,loadal(t)+ΔPp,loadal(t))=0ΔPi,load,minalΔPi,loadal(t)ΔPi,load,maxal (6)
Pout=CPinPin=[Pin,e  Pin,p  Pin,h]TPout=[Pout,e  Pout,p  Pout,h]TC=α1β1ηpeSOFC/GT-CHP0α2ηepP2G1-2β10(1-α1-α2)ηehGSHPβ1(ηphSOFC/GT-CHP+ηphGB)1 (7)

where subscript i=e, p, or h is the index of the load type, which denotes the electric, natural gas, or heat load; Pi,loadals0(t) is the load before substitution during period t; T is the total number of periods; Pi,loadal(t) is the load after substitution during period t; ΔPi,loadal(t) is the substitution amount for each load; Sialin(t) and Sialout are the substitution flags for Pi,loadalin(t) and Pi,loadalout(t) during period t, respectively; ΔPi,load,minal(t) and ΔPi,load,maxal(t) are the lower and upper limits of the substitution amount for each load during period t, respectively; Pin and Pout are the matrices of the power of energy input and output of the CIES, respectively; C is the matrix of the energy coupling; α1 and α2 are the energy distribution coefficients that determine the proportional distribution coefficients of the energy flow in the coupled device; β1 and β2 are the energy distribution coefficients that determine the proportional distribution coefficients of the energy flow in the coupled loads; η is the energy conversion factor, which determines the efficiency of CIES in the energy conversion process; the superscripts P2G, SOFC/GT-CHP, GB, and GSHP denote the device types; and the subscripts eh, ep, ph, and pe denote the electricity to heat, electricity to natural gas, natural gas to heat, and natural gas to electricity, respectively.

C. Transferable Loads

A transferable load is the demand response that directly compensates the user for the time adjustment of energy consumption in the form of incentive compensation. This compensation needs to be carried out in a way to satisfy the balance between the supply and demand of the IES and stable operation. The demand response of transferable loads should be negotiated with users in advance. When the supply and demand relationship is tense, economic means should be provided to compensate for the user adjustments in hours, alleviating the imbalance between the supply and demand while also ensuring the stability of the IES. Load shifting from one period to another period can have a peak-shaving and valley-filling effect on the load curve. Although the total amount shifted remains the same throughout the cycle, it can affect user comfort. Therefore, there is a comfort cost.

Pi,loadex(t)=Pi,loadex0(t)+ΔPi,loadex(t)ΔPi,loadex(t)=Siexin(t)Pi,loadexin(t)-Siexout(t)Pi,loadexin(t)Siexin(t)+Siexout(t)1t=1TΔPi,loadex(t)=0ΔPi,load,minex|ΔPi,loadex(t)|ΔPi,load,maxex (8)

where Pi,loadex0(t) is the predicted value of the load during period t; Pi,loadex(t) is the load after shifting during period t; Siexin and Siexout are the 0-1 variables for load shifting during period t (Siexin=1 denotes transfer in, Siexout=1 denotes transfer out, and Siexin=Siexout=0 denote no load shifting); and ΔPi,load,maxex and ΔPi,load,minex are the upper and lower transfer limits, respectively.

IV. Objective Function and Constraints

A. Objective Function

This paper aims to minimize the total cost of the CIES, which includes energy purchasing costs CP from the primary power grid and natural gas network, carbon trading expenses CCO2, comfort compensation costs CCOM, and wind abandonment penalty costs CW.

min C=min (CP+CCO2+CCOM+CW) (9)

where C is the total cost of the CIES.

1) Energy Purchasing Costs

CP=t=124(Pe,P(t)ce(t)+Pp,P(t)cp(t)) (10)

where Pe,P(t) is the amount of electricity purchased from the main power grid during period t; Pp,P(t) is the amount of natural gas purchased from the main natural gas network during period t; and ce(t) and cp(t) are the unit prices of purchased electricity and natural gas during period t, respectively.

2) Carbon Trading Expenses

The actual carbon emissions for the electricity and heat supply of the CIES are determined via the following equations.

EIES=t=1T(a1+b1Pe,P(t)+c1Pe,P2(t)+a2+b2Pptr(t)+c2Pe,P2(t)) (11)
Pptr(t)=Pp,GB(t)ηGB (12)

where EIES is the actual carbon emission of the CIES; a1, b1, and c1 are the carbon emission calculation coefficients of electricity purchased from the main power grid; a2, b2, and c2 are the carbon emission coefficients of heating; Pp,GB(t) is the input power of the GB during period t; and ηGB is the energy conversion efficiency of the GB.

Multiple carbon credits purchase bands are delineated by the stepped carbon trading model. As the CIES needs to purchase more carbon credit allowances, the purchase price of the corresponding band is higher, thus limiting the output of the high-emission device [

6]. The modified carbon trading model is given as:

CCO2=aEIES                                                       EIESda(1+ω)(EIES-d)+ad                        d<EIES2da(1+2ω)(EIES-2d)+(2+ω)ad       2d<EIES3la(1+3ω)(EIES-3d)+(3+3ω)ad     3d<EIES4la(1+3ω)(EIES-4d)+(4+6ω)ad     4d<EIES (13)

where a is the base price for carbon trading; d is the interval length of carbon emissions of the CIES; and ω is the rate of price increment.

3) Comfort Compensation Costs

Substitutable loads can choose different energy supply methods to meet their energy demand simultaneously, and since they do not change the energy demand of the users, they do not incur comfort compensation costs. Heat loads that consider the thermal inertia of the building also satisfy the comfort zone of the user and therefore incur no comfort compensation costs.

CCOM=t=124i=13ε|ΔPi,loadex(t)| (14)

where ε is the unit compensation coefficient for transferable loads.

4) Wind Abandonment Penalty Costs

CW=φt=1TPWP,cut(t) (15)

where φ is the penalty coefficient of the discarded wind power; and PWP,cut(t) is the abandoned wind power.

B. Constraints

1) Power Balance Constraints

Pe,P(t)+PWP(t)+Pe,SOFC/GT-CHP(t)=Pe,Load(t)+Pe,P2G(t)+                                                                  Pe,GSHP(t)+Pe,ES(t)Ph,GB(t)+Ph,SOFC/GT-CHP(t)+Ph,GSHP(t)=Ph,Load(t)+Ph,ES(t)Pg,P(t)+Pg,P2G(t)=Pg,Load(t)+Pg,SOFC/GT-CHP(t)+                                    Pg,GB(t)+Pg,ES(t) (16)

where PWP(t) is the actual wind power output during period t; Pi,Load(t) is the overall load after the demand response during period t; Pe,P2G(t) and Pp,P2G(t) are the electricity consumption and natural gas output of the P2G during period t, respectively; Pp,GB(t) and Ph,GB(t) are the natural gas consumption and heat output of the GB during period t, respectively; Pe,SOFC/GT-CHP(t), Ph,SOFC/GT-CHP(t), and Pg,SOFC/GT-CHP(t) are the the electricity output, heat power output, natural gas consumption of the SOFC/GT-CHP during period t, respectively; and Pi,ES(t) is the charging or discharging power of the CESS during period t.

2) Device Constraints

min PikPik(t)maxPikmin ΔPik(t)Pik(t)-Pik(t-1)maxΔPik(t) (17)

where minPik and maxPik are the lower and upper limits of the output of the device during period t, respectively; and min ΔPik(t) and max ΔPik(t) are the lower and upper limits of the ramping power of the device during period t, respectively.

3) Wind Power Output Constraints

0PWP(t)PWPer(t) (18)

where PWP(t) is the wind power output during period t; and PWPer(t) is the forecast of wind power output during period t.

4) Purchased Power Constraints

min Pe,PPe,P(t)max Pe,Pmin Pp,PPp,P(t)max Pp,P (19)

where min Pe,P and max Pe,P are the lower and upper limits of the purchased power, respectively; and min Pp,P and max Pp,P are the lower and upper limits of the purchased natural gas, respectively.

5) CESS Constraints

CESS constraints include those related to charging and discharging power limitations of the CESS, those associated with energy conversion devices, energy state constraints, and energy storage capacity limitations.

0Pi,chaCESS (t)Si,chaCESS(t)Pi,cha,maxCESS0Pi,disCESS (t)Si,disCESS(t)Pi,dis,maxCESS0Peme2p(t)Sme2pmin{0.2SeCESS(t-1),Pe,maxme2p}0Pemeh(t)Smehmin{0.2SeCESS(t-1),Pe,maxmeh}SeCESS(t)=SeCESS(t-1)+Pe,chaCESS(t-1)ηchaCESS-Pe,disCESS(t-1)ηdisCESSΔtShCESS(t)=ShCESS(t-1)+(Ph,chaCESS(t-1)+Phmeh(t-1))ηchaCESS-                 Ph,dis CESS(t-1)ηdisCESSΔtSpCESS(t)=SpCESS(t-1)+(Pp,chaCESS(t-1)+Ppme2p(t-1))ηchaCESS-                 Pp,disCESS(t-1)ηdisCESSΔtSi,chaCESS+Si,disCESS1Se,com,e2p+Sp,disCESS1Se,com,eh+Sh,disCESS1SiCESS(1)=SiCESS(24)Pi,cae,minSiCESS(t)ΔtPi,cae,max (20)

where Si,chaCESS and Si,disCESS are the indicators for CESS charging and discharging, respectively; Pi,cha,maxCESS and Pi,dis,maxCESS are the maximum charging and discharging power of the CESS, respectively; Sme2p and Smeh are the indicators of conversing electricity to natural gas and electricity to heat, respectively; Pe,maxme2p and Pe,maxmeh are the upper capacity limits of Peme2p and Pemeh, respectively; SiCESS(t) is the energy state of each storage unit in the CESS during period t; ηchaCESS and ηdisCESS are the charging and discharging efficiencies of the CESS, respectively; and Pi,cae,min and Pi,cae,max are the lower and upper capacity limits of each storage unit in the CESS, respectively.

C. Simulation Environment

This paper uses a PC with an Intel Core i5 processor and 8 GB of RAM to build an optimal scheduling model via MATLAB simulation with the YALMIP toolkit. By using the GUROBI commercial solver, we can optimize the controllable variables while satisfying the constraints and obtain an optimal solution to the objective function.

V. Case Study

This paper takes a CIES in a residential area in the northeastern of China as the research object, with N=300 heating users, assuming that the structure of each building in the residential area is the same and the typical room classification is shown in Fig. 8. However, the area of each room in the same building is different. The experimental rooms are divided into three layers and five categories according to the floors, orientations, and room areas. The parameters of typical room construction are shown in Appendix A Table AI. The parameter fitting results of the typical room are in Appendix A Table AII. The unit comfort compensation cost of load shifting and substitution is 0.1 ¥/kWh, the unit abandoned wind penalty cost is 0.07 ¥/kWh, the growth rate of the step carbon trading price is 25%, the carbon base price is 0.252 ¥/kg, the transferable and alternative loads of electricity account for 10% of the total load, and the transferable and alternative loads of heat and natural gas account for 5% of the total load. The day-ahead forecast curves of loads and wind power are shown in Appendix A Fig. A1.

Fig. 8  Typical room classification.

In this paper, five different scenarios are investigated to evaluate the performance of the CIES, and the information of the five scenarios is shown in Table II. Scenario 1 is the conventional method, which includes a conventional GT-CHP.

TABLE Ⅱ  Information of Five Scenarios
Scenario No.Conventional GT-CHPSOFC/GT-CHP simple controlSOFC/GT-CHP segmental controlLoad-side flexible resourceCESS
1 × × × ×
2 × × × ×
3 × × × ×
4 × × ×
5 × ×

A. Comparative Analysis of Two CHP Systems

To analyze the advantages of the SOFC/GT-CHP over conventional GT-CHP, this paper compares various scenarios. Compared with the previous analysis, the SOFC/GT-CHP has the characteristic of high efficiency and low-carbon operation. Figure 9 shows the comparison of carbon emissions between scenarios 1 and 2. According to Fig. 9, the SOFC/GT-CHP produces almost no carbon emissions during operation, especially during the period of high electricity prices.

Fig. 9  Comparison of carbon emissions between scenarios 1 and 2.

SOFC/GT-CHP becomes the leading heating equipment, which significantly decreases carbon emissions of the CIES with the increase of heating demand. Table III shows the costs and wind abandonment rates in various scenarios. According to Table III, after adopting the SOFC/GT-CHP, the total cost of the CIES is reduced by ¥1374.2, mainly due to the reduced carbon trading cost. Moreover, because the efficiency of SOFC/GT-CHP is higher than that of the conventional GT-CHP, to a certain extent, the conventional GT-CHP generated by the enormous heat load during the heating period is reduced, the thermoelectric coupling of the CHP is weakened, and some of the energy purchasing costs of the CIES are eliminated.

B. Comparative Analysis of Two Kinds of Controls of SOFCs

According to Table III, when we compare the results under the SOFC/GT-CHP simple control (scenario 2) and the SOFC/GT-CHP segmental control (scenario 3), we can find that scenario 3 reduces the total costs by nearly ¥344 compared with scenario 2. The energy purchasing cost is reduced from ¥5368.1 to ¥5062.2. As a result of the segmental control in scenario 3, the SOFC array is highly flexible when faced with different power generation demands. Thus, the SOFC/GT-CHP can operate more efficiently, producing more electricity and heat at the same energy accommodation level while minimizing the energy purchasing costs.

TABLE Ⅲ  Costs and Wind Abandonment Rates in Various Scenarios
Scenario No.Total cost (¥)Energy purchasing cost (¥)Carbon trading cost (¥)Wind abandonment penalty cost (¥)Comfort compensation cost (¥)CESS investment cost/day (¥)Wind curtailment rate (%)
1 7987.6 5728.4 1596.35 441.93 0 0 25.37
2 6613.4 5368.1 557.20 441.76 0 0 25.36
3 6269.2 5062.2 514.17 437.80 0 0 25.13
4 5904.7 4769.7 479.94 397.55 60.2 0 22.82
5 5804.0 4670.6 468.78 356.31 62.4 13.2 20.83

C. Flexibility Analysis of Load-side Resources

1) Load Shifting and Load Substitution

Scenario 4 considers load-side resources, which include transferable, substitutable, and heat loads. Figure 10 shows the load response in scenario 4. As observed from Fig. 10, for the transferable electric load, owing to the low wind power output during 11:00-20:00, the supply and demand relationship of electricity is tense. Moreover, since the SOFC/GT-CHP output reaches its capacity limit, the CIES incurs part of the energy purchasing cost. Thus, to reduce the energy purchasing cost of the CIES, the electric load can mainly transfer the load during the peak load period from 11:00 to 20:00 to the valley load period from 20:00 to 01:00. This measure can promote wind power accommodation, reduce the energy purchasing cost of the CIES, and also lower the heat load and natural gas load by adjusting the load distribution more reasonably. According to Fig. 11, for substitutable electric load, the equivalent natural gas price is lower than the equivalent electricity price during 00:00-08:00, on the basis of the equivalent electricity price calculated from the CIES. To maximize the user performance of the economy in terms of energy accommodation, natural gas and heat loads are substituted for electric loads as much as possible during this period. The same applies for other periods.

Fig. 10  Load response in scenario 4. (a) Electric load. (b) Heat load. (c) Natural gas load.

Fig. 11  Equivalent energy prices.

2) Flexible Heat Load Accounting for Thermal Inertia of Building

Temperature monitors are set up in the experimental rooms to analyze the changes in temperature and heat load in the building. The temperature curve of a random room in the community in the arithmetic example is shown in Fig. 12.

Fig. 12  Temperature curve of typical room

Figure 13 shows the power balance in scenario 4, which shows that as the wind power output increases between 22:00-24:00, the heat load is mainly met by the GSHP. To some extent, the change in the flexible heat load determines how much energy the CIES consumes. The flexible heat load should be increased as much as possible to dissipate excess wind power output during this period. Owing to the existence of the fuzzy intervals in the human body’s perception of temperature, the indoor temperature is increased. In this case, the room temperature does not affect the user comfort, so it reaches the upper comfort limit of the indoor temperature, at which point it stops increasing. In contrast, the flexible heat load increases significantly during 06:00-07:00. Owing to the change in the comfort zone, more heat is required indoors to meet the minimum temperature requirements. Afterward, to reduce the energy consumption and minimize the flexible heat load, the power increases during 10:00-11:00 when the equivalent electricity price is higher than the equivalent natural gas price. The CIES selects the more economical SOFC/GT-CHP to satisfy the heat load and recharges the batteries and heat storage tanks to cope with the period with high electricity price, and the CIES generates excess heat energy, which increases the indoor temperature. The CIES then cuts the flexible heat load to reduce the energy purchasing cost due to a gradual decrease in outdoor temperature. To reach the comfort zone, the CIES gradually increases the flexible heat load, and there is an air abandonment phenomenon during 20:00-21:00. Nevertheless, the indoor temperature does not change at this time because the climbing power of the GSHP is limited and the outdoor temperature decreases, which forces the GSHP to lift heat to counteract the change in outdoor temperature. The indoor temperature remains unchanged, causing the air abandonment phenomenon.

Fig. 13  Power balance in scenario 4. (a) Electric power balance. (b) Heat power balance. (c) Natural gas balance.

The scheduling process of flexible heat load participating in CIES is shown in Fig. 14. Combined with the analysis in Table III, scenario 4 reduces the total cost by ¥364.5 compared with scenario 3. It also reduces the wind abandonment rate by 2.31%, which shows that the flexibility of load-side resources considered in this paper can incentivize users to change their energy use habits to change the load profile. Therefore, wind power is accommodated more. Moreover, during the winter heating period, the CIES operator can control the overall heating for users, which can minimize part of the energy purchasing cost of the CIES while meeting the heat demands of the users, thus reducing the operating costs of the operator.

Fig. 14  Scheduling process of flexible heat load participating in CIES.

D. Comparative Analysis of CESS

Scenario 5 is the application of the CESS based on scenario 4. Figure 15 shows the energy storage comparison between scenarios 4 and 5. In scenario 4, the CIES may try to accommodate wind power due to the occurrence of the wind abandonment period. The electric energy storage space is fully utilized. In contrast, during the wind abandonment period, the demand for heat load is significant. Therefore, there is no excess heat to be stored in the heat storage tank, which results in a low utilization rate of the heat storage tank. The part of the space that needs to be fully utilized is for the scenarios in which wind power has a high utilization rate. As in scenario 4, the energy storage space in the natural gas storage tank needs to be fully utilized. In scenario 5, the charging and discharging pattern of the battery is unchanged in the CESS because its space is fully utilized. With the control system and micro-conversion device in the CESS, wind power output can be converted into heat and natural gas energy during wind abandonment period, and be stored in the heat and natural gas storage tanks. The control system stores heat and natural gas energy into heat and natural gas storage tanks and releases it during the subsequent periods when heat and natural gas energy is needed.

Fig. 15  Energy storage comparison between scenarios 4 and 5. (a) Battery energy state. (b) Energy state of heat storage tank. (c) Energy state of natural gas storage tank.

As shown in Fig. 15, since the amount of heat charged in the heat storage tank in the CESS is derived partly from the CIES heat flow and partly from the micro-conversion device in the CESS, the total amount of heat charged and released from the heat storage tank is increased in the overall scheduling process. The increased part is shown as the green area in Fig. 15, which is the additional heat energy stored in the heat storage tank. Similarly, the micro-conversion device increases the total charging and discharging power of the CESS. In scenario 5, compared with scenario 4, 12.34% more heat storage tanks are utilized, and 30.52% more natural gas storage tanks are utilized. Combined with the analysis of Table III, after the application of the CESS in the CIES, the total cost of the CIES is decreased by nearly ¥100. With the CESS, the wind abandonment rate is reduced by approximately 2%, so it can reduce the operating costs of the CIES and improve the system efficiency to a certain extent.

E. Impact of Thermal Resistance and Thermal Capacitance of Room

To explore the impact of different insulation materials on the total cost of the CIES, different thermal resistances and thermal capacitances are used in scenario 4 for a comparative analysis. In scenario 4, the thermal capacitance and thermal resistance are increased by 5% in turn. To analyze the effects of thermal capacitance and thermal resistance separately, the thermal capacitance is kept unchanged while the impacts of different thermal resistances are analyzed. Figure 16 shows the trends of total cost and abandonment rate under growth rate of different thermal resistances (R1 and R2 at the same rate at the same time as the growth rate).

Fig. 16  Trends of total cost and abandonment rate under growth rate of different thermal resistances.

With the increase of the growth rate of the thermal resistance, the total cost of the CIES shows a downward trend. The main reason is that the physical meaning of the thermal resistance is the performance of the insulation. Therfore, the greater the thermal resistance of the room, the better the performance of the thermal insulation. When the indoor temperature is low, the heating device releases heat to warm the room. The better the thermal insulation, the slower the indoor heat transfer to the outdoors, i.e., the slower the heat loss, which is reflected on the energy supply side. Thus, the demand for heat energy is reduced. The energy demand of the CIES is further reduced, i.e., the energy purchasing cost of the CIES is reduced. During high wind power output periods, the wind energy consumption is reduced, and the wind abandonment rate is increased. However, the reduction in the energy purchasing cost is more related to the increase in the wind abandonment penalty cost, so the total cost gradually decreases as the wind abandonment penalty cost increases.

Similarly, in the following analysis of the impact of different thermal capacities of the room on the operating cost of the CIES and the wind abandonment rate, the trends in total cost and wind abandonment rate under different thermal capacitances are shown in Fig. 17. The thermal capacity refers to the ability of the room to store heat energy, which depends on its mass, material, temperature, and other factors. The greater the thermal capacity of the room is, the more the heat which is needed to increase the indoor temperature. As the thermal capacity increases, the heat needed to increase the indoor temperature increases. Hence, the energy purchasing cost of the CIES is greater, and the amount of wind power accommodated during the wind abandonment period increases, resulting in a decrease in the wind abandonment rate. However, the higher unit energy purchasing cost results in an increase in the total cost of the CIES.

Fig. 17  Trends in total cost and wind abandonment rate under different room thermal resistances.

According to the above analysis, the parameters of the typical room have a particular impact on the economic operation of the CIES, so the thermal resistances and thermal capacitances are carefully considered in the selection of room insulation materials and CIES operation. This paper studies only a typical day in winter. Under high cooling demand on a typical summer day, distinct trends and their impacts become evident.

VI. Conclusion

This paper proposes an optimal scheduling model of an IES in low-carbon communities considering flexibility of resources and segmental control of SOFC, addressing the strong thermoelectric coupling of the CHP, the low utilization rate of energy storage, and underexploitation of flexibility of load-side resources. The conclusions are as follows.

1) When the SOFC/GT-CHP system is introduced, the CIES carbon trading cost is reduced by nearly ¥1020 relative to that of the conventional scenario. The SOFC array control strategy based on the segmental control method proposed in this paper improves the overall efficiency of the CHP, weakens the thermoelectricity coupling of the CHP, and reduces the energy purchasing cost of the CIES by nearly ¥360.

2) The CESS applied in this paper effectively improves the utilization of the storage capacity in natural gas and heat storage tanks, stores electric energy that cannot be utilized owing to system constraints, reduces the total cost, and lowers the wind abandonment rate.

3) In this paper, transferable, substitutable, and heat loads are considered, and the thermal inertia of buildings is considered on the load side. Community buildings during the heating period in northern cities in China are explored in detail as flexible resources. Rooms are classified by comparing the floors, orientations, and room areas of buildings to characterize the differences in user characteristics, which effectively improves the utilization rate of wind power of the CIES without affecting the user comfort, benefits both the user and the CIES operator, and reduces the total operating costs of the CIES. The total operating cost of the CIES increases inversely with the thermal capacitance and thermal resistance when different room thermal capacities and room thermal resistances are compared.

Appendix

Appendix A

TABLE AI  Parameters of Typical Room Construction
FloorFormPlacement

Room size

(m2)

Room volume

(m3)

External wall area

(m2)

External window area

(m2)

Top I Shaded-middle 9.9 24.0 17.8 2.7
II Sunny-center 1 15.6 37.7 23.5 3.1
III Sunny-center 2 11.3 27.4 20.1 3.5
IV Shaded-side 10.1 24.4 25.4 2.7
V Sunny-side 15.6 37.7 35.1 3.1
Middle I Shaded-middle 9.9 24.0 7.9 2.7
II Sunny-center 1 15.6 37.7 7.9 3.1
III Sunny-center 2 11.3 27.4 8.8 3.5
IV Shaded-side 10.1 24.4 15.4 2.7
V Sunny-side 15.6 37.7 19.5 3.1
Ground I Shaded-middle 9.9 24.0 17.8 2.7
II Sunny-center 1 15.6 37.7 23.5 3.1
III Sunny-center 2 11.3 27.4 20.1 3.5
IV Shaded-side 10.1 24.4 25.4 2.7
V Sunny-side 15.6 37.7 35.1 3.1

Fig. A1  Day-ahead forecast curves of loads and wind power.

TABLE AⅡ  Parameter Fitting Results of a Typical Room
FormFloorR1 (℃/W)R2 (℃/W)Cin (105 J/℃)Cwall (105 J/℃)
I Top 0.0038 0.044 2.1 113.0
Middle 0.0041 0.045 1.9 174.3
Ground 0.0042 0.038 1.7 213.8
II Top 0.0044 0.076 1.6 332.1
Middle 0.0036 0.077 2.2 186.9
Ground 0.0031 0.047 3.3 253.9
III Top 0.0033 0.067 2.7 246.0
Middle 0.0019 0.083 4.0 235.7
Ground 0.0039 0.043 1.9 179.8
IV Top 0.0037 0.052 2.6 185.8
Middle 0.0041 0.047 2.2 177.3
Ground 0.0023 0.028 4.0 220.2
V Top 0.0038 0.032 2.0 207.1
Middle 0.0036 0.034 3.0 371.0
Ground 0.0022 0.022 2.2 620.6

References

1

G. Chen, Y. Dong, Z. Liang et al., “Analysis and reflection on high-quality development of new energy with chinese characteristics in energy transition,” Proceedings of the CSEE, vol. 140, pp. 5493-5506, Aug. 2020. [Baidu Scholar] 

2

X. Guo, S. Lou, Z. Chen et al., “Stochastic optimal scheduling considering reserve characteristics of retrofitted combined heat and power plants,” International Journal of Electrical Power & Energy Systems, vol. 140, p. 108051, Sept. 2022. [Baidu Scholar] 

3

Y. Zhou, J. Wang, Y. Li et al., “A collaborative management strategy for multi-objective optimization of sustainable distributed energy system considering cloud energy storage,” Energy, vol. 280, p. 128183, Oct. 2023. [Baidu Scholar] 

4

Y. Qiu, S. Lu, H. Lu et al., “Flexibility of integrated energy system: basic connotation, mathematical model and research framework,” Automation of Electric Power Systems, vol. 46, no. 18, pp. 16-43, Sept. 2022. [Baidu Scholar] 

5

J. J. Fitzpatrick and R. Dunphy, “Comparison of the carbon footprints of combined heat and power (CHP) systems and boiler/grid for supplying process plant steam and electricity,” Chemical Engineering Science, vol. 282, p. 119303, Dec. 2023. [Baidu Scholar] 

6

J. Chen, Z. HU, J. Chen et al., “Optimal dispatch of integrated energy system considering ladder-type carbon trading and flexible double response of supply and demand,” High Voltage Engineering, vol. 47, pp. 3094-3106, Sept. 2022. [Baidu Scholar] 

7

K. Abouemara, M. Shahbaz, G. McKay et al., “The review of power generation from integrated biomass gasification and solid oxide fuel cells: current status and future directions,” Fuel, vol. 360, p. 130511, Mar. 2024. [Baidu Scholar] 

8

S. K. Effatpanah, H. R. Rahbari, M. H. Ahmadi et al., “Green hydrogen production and utilization in a novel SOFC/GT-based zero-carbon cogeneration system: a thermodynamic evaluation,” Renewable Energy, vol. 219, p. 119493, Dec. 2023. [Baidu Scholar] 

9

M. Bąkała, K. Martsinchyk, K. Motyliński et al., “Numerical analysis of natural gas, hydrogen and ammonia fueled solid oxide fuel cell based micro cogeneration units with anodic gas recirculation,” International Journal of Hydrogen Energy, vol. 52, pp. 952-964, Jan. 2024. [Baidu Scholar] 

10

J. Cai, Q. Xu, X. Yuan et al., “Configuration strategy of large-scale battery storage system orienting wind power consumption based on temporal scenarios,” High Voltage Engineering, vol. 45, pp. 993-1001, Mar. 2019. [Baidu Scholar] 

11

J. Li, Y. Fu, C. Li et al., “Economic optimal configuration of hybrid energy storage system for improving wind power consumption,” Power System Technology, vol. 44, no. 12, pp. 4547-4557, Dec. 2020. [Baidu Scholar] 

12

W. Huang, D. Liu, S. Li et al., “Two-level optimization of integrated energy system with P2G and demand response under dual carbon objective,” Electrical Measurement & Instrumentation, vol. 59, pp. 8-17, Nov. 2022. [Baidu Scholar] 

13

Y. Cui, F. Guo, X. Fu et al., “Source-load coordinated optimal dispatch of integrated energy system based on conversion of energy supply and use to promote wind power accommodation,” Power System Technology, vol. 46, no. 4, pp.1437-1447, Apr. 2022. [Baidu Scholar] 

14

P. Li, D. Wu, Y. Li et al., “Optimal dispatch of multi-microgrids integrated energy system based on integrated demand response and stackelberg game,” Proceedings of the CSEE, vol. 41, pp.1307-1321, Feb. 2021. [Baidu Scholar] 

15

Z. Yi and Z. Li, “Combined heat and power dispatching strategy considering heat storage characteristics of heating network and thermal inertia in heating area,” Power System Technology, vol. 42, no. 5, pp. 1378-1384, May 2018. [Baidu Scholar] 

16

Y. Cui, C. Gu, X. Fu et al., “Low-carbon economic dispatch of integrated energy system with carbon capture power plants considering generalized electric heating demand response,” Proceedings of the CSEE, vol. 42, pp. 8431-8446, Dec. 2022. [Baidu Scholar] 

17

Z. Liu. (2018, May). Research on the application of complementary heating system based on GSHP+ASHP+solar energy in northern towns. [Online]. Availbale: https://kns.cnki.net/kcms2/article/abstract?v=upp0mjVwiKUvY2j_RCRbE35CqiGyV17ngyNl8JttDckBfHcYWF64gGVpGB1FMnhkiLycVrGtIrTOeRWNVACSP7bFYrKt_ko56md0Z4n7_tPqCHPz4Rxw6BUfY14FFFiFzVze6zFYynQ73ung51KmSUDMDEGhBpR-jZfdi4fg5JclawaVVIYq0U31yYFBMyqLCBfdPTQxK40=&uniplatform=NZKPT&language=CHS [Baidu Scholar] 

18

H. Liu, Y. Wang, J. Li et al., “Coordinated heat and power dispatch of micro-energy network of countryside considering heat balance model of building and flexible indoor comfort constraints,” Automation of Electric Power Systems, vol. 43, no. 10, pp. 50-58, May 2019. [Baidu Scholar] 

19

H. Jia, Y. Lei, X. Jin et al., “Collaborative optimization of building aggregation and integrated community energy system considering difference in insulation performance,” Automation of Electric Power Systems, vol. 47, no. 24, pp.31-38, Dec. 2023. [Baidu Scholar]