An Annual Electric Energy Trade Scheduling Model under the Dual Track Mode

: The annual electricity trade scheduling is the basis of long-term power generation scheduling. In recent decades, the ratio of new energy generation in China has increased annually, and the electricity market has operated under the “market electricity” and “planned electricity” double track mode in recent years. However, the existing annual electricity trade scheduling methods are extensive and cannot adapt to the new situation of “market electricity” and large-scale new energy generation. The annual scheduled energy of the power units is set as a decision variable, and a novel annual energy scheduling optimization model based on Gini coefﬁcient of fairness is presented in this paper. In this model, “market electricity” capacity is conversed monthly, considering peaking reserve demand and monthly characteristics of new energy generation. The fairness constraint set based on Gini coefﬁcient is introduced into the optimization model to solve various fairness problems. Simulation results show that the introduction of the Gini coefﬁcient and the optimization model considering the monthly conversion of marketing electricity capacity can obtain more accurate and reasonable electricity distribution results, and the peaking demand can be considered more fairly and effectively. The proposed method provides a feasible solution to the annual electric energy scheduling for dual-track operation country such as China.


Introduction
The annual electricity trade scheduling is usually considered as the basement of the formulation and implementation of the electricity trade scheduling in China. It is the most important processes of the long-term power generation scheduling. The annual electricity generation demand should be reasonably allocated to each unit according to certain principles. The average allocation method according to the generating hours of each unit is the most primitive and traditional method for annual electricity allocation in China, which ensures the "fairness" in traditional sense and is the method of electricity energy scheduling continuously adopted in China in the past years.
At the beginning of this century, "energy-saving generation scheduling" was put forward in China with increasing attention to the energy environment. Under the background of energy-saving generation dispatching, provinces generally adopt the differential electric quantity model. All fossil energy generating units were classified or graded according to the unit type, energy consumption grade, capacity difference, environmental protection level, location and other factors. Then, according to the total electric quantity and the total number of "gears" or "stages", the difference of power generation utilization hours or basic generation utilization hours between each gear is determined [1][2][3][4]. The differences of unit power generation utilization hours further expanded between large capacity high efficiency income and consumption of equity, allocation fairness of resources, and any other things of the distribution of equal status. It is the internationally recognized measure of things distribution fairness in recent years. It has been widely applied in economics, humanities, environmental science, power system, and other fields [13][14][15][16][17][18].
If the annual electricity energy scheduling was modeled by introducing the Gini coefficient to measure the fairness of the energy scheduled for each unit, and the annual scheduled electricity energy for each unit is directly set as the decision variables, the balance between energy saving and emission reduction can be realized.
Considering the current "market electricity" and "planned electricity" double-track operation mode and the new situation of new energy generation ratio increasing year by year, a novel annual electricity energy scheduling optimization model was presented in this paper based on the Gini coefficient fairness index. In the model, the problem of traditional annual electricity planning method being too extensive and lack of theoretical basis for fairness was solved. Through introducing the monthly electricity characteristics considering load balance constraints, based on the user load reserve demand monthly capacity conversion "market electricity" and introducing the Gini coefficient of fairness constraint set, the balance among the fairness, the energy saving and emission reduction could be achieved within acceptable range of the generation producers under the doubletrack operation mode. The proposed method provides a feasible solution to the annual electric energy scheduling for dual-track operation country government, such as China.

The Concept of Gini Coefficient
The Gini coefficient was proposed on the basis of a Lorentz curve. It is used to reflect the social statistics of income level allocation fairness. Gini coefficient is an economic statistic index reflecting fairness on the whole. The Lorentz curve is a curve that reflects the fairness of a given income allocation which is shown in Figure 1. distribution of equal status. It is the internationally recognized measure of th bution fairness in recent years. It has been widely applied in economics, hum vironmental science, power system, and other fields [13][14][15][16][17][18].
If the annual electricity energy scheduling was modeled by introducing efficient to measure the fairness of the energy scheduled for each unit, and scheduled electricity energy for each unit is directly set as the decision variab ance between energy saving and emission reduction can be realized.
Considering the current "market electricity" and "planned electricity" d operation mode and the new situation of new energy generation ratio increas year, a novel annual electricity energy scheduling optimization model was p this paper based on the Gini coefficient fairness index. In the model, the probl tional annual electricity planning method being too extensive and lack of theo for fairness was solved. Through introducing the monthly electricity characte sidering load balance constraints, based on the user load reserve demand mon ity conversion "market electricity" and introducing the Gini coefficient of fa straint set, the balance among the fairness, the energy saving and emission could be achieved within acceptable range of the generation producers under track operation mode. The proposed method provides a feasible solution to electric energy scheduling for dual-track operation country government, such

