Preprint
Article

Optimized Dispatch of Regional Integrated Energy System Considering Wind Power Consumption in Low-temperature Environment

Altmetrics

Downloads

51

Views

9

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

11 November 2023

Posted:

13 November 2023

You are already at the latest version

Alerts
Abstract
To address the problems of wind abandonment in the regional integrated energy system (RIES), which occurs when cogeneration unit operates in a heat-determined electricity operation mode, and the low efficiency of the energy storage system in a low-temperature environment, we propose an optimal scheduling method for RIES based on fine energy storage and wind power dissipation. First, a fine energy storage model is established on the source side and the conditional value at risk(CvaR) theory is used to quantify the uncertainty of wind power; then a combined heat and power demand response mechanism is introduced on the load side to reduce the peak-to-valley difference between heat and power loads, and to promote the consumption of wind power. Finally, with the objective of minimizing the total cost of RIES optimal dispatch, the example is solved on the MATLAB platform. The simulation results show that compared with the traditional model, the proposed model is not only more adaptable to the low-temperature environment, but also can effectively improve the system economy and wind power consumption.
Keywords: 
Subject: Engineering  -   Electrical and Electronic Engineering

1. Introduction

With the intensification of the energy crisis, increasing the utilization rate of renewable energy and reducing fossil energy dependence has become an important research direction. Integrated regional energy systems have rapidly developed as an effective way to increase the utilization of renewable energy and reduce pollutant emissions [1]. Cogeneration units centered on Combined, Heating and Power (CHP) technology are capable of generating both electricity and heat and are used in a wide range of applications. As CHP units are often in heat-determined power operation in winter, the heat load is at its peak during the nighttime when the electric load is in the trough, and the system peaking capacity is smaller than that of conventional units, making it more prone to the phenomenon of wind abandonment [2].
In an integrated energy system, energy storage devices change load distribution in the space-time dimension and improve system stability and economy. Reference [3] based on the RIES optimal dispatch model including ground source heat pumps and energy storage, verified the effect of ground source heat pumps and energy storage devices on system economy. Reference [4] proposed a day-ahead optimization scheduling method based on park RIES based on the multi-energy coupling and energy conversion characteristics of cold, heat, electricity and gas, combined with the energy storage characteristics of the system. Reference [5] analyzes the energy storage capacity and economics of RIES from the perspective of integrated demand response and uncertainty. References [6,7] proposed the possibility of using shared energy storage instead of traditional energy storage to realize low-cost wind power consumption. However, none of the above literature takes into account the differences in the characteristics of the energy storage devices in RIES, such as the effect of charging and discharging power of the battery at low temperatures and the effect of the heat loss of the thermal storage tanks by the ambient temperature. Therefore, it is particularly important to consider the characteristics of energy storage devices under low temperature and other conditions for accurate scheduling of integrated energy systems. The RIES addresses the phenomena of wind power and photovoltaic uncertainty and large load peak-to-valley differences on the source and load side, which leads to the problems of impeded wind power and photovoltaic consumption and high integrated cost of dispatch. On the energy supply side, the current mainstream methods for dealing with the wind power and photovoltaic uncertainty problem include stochastic optimization and robust optimization. Reference [8] coordinates the robustness and economy of the system by flexibly adjusting the robust parameters to improve the system's consumption of wind power. Reference [9,10] established a stochastic optimization model of an integrated energy system containing electric vehicles, and simulated wind power output based on the scenario reduction technique with improved probabilistic distance, which improved the accuracy of the scheduling results. References [11,12] address the feature that robust optimization is only affected by the worst-case scenarios of uncertain input parameters, and the final results obtained are too conservative, and propose to deal with the uncertainty set with conditional value at risk (CVaR) theory to obtain a more economically efficient optimization scheme. On the energy demand side, reference [13] developed a tariff-based demand response optimization model based on RIES demand side characteristics. Reference [14] established an integrated energy system containing electricity-to-gas and energy storage devices based on gas-electricity joint demand response, and established a demand response optimization model based on electricity tariffs based on the characteristics of the energy demand side of RIES. Reference [15] developed an integrated energy system economic dispatch model for thermoelectricity integrated demand response (IDR) with the objective of minimizing the wind abandonment rate. However, these literatures are single consideration of the influence of source-side or load-side characteristics on the system operation, and do not consider the source and load sides comprehensively, so there is still room for further research in source-load coordinated optimal scheduling.
Aiming at the above problems, this paper proposes a coordinated optimization model of RIES based on fine energy storage and wind power consumption. Firstly, a fine energy storage model more in line with the low-temperature environment is established, and consideration is given to increasing the thermal load to safeguard the available capacity of the storage battery and the charge/discharge speed, and to improve the utilization rate of the energy storage system. Secondly, to address the uncertainty of wind power output on the energy supply side, the Latin hypercube-synchronous back generation elimination method is used to obtain more accurate wind power output based on the historical data of wind power output, and the CVaR theory is used to quantify the risk of revenue and enhance the system's ability to consume wind power. At the same time, the introduction of heat and power demand response mechanism on the load side cuts down the peak-to-valley difference of the system load and further improves the system's level of wind power consumption. After comprehensively considering the energy conversion characteristics of the source and load sides and the operation characteristics of the energy storage device, the model is solved using MATLAB to obtain the wind abandonment rate and operation scheme under the lowest operation cost. The example results show that compared with the traditional model, the RIES model based on fine energy storage and wind power consumption can effectively improve the level of wind power consumption and reduce the operating cost of the system.

