Next Article in Journal
The Schaake Shuffle Technique to Combine Solar and Wind Power Probabilistic Forecasting
Previous Article in Journal
Bridging Tools to Better Understand Environmental Performances and Raw Materials Supply of Traction Batteries in the Future EU Fleet
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Two-Stage Algorithm of Locational Marginal Price Calculation Subject to Carbon Emission Allowance

1
Guangdong Power Exchange Center Co., Guangzhou 510080, China
2
School of Electric Power Engineering, South China University of Technology, Guangzhou 510640, China
*
Author to whom correspondence should be addressed.
Energies 2020, 13(10), 2510; https://doi.org/10.3390/en13102510
Submission received: 9 March 2020 / Revised: 8 May 2020 / Accepted: 14 May 2020 / Published: 15 May 2020
(This article belongs to the Section C: Energy Economics and Policy)

Abstract

:
To analyze the effect of carbon emission quota allocation on the locational marginal price (LMP) of day-ahead electricity markets, this paper proposes a two-stage algorithm. For the first stage of the algorithm, a multi-objective optimization model is established to simultaneously minimize the total costs and carbon emission costs of power systems. Hence, an evenly distributed Pareto optimal solution can be solved effectively by means of the normalized normal constraint method. For the second stage, a tracing model is built with the goal of minimizing the total costs of power systems and satisfying the constraints generated based on the Pareto optimal solution obtained from the first stage. Furthermore, the influence of carbon emission quota allocation on the LMP of electricity markets is analyzed, and different schemes to allocate carbon emission quotas are evaluated on a real 1560-bus and 52-unit system.

Graphical Abstract

1. Introduction

With the rapid development of the modern industrial system and the overconsumption of resources and energy, huge amounts of greenhouse gases have been emitted into the atmosphere over the last 20 years, which has brought about environmental pollution and damage to the ecosystem [1,2,3]. Therefore, numerous industrialized countries have promoted international actions to reduce carbon emissions and achieve sustainable development [4,5]. In the European Union, many countries have committed to reducing carbon emissions by 20% by 2025, and 45% by 2050, in order to limit global temperature rise [6]. As a major contributor of carbon emissions, China set the goal of reducing carbon emissions by 40%–50% by 2020 compared to 2005 [7,8]. Within this context, in addition to vigorously developing the renewable energy industries, many countries also utilize financial incentives such as carbon emission trading (CET) to encourage emission reduction. China has already started CET pilots in seven provinces or cities, including Guangdong Province, and planned to further promote CET mechanisms across the country. Moreover, the government executed the pilot operation of provincial spot electricity markets in Guangdong, Mengxi, Zhejiang, Shanxi, Shandong, Fujian, Sichuan, and Gansu at the end of 2018.
According to studies on greenhouse gas emissions, the carbon emissions of the power generation industry by fossil fuel combustion account for a very high proportion of the total emissions in many countries, especially in developing countries. In China, over 40% of carbon emissions are produced by fossil fuel combustion during the process of power generation [9]. Therefore, the implementation of the CET mechanism has a significant impact on the development of the power industry, and it is significant to study how the CET mechanism and the carbon emission quota allocation mechanism will impact the electricity market in China at this stage, especially the marginal price of the spot market.
To date, numerous valuable studies have analyzed the impact of CET mechanisms on the electricity market. Wang et al. proposed a two-stage scheduling model to comprehensively investigate the environmental benefits of consumers participating in both electricity and CET markets [9]. Zhou et al. aimed to analyze these potential profit impacts and the possible compensation to affected generation companies through modeling the Australian National Electricity Market under a potential emission trading scheme [10]. In [11], a model was proposed to determine equilibrium fuel, power, and emission trades in an electric power supply chain framework with a load-based emission trading program. In [12], a novel model was proposed to optimize the decision-making process in the interrelated electricity, fuel, and carbon markets in order to maximize the GENCO’s profit. In [13,14], the short-term unit combination problem with carbon emission constraints was investigated, and a solution was developed in which the carbon emission costs will be added to the system operation costs when the total carbon emissions of a system exceed its allocated quota. Based on the classic Cournot equilibrium model, Li et al. modeled different types of units separately and established a bi-level programming model to study the impact of carbon emission costs on the equilibrium state of the market and the strategic behavior of market members [15]. Each market member can achieve the maximum benefit in this case. However, the marginal price of the electricity market can be obtained only when the independent system operator (ISO) figures out the inverse demand function of the market in advance, which is inconvenient in practice. To analyze the impact of different carbon emission quota allocation schemes on the unit combination results, system operation costs, and changes in the benefits structure, Wang et al. established a security-constrained unit combination model considering the allocation of carbon emission quotas [16]. However, it ignored the conflicts of interest between different entities of the market and failed to make a scheduling result that could take into account the different entities. Therefore, the model could not truly reflect the impact of carbon emission quotas on the marginal price of the electricity market.
This paper proposes a two-stage algorithm for calculating the locational marginal price (LMP) in the day-ahead market within the CET mechanism to analyze the impact of carbon emission allowance on the day-ahead electricity market. First, a multi-objective optimization model to minimize both system operation costs and carbon emission costs is proposed and solved by the normalized normal constraint (NNC) method. Next, a tracing model with the goal of minimizing the total costs of the power system and satisfying the constraints generated by the Pareto optimal solution is built to calculate the LMP. Finally, the proposed algorithm is applied to a real 1560-bus and 52-unit system. The impacts of different carbon emission allocation schemes on the LMP and other clearing results of the day-ahead market are analyzed. The results of the simulation show that the proposed algorithm can balance the conflicts between different generation units and reduce the total carbon emissions of the system.
The remainder of this paper is organized as follows. In Section 2, we introduce the framework and assumptions. In Section 3, we establish the two-stage model for calculating the LMP. The solution methodology is proposed in Section 4. Numerical results are presented in Section 5, and conclusions are drawn in Section 6.

2. Market Framework and Model Assumptions