The Concept of Gini Coefficient
The Gini coefficient was proposed on the basis of a Lorentz curve. It is us the social statistics of income level allocation fairness. Gini coefficient is an ec tistic index reflecting fairness on the whole. The Lorentz curve is a curve that fairness of a given income allocation which is shown in Figure 1. Based on the Lorentz curve, the Gini coefficient is: The Gini coefficient is equal to twice the area value of A S . For convenien lation, it is usually converted to the area value of B S . According to Formula (1 Based on the Lorentz curve, the Gini coefficient is: The Gini coefficient is equal to twice the area value of S A . For convenience of calculation, it is usually converted to the area value of S B . According to Formula (1), when the area of S A is zero, the allocation is absolutely fair. The larger the partial area is, the more uneven the distribution is.

Gini Coefficient Calculation Method and the Meaning of Value Range
The calculation methods of the Gini coefficient mainly include direct calculation method, regression curve method, radix bisection method and group decomposition method [19]. According to regulations of relevant United Nations organizations, the meanings of different value ranges of the Gini coefficient are shown in Table 1. If the value of the Gini coefficient is less than 0.4, the fairness is basically achieved. If the Gini coefficient value is less than 0.3, the fairness is better.

Applicability of the Gini Coefficient as Electric Energy Fairness Index
As mentioned above, the Gini coefficient is an internationally recognized indicator of economic income equity. It can be used to measure the fairness of income, consumption and resource allocation and the equal allocation of any other things [20]. Compared with the traditional standard difference index, the Gini coefficient can not only better measure the fairness of allocation of things and the system stability under different equal allocation conditions, but also the United Nations has given its index value to measure fairness. When the Gini coefficient value is between 0.3 and 0.4, the income gap is relatively fair and reasonable, which can maintain social stability [21]. In recent years, the Gini coefficient has been widely used in many fields such as economics, environmental science, plant ecology and power system. Under the current electricity market operation mode in China, the economic income of generator sets mainly comes from electricity income. Therefore, the fairness of electricity distribution actually represents the fairness of its economic income. As the Gini coefficient is a widely accepted indicator to measure the equity of economic income, it is suitable and easy to be accepted by all power producers as an indicator to measure the equity of unit electric energy distribution.

Modeling Ideas
As mentioned above, the double-track operation mode of "planned electricity" and "market electricity" will exist for a long time in my countries such as in China. The so-called "planned electricity" refers to the electricity generated, supplied and used by the relevant state organs in accordance with the corresponding policies and regulations. "Market electricity" refers to electricity generation and consumption whose price is determined through market competition. "Planned electricity" and "market electricity" dual track operation is the outcome of the current round of power system reform in China. It means that in all power generation and consumption scheduling, parts of them are determined by the government and fixed by the government, while the other scheduling is generated by market competition and settled at market prices. "Planned electricity" and "market electricity" dual track operation is one of the characteristics of electricity market operation in China, which does not exist in other countries.
Under the dual-track operation mode, generation units that participate in market transactions are called "market units". Market units usually sign a percentage of "market power" contracts. Under current market operation mode in China, all types of units have a share of "planned power". Therefore, the electric quantity of market unit includes "market electricity" and "planned electricity". Units that do not participate in market transactions are called "non-market units", and the generation capacity of such units only includes "planned electricity". Energy composition for different type of generation units is shown in Figure 2. a share of "planned power". Therefore, the electric quantity of market ket electricity" and "planned electricity". Units that do not participate tions are called "non-market units", and the generation capacity of cludes "planned electricity". Energy composition for different type of shown in Figure 2. According to the "Opinions on the Implementation of the Plan f and Development of Electricity", in principle, when arranging the plan power generation capacity should be deducted according to the direc tion. Direct transaction of large users has become an important factor the formulation of annual electricity energy scheduling. On the one ha of the corresponding capacity directly traded by large users involves istics of users and regions, and further involves the fairness of auxilia peak and frequency modulation provided by each generator set, and of the electricity energy schedule. On the other hand, in recent years, proportion of new energy generation integrated into the grid, such a photovoltaic, system regulation faces greater pressure [22][23][24]. Under strong seasonality of new energy, the traditional electric energy sched not carefully consider the specific electric quantity implementation in ing in the increase of electric energy scheduling implementation pres year by year. Therefore, the factors such as "market electricity" and n tion should be considered in the annual electricity energy scheduling m power generation energy scheduling more refined and fair. Among th also include the fairness of stable load rate and the fairness of encoura tricity", that is, fairness should be considered more refined. Therefore, ity energy scheduling optimization model based on the Gini coefficien "dual track" operation mode of planned electricity and market electri this paper. By introducing monthly electricity balance constraint, conve tricity" capacity by month and introducing the Gini coefficient fairness ious problems existing in current annual electricity planning can be be When the scale of direct transaction of large users increases, it is the volume conversion of direct transaction. Since there are great diffe acteristics, system load characteristics, and output characteristics of a ergy generation in different months, the capacity of generating units pa According to the "Opinions on the Implementation of the Plan for Orderly Release and Development of Electricity", in principle, when arranging the planned electricity, the power generation capacity should be deducted according to the direct transaction situation. Direct transaction of large users has become an important factor to be considered in the formulation of annual electricity energy scheduling. On the one hand, the conversion of the corresponding capacity directly traded by large users involves the load characteristics of users and regions, and further involves the fairness of auxiliary services, such as peak and frequency modulation provided by each generator set, and even the feasibility of the electricity energy schedule. On the other hand, in recent years, with the increasing proportion of new energy generation integrated into the grid, such as wind power and photovoltaic, system regulation faces greater pressure [22][23][24]. Under the background of strong seasonality of new energy, the traditional electric energy scheduling method does not carefully consider the specific electric quantity implementation in each month, resulting in the increase of electric energy scheduling implementation pressure and deviation year by year. Therefore, the factors such as "market electricity" and new energy generation should be considered in the annual electricity energy scheduling model, so as to make power generation energy scheduling more refined and fair. Among them, fairness should also include the fairness of stable load rate and the fairness of encouraging "market electricity", that is, fairness should be considered more refined. Therefore, an annual electricity energy scheduling optimization model based on the Gini coefficient fairness under the "dual track" operation mode of planned electricity and market electricity is presented in this paper. By introducing monthly electricity balance constraint, converting "market electricity" capacity by month and introducing the Gini coefficient fairness constraint set, various problems existing in current annual electricity planning can be better solved.
When the scale of direct transaction of large users increases, it is directly affected by the volume conversion of direct transaction. Since there are great differences in load characteristics, system load characteristics, and output characteristics of all kinds of new energy generation in different months, the capacity of generating units participating in large user direct transaction should be converted at least in the unit of month. However, if the capacity is converted on a weekly or daily basis, on the one hand, the calculation scale of 365 days in a year is too large; on the other hand, the consideration is too detailed, which is not good for the generation manufacturers to understand. Moreover, accurate information in daily basis is not available at the beginning of the year, so such fine granularity is practically meaningless. Therefore, in this paper, the participating trading capacity of each unit participating in direct trading is converted on a monthly basis, and the electric quantity balance model is established on a monthly basis. For the electricity capacity conversion method, the fairness of the auxiliary services such as peaking undertaken by each unit is taken into account [25,26]. From the point of view of direct transaction, the peak load of electricity users needs the corresponding installed capacity of power generation enterprises to match. Users with high peak load occupy more installed capacity resources, so the calculation method of power users' direct transaction capacity considering load characteristics is adopted, that is, the power generation enterprises should undertake the task of peak load regulation with direct transaction load.
In view of the fairness of scheduled electricity decomposition, the principle in this paper is to encourage power generation enterprises to participate in direct market transaction as much as possible, so the fairness constraint method of a double Gini coefficient is adopted for the generation utilization hours of planned electricity after deducting market capacity and the generation utilization hours of planned electricity without deducting market capacity. The Gini coefficient fairness constraint of planned electricity generation utilization hours after deducting market capacity can be relaxed by setting a larger Gini coefficient threshold.

The Objective Function
The mathematical model in this paper aims to minimize the sum of energy consumption cost and pollutant emission cost of all units throughout the year: where F g is the sum fuel consumption value of all thermal power units including the plants using gas, fuel oil and biomass and also other hydrocarbons, which is converted into standard coal.
where N k is the number of unit types, N T k is the number of unit type k, f k ( ) is the fuel consumption function of unit type k, and Q i is the annual scheduled energy of unit i.
For a market-oriented generator unit, the annual total electric energy is the sum of the planned electricity energy finally distributed and the medium and long term electricity energy already signed: where Q Planed,i is the planned electricity energy distributed of unit, Q Market,i is the trading energy, and N Market is the number of marketing unit. For non-market generators, the annual total electric energy is the planned electric energy finally allocated: where Q m i is the scheduled energy for unit i in month m.
where S g is the total sulfur dioxide emission, λ is the Coefficient of conversion of standard coal to raw coal according to average calorific value, η is average sulfur content of raw coal, and D i is the desulfurization rate of unit i.

The Constraint Conditions
In order to consider the generation characteristics of new energy and the capacity conversion of "market electricity" in a more refined way, the annual overall constraint and the monthly specific constraint are considered.

Annual Overall Constraints
Fairness is a very important index of annual electricity allocation under the principle of impartial and open dispatching. Compared with shorter time scale generation scheduling and economic dispatching, the annual electricity energy allocation fairness should be more strict. It should not only consider the difference of electricity energy schedule of all units within a certain fair range, but also consider the difference of electricity energy schedule of the same type of units and the same capacity units within the required fair range. The Gini coefficient is introduced in this paper as the power of all the units and various classification units set plan fairness index, the main constraint corresponding electricity power utilization hours of fairness, namely market schedule of the unit electric part of power utilization hours and non-market unit's overall power generation using the hours of fairness, fairness constraint set power plan is established.
There are four methods to calculate the Gini coefficient value: direct method, regression curve method, radix bisection method and group decomposition method. In consideration of accuracy and operability, the direct calculation method without error factors is applied in this paper. According to the direct calculation method, the Gini coefficient is expressed as follows: where u is the mean income, N is the sample size, and X i − X j is the modulus of the difference between any two samples in the total sample.
Generally, the Gini coefficient is required to be less than or equal to the specified Gini coefficient threshold when it is used as a constraint condition.
where G 0 is a specified Gini coefficient threshold. The Gini coefficient is used as the fairness evaluation index of generating hours. Therefore, the original revenue X is replaced by the generation utilization hours T of each unit: Specifically, the fairness constraints considered include: (1) Overall fairness constraint of the planned annual generation utilization hours of thermal power units: Specifically, the fairness constraints considered include: where G Planed is the Gini coefficient of planned generating hours of all thermal power units, G 0 is the Gini coefficient index threshold of planned generating hours of thermal power units set by the system, N is the number of thermal power units, and T Planed,i and T Planed,j are the generation hours of annual planned electricity energy of unit i and unit j respectively, and u is the average annual generating hours for all thermal power units.
where the power generation utilization hours of part of planned electricity of marketoriented units are equal to the sum of the planned power generation utilization hours of each month: where T m Planed,i is the generation hour of planned electricity for unit i in month m: (2) Fairness constraints on annual generation utilization hours of the same type of units: where k = 1, 2, . . . , N k , N k is the Gini coefficient of planned generating hour of type k, G 0 is the Gini coefficient index value of planned generating hours of the same type unit set by the system, and u k is the generation hour of planned electricity for unit type k.
(3) Fairness constraints on annual generation utilization hours of units with the same capacity range: where k = 1, 2, . . . , N C , N C is the number of unit capacity zones, N C k is the total number of units in the capacity zone k, G C Planed, k is the Gini coefficient of annual planned electric energy of units in the capacity zone k, G C 0 is the Gini coefficient index of annual planned electricity of units in the same capacity partition, and u C k is the average annual generation utilization hours of planned electricity of units in the capacity zone k: (4) Fairness constraint of planned electricity generation utilization hours after deducting market capacity: where, G A conversion is the Gini coefficient of power utilization hours of each generator set after deducting market capacity, T conversion, i is the generating utilization hours of each generating set after deducting the market capacity of generating set i, G A conversion,0 is the Gini coefficient threshold. When setting the actual Gini coefficient threshold, G A conversion,0 can be set to a larger value to encourage generating units to participate in market transactions.
The market capacity deducted by the market unit shall be considered as the average of the market capacity converted into each month.
Among the above four constraints in the Gini Coefficient fairness constraint set, the overall fairness constraint for all units is mandatory. The fairness constraint for units of the same type and of the same capacity partition is optional. In general, the Gini coefficient of the same type or capacity partitioned units is set at a lower value than the overall Gini coefficient, i.e., its fairness requirement is higher than the overall fairness requirement. In practice, the fairness constraint can be increased or decreased for each classification of units as needed.

Monthly Constraints
(1) Electricity balance constraint where, Q m sum is the total electricity demand of the system in month m, Q m w , Q m H , Q m T and Q m O are the forecast month m for wind, hydro, nuclear and other types of power generation according to the operating arrangements.
(2) Market-based unit capacity commutation constraint where, C m Market,i is the converted market electricity capacity of marketed unit i in month m, D m Market,max,d is the maximum load in month m of directly traded user d of market-based unit i, Nd i is the number of users who carry out direct transactions with market-based unit i, k d peak,i is the matching factor between the maximum load of the directly traded user d of the marketed unit i and the maximum load of the system, and k adj is the adjustment factor to reflect the electricity consumption characteristics of the power system and the need for peaking reserve. Usually, the k adj setting is the same for all direct transaction users. k Market is the market participation incentive factor, which takes a value of 1 to indicate neither encouragement nor opposition, while a value greater than 1 indicates a higher degree of encouragement to participate in the market.
(3) Marketed unit market capacity in relation to planned capacity constraints where, C m Planed,i is the planned electrical generation capacity of unit i in month m. (4) Constraints on the relationship between market-based unit plans and market power where, Q m Planed,i is the planned electricity consumption of unit i in month m, Q m Market,i is the amount of electricity traded directly by unit i in month m.
(5) Market-based unit market power balance constraints where, Q d Market,i is the annual direct trading of electricity between unit i and customer d.
where, T m Market,i is the number of hours of market electricity generation utilized by marketed unit i in month m.
When all market power is not liberalized, i.e., under the planned power system, it is unnecessary to distinguish between market and non-market units. Moreover, it is unnecessary to distinguish between market and planned power components. Fairness constraint concentration does not need to consider constraint (4), that is, the fairness constraint of planned electricity generation utilization hours after deducting market electricity capacity, and the rest is not different from the mathematical model under "planned electricity" and "market electricity" dual track operation.

Model Solving Methods
This chapter models a non-linear programming problem with only finite non-linear constraints due to the addition of the Gini coefficient non-linear constraint on the fairness of the electricity schedule. The nonlinear constraint is a nonlinear constraint with absolute values, and can be solved directly with branch-and-bound and Lagrangian dual methods. In recent years, commercial optimization software, such as IBM's CPLEX optimization package, has been widely used to solve various linear and nonlinear programming problems in an efficient and simple manner.

System Description
A system with 20 thermal power units is used as an example for the simulation analysis. The parameters are derived from the actual operating data of a provincial network, with necessary modifications. The parameters of each unit are shown in Table A1 in the Appendix A.
According to the compilation principle, the sulfur content is assumed to be 2%, the total annual electricity demand of thermal power units is 14,950 GWh, the standard coal to raw coal factor is 1.4017, and the average sulfur content of raw coal is 0.02. Considering the fairness constraint of annual power generation utilization hours of units in the same capacity range, the capacity range of units is divided into three bands: 0-100 MW, 101-200 MW, and 201-500 MW. The capacity range is divided into three zones: 0-100 MW, 101-200 MW, and 201-500 MW.

Analysis of Calculation Results in Fully Scheduled Electricity Mode
Assuming that all units do not participate in market trading under the full planed electricity mode, two sets of fairness calculation scenarios are considered and simulated: (1) the overall Gini coefficient, the Gini coefficient of the same type of unit, and the Gini coefficient of the same capacity range unit are all taken as 0.45; and (2) the overall Gini coefficient is taken as 0.3, the Gini coefficient of the same type of unit is taken as 0.2, and the Gini coefficient of the same capacity range unit is taken as 0.1. The results of the annual electricity plan for each unit under the two scenarios are shown in Table A2 in the Appendix A. The simulations were programmed, run, and data processed on a PC with an Intel(R) Core(TM) i3-4160 CPU @ 3.60 GHz and 8.0 GB of running memory, with two sets of simulations taking 11.84 s and 12.35 s respectively.
As can be seen from the results in Table A2, the difference in generation hours utilized by each unit in Scenario 1 is greater than in Scenario 2 because the Gini coefficient for Scenario 1 is larger and the requirement for fairness in generation hours utilized by each unit is more relaxed. Therefore, units with low coal consumption and low emissions can be allocated more power generation hours and electricity, while units with high coal con-sumption and high emissions can be allocated less power generation hours and electricity, resulting in better overall power generation efficiency of thermal power units.
The annual generation scheduling of the calculation system is developed using the dynamic polarization method [23], which is based on the following principles.
Thermal power units are classified according to the following principles, with each class having the same annual utilization hours. The difference in annual utilization hours is set between the different classes of units. The difference in utilization hours is 50 h for each priority level. (1) Classification according to the coal consumption of the supply: 1st Class, above 360 g/kWh; and 2nd Class, below or equal to 360 g/kWh. Each class has a priority difference of 1. (2) Classification according to installed capacity: 1st Class, units with a single capacity of up to 100 MW; 2nd Class, units with a single capacity of 100 to 300 MW; and 3rd Class, units with a single capacity of 300 to 600 MW. Each class has a priority difference of 2. (3) Divided according to environmental differences: 1st class, FGDs installed; and 2nd class, no FGDs installed. The priority difference for each class is 1. The base utilization hours are set at 1200 h. The results of the annual power plan, total coal consumption and total emissions for each unit calculated by the dynamic polarization method are shown in Table A3 in the Appendix A.
As can be seen from the results in Table A3, the application of the dynamic polarization method basically divides all units into a number of classes, with the number of hours utilized by the units within the same class being identical and the difference in hours utilized between classes increasing in equal proportions, which is a relatively sloppy allocation. The results of two simulation scenarios with the method presented in this paper and the dynamic polarization method are shown in Table 2. Under the three scenarios, the total coal consumption of all thermal power units was 4.205 million tonnes, 4.2611 million tonnes, and 4.3279 million tonnes, respectively; the total SO 2 emissions were 685.15 tonnes, 709.63 tonnes, and 784.63 tonnes, respectively; and the overall Gini coefficients for each unit were 0.45, 0.3, and 0.54, respectively. It can be seen that the total coal consumption and emissions of power generation using the dynamic polarization method are higher compared to the results of the two methods in this paper. In the case of high coal consumption, the overall Gini coefficient value generating hours obtained by the dynamic polarization method is higher than that by the method presented in this paper. This means that due to the total electric energy allocation is extensive, poorer energy saving and emission reduction effect has been achieved with the dynamic polarization method compared with the method presented in this paper. This would explain why there are always many doubts about the dispatching fairness from most of the energy producers. By using the method proposed in this paper, the overall Gini coefficient is varied to obtain the total number of hours utilized by each unit for various values of the overall Gini coefficient as shown in Figure 3. When the overall Gini coefficient is 0.1, the requirement for fairness in the number of hours generated by each thermal unit is very strict, and the difference in the number of hours generated by each thermal unit is small, and the number of hours generated by Units 3 and 20 is equal. As the Gini coefficient increases, the equity requirement is gradually relaxed and the fluctuation in the number of hours of generation use between units increases. When the Gini coefficient is greater than 0.45, the number of hours utilized by each unit remains more or less the same, which means that the Gini coefficient of each unit is around 0.45 when the fairness constraint of each unit is not considered.
By using the method proposed in this paper, the overall Gini coefficient is varied to obtain the total number of hours utilized by each unit for various values of the overall Gini coefficient as shown in Figure 3. When the overall Gini coefficient is 0.1, the requirement for fairness in the number of hours generated by each thermal unit is very strict, and the difference in the number of hours generated by each thermal unit is small, and the number of hours generated by Units 3 and 20 is equal. As the Gini coefficient increases, the equity requirement is gradually relaxed and the fluctuation in the number of hours of generation use between units increases. When the Gini coefficient is greater than 0.45, the number of hours utilized by each unit remains more or less the same, which means that the Gini coefficient of each unit is around 0.45 when the fairness constraint of each unit is not considered.  Figure 4 shows the total cost of coal consumption for various values of the overall Gini coefficient. As can be seen, the total coal cost decreases as the Gini coefficient increases. When the Gini coefficient value is small, the total cost of coal consumption decreases more quickly as the Gini coefficient increases, and when the Gini coefficient value is greater than 0.4, the total cost of coal consumption decreases slowly as the Gini coefficient increases. Therefore, in this example, 0.4 is the recommended value for the overall Gini coefficient (0.3-0.4 is a fair range as determined by the UN).  Figure 4 shows the total cost of coal consumption for various values of the overall Gini coefficient. As can be seen, the total coal cost decreases as the Gini coefficient increases. When the Gini coefficient value is small, the total cost of coal consumption decreases more quickly as the Gini coefficient increases, and when the Gini coefficient value is greater than 0.4, the total cost of coal consumption decreases slowly as the Gini coefficient increases. Therefore, in this example, 0.4 is the recommended value for the overall Gini coefficient (0.3-0.4 is a fair range as determined by the UN).  It can be seen from the simulation results that the use of this method can reasonably consider the fairness of the annual generation hours of each thermal power unit, while keeping the system parameters constant, which is acceptable to the major power produc ers. At the same time, the method can better control the coal consumption and pollutan emissions of power generation, which is important for the implementation of the impar tial and open dispatching that takes into account coal consumption and environmenta benefits. It can be seen from the simulation results that the use of this method can reasonably consider the fairness of the annual generation hours of each thermal power unit, while keeping the system parameters constant, which is acceptable to the major power producers. At the same time, the method can better control the coal consumption and pollutant emissions of power generation, which is important for the implementation of the impartial and open dispatching that takes into account coal consumption and environmental benefits.

Analysis of the Results of the "Planned Electricity" and "Market Electricity" Dual-Track System
Assuming that the odd-numbered units are marketed, the amount of electricity traded directly by each unit is shown in Table 3. In order to simplify the calculation and analysis, it is assumed that a customer with the same market-based unit is equivalent to one customer, i.e., one market-based unit for one directly traded customer. The maximum monthly load for each direct trading customer is shown in Figure 5, and the forecast for other types of units for each month is shown in Appendix A Table A4.  It is assumed that the regulation coefficient kadj is 0.9, the participation incentive coefficient kMarket is 1.1, and the overall Gini coefficient is 0.3; the Gini coefficient of units of the same type is 0.2, and the Gini coefficient of units of the same capacity range is 0.1. The Gini coefficient index of power generation utilization hours of planned electricity after deducting market capacity is 0.6. The converted capacity of each market-based unit for direct trading is shown in Table A5 in the Appendix A.
The total annual electricity consumption, planned electricity consumption, planned electricity generation hours, and total electricity generation hours utilized by each unit are shown in Table A6 in the Appendix A. The overall Gini coefficient value was changed to obtain the total generation hours utilized by each unit for various Gini coefficient values as shown in Figure 6. When the overall Gini coefficient for planned electricity is 0.1, there is a very strict requirement for fairness in the number of generation hours used by each thermal unit. The difference in generation hours between thermal units is relatively small, but there is still some difference between marketed and non-marketed units. As the Gini coefficient value increases, the equity requirement is gradually relaxed and the fluctuation in the number of hours of generation use between units increases. Similar to the results of the conventional scheduling model, the Gini coefficient value is greater than 0.45, and the number of hours utilized by each unit remains the same. This means that the Gini coefficient value of each unit is around 0.45 when the fairness constraint of each unit is not considered. It is assumed that the regulation coefficient k adj is 0.9, the participation incentive coefficient k Market is 1.1, and the overall Gini coefficient is 0.3; the Gini coefficient of units of the same type is 0.2, and the Gini coefficient of units of the same capacity range is 0.1. The Gini coefficient index of power generation utilization hours of planned electricity after deducting market capacity is 0.6. The converted capacity of each market-based unit for direct trading is shown in Table A5 in the Appendix A.
The total annual electricity consumption, planned electricity consumption, planned electricity generation hours, and total electricity generation hours utilized by each unit are shown in Table A6 in the Appendix A. The overall Gini coefficient value was changed to obtain the total generation hours utilized by each unit for various Gini coefficient values as shown in Figure 6. When the overall Gini coefficient for planned electricity is 0.1, there is a very strict requirement for fairness in the number of generation hours used by each thermal unit. The difference in generation hours between thermal units is relatively small, but there is still some difference between marketed and non-marketed units. As the Gini coefficient value increases, the equity requirement is gradually relaxed and the fluctuation in the number of hours of generation use between units increases. Similar to the results of the conventional scheduling model, the Gini coefficient value is greater than 0.45, and the number of hours utilized by each unit remains the same. This means that the Gini is a very strict requirement for fairness in the number of generation hours used by each thermal unit. The difference in generation hours between thermal units is relatively small, but there is still some difference between marketed and non-marketed units. As the Gini coefficient value increases, the equity requirement is gradually relaxed and the fluctuation in the number of hours of generation use between units increases. Similar to the results of the conventional scheduling model, the Gini coefficient value is greater than 0.45, and the number of hours utilized by each unit remains the same. This means that the Gini coefficient value of each unit is around 0.45 when the fairness constraint of each unit is not considered.  Changing the market participation incentive factor k Market from 1 to 1.5, the planned electricity and total generation hours of each unit are shown in Tables A7 and A8 in the Appendix A, respectively. The results in the table show that as k Market increases, the amount of planned electricity and total generation hours allocated to market-based units increases, while the amount of planned electricity and total generation hours allocated to non-marketbased units decreases, which achieves the purpose of adjusting the degree of incentive to participate in market trading by adjusting k Market .

Conclusions
This paper establishes an optimal model for annual electricity scheduling based on the fairness constraint set of the Gini coefficient. The annual electricity of each unit is the decision variable. It can realize the dual-track operation mode of "planned electricity" and "market electricity" for the power grid of high proportion of new energy and fully planned electricity. The annual electricity plan can be formulated as a reasonable manner under the dual-track operation mode of "planned electricity" and "market electricity", as well as the fully planned electricity mode. The theoretical and numerical analysis shows the following conclusions.
(1) Under traditional planned power dispatch model, the optimal power allocation result is obtained by solving a continuous optimization model. The Gini coefficient, a recognized fairness indicator, is used to constrain the overall power utilization hours of each unit. It can reduce energy, resource consumption, and pollutant emissions. A more fair annual power plan allocation scheme can be obtained. The Annual Electricity Plan (AEP) is designed to focus on different levels of equity by adjusting the Gini coefficient.
(2) Under the dual-track operation mode of "planned electricity" and "market electricity", the proposed model can take into account the uncertainties of new energy generation such as wind and photovoltaic in a more detailed manner by adding a monthly power balance constraint. It converts the capacity of "market electricity" to take into account the peak and frequency regulation reserve demand on a monthly basis and introducing a fairness constraint set of the Gini coefficient. The optimal allocation of annual power plan is achieved under the dual-track operation mode. By adjusting the Gini coefficient, a balance can be struck between fairness and affordability. By adjusting the market participation incentive coefficient, the level of incentive for market transactions can be effectively controlled.
In this paper, fairness problems for the power generation company, including annual power generation and peaking responsibility are considered. However, fairness of undertaking frequency regulation responsibility for all the generation units should also be considered, which will influence the market capacity conversion results. The annual