2. Integrated Energy System Architecture

Due to the cold winter climate in the north, there is a high demand for heating. At the same time, the region is rich in geothermal energy, wind energy, solar energy ,and other renewable energy sources, which provides a good foundation for the application of ground source heat pumps and CHP in the region [16,17]. Compared with the traditional RIES which is limited to the independent operation of a single form of energy, this paper establishes a RIES containing integrated energy storage (IES) equipment, wind turbines (WT), photovoltaic (PV), ground source heat pumps (GSHP) and CHP units, with the structure shown in Figure 1, which is able to utilize and synergize the outputs of different forms of energy in a more effective way. The power load demand is met by the superior grid, wind turbines, photovoltaic power plants and cogeneration units; the heat load demand is met by the cogeneration units and ground-source heat pumps, and the energy storage equipment includes heat storage tank (HST) and IES to provide support for the heat and power loads. Considering that the actual available capacity of the IES is greatly affected by the temperature in winter, the working efficiency of the IES is guaranteed by increasing part of the thermal load.

2.1. CHP System

CHP systems are capable of delivering both heat and power at the same time. It converts high-quality heat from the combustion of natural gas into electricity and absorbs the waste heat to meet the heat load by using a bromine chiller, which usually consists of a bromine chiller(BC) and a gas turbine(GT)[18], and is mathematically modeled as follows:
P B , h ( t ) = Q G ( t ) δ B λ B , h
Q G ( t ) = P G   T ( t ) ( 1 β G β l ) / β G
Where: P B , h ( t ) , P G T ( t ) , Q G ( t ) represent the heating power of BC and the electric power output and waste heat generated by GT at the moment t, respectively; δ B , λ B , h , βG, βl are, in order, the recovery rate of the flue gas of BC, the heat production coefficient of the flue gas of BC, the power generation efficiency of GT, and the heat loss coefficient of GT.

2.2. Ground Source Heat Pump

Compared with traditional heating equipment such as electric boilers, ground source heat pumps convert a large number of low-temperature heat sources in the shallow ground into high-quality energy to meet heating and cooling needs with a small electrical energy input. The mathematical model for its heat supply is as follows:
P lh , l ( t ) = λ l P lh ( t )
Where: Plh,l(t) and Plh(t) are the heating power and power consumption of GSHP at the moment t, respectively; λl is its heating efficiency.

2.3. Fine Energy Storage Equipment

As an important component of RIES, energy storage equipment mainly includes thermal storage tanks and batteries. These devices absorb and release energy, thus reducing system fluctuations, especially suppressing fluctuations due to the uncertainty of new energy output. The charging and discharging mathematical model of the energy storage system is as follows:
S IES ( t ) = ( 1 η IES ) S IES ( t 1 ) + ( M P in ( t ) λ in Δ t N P out ( t ) λ out Δ t )
M + N 1 ,   M { 0 , 1 } ,   N { 0 , 1 }
Where: SIES (t) is the remaining energy of the energy storage device at time t; P in ( t ) , P out ( t ) represent the energy input and output power of the energy storage device at time t, respectively; λ in , λ out represent the conversion efficiencies of the energy input and output of the energy storage device, respectively; M, N are Boolean variables that represent the state of the energy storage device, which is in a state of charging at M=1 and in a state of energy release at N=1; and η IES is the energy loss coefficient of the energy storage device.

2.3.1. Battery Storage Model

The traditional energy storage model does not consider the effect of ambient temperature on the battery, but the actual usable capacity is related to the surface temperature of the battery [19,20], which is mathematically expressed as follows:
S IESN ( t ) = S EIESN [ 1 + ξ a ( T a ( t ) T N ) ]
Where: for the standard state of the rated capacity of the battery; for the battery surface temperature at the moment of t; for the standard state of the temperature, taking the value of 25 ℃.
The surface temperature of the battery is determined by its heat generation and heat dissipation during the working period, and its mathematical expression is as follows
T a ( t ) = T a ( t 1 ) + Φ in ( t ) Φ out ( t ) K tep S ele Φ in , e ( t ) = λ IES [ M ele · P IES , in ( t ) + N ele · P IES , out ( t ) ] Δ t   Φ out , e ( t ) = K tep S ele ( T a ( t 1 ) T evr ( t ) )
Where: Φ in , e ( t ) , Φ out , e ( t ) are the heat generation and dissipation of the battery at the moment of t; P IES , in ( t ) , P IES , out ( t ) are the charging and discharging power of the battery at the moment of t; M ele , N ele indicates the working state of the battery, M ele =1 when the battery is in the charging state, N ele =1 when the battery is in the discharging state; T evr ( t ) for the ambient temperature at the moment of t; ξ a for the capacity temperature coefficient; S ele for the surface area of the battery; K tep for the heat convection coefficient; λ IES for the battery heat generation coefficient.
In addition, the service life of the battery is also affected by its charging and discharging rate. If the charge/discharge rate is too high, the service life will be seriously damaged. Therefore, the maximum charging and discharging power per hour should be limited to not more than 20% of the actual usable capacity [21], and its mathematical model is as follows:
P IES , in ( t ) 0.2 S IESN ( t ) / Δ t P IES , out ( t ) 0.2 S IESN ( t ) / Δ t
In order to improve the utilization of the battery in practical applications, this paper optimizes the problems that may occur in the energy storage process. Specific measures include the rational design of battery charging and discharging programs and artificially increasing the ambient temperature to increase the actual available capacity.