2.1. Market Framework

Figure 1 shows the transaction framework in which generation units (i.e., conventional units, renewable units, and pumped storage units) participate in both the CET market and the day-ahead electricity market. During electricity generation, conventional units emit greenhouse gases, which are monitored by the government. Additionally, the government is responsible for both the allocation of annual carbon emission quotas to generation units at the beginning of the year and the verification of total carbon emissions to those units at the end of the year.
The interaction between the CET market and day-ahead electricity market can be summarized as follows. Each generation unit offers price and capacity bidding to ISO daily within the day-ahead electricity market. Conventional units trade emission quotas in the CET market quarterly. Hence, the carbon price is relatively fixed for a certain operational day of the day-ahead electricity market. Conventional units are able to trade their surplus quotas with each other to obtain profits; however, they need to pay for carbon emission costs when their emissions exceed their prescribed free allocation quota, which will increase their generation costs. Therefore, they will evaluate the predicted carbon emission costs according to the allocated quotas obtained from the government before bidding in the day-ahead electricity market.
According to the interaction mentioned above, the CET market does have a unilateral influence on the day-ahead market within the CET mechanism, and the details are given in Section 3.

2.2. Model Assumptions

2.2.1. Carbon Emission Quota Allocation Approaches

Currently, the government allocates carbon emission quotas for various industries annually. In order to study the impact of carbon emission quotas on the marginal price of day-ahead electricity markets, this paper modifies the time scales of the quota allocation from years to days and assumes that each conventional unit can only obtain quotas from the government. The quotas are unable to deal between different units.
Before allocating carbon emission quotas for conventional units, it is necessary to specify the total amount of carbon emissions available for quota allocation. Generally, we set the benchmark by solving the unit commitment (UC) problem with CET and calculating the total carbon emissions of conventional units. Then, the total emission quota with CET can be expressed as:
Q = ( 1 α ) E b ,
where Q is the total daily quota available to allocate; Eb represents the total daily emission of conventional units without CET; and α denotes the emission reduction factor.
To allocate the quotas to generating units, two concrete approaches are proposed [10], including the historical carbon emission-based method and the power generation performance-based method.
(1) The historical carbon emission-based method:
The quota Eq,i of unit i is allocated based on the historical emissions quantities without CET, which can be expressed as:
E q , i = E i b E b Q ,
where E i b is the carbon emissions of unit i without CET.
(2) The power generation performance-based method:
Power generation performance G is the amount of CO2 emitted by conventional unit i for producing one unit of electricity at time t, which is expressed as:
G = Q i = 1 N 1 t = 1 24 P i , t ,
where Pi,t is the output of unit i at time t, and N1 represents the total number of conventional generating units.
Hence, the quota Eq,i of the unit i can be defined as:
E q , i = G t = 1 24 P i , t .

2.2.2. Stepwise Bidding of Conventional Units

In this paper, we assume that the generating units bid and offer in the day-ahead electricity market, while the consumers only declare quantity in the demand side. Without CET, the quotation curve of a conventional unit is segmented, as shown in Figure 2, where Pi,min and Pi,max represent the lower and upper output limits of unit i, respectively. The quotation of unit i in the kth step at time t can be written as:
C i . t k = a i ( P i , t k ) 2 + b i P i , t k + c i ,
where ai, bi, and ci denote the production cost coefficients, and P i , t k represents the bidding quantities corresponding to the quotation of unit i in the kth step at time t, which is written as:
P i , t k = { P i , min + L i P i , min + 2 L i P i , min + 3 L i P i , min + 4 L i k = 1 , 2 k = 3 k = 4 k = 5 ,
where Li represents the length of the multi-section quotation quantity of conventional unit i.
In general, the sum of the winning quantity for each segment must be equal to the total output of the unit i in Figure 1, that is,
P i , t = k = 1 K P i , t k ,
where P i , t k represents the capacity offer in the kth step by conventional unit i at time t. Each segment of quantity needs to meet the upper and lower limits of the segment, that is,
{ P i , t k = P i , min I i , t , k = 1 0 P i , t k L i I i , t , 1 < k K .
Considering the mechanism of CET, conventional unit i has to pay its carbon emission costs when its emissions exceed its allocation quota. Therefore, the total costs of generating unit i will rise, and then unit i may adjust its online quotation strategy by adding the predicted carbon emission costs, corresponding to its quotas, in each segment of the declared prices, that is,
C i . t k = a i ( P i , t k ) 2 + b i P i , t k + c i + ( t = 1 24 β c i P i , max η E q , i ) × p r e t = 1 24 P i , max ,
where βci represents the emission coefficient of conventional unit i, and pre is the carbon emission trading price, which we set to 15 yuan/ton.

3. Two-Stage Model of Locational Marginal Price Calculation with Carbon Emission Trading (CET)

To analyze the effect of carbon emission quota allocation on the LMP of the day-ahead electricity market, a two-stage model within the CET mechanism is proposed in this paper. The relationship between the two stages of the model is shown in Figure 3.
In the first stage, there is an optimization for both system operation costs and carbon emission costs to determine the outputs and statuses of the generation units. The solutions obtained from the first stage can balance the conflicts between different objectives but fail to analyze the impact on the LMP. We utilize these solutions to limit the variables in the second stage to make sure that the outputs and statuses of generation units are the same as those obtained from the first stage. Then, the LMP calculated in the second stage is used to analyze the impact of the CET mechanism.

3.1. The First Stage of the Model: Mutli-Objective Optimization

In the first stage, a multi-objective optimization model is established with the objective of simultaneously minimizing (1) the operation costs of the power system and (2) the total costs of carbon emissions paid by traditional generator units. Meanwhile, this model should satisfy the following constraints: (1) the unit’s startup and shutdown cost constraints, (2) power balance constraints, (3) lower and upper output limits of conventional units, (4) ramp rate constraints, (5) minimum online and offline duration time constraints, (6) pumped storage unit’s operational constraints, (7) spinning reserve constraints, (8) consistent start–stop status of gas-steam combined cycle unit, (9) lower and upper output limits of renewable units, and so on.
In this stage, we establish day-ahead scheduling to determine start–stop statuses and outputs of conventional units, outputs of renewable units, and discharging/charging statuses and outputs of pumped storage units.

3.1.1. Objectives

The first objective function J1(⋅) is defined as the operation costs of the power system, including the costs of electricity purchased from generating units, the costs of spinning reserve, the costs of the unit’s startup and shutdown, and the unit no-load costs, which can be expressed as follows:
min J 1 = t = 1 24 { i = 1 N 1 [ k = 1 K C i , t k P i , t k + C i , t U s u i , t + C i , t D s d i , t + C i N L I i , t + C i U , t + C i D , t ] + s = 1 N 2 ( C s , t U s u s , t + C s , t D s d s , t + C s U , t + C s D , t ) } ,
where C i , t U and C i , t D , respectively, denote the upward and downward reserve prices bid by conventional unit i at time t; sui,t and sdi,t, respectively, represent the upward and downward reserve capacities offered by conventional unit i at time t; CiU,t and CiD,t, respectively, denote the startup and shutdown costs of the conventional unit i at time t; C s , t U and C s , t D , respectively, represent the upward and downward reserve prices bid by pumped storage unit s at time t; sust and sds,t, respectively, represent the upward and downward reserve capacities offer by pumped storage unit s at time t; CsU,t and CsD,t, respectively, denote the startup and shutdown costs of the pumped storage unit s at time t; K is the number of quotation segments for conventional units; and C i , t k represents the price bid in the kth step by conventional unit i at time t. In this model, we mandate that the pumped storage unit and renewable unit take a zero quotation in the day-ahead market to ensure they win the bid first.
The second objective function J2(⋅) is defined as the total cost of carbon emissions paid by traditional generator units as follows:
min J 2 = i = 1 N 1 E i , c × p r e ,
where Ei,c represents the overall emissions exceeding the quota of conventional unit i in a day, which can be expressed as follows:
E i , c = t = 1 24 E i , t η E q , i ,
where η denotes the rate of free quota; and Ei,t is the emissions of conventional unit i at time t, which can be expressed as follows:
E i , t = β c i P i , t .

3.1.2. Unit’s Startup and Shutdown Cost Constraints

{ C i U , t K i ( I i , t I i , t 1 ) , C i U , t 0 C s U , t K s ( Z s , t Z s , t 1 ) , C s U , t 0 ,
{ C i D , t J i ( I i , t 1 I i , t ) , C i D , t 0 C s D , t J s ( Z s , t 1 Z s , t ) , C s D , t 0 ,
where Ki and Ji represent the startup and shutdown cost coefficients of the conventional unit i, respectively; Ks and Js represent the startup and shutdown cost coefficients of the pumped storage unit s, respectively; and Ii,t and Zs,t denote the on/off statuses of the conventional unit i and the discharging/charging statuses of the pumped storage unit s at time t, respectively.

3.1.3. Lower and Upper Limits of Conventional Units

I i , t P i , min P i , t I i , t P i , max .

3.1.4. Ramp Rate Constraints of Conventional Units

{ P i , t P i , t 1 r u i T 60 I i , t 1 + P i , min ( I i , t I i , t 1 ) P i , t 1 P i , t r d i T 60 I i , t + P i , min ( I i , t 1 I i , t ) ,
where rui and rdi represent the upward and downward ramping rate limits of unit i, respectively.

3.1.5. Pumped Storage Unit’s Operational Constraints

{ 0 P p g , s , t P p g , s , max Z p g , s , t P p p , s , max Z p p , s , t P p p , s , t 0 ,
Z s , t = Z p g , s , t + Z p p , s , t 1 ,
t = 1 T P p g , s , t + ξ t = 1 T P p p , s , t = 0 ,
where Ppg,s,t and Ppp,s,t denote the discharging and charging power of the pumped storage unit s at time t, respectively; Ppg,s,max and Ppp,s,max represent the upper limits of discharging and charging power of unit s, respectively; Zpg,s,t and Zpp,s,t indicate the discharging/charging status of the pumped storage unit s, respectively; and ξ is the conversion efficiency of the pumped storage unit, which is set to 75%.

3.1.6. Spinning Reserve Constraints

{ S u , t = i = 1 N 1 + N 2 s u i , t P L o a d , t L u % S d , t = i = 1 N 1 + N 2 s d i , t P L o a d , t L d % ,
{ 0 s u s , t P p g , s , max Z p g , s , t P p g , s , t P p p , s , t 0 s d s , t P p g , s , t + P p p , s , t P p p , s , max Z p p , s , t ,
{ 0 s u i , t min ( P i , max I i , t P i , t , r u i T 60 I i , t ) 0 s d i , t min ( P i , t P i , min I i , t , r d i T 60 I i , t ) ,
where Su,t and Sd,t indicate the capacity of the upward and downward spinning reserve requirements at time t, respectively; and Lu% and Ld% denote the coefficients of load forecast deviation for the upward and downward reserve capacity of the system, respectively.

3.1.7. Lower and Upper Output Limits of Renewable Units

{ 0 P v , t P v , t , max 0 P w , t P w , t , max ,
where Pw,t and Pv,t represent the output of the wind farm unit w and the photovoltaic power generation unit v at time t, respectively.

3.1.8. Power Balance Constraints

i = 1 N 1 P i , t + s = 1 N 2 ( P p g , s , t + P p p , s , t ) + j = 1 N 3 P j , t + w = 1 N 4 P w , t + v = 1 N 5 P v , t = d = 1 N D P d , t + P L o s s , t ,
where Pj,t is the injected power of the tie line unit j at time t; Pd,t denotes the demand of load bus d at time t; and PLoss,t indicates the power loss of the system at time t, which is set to 0.075% of the total load demand.

3.1.9. Network Transmission Constraints

P l , t = i = 1 N 1 G l , i P i , t + s = 1 N 2 G l , s ( P p g , s , t + P p p , s , t ) + j = 1 N 3 G l , j P j , t + w = 1 N 4 G l , w P w , t + v = 1 N 5 G l , v P v , t d = 1 N D D l , d P d , t .
P l max P l , t P l max l = 1 , 2 , , N L ; t = 1 , 2 , , T .

3.1.10. Minimum Online and Offline Duration Time Constraints

{ I i , t = 1 , t [ 1 , U i ] , U i = min { T , ( T o n _ i X o n _ i , 0 ) I i , 0 } n = t t + T o n _ i 1 I i , n T o n _ i ( I i , t I i , t 1 ) , t [ U i + 1 , T T o n _ i + 1 ] n = t T [ I i , n ( I i , t I i , t 1 ) ] 0 , t [ T T o n _ i + 2 , T ]
{ I i , t = 0 , t [ 1 , D i ] , D i = min { T , ( T o f f _ i X o f f _ i , 0 ) ( 1 I i , 0 ) } n = t t + T o f f _ i 1 ( 1 I i , n ) T o f f _ i ( I i , t 1 I i , t ) , t [ D i + 1 , T T o f f _ i + 1 ] n = t T [ 1 I i , n ( I i , t 1 I i , t ) ] 0 , t [ T T o f f _ i + 2 , T ] ,
where Ui/Di represents the number of hours unit i must be initially on/off due to its minimum on/off time limit; Ton_i/Toff_i represents the minimum on/off time of unit i; and Xon_i,0/Xoff_i,0 represents the initial on/off hours of unit i.

3.1.11. Consistent Start–Stop Status Constraint of Gas-Steam Combined Cycle Units

I i , t = I j , t ,
where Ii,t and Ij,t denote the on/off status of a large-capacity unit i and a small-capacity unit j, respectively. In the proposed model, we separately dispatch two different units of a gas-steam combined cycle unit and control their start–stop states to be consistent.

3.1.12. Simplified Formation of the First Stage of the Model

For the convenience of later discussion, we rewrite the multi-objective optimization model (10)–(30) as the following compact form:
min ( J 1 , J 2 ) ,
s . t . f ( x 1 , x 2 ) = 0 ,
g ( x 1 , x 2 ) 0 ,
where x1 represents the vector of outputs of conventional units, pumped storage units, and renewable units; and x2 denotes the vector of the on/off status of the conventional units and the discharging/charging status of pumped storage units.

3.2. The Second Stage of the Model: Tracking Pareto Optimal Solution

The day-ahead market-clearing model is essentially a security-constrained unit commitment problem [17]. Therefore, the analytic expression of the marginal electricity price in the day-ahead market can be obtained using information such as the Lagrangian multipliers of the power balance constraints. In the first stage, we strike a better balance among the conflicts between different objectives but fail to obtain the marginal electricity price information directly. In the second stage, we propose a typical economic dispatch model to avoid reducing the convexity of the clearing problem for calculating the LMP [18]. We utilize the solutions obtained from the first stage of the model to fix the on/off status of the conventional units and the discharging/charging status of the pumped storage units as well as to limit the outputs of those generation units.

3.2.1. Objective

The objective of the second stage is to minimize the costs of system operation with the fixed on/off status of the conventional units and the fixed discharging/charging status of the pumped storage units, which can be written as Equation (34).
min J = t = 1 24 { i = 1 N 1 [ k = 1 K C i , t k P i , t k + C i , t U s u i , t + C i , t D s d i , t ] + s = 1 N 2 ( C s , t U s u s , t + C s , t D s d s , t ) } .

3.2.2. Constraints

The second stage of the model should satisfy not only the constraints (16)–(18) and (20)–(27) mentioned in the first stage but also the additional constraints (35)–(39) listed below, which can ensure the solution obtained from the second stage of the model is the same as that obtained from the first stage:
P i , t ε P i , t P i , t + ε ,
P w , t ε P w , t P w , t + ε ,
P v , t ε P v , t P v , t + ε ,
( P p g , s , t + P p p , s , t ) ε P p g , s , t + P p p , s , t ,
P p g , s , t + P p p , s , t ( P p g , s , t + P p p , s , t ) + ε ,
where P i , t , P w , t , P v , t , and ( P p g , s , t + P p p , s , t ) represent the output of the conventional unit i, the wind farm w, the PV power generation v, and the pumped storage unit s at time t, obtained from the first stage model, respectively; ε denotes the relaxation error and is set to 0.0000001.
In addition, the LMP for the settlement [19] is expressed as follows:
L M P k , t = λ t l = 1 N L ( τ l , t max τ l , t min ) G l , k ,
where LMPk,t is the marginal price of the node k at time t; λt represents the Lagrangian multipliers of the power balance constraints at time t; τ l , t max and τ l , t min denote the Lagrangian multipliers of the maximum forward power flow constraint and the maximum backward power flow constraint for line l at time t, respectively; and Gl,k indicates the transfer distribution factor of generator output from generator node k to line l.
According to Equation (40), if transmission is constraining, the LMPs at the two ends of the congested line are different.

4. Solution and Methodology

In this section, we discuss the optimization algorithm to solve the two-stage model. Obviously, the second stage of the model is a linear programming problem and can be easily solved using the GAMS/CPLEX solver [20], whereas the first stage of the model in (31)–(33) is a multi-objective mixed-integer linear programming problem.
We apply the NNC method [21] to solve the multi-objective mixed-integer linear programming problem, which is elaborated as follows.
Step 1: Solve the single objective model, which aims to minimize J1 and is subject to (32) and (33), to obtain the minimum value of objective J1, written as J1 min, and the maximum value of objective J2, written as J2 max.
Step 2: Solve the single objective model, which aims to minimize J2 and is subject to (32) and (33), to obtain the minimum value of objective J2, written as J2 min, and the maximum value of objective J1, written as J1 max.
Step 3: Normalize the objective (31) by J1 min, J1 max, J2 min, and J2 max obtained in Steps 1 and 2, that is,
J ¯ = [ J ¯ 1 , J ¯ 2 ] = [ J 1 J 1 min J 1 max J 1 min , J 2 J 2 min J 2 max J 2 min ] .
Step 4: Take J ¯ 1 as the abscissa, J ¯ 2 as the ordinate, and A1 (0, 1) and A2 (1, 0) as the anchor points, and depict the Pareto front, Utopia line, and normalized objective space in Figure 4.
Step 5: Divide the Utopia line into m equal segments, and generate m + 1 evenly distributed segmentation points in the Utopia line, that is,
X ¯ p j = ( 1 j m ) ( 0 , 1 ) + j m ( 1 , 0 ) = ( j m , 1 j m ) j = 0 , 1 , , m .
Step 6: Draw the normal P j of the Utopia line across X ¯ p j , and generate a corresponding Pareto point Bj by intersecting this normal with the Pareto front. In order to find the optimal point Bj in the Pareto front, minimize J ¯ 2 , which is subject to Equations (32) and (33). This optimization process should be forced to stay within the feasible region of the normalized objective space, i.e., the region labeled with the dot in Figure 4. Mathematically, the single objective optimization problem can be expressed as follows:
min J ¯ 2 ,
s . t . f ( x 1 , x 2 ) = 0 ,
g ( x 1 , x 2 ) 0 ,
N 1 ( J ¯ X ¯ p j ) T 0 ,
where J ¯ = ( J ¯ 1 , J ¯ 2 ) represents the point in normalized objective space; Equation (45) is an additional inequality to force the upper left region of P j to be feasible, as denoted with the dot in Figure 4; and the vector N 1 = 1 45 denotes the direction from A1 to A2.
In this paper, we use the GAMS/CPLEX solver to solve the single objective optimization model in Equations (42)–(45), which is a mixed-integer linear programming problem. When P j moves along the Utopia line from A1 to A2 by setting j = 1 , 2 , , m 1 , the Pareto points except for the endpoints are obtained by solving the above single objective optimization problems.

5. Case Studies and Numerical Results

5.1. Overview of Test System

This study proceeds with calculations using the operation data of a real 1560-bus system in China to verify the performance of the proposed methodology. This system consists of 729 branches, 1028 transformer branches, and 52 generators, including four nuclear power plants with total capacities of 4152 MW, whose daily outputs are constant; two pumped storage units with a capacity of 350 MW; two wind farms with capacities of 360 MW and 280 MW, respectively; two PV units with capacities of 250 MW and 220 MW, respectively; four refuse incineration power plants; six coal-fired units; twenty-eight gas-fired units; and four combined heat and power units (all gas-fired). The total installed capacity of the system is 11,543 MW. This system has peak and valley loads of 17,434 MW and 9370 MW, respectively. Figure 5 presents the forecast total daily load, the outputs of two wind farms, and the outputs of two PV units. For this system, Lu% and Ld% are set to 3% and 1%, respectively; η is set to 95%; and the emission coefficients of coal-fired and gas-fired units are set to 0.85 t/MWh and 0.35 t/MWh [22], respectively.

5.2. Normalized Normal Constraint (NNC) Method for Solving the First Stage of the Model

First, we use the historical carbon emission-based method to allocate carbon emission quotas for the conventional units. According to the NNC method mentioned above, we can obtain several optimal solutions from the first stage of the model with minimum system operation costs and minimum carbon emission costs. In the simulation, we set m as 10, which indicates that the Utopia line is divided into 10 segments. Then, we can use the GAMS/CPLEX solver to obtain the Pareto front. Table 1 and Figure 6 illustrate the 11 optimal solutions of the Pareto front, representing the costs of system operation and the costs of carbon emission.
According to Table 1, the value of j varies from 0 to 10, indicating that the weights of the objectives J ¯ 1 and J ¯ 2 change accordingly. After obtaining the optimal solution set, we need to select the compromise solution to fix the on/off statuses and outputs of generating units. Various techniques to select the compromise solution were proposed in [23,24]. In this paper, we expect to find a compromised optimal solution that minimizes the sacrifice to J1 while greatly reducing the value of J2. Therefore, we can find that the most suitable solution is at j = 4, according to the comparison between the nodes of the Pareto front and their values corresponding to objectives J ¯ 1 and J ¯ 2 in Table 1 and Figure 4. With this solution, the values of J ¯ 1 and J ¯ 2 are 0.01083 and 0.21083, respectively, indicating that the system operation costs are very close to the minimum, and the total carbon emission costs paid by conventional units are small. Hence, we choose the point ( J ¯ 1 , J ¯ 2 ) = ( 0.01083 , 0.21083 ) as the compromise solution.
The traditional clearing model for calculating the LMP only minimizes the system operation costs, which will greatly reduce the control of carbon emissions by the CET mechanism. To explain the advantage of using the proposed two-stage model, we first conduct a simulation in the real 1560-bus and 52-unit system with two scenarios, including clearing without the CET mechanism and clearing in the CET mechanism with the traditional model. As shown in Table 2, we can see that the total carbon emissions remain high in the CET mechanism with the traditional clearing model compared with the scenario without the CET mechanism, meaning it does little to reduce emissions.
In order to illustrate the performance of the proposed methodology, we compare the system operation costs and carbon emission costs obtained from the two-stage model with those of the traditional clearing model that does not consider carbon emission costs in the CET mechanism. According to Table 3, the system operation costs of the two-stage model are only 0.68% higher than those of the traditional clearing model, but the total carbon emission costs paid by conventional units are reduced by 67.25%, and the carbon emissions are reduced by 20.82%. Therefore, although the proposed method does not reduce the system operation costs, it greatly reduces the carbon emission costs paid by conventional units, thus alleviating conflicts between different objectives and effectively reducing system carbon emissions.

5.3. Calculation of Locational Marginal Price (LMP) Considering the Impact of Carbon Emission Quotas

In the simulation, we use the two-stage model to solve the LMP of the day-ahead electricity market considering the impact of carbon emission quotas. In order to ensure the accuracy of the LMP calculated by the two-stage model, we compare the solutions obtained from the first stage of the model and the second stage of the model. Figure 7 shows the comparison of outputs of renewable energy units and pumped storage units. Figure 8 illustrates the comparison of outputs of the coal-fired unit with a capacity of 350 MW.
Figure 7 and Figure 8 indicate that the output of each unit in the second stage is consistent with that in the first stage. Simultaneously, the system operation costs and carbon emission costs obtained from the two stages in Table 4 are considerably close to each other. In addition, the relative errors do not exceed 10−7. Therefore, we can ensure the accuracy of the LMP obtained from the two-stage model.
We compare the LMP calculated by the two-stage model in the CET mechanism with that obtained from the traditional clearing model without CET in the past and evaluate the impact of carbon emission quotas on the clearing of the day-ahead electricity market. In this section, we take the calculation results without CET as the benchmark. Figure 9 and Figure 10 illustrate the comparison of the LMP.
According to Figure 9 and Figure 10, the LMP has changed after considering the impacts of the CET mechanism. The LMP in the CET mechanism during the peak load period is higher than that without the impact of carbon emission quota allocation, with a maximum increase of 27.65%, while it drops during the valley load period with a maximum decrease of 25.54%. Both daily LMP variation trends are the same as the load curve.

5.4. Clearing Results of Different Carbon Emission Quota Allocation Methods

Figure 11 illustrates the comparison of carbon emissions of conventional units between two different carbon emission quota allocation methods. The on/off statuses of conventional units are the same under two different allocation methods. As shown in Figure 9, units 1–6 are coal-fired units, and units 7–26 are gas-fired units. When the quota is allocated based on the historical carbon emissions, all of the carbon emissions of the coal-fired units do not exceed their free quotas, and only three of the twenty gas-fired units are below the allocated quotas. In the case of using the power generation performance-based method to allocate quotas, the carbon emissions of some coal-fired units do exceed their reduced quotas due to their lower power generation performance. On the contrary, only six of the twenty gas-fired units have exceeded their quotas because of their increased quotas and their low carbon emission coefficients. The possibility that the carbon emissions of gas-fired units will exceed their quotas has been greatly reduced.
Figure 12 displays the comparison of outputs of coal-fired units between two different carbon emission quota allocation methods. Due to the high carbon emission coefficient, the output of coal-fired units will be cut down to reduce their total carbon emission costs in the CET mechanism. The reduction of output is more obvious during the valley load periods. The output of coal-fired units will be further reduced (up to 59.62%) if their quotas are allocated based on their power generation performance.
Figure 13 shows the comparison of outputs of gas-fired units between two different carbon emission quota allocation methods. Although the marginal costs of gas-fired units are higher than those of coal-fired units, the carbon emission costs are lower due to their lower carbon emission coefficients. In order to meet the power balance constraints, the outputs of the gas-fired units will increase in the CET mechanism. The increase of gas-fired units’ outputs can reach a maximum of 58.4% when the emission quotas are allocated based on power generation performance.
Figure 14 demonstrates the comparison of outputs of renewable energy units between two different carbon emission quota allocation methods. Renewable units do not emit carbon during power generation. Due to the sufficient peak shaving capacity of the proposed system, the outputs of renewable units are not affected by the CET.
Figure 15 presents the comparison of outputs of pumped storage power stations. The pumped storage units are all in charging mode during periods of valley load, and their outputs are exactly the same as the outputs without CET. When the emission quotas are allocated based on historical carbon emissions, the changes of outputs of conventional units are small during the peak load period, so the outputs of pumped storage units do not change significantly. However, the outputs of conventional units do change greatly due to the large difference in power generation performance with the power generation performance-based method. Therefore, the behavior of pumped storage units’ output will also change significantly during the peak load period.
Figure 16 presents the comparison of the LMP with different carbon emission quota allocation methods. When the emission quotas for conventional units are allocated by the power generation performance-based method, coal-fired units with lower generating costs obtain lower carbon quotas due to their lower power generation performance; hence, their quotation prices will significantly increase in the day-ahead market. The quotation prices of all conventional units in the market are at a relatively high level. Compared with the historical carbon emission-based method, the LMP of the electricity market has been increased, most significantly (>62%) during peak load periods.

6. Conclusions

In this study, a two-stage model is formulated for the LMP calculation with CET. For the first stage of the model, a multi-objective optimization model to minimize both system operation costs and carbon emission costs is proposed and solved by the NNC method. For the second stage of the model, a tracing model with the goal of minimizing the total costs of the power system and satisfying the constraints generated by the Pareto optimal solution is built to calculate the LMP. Based on the data of a real-world 52-unit and 1560-bus system, the performance of the proposed methodology is studied in detail to analyze the impact of the CET.
According to the simulation, the outputs of coal-fired units are reduced, and the outputs of gas-fired units are increased. These results are even more pronounced in the allocation scheme based on the power generation performance-based method. The total operation costs of conventional units increase in the CET mechanism, and the quotation prices of conventional units increase in the day-ahead market, thereby raising the LMP sharply, especially during peak load periods. Furthermore, the LMP may rise up to 62% in the case of quotas allocated by the power generation performance-based method, which will significantly increase electricity costs and reduce social welfare. Therefore, this paper can provide a reference for the study of market clearing in consideration of the CET mechanism.

Author Contributions

Methodology, M.W., Z.L. and M.L.; Supervision, Q.C., T.Z., E.L. and W.L.; Writing—original draft, M.W., Z.L. and M.L.; Writing—review & editing, M.W., Z.L. and M.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research is funded by the Science and Technology Project of Guangdong Power Exchange Center Co. under Grant GDKJXM20172986.

Acknowledgments

The authors would like to thank the editor and reviewers for their sincere suggestions on improving the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, J.P.; Zhang, X.B.; Song, X.H. Regional carbon emission evolution mechanism and its prediction approach driven by carbon trading—A case study of Beijing. J. Clean. Prod. 2018, 172, 2793–2810. [Google Scholar]
  2. Chen, X.; Lv, J.; McElroy, M.B.; Han, X.; Nielsen, C.P.; Wen, J. Power system capacity expansion under higher penetration of renewables considering flexibility constraints and low carbon policies. IEEE Trans. Power Syst. 2018, 33, 6240–6253. [Google Scholar]
  3. Lu, S.; Wu, Y.; Lou, S.; Yin, X. A model for optimizing spinning reserve requirement of power system under low-carbon economy. IEEE Trans. Sustain. Energy 2014, 5, 1048–1055. [Google Scholar]
  4. Yuchi, Q.; Wang, N.; Li, S.; Yang, Z.; Jiang, B. A bi-objective reverse logistics network design under the emission trading scheme. IEEE Access 2019, 7, 105072–105085. [Google Scholar]
  5. Huang, J.; Xue, Y.; Jiang, C.; Wen, F.; Xue, F.; Meng, K.; Dong, Z.Y. An experimental study on emission trading behaviors of generation companies. IEEE Trans. Power Syst. 2015, 30, 1076–1083. [Google Scholar]
  6. Zhang, Z.; Zhu, K. Border carbon adjustments for exports of the United States and the European Union: Taking border-crossing frequency into account. Appl. Energy 2017, 201, 188–199. [Google Scholar]
  7. Ding, G.; Deng, Y.; Lin, S. A study on the classification of China’s provincial carbon emissions trading policy instruments: Taking Fujian province as an example. Energy Rep. 2019, 5, 1543–1550. [Google Scholar]
  8. Lin, B.; Jia, Z. Is emission trading scheme an opportunity for renewable energy in China? A perspective of ETS revenue redistributions. Appl. Energy 2020, 263, 114605. [Google Scholar]
  9. Wang, Y.; Qiu, J.; Tao, Y.; Zhao, J. Carbon-oriented operational planning in coupled electricity and emission trading markets. IEEE Trans. Power Syst. 2020. [Google Scholar] [CrossRef]
  10. Zhou, X.; James, G.; Liebman, A.; Dong, Z.Y.; Ziser, C. Partial carbon permits allocation of potential emission trading scheme in Australian electricity market. IEEE Trans. Power Syst. 2010, 25, 543–553. [Google Scholar]
  11. Yang, H.; Chung, C.Y.; Wong, K.P. Optimal fuel, power and load-based emissions trades for electric power supply chain equilibrium. IEEE Trans. Power Syst. 2012, 27, 1147–1157. [Google Scholar]
  12. Li, X.R.; Yu, C.W.; Xu, Z.; Luo, F.J.; Dong, Z.Y.; Wong, K.P. A multimarket decision-making framework for genco considering emission trading scheme. IEEE Trans. Power Syst. 2013, 28, 4099–4108. [Google Scholar]
  13. Kuo, M.T.; Lu, S.D.; Tsou, M.C. Considering carbon emissions in economic dispatch planning for isolated power systems: A case study of the taiwan power system. IEEE Trans. Ind. Appl. 2018, 54, 987–997. [Google Scholar]
  14. D Olsen, D.J.; Dvorkin, Y.; Fernandez-Blanco, R.; Ortega-Vazquez, M.A. Optimal carbon taxes for emissions targets in the electricity sector. IEEE Trans. Power Syst. 2018, 33, 5892–5901. [Google Scholar]
  15. Li, J.L.; Chen, Y.G.; Liu, S.j.; Wang, N.; Zou, P.; Chen, Q.X. Electricity market equilibrium analysis considering carbon emission cost. Power Syst. Technol. 2016, 40, 1558–1563. [Google Scholar]
  16. Wang, T.; Jiang, C.W.; Li, L.; Zhao, Y. Security constrained unit commitment problem considering carbon emission allocation and demand side resources. Power Syst. Technol. 2016, 40, 3355–3361. [Google Scholar]
  17. Fernández-Blanco, R.; Arroyo, J.M.; Alguacil, N. On the solution of revenue- and network-constrained day-ahead market clearing under marginal pricing—Part I: An exact bilevel programming approach. IEEE Trans. Power Syst. 2017, 32, 208–219. [Google Scholar]
  18. Knueven, B.; Ostrowski, J.; Castillo, A.; Watson, J.P. A Computationally Efficient Algorithm for Computing Convex Hull Prices. Available online: http://www.optimization-online.org/DB_FILE/2019/09/7370.pdf (accessed on 9 January 2020).
  19. Zhao, Z.; Wu, L. Impacts of high penetration wind generation and demand response on LMPs in day-ahead market. IEEE Trans. Smart Grid 2014, 5, 220–229. [Google Scholar]
  20. Rosenthal, R.E. GAMS, A User’s Guide; GAMS Development Corp.: Washington, DC, USA, 2014. [Google Scholar]
  21. Li, Q.; Liu, M.; Liu, H. Piecewise normalized normal constraint method applied to minimization of voltage deviation and active power loss in an AC–DC hybrid power system. IEEE Trans. Power Syst. 2015, 30, 1243–1251. [Google Scholar]
  22. Yu, S.; Wei, Y.M.; Guo, H.; Ding, L. Carbon emission coefficient measurement of the coal-to-power energy chain in China. Appl. Energy 2014, 114, 290–300. [Google Scholar]
  23. Marler, R.T.; Arora, J.S. The weighted sum method for multi-objective optimization: New insights. Struct. Multidiscip. Optim. 2010, 41, 853–862. [Google Scholar]
  24. Messac, A.; Ismail-Yahaya, A.; Mattson, C.A. The normalized normal constraint method for generating the Pareto frontier. Struct. Multidiscip. Optim. 2003, 25, 86–98. [Google Scholar]
Figure 1. Market framework considering the interaction between carbon emission trading (CET) market and day-ahead electricity market.
Figure 1. Market framework considering the interaction between carbon emission trading (CET) market and day-ahead electricity market.
Energies 13 02510 g001
Figure 2. Multi-segment bid curve for conventional units (K = 5).
Figure 2. Multi-segment bid curve for conventional units (K = 5).
Energies 13 02510 g002
Figure 3. Relationship between the two stages of the model.
Figure 3. Relationship between the two stages of the model.
Energies 13 02510 g003
Figure 4. Pareto front, Utopia line, and normalized objective space.
Figure 4. Pareto front, Utopia line, and normalized objective space.
Energies 13 02510 g004
Figure 5. Forecast total daily load, the outputs of two wind farms, and the outputs of two PV units.
Figure 5. Forecast total daily load, the outputs of two wind farms, and the outputs of two PV units.
Energies 13 02510 g005
Figure 6. Pareto front and optimal compromise solution.
Figure 6. Pareto front and optimal compromise solution.
Energies 13 02510 g006
Figure 7. Comparison of outputs of renewable energy units and pumped storage units. (a) Comparison of outputs of PV units; (b) Comparison of outputs of wind farms; (c) Comparison of outputs of pumped storage units.
Figure 7. Comparison of outputs of renewable energy units and pumped storage units. (a) Comparison of outputs of PV units; (b) Comparison of outputs of wind farms; (c) Comparison of outputs of pumped storage units.
Energies 13 02510 g007
Figure 8. Comparison of outputs of the coal-fired unit with capacity of 350 MW.
Figure 8. Comparison of outputs of the coal-fired unit with capacity of 350 MW.
Energies 13 02510 g008
Figure 9. Comparison of the marginal price of the day-ahead market.
Figure 9. Comparison of the marginal price of the day-ahead market.
Energies 13 02510 g009
Figure 10. Variation magnitude of marginal price of the day-ahead electricity market.
Figure 10. Variation magnitude of marginal price of the day-ahead electricity market.
Energies 13 02510 g010
Figure 11. Comparison of carbon emissions of conventional units.
Figure 11. Comparison of carbon emissions of conventional units.
Energies 13 02510 g011
Figure 12. Comparison of outputs of coal-fired units.
Figure 12. Comparison of outputs of coal-fired units.
Energies 13 02510 g012
Figure 13. Comparison of outputs of gas-fired units.
Figure 13. Comparison of outputs of gas-fired units.
Energies 13 02510 g013
Figure 14. Comparison of outputs of renewable units.
Figure 14. Comparison of outputs of renewable units.
Energies 13 02510 g014
Figure 15. Comparison of outputs of pumped storage units.
Figure 15. Comparison of outputs of pumped storage units.
Energies 13 02510 g015
Figure 16. Comparison of locational marginal price (LMP) with different carbon emission quota allocation methods.
Figure 16. Comparison of locational marginal price (LMP) with different carbon emission quota allocation methods.
Energies 13 02510 g016
Table 1. Total system operation costs and total carbon emission costs represented by 11 Pareto frontier points.
Table 1. Total system operation costs and total carbon emission costs represented by 11 Pareto frontier points.
j J ¯ 1 J ¯ 2
0−2.38 × 10−151
11.37 × 10−40.80016
24.09 × 10−40.60041
30.001160.40216
40.010830.21083
50.052530.05253
60.25.29 × 10−15
70.43.66 × 10−15
80.64.88 × 10−15
90.83.66 × 10−15
1011.22 × 10−15
Table 2. Comparison of total system carbon emissions under different clearing methods.
Table 2. Comparison of total system carbon emissions under different clearing methods.
ScenarioTotal Carbon Emissions (tons)
Clearing without CET mechanism198,259.35
Clearing in CET mechanism with traditional model191,733.25
Table 3. Comparison of results between two optimization models.
Table 3. Comparison of results between two optimization models.
Traditional Clearing Model Two-Stage ModelRangeability
System operation costs (¥)28,650,271.9228,845,093.770.68%
Carbon emission costs (¥)3,045,165.34997,291.65−67.25%
Carbon emissions (tons)191,733.25151,814.39−20.82%
Table 4. Comparison of the system operation costs and carbon emission costs obtained from the two stages of the proposed model.
Table 4. Comparison of the system operation costs and carbon emission costs obtained from the two stages of the proposed model.
System Operation Costs (¥)Carbon Emission Costs (¥)
First stage of model28,845,093.77997,291.65
Second stage of model28,845,093.44997,291.78
Error−0.000001141%0.000012677%

Share and Cite

MDPI and ACS Style

Wu, M.; Lu, Z.; Chen, Q.; Zhu, T.; Lu, E.; Lu, W.; Liu, M. A Two-Stage Algorithm of Locational Marginal Price Calculation Subject to Carbon Emission Allowance. Energies 2020, 13, 2510. https://doi.org/10.3390/en13102510

AMA Style

Wu M, Lu Z, Chen Q, Zhu T, Lu E, Lu W, Liu M. A Two-Stage Algorithm of Locational Marginal Price Calculation Subject to Carbon Emission Allowance. Energies. 2020; 13(10):2510. https://doi.org/10.3390/en13102510

Chicago/Turabian Style

Wu, Mingxing, Zhilin Lu, Qing Chen, Tao Zhu, En Lu, Wentian Lu, and Mingbo Liu. 2020. "A Two-Stage Algorithm of Locational Marginal Price Calculation Subject to Carbon Emission Allowance" Energies 13, no. 10: 2510. https://doi.org/10.3390/en13102510

APA Style

Wu, M., Lu, Z., Chen, Q., Zhu, T., Lu, E., Lu, W., & Liu, M. (2020). A Two-Stage Algorithm of Locational Marginal Price Calculation Subject to Carbon Emission Allowance. Energies, 13(10), 2510. https://doi.org/10.3390/en13102510

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop