Multiobjective TOU Pricing Optimization Based on NSGA2

. Fast and elitist nondominated sorting generic algorithm (NSGA2) is an improved multiobjective genetic algorithm with good convergence and robustness. The Pareto optimal solution set using NSGA2 has the character of uniform distribution. This paper builds a time-of-use (TOU) pricing mathematical model considering actual constraint conditions and puts forward a new method which realizes multiobjective TOU pricing optimization using NSGA2. Avariety ofobjective TOU pricingschemes can beprovided for decision makers compared with traditional method. Furthermore, the multiple attribute decision making theory is applied in processing the Pareto optimal solution set to calculate the optimal compromise price scheme. The simulation results have shown that the TOU pricing scheme determined by the method proposed above can achieve a better effect of clipping the peak load to fill the valley load. Consequently, the study in this paper is innovative and is a successful exploration of coordinating the relation of various objective functions concerned in TOU pricing optimization problem.


Introduction
Time-of-use price is an important demand-side management (DSM) method, which could stimulate and encourage the consumers to change the consumption behaviors and achieve the purpose of clipping the peak load to fill the valley load as well as changing the load curve, so that the economy, security, and reliability of the power system are improved.So far, three kinds of TOU pricing are studied: the first is based on the relation between sale and retail prices [1], the second is based on the consumer response curve or price elasticity matrix [2,3], and the third is based on game control [4].However, the research on the optimization methods, which can decrease peak load and peak-valley difference at the same time, has not been studied enough.It is a typical multivariable constrained multiobjective optimization problem.
The traditional algorithms of multiobjective TOU pricing optimization are weighted method [5], constraint method [6], and so on.The multiobjective optimization problem is transferred into a single objective optimization problem.Then, the problem is solved by the relatively mature single objective optimization algorithm in which only one solution is obtained.However, there is a solution set in the multiobjective optimization problem.Moreover, the weights of the objectives are defined subjectively rather than objectively.So, the optimal extent is affected.Multiobjective genetic algorithm is an evolutionary algorithm.Its core thought is to coordinate the relation between the objective functions to search for optimal solution set [7].NSGA2, an improved multiobjective generic algorithm, overcomes the shortcomings mentioned above and has the characters of good convergence and excellent robustness.Its noninferior solutions are uniformly distributed in objective space.Considering TOU pricing is a multivariable constrained multiobjective optimization problem; this paper uses NSGA2 to solve the problem and puts forward a new kind of TOU pricing optimization mechanism which can achieve a better effect in balancing the relation between the objective functions.Case analysis shows that the proposed TOU pricing multiobjective optimization based on NSAG2 takes on better effectiveness and superiority.

Multiobjective Optimization Pareto Optimal
Solution Set and Pareto Genetic Algorithm

The Concept of Nondominated Solution (Pareto Solution).
Given two solutions, X  and X V , if all the objectives of X  are superior to the objectives of X V , we say that X  dominates X V or X V is dominated by X  .If X  is not dominated by the others, it is called nondominated solution or Pareto solution.
The set of nondominated solutions is called Pareto front.The solution of Pareto front is not dominated by the ones outside the Pareto front.Compared with other solutions, Pareto solution has less objective conflicts; thus, a better selection space is supplied for the decision maker.

The Advantages of Pareto
Front and NSGA2.For a multiobjective optimization problem, the optimal solution is not unique, and it is a set of many satisfied solutions, which is called Pareto optimal set.Pareto optimal set can provide more than one optional solution for decision makers.According to the practical needs and combined with the theory of multiple attribute decision making, decision maker can choose an optimal compromise solution [8], which owns the highest satisfaction.
In the multiobjective optimization, there may be conflicts between different objectives and it is hard to make a judicious decision by comparing all these objectives.If a solution is the best one based on a certain objective, it may be weakly worse for the other objectives.The core of multiobjective genetic algorithms is to coordinate the relation between the objective functions.Nondomination rank and selection operator are used to make the population evolve towards Pareto optimal set in the formulation of the problem.Then, the Pareto optimal set is obtained, which can make the objective function relatively maximized (or relatively minimized).Thus, several kinds of multiobjective genetic algorithms based on Pareto optimal set come into being, in which nondominated sorting generic algorithm (NSGA) [9] is a typical example.NSGA2 [10] is the improvement of NSGA which is proposed by Deb and others.Its prominent characteristics can be listed as follows: (1) NSGA2 is based on nondomination rank, which makes the complexity decrease from ( 3 ) of NSGA to ( 2 ); (2) crowding distance comparison operator is adopted for optimization.So that individuals in the population can extend to the entire Pareto domain and maintain the diversity of the population to avoid local convergence; (3) elitist strategy is used in NSGA2.The next generation is got from the sorting and competition of parent generation and offspring generation.Thus, the quality of the next generation is guaranteed.