2.3.2. Model of Heat Storage Tank

The heat storage tank model consists of an atmospheric pressure hot water storage tank, a circulating pump and additional equipment such as a heat exchanger. In practice, due to the influence of the ambient temperature, there is a certain amount of heat loss in the heat storage tank, mainly the inner wall of the tank to transfer heat and the outer surface of the tank to the surrounding air to emit heat [22,23]. The mathematical model of heat storage tank storage and heat release is as follows:
S LES ( t ) = η wat m T LES ( t ) T LES ( t ) = T LES ( t 1 ) + Φ in , w ( t ) Φ out , w ( t ) K wat S wat Φ in , w ( t ) = H dis ( t ) ω dis H ch ( t ) ω ch Δ t Φ out , w ( t ) = K wat S wat ( T LES ( t 1 ) T evr ( t ) )
Where: S LES ( t ) is the heat storage capacity of the heat storage tank at time t; η wat is the specific heat capacity of water; m is the mass of hot water in the heat storage tank; K wat and S wat are the integrated thermal convection coefficient and surface area of the heat storage tank, respectively; Φ in , w ( t ) and Φ out , w ( t ) are the heat input and heat loss from the heat storage tank, respectively; H ch ( t ) and H dis ( t ) are the storage and discharge thermal power of the heat storage tank, respectively; ω ch and ω dis are the storage and discharge efficiencies of the heat storage tank, respectively.

3. Combined Heat and Power Demand Response Mechanism

3.1. Electric Load Demand Response

A price-based DR (PBDR) model is established based on the price-based demand response mechanism. Through PBDR, the demand side electricity consumption is regulated so as to effectively realize peak shifting and valley filling, and to reduce the peak pressure of the system [24,25]. The demand elasticity of consumer load is as follows:
M t 1 , t 2 = Δ K t 1 K 0 t 1 · C 0 t 2 Δ C t 2 M t 1 , t 2 0 , t 1 t 2 M t 1   , t 2 0 , t 1 = t 2
Where: is the amount of electric load at t1 before taking demand response; is the change of electric load at t1 after taking demand response; is the electric price at t2 before taking demand response; is the change of electric price at t2 after taking demand response; is the cross-elasticity coefficient when t1 ≠ t2; and is the self-elasticity coefficient when t1=t2. After the introduction of the PBDR model, the change amount of electric load is as follows:
  Δ K 1 / K 0 1   Δ K 2 / K 0 2 Δ K 24 / K 0 24 = M 11   M 124 M 21   M 224     M 241   M 2424 Δ C 1 / C 0 1 Δ C 2 / C 0 2   Δ C 24 / C 0 24
With the introduction of the PBDR model, the magnitude of the electrical load at time t is as follows:
K load , e t = K 0 t ( 1 + Δ K t / K 0 t )

3.2. Heat Load Demand Response

The heat load of a residential area consists of two parts: domestic hot water and heating heat load. The heating heat load is more elastic and accounts for a higher proportion in the total heat load, while domestic hot water is a rigid load with little change and accounts for a smaller proportion in the total heat load. In order to accurately assess the impact of heat load on wind power consumption in RIES, only analyzes and calculates the heating heat load [26]. The current outdoor temperature and the expression for the relationship between the required thermal power and the amount of change in the indoor temperature of the house are as follows:
P l , h ( t ) = N 0 R ( T i ( t + 1 ) e Δ t / τ T i ( t ) 1 e Δ t / η T o ( t ) ) η = R C 0
Where: P l , h ( t ) represents the thermal power required in the room at the moment t; N 0 is the number of heating households in the residential area; C 0 is the air heat capacity in the room; R is the equivalent thermal impedance of the house; T i ( t ) , T o ( t ) are the indoor and outdoor temperatures of the house at the moment t.
Due to the human body's ambiguity in perceiving temperature changes in the surrounding environment, when the indoor temperature changes in a slow and small range, it does not cause discomfort. In order to assess the degree of living comfort, reference [27] used the Predicted Mean Vote (PMV) as a quantitative indicator. This literature states that when γ PMV is between -0.5 and 0.5, the comfort level of the inhabitants is not affected. Based on the above, the constraints on indoor temperature are as follows:
T i , min T i ( t ) T i , max
Where: T i , max and T i , min are the maximum and minimum indoor temperatures, respectively.

4. RIES Coordinated Optimization Model Considering CVaR

4.1. Uncertain Scene Selection

In this paper, using historical wind power output data, the kernel density estimation method is used to establish the actual probability density distribution function of the prediction error of wind power output at each scheduling moment, and Latin hypercubic sampling is applied to obtain the set of prediction error samples at each scheduling time point of the system [28]. In order to improve the computational efficiency, this paper adopts the synchronized back generation elimination method to cut down the scenes and get the corresponding probability values. Let the set of scenes be W s = W s , 0 , W s , 1 , W s , t , W s , T , where W s , t is the outgoing value of scene s at time t, and the sum of the probability of occurrence of each scene is 1. The expressions for the distances between the different scenarios are as follows:
D ( W s , W j ) = [ t = 1 T ( W s ( t ) W j ( t ) ) 2 ] 1 2
In order to make the sample set better approximate the original scene set, the paper adopts the scene cutting technique to obtain multiple scene sets with corresponding probabilities by clustering similar scenes. At the same time, it is also necessary to satisfy the minimization of the probability distance between the cut scene sets. The mathematical expression is as follows:
min s J , j J T D ( W s , W j )

4.2. CVaR Overview

CVaR is developed from Value at Risk (VaR), which represents the worst possible loss that a portfolio can suffer at a certain confidence level [29]. Specifically, assuming that the decision variable is u, the uncertainty variable is v and its probability density function is R(v), the loss function of u is S(u, v); then when u is fixed, the loss function S(u, v) caused by the loss function S(u, v) does not exceed the threshold θ of the cumulative distribution function is as follows:
π ( u , θ ) = S ( u , v ) θ   R ( v ) d v
At a specific confidence level α, VaR represents the maximum loss of revenue that RIES may face in a future dispatch cycle, which is mathematically modeled as follows:
f VaR ( u ) = min { θ : π ( u , v ) α }
CVaR is a measure of the average return loss of RIES in the event that a given VaR value is exceeded in the system return. Its mathematical model can be expressed as follows:
f CVaR ( u ) = 1 1 α S u , v f VAR ( u )   S ( u , v ) R ( v ) d v
Considering that Equation (18) is more difficult to calculate CVaR directly, a simple calculation method is proposed as follows:
f CVaR ( u ) = min θ + 1 1 α S = 1 S max { S u , v θ , 0

4.3. Objective Function

Due to the uncertainty of wind turbine output, there is a potential risk in the economic dispatch scheme. In order to comprehensively consider the performance of the model in different operating environments, this paper adopts the conditional value at risk (CVaR) to quantify the risk of returns in various situations [30,31]. The RIES coordination optimization model is established with the objective function of minimizing the total cost of RIES operation, which is shown as follows:
min f = β f EX + ( 1 β ) f CVaR f EX = S = 1 S P S · f S f CVaR ( x ) = θ + 1 1 α S = 1 S P S · max { f S θ , 0 } f S = t = 1 T C w S ( t ) + C grid S ( t ) + C e S ( t ) + C fu S ( t ) + C wt S ( t )
Where: f is the total economic cost of system operation; fCVaR is the CVaR cost; fEX is the expected cost; θ is the auxiliary variable; α is the confidence level; β is the weighting coefficient; PS is the corresponding probability of scenario S; S is the number of scenarios; T is the period of scheduling; and C w S t , C grid S t , C e S t , C fu S t , C wt S t are the Operation and maintenance (O&M) cost, the cost of power interaction with the superior grid, the cost of environmental cost with the superior grid, the cost of gas purchased with the superior grid, and the cost of wind abandonment with the superior grid, at the time t under scenario S, respectively.
  • O&M costs;
  • The O&M costs consist of the operation and maintenance costs of controllable units, new energy units and energy storage equipment, which are expressed as follows:
    C w S ( t ) = i = 1 K P a   S ( t ) W a + j = 1 N 1 P b   S ( t ) W b + n = 1 N 2 P ES , n   S ( t ) W n
  • Where: P a   S ( t ) , P b   S ( t ) , P ES , n   S ( t ) are the output power per unit O&M cost of controllable unit a, renewable energy machine b and energy storage device n at moment t under scenario S, and W a , W b , W n are the unit O&M cost per unit of controllable unit a, renewable energy machine b and energy storage device n at moment t under scenario S, respectively.
  • Electricity Interaction Costs;
3.
The cost of electricity interaction represents the sum of the cost of electricity purchased by the RIES from the main grid and the benefit of electricity sold to the main grid, expressed as follows:
4.
C grid S ( t ) = L   g   ( t ) · max { P   d S ( t ) , 0 } L s   ( t ) · max { P d S ( t ) , 0 }
5.
Where: L   g   ( t ) and L s   ( t ) represent the prices at which the grid makes power purchases and sells power at moment t, respectively; P   d S ( t ) denotes the power on the contact line of the grid at moment t under scenario S, which is greater than 0 for power purchases and less than 0 for power sales.
  • Environmental costs;
6.
The environmental costs represent the GHG emissions emitted during the operation of the controllable units and the GHG treatment costs caused by purchasing electricity from the main grid, which are specified as follows:
C e S ( t ) = m = 1 M λ m [ E m P a S ( t ) + E m P d , b S ( t ) ]
7.
Where: P d , b S ( t ) is the power purchased by the system at moment t under scenario S; Em and λm are the share coefficient of the mth greenhouse gas and the treatment cost, respectively; and M is the type of greenhouse gas.
  • Wind abandonment costs;
8.
The cost of maintaining the wind turbine system during the period of wind abandonment is the cost of wind abandonment, which is expressed as follows:
C wt S ( t ) = P wt S ( t ) D wt
9.
Where: P wt S ( t ) and D wt are the abandoned wind power and unit abandoned wind cost at moment t under scenario S, respectively.
  • Cost of gas purchases;
10.
The cost of purchasing gas for the CHP system is the cost of purchasing gas, which is specifically expressed as:
11.
C fu S ( t ) = P a S ( t ) β G · Δ t L CH 4 D CH 4
12.
Where: L CH 4 and D CH 4 are the average calorific value and gas price of the gas, respectively.

4.4. Restrictive Condition

Considering the operating parameters of the various devices, the RIES needs to satisfy the following constraints.
  • Electrical and thermal power balance constraints;
    K   load , e   ( t ) + P   lh S ( t ) = P   GT S ( t ) + P WT S ( t ) + P   d S ( t ) + P   PV S ( t ) + P   IES S ( t ) P   l , h S ( t ) = P BC , h S ( t ) + P   lp , h S ( t ) + H in S ( t )
Where: K   load , e   ( t ) and P   l , h S ( t ) are the power of the system electrical and thermal loads, respectively.
  • Controllable unit operating constraints;
    P a , min P   a S ( t ) P a , max λ i , min Δ t P   a S ( t ) P   a S ( t 1 ) λ i , max Δ t
Where: λ i , max , λ i , min are the upper and lower limits of the climbing rate of controllable unit i, respectively; P a , max , P a , min are the upper and lower limits of the output force of controllable unit a, respectively.
  • Interactive power constraints with the higher grid;
    P d , min P   d S ( t ) P d , max
Where: P d , max , P d , min are the maximum and minimum values of grid interaction power.
  • Wind power output constraints
    0 P WT S ( t ) P ut S ( t ) P WT S ( t ) P ut S ( t ) h ut
Where: P WT S ( t ) and P ut S ( t ) denote the actual and predicted power of wind power at moment t under scenario S, respectively, and h ut is the minimum rate of wind power consumption.
  • Energy storage device constraints
    P   min S ( t ) P   IES S ( t ) P     max S ( t ) E   min S ( t ) E   IES S ( t ) E max S ( t ) E   IES S ( 0 ) = E   IES S ( 24 )
Where: P     max S ( t ) and P   min S ( t ) are the maximum and minimum values of the transmission power of the energy storage device at moment t under scenario S within the constraints; E max S ( t ) and E   min S ( t ) are the maximum and minimum available energy storage capacity of the energy storage device at moment t under scenario S.
  • Electric load demand response constraints
    t = 1 24 K up S ( t ) = t = 1 24 K down S ( t ) 0 K up S ( t ) δ up · K   load , e   ( t ) 0 K up S ( t ) δ down · K   load , e   ( t )
Where: K up S ( t ) and K down S ( t ) denote the amount of load transfer in and out at moment t under scenario S, respectively; δ up and δ down are the proportions of the maximum upward and maximum downward load transfers, respectively.

5. Calculus Analysis

5.1. Parameterization

In order to verify the RIES model in the paper in terms of the economy of system operation and scheduling and the timeliness of wind power consumption, the paper sets up the simulation of the RIES coordinated optimization model in the following four scenarios for verification. Scenario 1 is the traditional RIES model without considering demand response and battery heating; Scenario 2 is the RIES model with fine energy storage, without considering demand response and battery heating; Scenario 3 is the RIES model with fine energy storage, considering demand response and battery heating; Scenario 4 is the RIES model with fine energy storage, considering demand response and battery heating; Scenario 4 is the RIES model with fine energy storage, considering demand response and battery heating; Scenario 4 is the RIES model with fine energy storage, considering demand response and battery heating; Scenario 4 is the RIES model with fine energy storage, considering demand response and battery heating; where scenario 4 is the model proposed in the paper. The model parameters are taken as confidence level α=0.85 and weighting factor β=0.9. The parameters of the energy storage system are shown in Table 1, and the pollutant treatment cost and unit operation parameters are shown in Table 2 and Table 3. The forecast data of PV, wind power and electric heat loads are shown in Figure 2, and the temperature and the purchase and sale price of electricity are shown in Figure 3. Taking the economy of system optimization operation and the timeliness of wind power consumption as the goal, the MATLAB platform is used to build the RIES coordinated optimization operation model and call the CPLEX solver to solve the model.
We introduce the Latin Hypercubic Sampling method to deal with the uncertainty introduced by wind power output. The method can quickly generate 1000 typical wind power scenarios and reduce them to 5 typical scenarios by synchronized back generation elimination. Its sample space and reduced scenarios are shown in Figure 4 and Figure 5, and the probability sum after the reduced scenarios is 1. It is assumed that the wind turbine output day-ahead prediction error is within 20%.

5.2. Optimized Operation Results and Analysis

The optimal scheduling results of RIES are shown in Figure 6, Figure 7, Figure 8 and Figure 9: The scheduling results of Scenario 1 are shown in Figure 6, at this time, the electric load of RIES mainly consists of electricity purchased and sold between CHP, wind power, PV, storage battery and the main grid. During the hours of 1:00-7:00 and 23:00-24:00, the electricity price is in the trough, so the CHP units and the ground source heat pumps are operating in the heat-determined electricity operation mode. However, limited by the transmission power of the contact line, there is wind abandonment in the system, resulting in energy waste.
Scenario 2 considers the influence of low temperature environment on the available capacity and energy transfer rate of the energy storage system on the basis of Scenario 1, and its scheduling results are shown in Figure 7. As the available capacity and charging/discharging rate of the storage battery are affected, the actual available capacity of the storage system is limited to be insufficient in the 1:00-7:00 time period, and it is difficult to consume the wind power in this time period, which leads to a significant increase in wind abandonment phenomenon, the peak-to-valley difference of the system increases by 2.33% compared with Scenario 1, and the wind abandonment rate increases by 4.12%.
Scenario 3 introduces a combined heat and power demand response mechanism on the basis of Scenario 2. The use of electric load demand response mechanism regulates the user's power consumption habits, optimizes the distribution of power consumption, makes the power consumption increase in the price of electricity in the trough time, and reduces the power consumption in the price of electricity in the peak time; effectively reduces the peak-valley difference of the electric load; the scheduling results of the system are shown in 8, in the time period of 1:00-7:00, the wind power consumption rate of the system increases by 5.95% compared to Scenario 2, and the peak-valley difference decreases by 5.04%.
Scenario 4 adds a battery insulation device on the basis of Scenario 3. The battery is always in the best working condition by adding part of the thermal load. The scheduling results are shown in Figure 9, the system achieves complete consumption of wind power from 1:00 to 7:00, and the peak-to-valley difference of the system is reduced by 10.03% compared with Scenario 3.
The results of the optimized scheduling of heat load and the comprehensive cost of each scenario in this process are shown in Figure 10 and Table 4. Scenario 2 takes into account that the heat dissipation of the heat storage tank in the low-temperature environment increases significantly, and compared with Scenario 1 Scenario 2 the heat load increases by 4.71% due to heat loss, and the total cost increases by 1.53%.
Scenario 3 introduces a heat load response mechanism to regulate the indoor temperature on the basis of scenario 2, which reduces the demand for heat load without affecting the comfort of the residents. Its total cost is reduced by 1.74% compared to Scenario 2, but the heat load is only reduced by 0.54%, which is due to the fact that the power purchase price is higher than the generation cost of the CHP unit from 7:00-11:00 and 19:00-22:00, and the room temperature is increased in order to reduce the amount of purchased power from the grid and increase the generation of the CHP unit, which leads to an insignificant effect of the optimization of the heat load in this scenario.
Scenario 4 increases the thermal load of the system compared to Scenario 3 due to the additional consumption of part of the thermal load by the storage heater, but the total cost is reduced by 2.58% compared to Scenario 2 by realizing the complete consumption of wind power, which effectively improves the comprehensive benefit of the system.

5.3. Analysis of the Effect of Weighting Coefficients on the Results of RIES Scheduling

The weighting coefficient can reflect the degree of risk preference of the decision maker, in order to explore the influence of the weighting coefficient β on the scheduling results of RIES, the paper derives the scheduling results under different weighting coefficients under scenario 4 according to the change of weighting coefficients from 0.1 to 0.9 as shown in Figure 11. The results show that in the process of gradually increasing weight coefficients, the total system cost gradually increases, the expected cost share gradually increases, and the CVaR share gradually decreases. This indicates that the greater the degree of risk preference of the decision maker, the smaller the CVaR percentage, the more it can improve the operational safety of the system, but at the same time, it will reduce the unity economy. In actual scheduling, operators should choose appropriate weighting coefficients to balance the benefits and risks of the system according to specific needs.

6. Conclusions

By establishing a coordinated and optimized operation model of RIES considering wind power consumption in low-temperature environment, the following conclusions can be drawn based on the comparative analysis of the four scenarios:
(1)
The fine energy storage model takes into account the influence of ambient temperature, capacity and other constraints, and the way of equipment output is more in line with the actual situation than the traditional energy storage model, and the reliability of system optimization and scheduling is higher.
(2)
Comprehensive consideration of electricity and heat flexible load demand response, can effectively reduce the peak and valley difference of the user load, while taking into account the wind power consumption capacity and economic benefits.
(3)
Compared to traditional deterministic models, it is more reasonable to utilize CVaR theory to describe the risk of returns from wind power uncertainty.
(4)
The adjustment of the weight coefficients in the CVaR theory according to the historical operation information and the actual situation during the actual scheduling process can further improve the comprehensive efficiency of the system.

Author Contributions

Conceptualization, L.L. and J.H.; methodology, L.L.; software, L.L.; validation, J.H. and Z.L.; formal analysis, L.L.; investigation, H.Q.; resources, H.Q.; data curation, H.Q.; writing—original draft preparation, L.L.; writing—review and editing, L.L.; visualization, L.L.; supervision, J.H.; project administration, Z.L.; funding acquisition, Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China, grant number 52077120.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank the hard work of all anonymous reviewers to improve the quality of this submission.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, Y.; Wang, Y.; Huang, Y.; Li, F.; Zeng, M.; Li, J.; Wang, X.; Zhang, F. Planning and operation method of the regional integrated energy system considering economy and environment. Energy 2019, 171, 731–750. [Google Scholar] [CrossRef]
  2. Wang, Y.; Lu, Y.; Ju, L.; Wang, T.; Tan, Q.; Wang, J.; Tan, Z. A multi-objective scheduling optimization model for hybrid energy system connected with wind-photovoltaic-conventional gas turbines, CHP considering heating storage mechanism. Energies 2019, 12, 425. [Google Scholar] [CrossRef]
  3. Wang, J.; Deng, H.; Qi, X. Cost-based site and capacity optimization of multi-energy storage system in the regional integrated energy networks. Energy 2022, 261, 125240. [Google Scholar] [CrossRef]
  4. Xu, J.; Wang, X.; Gu, Y.; Ma, S. A data-based day-ahead scheduling optimization approach for regional integrated energy systems with varying operating conditions. Energy 2023, 283, 128534. [Google Scholar] [CrossRef]
  5. Guo, Z.; Zhang, R.; Wang, L.; Zeng, S.; Li, Y. Optimal operation of regional integrated energy system considering demand response. Appl. Therm. Eng. 2021, 191, 116860. [Google Scholar] [CrossRef]
  6. Zhang, S.; Miao, S.; Li, Y.; Yin, B.; Li, C. Regional integrated energy system dispatch strategy considering advanced adiabatic compressed air energy storage device. Int. J. Electr. Power Energy Syst. 2021, 125, 106519. [Google Scholar] [CrossRef]
  7. Zhu, H.; Li, H.; Liu, G.; Ge, Y.; Shi, J.; Li, H.; Zhang, N. Energy storage in high renewable penetration power systems: Technologies, applications, supporting policies and suggestions. CSEE J. Power Energy Syst. 2020. [Google Scholar]
  8. Zhang, Y.; Huang, Z.; Zheng, F.; Zhou, R.; Le, J.; An, X. Cooperative optimization scheduling of the electricity-gas coupled system considering wind power uncertainty via a decomposition-coordination framework. Energy 2020, 194, 116827. [Google Scholar] [CrossRef]
  9. Sun, G.; Li, Y.; Chen, S.; Wei, Z.; Chen, S.; Zang, H. Dynamic stochastic optimal power flow of wind power and the electric vehicle integrated power system considering temporal-spatial characteristics. J. Renew. Sustain. Energy 2016, 8. [Google Scholar] [CrossRef]
  10. Ju, L.; Li, H.; Zhao, J.; Chen, K.; Tan, Q.; Tan, Z. Multi-objective stochastic scheduling optimization model for connecting a virtual power plant to wind-photovoltaic-electric vehicles considering uncertainties and demand response. Energy Convers. Manag. 2016, 128, 160–177. [Google Scholar] [CrossRef]
  11. Li, P.; Yu, D.; Yang, M.; Wang, J. Flexible look-ahead dispatch realized by robust optimization considering CVaR of wind power. IEEE Trans. Power Syst. 2018, 33, 5330–5340. [Google Scholar] [CrossRef]
  12. Levy, D.; Carmon, Y.; Duchi, J.C.; Sidford, A. Large-scale methods for distributionally robust optimization. Adv. Neural Inf. Process. Syst. 2020, 33, 8847–8860. [Google Scholar]
  13. Wang, Y.; Li, F.; Yang, J.; Zhou, M.; Song, F.; Zhang, D.; Xue, L.; Zhu, J. Demand response evaluation of RIES based on improved matter-element extension model. Energy 2020, 212, 118121. [Google Scholar] [CrossRef]
  14. Zhao, D.; Xia, X.; Tao, R. Optimal configuration of electric-gas-thermal multi-energy storage system for regional integrated energy system. Energies 2019, 12, 2586. [Google Scholar] [CrossRef]
  15. Yang, L.; Han, Q.; Li, X. Economic Dispatch of Multi Region Electric and Heat Energy System Using Two Stage Demand Response for Better Integration of Wind Power. J. Electr. Eng. Technol. 2022, 17, 2553–2563. [Google Scholar] [CrossRef]
  16. Zhang, Q.; Zhang, L.; Nie, J.; Li, Y. Techno-economic analysis of air source heat pump applied for space heating in northern China. Appl. Energy 2017, 207, 533–542. [Google Scholar] [CrossRef]
  17. Cai, B.; Li, H.; Hu, Y.; Liu, L.; Huang, J.; Lazzaretto, A.; Zhang, G. Theoretical and experimental study of combined heat and power (CHP) system integrated with ground source heat pump (GSHP). Appl. Therm. Eng. 2017, 127, 16–27. [Google Scholar] [CrossRef]
  18. Matjanov, E. Gas turbine efficiency enhancement using absorption chiller. Case study for Tashkent CHP. Energy 2020, 192, 116625. [Google Scholar] [CrossRef]
  19. Fatullah, M.A.; Rahardjo, A.; Husnayain, F. Analysis of discharge rate and ambient temperature effects on lead acid battery capacity. In Proceedings of the 2019 IEEE International Conference on Innovative Research and Development (ICIRD); 2019; pp. 1–5. [Google Scholar]
  20. Low, W.Y.; Abdul Aziz, M.J.; Idris, N.R.N. Modelling of lithium-titanate battery with ambient temperature effect for charger design. IET Power Electron. 2016, 9, 1204–1212. [Google Scholar] [CrossRef]
  21. Ruetschi, P. Aging mechanisms and service life of lead–acid batteries. J. Power Sources 2004, 127, 33–44. [Google Scholar] [CrossRef]
  22. Arslan, M.; Igci, A.A. Thermal performance of a vertical solar hot water storage tank with a mantle heat exchanger depending on the discharging operation parameters. Sol. Energy 2015, 116, 184–204. [Google Scholar] [CrossRef]
  23. Cabeza, L.F. Advances in thermal energy storage systems: Methods and applications. In Advances in Thermal Energy Storage Systems; Elsevier: 2021; pp. 37–54.
  24. Yang, C.; Meng, C.; Zhou, K. Residential electricity pricing in China: The context of price-based demand response. Renew. Sustain. Energy Rev. 2018, 81, 2870–2878. [Google Scholar] [CrossRef]
  25. Shewale, A.; Mokhade, A.; Funde, N.; Bokde, N.D. An overview of demand response in smart grid and optimization techniques for efficient residential appliance scheduling problem. Energies 2020, 13, 4266. [Google Scholar] [CrossRef]
  26. Wang, Y.; Huang, Y.; Wang, Y.; Zeng, M.; Yu, H.; Li, F.; Zhang, F. Optimal scheduling of the RIES considering time-based demand response programs with energy price. Energy 2018, 164, 773–793. [Google Scholar] [CrossRef]
  27. Yau, Y.; Chew, B. A review on predicted mean vote and adaptive thermal comfort models. Build. Serv. Eng. Res. Technol. 2014, 35, 23–35. [Google Scholar] [CrossRef]
  28. Gao, B.; Pan, Y.; Chen, Z.; Wu, F.; Ren, X.; Hu, M. A spatial conditioned Latin hypercube sampling method for mapping using ancillary data. Trans. GIS 2016, 20, 735–754. [Google Scholar] [CrossRef]
  29. Stoyanov, S.V.; Rachev, S.T.; Fabozzi, F.J. Sensitivity of portfolio VaR and CVaR to portfolio return characteristics. Ann. Oper. Res. 2013, 205, 169–187. [Google Scholar] [CrossRef]
  30. Dong, L.; Li, M.; Hu, J.; Chen, S.; Zhang, T.; Wang, X.; Pu, T. A hierarchical game approach for optimization of regional integrated energy system clusters considering bounded rationality. CSEE J. Power Energy Syst. 2023. [Google Scholar]
  31. Li, P.; Yang, M.; Wu, Q. Confidence interval based distributionally robust real-time economic dispatch approach considering wind power accommodation risk. IEEE Trans. Sustain. Energy 2020, 12, 58–69. [Google Scholar] [CrossRef]
Figure 1. Basic structure of RIES.
Figure 1. Basic structure of RIES.
Preprints 90251 g001
Figure 2. Predictive output of wind, photovoltaic, and electrical heating loads.
Figure 2. Predictive output of wind, photovoltaic, and electrical heating loads.
Preprints 90251 g002
Figure 3. Temperature and purchase-sale price.
Figure 3. Temperature and purchase-sale price.
Preprints 90251 g003
Figure 4. Wind Power Sample Space.
Figure 4. Wind Power Sample Space.
Preprints 90251 g004
Figure 5. Wind Powe landscape reduction.
Figure 5. Wind Powe landscape reduction.
Preprints 90251 g005
Figure 6. Scenario 1 Scheduling Results.
Figure 6. Scenario 1 Scheduling Results.
Preprints 90251 g006
Figure 7. Scenario 2 Scheduling Results.
Figure 7. Scenario 2 Scheduling Results.
Preprints 90251 g007
Figure 8. Scenario 3 Scheduling Result.
Figure 8. Scenario 3 Scheduling Result.
Preprints 90251 g008
Figure 9. Scenario 4 Scheduling Result.
Figure 9. Scenario 4 Scheduling Result.
Preprints 90251 g009
Figure 10. Heat load scheduling results for four scenarios.
Figure 10. Heat load scheduling results for four scenarios.
Preprints 90251 g010
Figure 11. Expected Cost and CVaR under Different Weighting Coefficients.
Figure 11. Expected Cost and CVaR under Different Weighting Coefficients.
Preprints 90251 g011
Table 1. Parameters of energy storage device.
Table 1. Parameters of energy storage device.
Parameter Batteries Heat storage tanks Parameter Heat storage tanks Batteries
Capacity/(kW·h) 200 300 Initial energy storage state 0.2 0.2
Charging and discharging rate 0.9 0.88 Maximum energy storage state 0.9 0.9
attrition rate 0.001 0.01 Minimum energy storage state 0 0.2
O&M unit price (USD/kW·h) 0.051 0.045 Maximum Charge and Discharge Power/kW 50 50
Table 2. Pollutant discharge and treatment costs.
Table 2. Pollutant discharge and treatment costs.
Type SO2 NOX CO2
Gas turbine emission standard (g/kWh) 0.023 4.795 170.16
Emission standard of purchased power (g/kWh) 6.4 2.32 696
Treatment cost of each pollutant(USD/t) 1000 1950 9.75
Table 3. Pollutant discharge and treatment costs.
Table 3. Pollutant discharge and treatment costs.
Parameter GT GSHP PW PV grid
Power upper limit/kW 500 30 400 100 100
Lower power limit/kW 50 0 0 0 0
Climbing rate upper limit/(kW/min) 6 4
Climbing rate lower limit/(kW/min) 5 3
Efficiency 0.24 3
O&M unit price (USD/kW·h) 0.053 0.026 0.029 0.025
Table 4. Four Scenario Scheduling Results.
Table 4. Four Scenario Scheduling Results.
Cost Category(USD/Day) Scenario1 Scenario2 Scenario3 Scenario4
Total Cost 5620.5 5706.2 5646.7 5502.5
Expected Cost 5558.8 5648.5 5548.3 5442.4
CVaR 6175.8 6234.4 6208.2 6043.5
Fuel Cost 4292.3 4319.3 4275 4266.9
Maintenance Cost 446.7 441.1 443.7 455.3
Environmental Costs 56.57 58.67 56.64 54.41
Purchase and Sale Costs 766.5 833.8 812.5 665.8
Wind Abandonment Cost 3.96 6.35 3.15 0
Wind Power Consumption Rate 95.54% 91.42% 97.37% 0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated