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

: 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.


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] and [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. 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.

Market Framework
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.  Figure 1. Market framework considering the interaction between carbon emission trading (CET) market and day-ahead electricity market.

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: where Q is the total daily quota available to allocate; E b 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: where E b i 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: 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:

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: where ai, bi, and ci denote the production cost coefficients, and P k' i,t represents the bidding quantities corresponding to the quotation of unit i in the kth step at time t, which is written as: 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, where P k i,t 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, 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, 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. Figure 2. Multi-segment bid curve for conventional units (K = 5).

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.  Figure 3. Relationship between the two stages of the model.

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.

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: where C U i,t and C D i,t , 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 U s,t and C D s,t , 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 k i,t 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: where Ei,c represents the overall emissions exceeding the quota of conventional unit i in a day, which can be expressed as follows: 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:

Unit's Startup and Shutdown Cost Constraints
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.

Ramp Rate Constraints of Conventional Units
where rui and rdi represent the upward and downward ramping rate limits of unit i, respectively.

Pumped Storage Unit's Operational Constraints
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%.

Spinning Reserve Constraints
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.

Lower and Upper Output Limits of Renewable Units
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.

Power Balance Constraints
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.    I  t  D D  T T  X  I   I  T  I  I   t D  T T   I  I  I  t T T  T

Network Transmission Constraints
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.

Consistent Start-Stop Status Constraint of Gas-Steam Combined Cycle Units
, , 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.

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: 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.

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.

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).

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: 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; τ max l,t and τ min l,t 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.

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 mixedinteger 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 J1min, and the maximum value of objective J2, written as J2max.
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 J2min, and the maximum value of objective J1, written as J1max.
Step 3: Normalize the objective (31) by J1min, J1max, J2min, and J2max obtained in Steps 1 and 2, that is, Step 4: Take as the abscissa, 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, Step 6: Draw the normal where  

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, 1,028 transformer branches, and 52 generators, including four nuclear power plants with total capacities of 4,152 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 9,370 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 tons/MWh and 0.35 tons/MWh [22], respectively.  Figure 6. Pareto front and optimal compromise solution.

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 1 J and 2 J 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] and [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 1 J and 2 J in Table 1 and 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 1560bus 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.

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 coalfired unit with a capacity of 350 MW.   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. Figures 9 and 10 illustrate the comparison of the LMP. According to Figures 9 and 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.  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 gasfired units will exceed their quotas has been greatly reduced.  Figure 11. Comparison of carbon emissions of conventional units. 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 coalfired 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.     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.

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 gasfired 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 performancebased 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.