The Mathematic Multiobjective Optimization Model
According to the typical daily load curve, peak, flat, and valley are defined as   ,   , and  V , respectively, and their corresponding prices are   ,   , and  V .Consumer response curve is  = ().The average price or market clearing price before the implementation of TOU pricing is  0 which is assumed as a known constant.

The Selection of Multiobjective Functions.
In this paper, minimizing the maximum load and minimizing the peakvalley difference are considered simultaneously as TOU pricing objective functions.
Objective function 1 is as follows: minimizing the maximum load min ( max ) . ( Objective function 2 is as follows: minimizing the peakvalley difference where where   ,   , and  V are the peak, flat, and valley loads before the implementation of TOU pricing and    ,    , and   V are the peak, flat, and valley loads after the implementation of TOU pricing.
Considering the production cost and the consumers' endurance, the prices of peak, flat, and valley loads are limited by the price supervision department.The limitations are as follows: where   min ,   max ,   min ,   max ,  V min , and  V max are the upper and lower limit prices of peak, flat, and valley loads.
The price of peak, flat, and valley loads varies within limits, and the effect of peak load shifting varies at the same time.This paper mainly studies the optimization of the effect of peak load shifting.  ,   , and  V are obtained according to the optimization of the mathematical model.They are substituted into the following inequality: where  =   / 0 is interest concessions coefficient,  0 is the total electricity cost before the implementation of TOU pricing, and   is the reduced electric power cost of the power company after the implementation of TOU pricing.The inequality is used to verify the model and  is obtained through the inequality, whose upper limitation is 10%.It is a reference for decision makers to optimize the economic indicator.

Start
Set the initial parameters of TOU pricing model.

Generate the Initial Population.
When generating a random initial population of individuals, the constraint conditions of price are used to limit the random number, to ensure that the random number is generated within the constraints.Then, set the appropriate number of initial population of individuals.It cannot be too small to converge or too big to prolong search time.It also should be selected randomly to ensure the diversity and the convergence to the optimal solution.

Fast Nondomination Sorting.
Before the process of selection operator for individual populations, according to the level of individual nondomination solutions, the population is sorted by fast nondomination sorting of NSGA2.The specific methods are as follows.
In the algorithm, two parameters,   and   , of individual  should be calculated.  is the number of individuals that dominate individual .  is the set of individuals that are dominated by individual .The specific steps of fast nondomination sorting are as follows.
(1) Find all the individuals whose parameter   is 0. The individuals are preserved in the set , which is the first nondomination layer.
( (3) Repeat the previous step until the entire population is sorted.
For a population, whose number of objective functions is K ( > 1) and size is N, two parameters,   and   , are calculated for each individual by NSGA2.Traversing the entire population, the total computational complexity is ( 2 ).Nondomination sorting is used in NSGA, and its computational complexity is ( 3 ).So, the calculation efficiency of fast nondomination is higher.

Virtual Fitness.
In NSGA2, the individual fitness includes nondomination rank and virtual fitness.Crowding distance is used by NSGA2 to form individual fitness to maintain the diversity of the individual and prevent the accumulation of the individual in the local.It is defined as the sum distance of the two points on either side of this point along each of the objectives.In Figure 2, the crowding distance of the th solution in its front (marked with solid circles) is the sum of sidelength of the cuboid (shown with a dashed box), which is [ 1 ( + 1) −  1 ( − 1)] + [ 2 ( − 1) −  2 ( + 1)].This approach makes the selected individuals distribute evenly in the objective space with good robustness.

Selection
Operator.Genetic selector operator is for optimization towards Pareto optimal solution set and to make the solutions distribute evenly.The best individual with the greatest probability of survival is preserved by selection operator to avoid the loss of good genes.Thereby, global convergence and computational efficiency are improved.Each individual has two parameters,  rank (nondomination rank) and   (crowding distance).Binary tournament selection operator is used in this paper.Two individuals are chosen randomly for comparison.If  rank <  rank or ( rank =  rank and   >   ), individual  is better than individual .That is to say, if two individuals are in different layers, the individual with the lowest nondomination rank is selected, which is firstsorted in the process of nondomination sorting with high level of nondomination.Or if two individuals are in the same layer, the individual with the largest crowding distance is selected, which has advantages in selection.

Crossover and Mutation.
The cooperation of crossover and mutation makes genetic algorithm perform well in local and global search.SBX crossover operator and polynomial mutation operator are adopted.SBX operator simulates binary crossover process for real-coded parent individuals.For a given random cross-point, switch two sides of the crosspoint of the parent of the individual parts.Polynomial mutation operator is based on polynomial distribution.Offspring generation   is generated through crossover and mutation operator.

Elitist Strategy.
Elitist strategy makes excellent individuals retain to the offspring generation directly.It is a necessary condition for the genetic algorithm converging with probability 1.The method is as follows.
(1) Parent generation   and offspring generation   are combined into a population   =   ∪   .
(2) Sort   by fast nondomination sorting and calculate crowding distance of every individual.Select individuals one by one based on nondomination rank, until the number of individuals reaches .Then, the new parent generation  +1 is generated.
Iteration is done again and again until convergence is reached.Then, fast nondomination sorting is used to obtain the Pareto optimal solution set, that is, Pareto front, which consists of the individuals in the first nondomination layer.

Optimal Compromise Solution.
The optimal compromise solution is the solution in the Pareto front with the greatest satisfaction.
In Pareto optimal solution set, the satisfaction of each objective function of each Pareto solution is represented by fuzzy membership function as follows: where   is the objective function and  = 1, 2, . . .,  obj . obj is the number of object functions.  min and   max are the maximum and minimum of the th objective function.
According to formula (8), satisfaction of each Pareto solution is calculated.Then, the solution with the greatest ℎ value is chosen as the optimal compromise solution

Case and Analysis
In order to verify the effectiveness and feasibility of NSGA2 in multiobjective TOU pricing optimization, large industrial consumers in a specific region are chosen to be optimized by NSGA2.The result is analyzed subsequently.

The Data before TOU Pricing Optimization
(1) Typical Daily Load Data.The typical daily load data is shown in Table 1.
(2) The Price and the Cost of Power Consumers.In the region, the average price or market clearing price is  0 = ¥0.5647/kW⋅ h.The cost of large industrial power consumers is  0 = ¥8.80 × 10 6 . (
Through investigation, consumer response curve is obtained, whose expression is  = () = −0.2807+ 1.2952, where  is the per unit value of price and y is the per unit load value of consumer.
(4) Constraint Conditions of Peak, Flat, and Valley Price.In the region, according to the power cost and consumers' endurance, the constraint conditions of peak, flat, and valley price are regulated as follows (shown in per unit value): 0.3 ≤  V ≤ 0.8, 0.8 ≤   ≤ 1.2, 1.2 ≤   ≤ 1.8.

The TOU Pricing Optimization by NSGA2.
In the calculation case, the parameters settings of NSGA2 are as follows: the population size is 500; the iterative time is 500; the crossover rate is 0.9; the mutation rate is 0.1.In NSGA2, the selection basis of individuals is the indexes of nondomination rank and crowding distance.For the individuals in parent generation, objective functions are calculated by the mathematic multiobjective optimization model.Then, 250 individuals (/2 = 250) are selected according to nondomination rank and crowding distance.They are used to crossover and mutate to generate the offspring generation.The parent generation and the offspring generation are combined together to generate the new parent generation by elitist strategy.
The case takes the 10th generation for example.Table 2 shows the individuals in the first 3 layers sorted by fast nondomination sorting.The distribution of these individuals in the objective space is shown in Figure 3.The crowding distance of the individuals in the Pareto front is shown in Table 3. From Figure 3, it is known that the accumulation degree of each nondomination layer is different, so that the global optimization is ensured.Crowding distance is taken into consideration when selecting individuals to ensure the diversity and avoid the local convergence.
When the convergent conditions are met, the iteration goes to an end.Then, the final Pareto optimal solution set (Pareto front) is obtained.A part of the representative Pareto optimal solutions is shown in Table 4.The distribution of Pareto front and the optimal compromise solution are shown in Figure 4.The optimal compromise price scheme  4 that the two objectives, the minimum peak load and the minimum peak-valley difference, have a contradiction.It is impossible to find a solution which is optimal for two objectives at the same time.Therefore, decision makers can make choices by objective needs in Pareto optimal solutions.If the minimum peak load is used as the main objective, the solution is selected from the upper left part of the Pareto solution set, such as the first solution in Table 4.If the minimum peak-valley difference is used as the main objective, the solution is selected from the   4. If the preference objective function is not specified, the optimal compromise solution is selected, just at the marked point in Figure 4.The load curve of optimal compromise solution is shown in Figure 5. Multiobjective optimization is transferred into single objective in traditional algorithm.However, the weights of each objective function are artificially regulated.The blindness and subjectivity of traditional optimization algorithm are avoided by NSGA2, thus realizing multiobjective TOU pricing optimization.
(2) Comparison with Traditional Method.The weighted method is a traditional single objective optimization method.When used for solving the multiobjectives optimization problem, the weighted method would transfer the problem into a single objective optimization problem by assigning a set of weight value to each objective.In this paper, The weighted method is also used in optimizing TOU pricing problem considering setting different weight values ( 1 and  2 ). 1 is the weight value of the objective function 1, the maximum load. 2 is the weight value of the objective function 2, the peak-valley difference.Besides, 1 +  2 = 1.Compared with NSGA2, the simulation results of these two methods are shown in Table 5.
It can be seen in Table 5 that different weight value distribution in the weighted method would lead to different solutions.Generally, the bigger the weight value of an objective is, the better the solution of corresponding objective is, since this objective is emphasized.When an objective is The load (MW) The original load The load of NSGA2 pricing overemphasized, the corresponding objective function value is superior to that of NSGA2, while the other objective function values would be much worse, just as the case of  1 = 0.1  2 = 0.9 shown in Table 5.When the weight values are chosen improperly, all the objective function values are worse than those of NSGA2, just as the case of  1 = 0.5  2 = 0.5 shows.For the weighted method, it is important but difficult to choose suitable weight of objectives because the weights are determined subjectively.By contrast, the optimal compromise solution of NSGA2 is obtained objectively.In conclusion, the TOU pricing scheme of NSGA2 has a significant advantage on clipping the peak load to fill the valley load.

Conclusions
NSGA2 is a newly developed algorithm.For multiobjective TOU pricing optimization, a new thought is put forward in this paper.NSGA2 is applied in TOU pricing, and more than one Pareto solutions are obtained.Choices are supplied for the decision makers for selection, which overcomes the shortcomings of the traditional algorithm and realizes the multiobjective TOU pricing optimization.Example analysis demonstrates that the remarkable effect of clipping the peak load to fill the valley load is achieved by NSGA2.Good technique support is supplied for power market operation which is affected by the industrial and economic structure.
) For each individual  of set , traverse each individual  in the set   and execute   =   −1.If   = 0, preserve individual  in the next nondomination layer.

Figure 2 :
Figure 2: Chart of local crowding distance.

Figure 3 :
Figure 3: The distribution of individuals in the first three layers in objective space.

Figure 5 :
Figure 5: Contrast of load curve between optimization by NSGA2 and the original one.
max = max[(  ,   ,   ,   ,   ,   )] is the maximum of the peak load of the day,  min = min[(  ,   ,   ,   ,   ,   )] is the minimum of valley loads,   ,   , and  V are the time price of peak, flat, and valley loads,   ,   , and  V are the time period of peak, flat, and valley loads, and (  ,   ,  V ,   ,   ,  V ) is the new load after the implement of TOU pricing according to the consumer response curve.

Multiobjective TOU Pricing Optimization Based on NSGA2
0 is used as the reference.  ,   , and  V are calculated into per unit value   ,   , and  V , which are used as control variables.The coding format in real-coded pattern is [  ,   ,  V ].The real decoding value of peak, flat, and valley price is as follows: t Figure 1: Flow chart of optimizing time-of-use price based on NSGA2.4.

Table 1 :
Data of the large industrial load.

Table 2 :
The first 3 layers of the 10th generation (accurate to 0.001).

Table 3 :
The crowding distance of the individuals in Pareto front.

Table 4 :
A set of Pareto optimal solutions.

Table 5 :
Contrast of NSGA2 and weighted method.