Next Article in Journal
Turbulence-Model Comparison for Aerodynamic-Performance Prediction of a Typical Vertical-Axis Wind-Turbine Airfoil
Next Article in Special Issue
A Two-Stage Feature Selection Method for Power System Transient Stability Status Prediction
Previous Article in Journal
Velocity-Controlled Particle Swarm Optimization (PSO) and Its Application to the Optimization of Transverse Flux Induction Heating Apparatus
Previous Article in Special Issue
Optimal Decision-Making to Charge Electric Vehicles in Heterogeneous Networks: Stackelberg Game Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Cross-Entropy-Based Hybrid Membrane Computing Method for Power System Unit Commitment Problems

1
School of Electric Power, South China University of Technology, Guangzhou 510640, China
2
State Grid Ganzhou Electric Power Supply Company, Ganzhou 341000, China
3
Guangzhou Power Supply Bureau Co., Ltd., Guangzhou 510620, China
*
Author to whom correspondence should be addressed.
Energies 2019, 12(3), 486; https://doi.org/10.3390/en12030486
Submission received: 7 January 2019 / Revised: 30 January 2019 / Accepted: 31 January 2019 / Published: 2 February 2019
(This article belongs to the Special Issue Optimization Methods Applied to Power Systems)

Abstract

:
The cross-entropy based hybrid membrane computing method is proposed in this paper to solve the power system unit commitment problem. The traditional unit commitment problem can be usually decomposed into a bi-level optimization problem including unit start-stop scheduling problem and dynamic economic dispatch problem. In this paper, the genetic algorithm-based P system is proposed to schedule the unit start-stop plan, and the biomimetic membrane computing method combined with the cross-entropy is proposed to solve the dynamic economic dispatch problem with a unit start-stop plan given. The simulation results of 10–100 unit systems for 24 h day-ahead dispatching show that the unit commitment problem can be solved effectively by the proposed cross-entropy based hybrid membrane computing method and obtain a good and stable solution.

1. Introduction

The unit commitment (UC) problem is a typical optimization problem for power systems. The main goal of UC is to schedule the start-stop state of units and generate power according to the load forecasting curve during the dispatch period, with the corresponding constraints so that the cost is minimized [1]. Usually the UC problem can be broken down into two sub-problems: the unit start-stop plan and economic dispatch [2].
Mathematically, the UC problem is a high-dimensional, non-convex and mixed-integer nonlinear programming problem. Its discrete and continuous variables, non-convex objective function and network constraints enhance its non-convexity and complex [3]. Moreover, with the increase in unit and calculation scale, it is difficult to obtain an accurate feasible solution in a reasonable time frame. Therefore, many methods have been proposed by scholars to solve the UC problem, which can be roughly divided into three categories: heuristic methods, mathematical optimization methods, and intelligent optimization methods.
Heuristic methods are represented by the priority list method [4], the earliest method applied to solve the UC problem, which generally sorts by some economic indicators with small and simple calculations, usually relies on the actual scheduling experience.
Mathematical optimization methods include mixed integer nonlinear programming, the Lagrangian relaxation (LR) method [5], etc. Mixed integer nonlinear programming methods include the branch-and-bound (BB) [6], Benders decomposition [7], and other methods, with the decomposition technique generally used to simplify the problem; the solving efficiency has been rapidly improved with improvements in mathematical optimization software [8]. For the dynamic programming method, the global optimal solution can be obtained with no special requirements on the behavior of the objective function; however, there will be a “dimensionality disaster” [9] when the number of units is large, and to simplify the problem, the optimal solution will be lost when the approximation method adopted. Moreover, it is difficult to consider time-dependent constraints such as the ramp rate constraint. Compared with the dynamic programming method, the LR method has an advantage for large-scale problems since the calculation complexity is approximately proportional to the unit scale; in addition, the Lagrangian multiplier is of practical economic significance, but it cannot prove whether the solution is optimal due to the dual gap. It is also inflexible when considering some kinds of constraints, such as the ramp rate limit, and the possible oscillation and singularity during iteration may lead to convergence difficulties [10].
Bioinspired optimization methods are algorithms that simulate biological evolution or its behavior, and include genetic algorithms (GA) [11], particle swarm algorithm (PSO) [12], and memetic algorithm (MA) [13]. An approximate optimal solution can be obtained with no special requirements on the behavior of the objective function, and feasible solutions can be obtained even when the unit scale is large where no feasible solution can be obtained by other methods; however, the solving efficiency is affected by how the constraints are processed since these methods are essentially unconstrained.
Membrane computing is a computational framework inspired by the living cell and its organization in tissues and other higher order structures, and was first proposed by Gheorghe Pǎun, an academic at the European Academy of Sciences. In recent years, optimization problems such as image processing [14], robot path planning [15], DNA sequence design [16], gasoline blending scheduling [17], the travelling salesman problem [18], and the minimum storage problem [19] have been successfully solved by this framework. In this paper, we propose a cross-entropy-based hybrid membrane computing (CEHMC) method to solve the UC problem by combining the genetic algorithm-based P system (GAPS) with the biomimetic membrane computing (BMC) method. The GAPS is based on the binary genetic algorithm. It is nested by multiple membranes with the unit start-stop state (0,1 binary variables) as evolution objects. The evolution rules adopt the crossover and mutation rules of binary coded genetic algorithms. The difference between GAPS and GA is the communication of optimal objects in different membranes, which means the optimization results of the outer membrane can be transmitted into the inner membranes constantly when genetic rules are executed in each membrane. The calculation method of biomimetic membrane is inspired by the important role of Golgi apparatus in living cells. With the inner membrane system as the calculation carrier, the required regional structure is designed to simulate the static or dynamic membrane structure, the biochemical or physical reactions in the inner membrane system are simulated by various evolutionary rules, so that the top-ranking evolution objects after evaluation are selected for transmission according to the permeability of substances in the transporting through the membrane. Based on the combination of GAPS with the biomimetic membrane computing method, the cross-entropy (CE) optimization method is introduced to strengthen the searching ability during the optimization.
Because the UC problem is not fully standardized and specified, we will not attempt to give a definitive algorithmic solution to the problem. Instead, our goal in this paper is to demonstrate that the proposed CEHMC approach is viable. We include a case study showing than our method outperforms other nonlinear optimization alternatives for a case study. Therefore, our formulation of the UC problem and its solution algorithm are for illustration purposes only. For that purpose and for the ease of exposition and understanding, the classical UC formulation which has been commonly used in the past decades is used in this paper. The simulation results of 10–100 unit systems for 24 h day-ahead dispatching showed that the UC problem could be solved effectively by the proposed method, and a satisfactory and stable solution was obtained.

2. The Mathematical Model for UC Problems

2.1. Objective Function

The total generation cost FC during the dispatching period is taken as the optimization goal, including fuel cost and startup cost of all units. Thus, the objective function is:
min F C = t = 1 T i = 1 N [ I i , t F i , t + ( I i , t I i , t 1 ) 2 C i , t ]
where T is the total dispatch period; N is the number of generators; Ii,t is the state of unit i at tth hour; value 1 represents startup, while 0 represents shutdown; Fi,t and Ci,t are the fuel cost and startup cost, respectively, of unit i at t-th hour; and the fuel cost Fi,t can be expressed by the quadratic function as follows:
F i , 1 = A i , 2 P i , t 2 + A i , 1 P i , t + A i , 0
where Pi,t is the generation of unit i at t-th hour; and Ai,2, Ai,1, and Ai,0 are the unit cost coefficients of unit i.
Startup cost Ci,t is related to the unit OFF time, which can be expressed by the following step function:
C i = { C h o t , i T _ i o f f T i , t o f f T _ i o f f + T i c o l d C c o l d , i T i , t o f f > T _ i o f f + T i c o l d
where Chot,i and Ccold,i are the hot start cost and cold start cost, respectively, of unit i; T _ i o f f is the minimum downtime of unit i; T i , t o f f is the OFF time of unit i at t-th hour; and T i c o l d is the cold start time of unit i.

2.2. Constraints

The constraints of the UC problem usually include the system power balance, system spinning reserve requirement, generation limits, ramp rate limits, and minimum up and down time limits.
  • System power balance:
    i = 1 N P i , t = D t
  • System spinning reserve requirement:
    i = 1 N P ¯ i I i , t D t + R t
  • Generation limits:
    P _ i I i , t P i , t P ¯ i I i , t
  • Ramp rate limits:
    { P i , t P i , t 1 P i u p I i , t 1 + P i s t a r t ( I i , t I i , t 1 ) P ¯ i ( 1 I i , t ) P i , t 1 P i , t P i d o w n I i , t + P i s h u t ( I i , t 1 I i , t ) P ¯ i ( 1 I i , t 1 )
  • Minimum up and down time limits:
    { ( T i , t 1 o n T _ i o n ) ( I i , t 1 I i , t ) 0 ( T i , t 1 o f f T _ i o f f ) ( I i , t I i , t 1 ) 0
Here, Dt is the system load demand at t-th hour; P ¯ i and P _ i are the maximum and minimum generation limits, respectively, of unit i; Rt is the system reserve demand, typically 5–10% of the system load demand; P i u p and P i d o w n are the up and down ramp rate limit, respectively; P i s t a r t and P i s h u t are the start and stop rate limit, respectively; T i , t 1 o f f is the OFF time of unit i at the (t−1)-th hour; and T _ i o n is the minimum up time.

3. Unconstrained UC Model for CEHMC Method

As mentioned above, the traditional UC model is a multi-constrained nonlinear optimization problem, yet natural computing frameworks such as membrane computing are essentially unconstrained optimization methods, which are commonly used to solve unconstrained optimization problems. As for the complex multi-constrained optimization problem, using the penalty function is usually a good choice [20]. To facilitate the proposed method, the quadratic penalty function was used to deal with the constraints, transforming the UC problem into an unconstrained model without providing any initial feasible solution. Since GAPS is for the unit start-stop plan with embedded generic rules. These rules are design to screen the unit start-stop plans, i.e., testing whether the system spinning reserve requirement and minimum up and down time limits are satisfied, which determines the feasible unit state combination. Consequently, for the unconstrained UC model, only system power balance and ramp rate limits need to be transformed. The objective function can be rewritten as a penalty function as follows:
min Q = F C + μ 1 h max ( | i = 1 N P i , t D t | δ , 0 ) 2 + μ 2 g max { [ P i , t P i , t 1 P i u p I i , t 1 P i s t a r t ( I i , t I i , t 1 ) P i ¯ ( 1 I i , t ) , 0 ] 2 + [ P i , t 1 P i , t P i d o w n I i , t P i s h u t ( I i , t 1 I i , t ) P i ¯ ( 1 I i , t 1 ) , 0 ] 2 }
Here, the penalty factors μ1, μ1 > 0; δ > 0 is the acceptable violation domain when converting the equality constraint into an inequality constraint. In fact, the value of penalty factors is related to the magnitude of objective function and constraints. If the value is too large, it is likely to get bad solutions, and if the value is too small, the searching direction might be far away from the feasible region. Therefore, the method of testing is usually used to make the penalty item have the same or larger magnitude as the objective function so as to determine the penalty factors.

4. CEHMC Method Applied to Solve UC Problem

Three basic elements of membrane computing are membrane structure, evolution object, and evolution rules. The membrane structure and evolution rules can be designed according to requirements [21]. Taking an intracellular membrane system as the computational framework, the main procedures of the membrane computing method are as follow: first design the membrane structure, then evolve objects according to evolution rules, and lastly select the optimal objects for transmission and communication.
When using the CEHMC method for solving the UC problem, we utilize GAPS to schedule the unit start-stop plan since the start-stop states of generators are discrete variables. While for the dynamic economic dispatch problem, the generation of units is a continuous variable. Since the evolutionary objects (solution vectors) of biomimetic membrane computing is only required to be a real number, not limited to discrete variables. If the evolutionary rules are set to apply to continuous variables, biomimetic membrane computing can used to solve continuous variable optimization problems, such as the dynamic economic dispatch problem in UC. Then, the cross-entropy optimization method is introduced combined with biomimetic membrane computing to improve the searching ability during iteration. Finally, the minimum generation cost is taken as the optimal result of the UC problem.

4.1. GAPS for Unit Start-Stop Plan

Figure 1 shows the membrane structure of GAPS; it is a nested structure [22], where the symbol N represents the N-th nested genetic membrane (i.e., the variable N in Table 1). In practice, it needs to be determined synthetically according to unit size, number of constraints, calculation time and convergence effect, etc. The start-stop states of generators are taken as the evolution objects I, and the crossover and mutation rules of the binary encoding GA are the evolution rules.

4.2. Biomimetic Membrane Computing Method for Dynamic Economic Dispatch

Figure 2 shows the membrane structure of the biomimetic membrane computing method used in this paper. It is a reticular structure including an outermost membrane, m basic membranes, and a quasi-Golgi membrane [23] represented by G, and the output of generators is taken as the evolution objects P, that is, P = ( p 1 , p 2 , , p l ) , p k { P i , t ,   i = 1 , 2 , , N ;   t = 1 , 2 , , T } . In addition, m is corresponding to the variable Nb in Table 1, which is usually determined synthetically according to unit size, number of constraints, calculation time and convergence effect, etc.
The evolution rules can be divided into two categories according to the environment, that is, rules in the basic membranes, and rules in the quasi-Golgi membrane. Besides, correction rules can be used in all conditions. All the basic rules are set with the calculation probability as the execution condition. When the calculation probability is satisfied, the rule will be executed immediately. The evolution rules are as follows:
1) Crossover rules
  • Numerical crossover rule. Each element on both evolution object is numerically crossed, and the cross ratio of each element is different, which can be described as follows:
    P a , P b P a , P b { P a = η × P a + ( 1 η ) × P b P b = ( 1 η ) × P a + η × P b
    where P and P are the original and new object, respectively, generated after executing evolution rule; and η is a vector whose value is a random number uniformly distributed on [0,1].
  • Interval crossover rule. First, select the interval to be randomly exchanged, and then swap the elements in the same interval for two objects. This is described as follows:
    P a , P b P a , P b { P a = ( p a 1 , , p a p , , p a q , , p a l ) P b = ( p b 1 , , p b p , , p b q , , p b l ) P a = ( p a 1 , , p b p , , p b q , , p a l ) P b = ( p b 1 , , p a p , , p a q , , p b l )
    where p , q [ 1 , l ] , [ p , q ] is the interval to be swapped; and l is the length of object. The procedure of crossover can be divided into two or three steps according to the length of the object. When the length is short, two-step crossover is preferred: first, execute the numerical crossover rule on two selected adjacent objects, and then execute the interval crossover rule on two randomly selected objects. Conversely, when the length of object is long, three-step crossover is preferred: first, execute the numerical crossover rule on two randomly selected objects, then execute the interval crossover rule on the two objects, and lastly execute the numerical crossover rule on two adjacent objects.
2) Mutation rule
A random increment vector is added to the original object, which can be described as follows:
P P { P = P + φ   φ = h × r × ( P ¯ P _ )
where φ is the mutation vector; P ¯ and P _ are the upper and lower limits of P , respectively; h is the mutation coefficient; and r is the random vector that follows the standard normal distribution.
3) Correction rule
After some certain rules, such as the mutation rule, some elements of the new object may exceed their limits and must be modified, that is:
p = { p ¯ p > p ¯ p _ p < p _
For the quasi-Golgi membrane, activation conditions are set, which means the following rules are executed only if the quasi-Golgi membrane is activated:
• Target indication rule
The communication objects sent into the quasi-Golgi are sorted by their function value, and then each communication object is subtracted from the former stored object by this order. Thus, a direction vector is generated, and the new target indication vector is the sum of this direction vector and last target indication vector. Regardless of whether the quasi-Golgi is activated, the direction vector and target indication vector should be calculated, and the target indication vector should be reserved. Furthermore, the target indication object can be generated only if the quasi-Golgi is activated. This is described as follows:
P P S 1 , P S 2 { λ = λ + ( P s o r t n e w P s o r t o l d ) P S 1 = P + w × λ P S 2 = P w × λ
where λ is the target indication vector; P s o r t n e w and P s o r t o l d are the sorted communication objects stored in the quasi-Golgi this time and the last time; P S 2 and P S 1 both are the new target indication objects; and w is the indication coefficient.
• Transition rule
First select two elements of the selected object randomly, and then swap these elements if they are in the range of the other’s value, that is:
P P { P = ( p 1 , , p p , , p q , , p l ) P = ( p 1 , , p q , , p p , , p l ) p q p p p q ¯ p p p q p p ¯
• Abstraction rule
The abstraction rule is designed only for optimal and suboptimal objects, where each element of the optimal object is replaced by the element on the same position of the suboptimal object, one after another. Then, the element from the suboptimal object is reserved if the new object is better than the old optimal one; otherwise, the old optimal object remains:
P a , P b P { P a = ( p a 1 , , p a j , , p a l ) P b = ( p b 1 , , p b j , , p b l ) P = ( p 1 , , p j , , p l )     p j { p a j , p b j }
When using the biomimetic membrane computing method, there is an iterative calculation and communication for each basic membrane, and the good objects in the basic membrane can be reserved and sent to the quasi-Golgi for local optimization and evolution. However, with the increment of unit scale, the searching ability and stability of the method is reduced, and thus we added the CE optimization method to the membrane computing method to enhance its searching ability and stability.

4.3. CE Optimization Method

The CE method is an optimization method proposed by Rubinstein and Kroese, who estimated probabilities of rare events for stochastic networks [24]. It was first used to solve combinatorial optimization problems, and then also applied to solve continuous optimization problems [25]. Cross-entropy is a measurement for the similarity degree of two probability distributions. The main idea of CE is to get a probability distribution whose difference with the unknown optimal probability distribution is minimal [26].
For the dynamic dispatch problem:
γ * = min P ρ F ( P )
where γ* represents the optimization result of γ for economic dispatching; ρ is the probability space. Obviously, Equation (17) is a minimization problem. When generation P obeys the distribution of f ( ; u ) , the original optimization problem can be transformed into an optimization problem of finding the optimal probability density function f ( ; u ) . Moreover, when we define S { F ( P ) γ } as a set of different indicator functions on ρ with value γ , the optimization problem of f ( ; u ) can be further transformed into a corresponding probability estimation problem as follows. The indicator functions in S { F ( P ) γ } is used to describe the characteristics of a stochastic process, usually formulated as mean function, variance function and correlation function, etc.:
l ( γ ) = H u ( F ( P ) γ ) = P S { F ( P ) γ } f ( P ; u ) = E u S { F ( P ) γ }
where H u and E u are the probability measure and expectation, respectively, of the optimal probability distribution f ( ; u ) .
Usually the importance sampling method is used to calculate l ( γ * ) : sampling the generation P i ( i = 1 , 2 , , N ) based on the probability distribution g on P:
l ^ = 1 N i = 1 N S { F ( P i ) γ } f ( P i ; u ) g ( P i )
when the probability distribution g is:
g * ( P ) = S { F ( P ) γ } f ( P ; u ) l
where l ^ is the value of l after the importance sampling method; g* is the assumed value of g for the importance sampling method.
There is an unbiased estimated zero variance with only one sample needed. Therefore, it is difficult to calculate g * since its value is related to l , and the probability distribution g can usually be selected from the probability distribution cluster { f ( ; v ) } . In this way, the original optimization problem is finally transformed into a determination of parameter v that minimizes the difference between the probability distribution of generation f ( ; v ) and optimal g * . In this paper, parameter v includes the mean value μ and standard deviation σ of samples.
Relative-entropy (i.e., Kullback-Leibler distance) and cross-entropy are commonly used measures for the similarity degree of two probability distributions [27]. In this paper, we used the CE method, as the following formula shows:
min v g * ( P ) ln f ( P ; v ) d P
Combining Equations (18) and (20), v * can be calculated by Equation (22):
v ^ * = arg max v 1 N i = 1 N S { F ( P i ) γ } ln f ( P i ; v )
Moreover, through the smoothing technique, the parameter estimation form v ^ of v * can be expressed as follows:
v ^ k = α v ~ k + ( 1 α ) v ^ k 1
where “~” represents the parameter of elite samples. The top-ranking dispatching solutions after evaluation are chosen as the elite samples. “^” represents the parameter of total samples. α is the smoothing factor (typically between 0.7 and 1); and v ˜ k is the value v of the elite sample after the kth smoothing. Therefore, it is easier to approach the optimal solution with the correction of v ^ through the elite sample.
In short, when using the CE optimization method, first predefine the parameter v and generate the candidate solution set according to the probability density function f ( ; v ) , and then update the value of v through the elite sample. Thus, the searching direction continues to approximate the optimal solution in the iteration.

4.4. Procedures of CEHMC Method for UC Problem

The main steps of the CEHMC method are in detail as follows, and illustrated in Figure 3.
Step 1: Initialization. Set the initial parameters, construct the genetic membrane structure, and generate the initial binary object Io in each membrane (i.e., generate the start-stop states).
Step 2: Evolution. Execute the evolution rules for all objects I in each membrane.
Step 3: Evaluation. For all start-stop states, solve the economic dispatch part with the CEHMC method. The numerical result F(P) is taken as the evaluation indicator to replace the old objects. In this dispatch part, the evolution object P is the generation of units, and the main steps are as follows:
3.1 
Initialization. Set the initial parameters, construct biomimetic membrane structure, and create Nco initial communication objects Pco in the outermost membrane and send them to the first basic membrane.
3.2 
Computation in the basic membrane. Create No initial objects Po in the current membrane, and execute the basic evolution rules in order for No optimal objects to be selected in this membrane. Then, select Nco optimal objects Pbest as the communication objects Pco and Ns suboptimal objects as the reservation objects. Finally, remove the remaining objects, and send the communication objects Pco into the quasi-Golgi membrane.
3.3 
Computation in the quasi-Golgi membrane. Update the target indicator vector, and check whether the activate condition of quasi-Golgi is satisfied. If not, the communication objects Pco are sent to the next basic membrane directly, and return to step 3.2; if satisfied, execute evolution rules in the quasi-Golgi for communication objects Pco, and then select the new communication objects Pco.
3.4 
Check whether the current computation cycle is completed. If not, send the communication objects Pco into the next basic membrane and return to step 3.2; if completed, start the CE optimization steps:
  • Initialization. Calculate the mean and standard deviation for all communication objects Pco in the current cycle as the initial sample.
  • Sampling. Take generation samples P 1 , P 2 , , P N based on the P ~ N ( μ ^ k , σ ^ k 2 ) .
  • Evaluation. Select the elite sample set, and then calculate its mean and standard deviation.
  • Update parameters. Update the mean μ ˜ k and standard deviation σ ˜ k according to Equations (24) and (25):
    μ ^ k = α μ ~ k + ( 1 α ) μ ^ k 1 ,
    σ ^ k = β k σ ~ k + ( 1 β k ) σ ^ k 1
    The standard deviation is usually updated by dynamic smoothing, that is, β k = β 0 β 0 ( 1 1 k ) r , where β 0 is the smoothing factor (typically between 0.8 and 0.99), k is the iteration number, and r is an integer (typically between 5 and 10).
  • Judgment of the termination condition. If the iteration was over, output the optimal object Pbest; if not, return to Step 3.2.
3.5 
Judgment of termination conditions. Check whether all computation cycles are completed. If not, send the communication objects Pco to the first basic membrane and return to step 3.2; if completed, output the function value of the optimal object Pbest.
Step 4: Communication. Select Ne optimal objects Ibest, and send them to the adjacent (sub-outer) membrane. At the same time, the outermost membrane is dissolved, and thus the previous sub-outer membrane becomes the new outermost membrane.
Step 5: Judgment of termination conditions. If all membranes are dissolved, output the best result Ibest of the UC problem; if not, return to Step 3.2.

5. Case Study

The 10-unit 24 h standard thermal power test system was taken as the test example to verify the effectiveness of the proposed method. The unit characteristics and load demand are detailed in the literature [28], and the ramp rate limits of units 1, 3, and 4 are all as follows [29]: Piup = Pidown = 40 MW/h, Pistart = Pishut = 2Pi. The initial parameters of the proposed method are list in Table 1. The quasi-Golgi activate condition is that the multiplication of current computation cycle and current basic membrane can be divided by 3 after the second computation cycle.
In order to analyze the performance of the proposed method, 20 complete independent simulations are conducted on each system. Two cases of the UC problem are simulated: (1) UC problem without ramp constraints; (2) UC problem with ramp constraints. The characteristics, including the number of continues and integer variables and number of constraints for these two cases are list in Table 2, which contribute to depict the magnitude and complexity of the investigated UC problem.
In the following sections, the corresponding simulation results for these two cases are discussed and compared. Moreover, the simulation of the proposed method is programmed by MATLAB 2014a (The MathWorks, Inc, Natick, MA, USA) on the PC with Intel Core i5-3470 CPU, 3.20 GHz, 4 GB Ram (Intel, Santa Clara, CA, USA).

5.1. Simulation Results of UC Problem without Ramp Constraints

The generation of a 10–100-unit 24 h system is shown in Figure 4, where the area of each unit is the generation value.
Comparison of the proposed CEHMC with six bioinspired optimization methods (EPSO, MRCGA, GA, MA, EP, HMC) and three mathematical optimization methods (SOCP, IPL, C&B) is shown in Table 3. To be more illustrative, the normalized comparison results of different algorithms are plotted in Figure 5. For the bioinspired optimization methods, the best and mean solutions of them are plotted as Figure 5a,b, respectively.
From the best solutions, it can be seen that CEHMC has the best performance in all the seven intelligent methods. Compared with the three mathematical optimization methods, the best solutions of CEHMC are lightly larger than IPL (with the smallest objective cost) when unit size is 60 and 80. In other cases of unit size, CEHMC gets the best results. That means, the proposed CEHMC has stable optimization process when unit size increases. Comparing the results of CEHMC and HMC, the solution of CEHMC is much better than HMC. It is shown that the hybrition of cross-entropy and membrane computing theory can significantly improve the optima searching ability and optimization efficiency.

5.2. Simulation Results of UC Problem with Ramp Constraints

The generation of a 10-unit 24 h system is shown in Figure 6, where the area of each unit is the generation value. Comparison of the proposed CEHMC with two intelligent optimization methods (PSO, HMC) and six mathematical optimization methods (MISOCP, OO, BB, HCMIP, SHCMIP, MILP) is shown in Table 4. To be more illustrative, the normalized comparison results of different algorithms are also plotted in Figure 7. For PSO, HMC and CEHMC, the best solutions of them are plotted as Figure 7a,b, respectively.
It is found that, with the unit size larger than 40, CEHMC has the best performance in all the nine optimization methods. In other cases of unit size, the best result of CEHMC is lightly larger than the optima of all the nine optimization algorithms. The introduction of cross-entropy theory greatly improves the optima searching ability, which makes the results of CEHMC generally much better than HMC.

5.3. Influences of the Membrane Number to Convergence

For ease of explanation, the case of 10-unit without ramp constraints is taken as an example to analyze the influences of membrane number to approximation of optimal solution for the proposed CEHMC method.

5.3.1. GAPS

The effect of membrane number in GAPS on the objective value is simulated and illustrated in Figure 8. It is found that No and Nco may have little influence on the convergence with all those curves have similar shapes. Considering the computation cost and matching with other parameters, No is set as 10, and Nco is 10.

5.3.2. BMC

BMC is for the dynamic dispatching problem with continuous variables. The influence of membrane number to convergence for BMC is based on a certain unit start-stop plan. The number of membranes is the product of Nc′ and Nb′, which can also represent the iteration number. Figure 9 shows the penalty value curves changing with membrane number of BMC. It is found that when the abscissa is 100, the curve is almost gentle. Therefore, Nc′ is set as 10, and Nb′ is 10.

5.4. Analysis of Calculation Efficiency

Figure 10 shows how Nco′ influences the approaching to the optimal solution. It is found that Nco′ have effect on convergence. When Nco′ is set 6, the convergence rate is fastest. when Nco′ is set 2, the convergence rate is the slowest. After comprehensive consideration, the value of Nco′ is set 4.
In order to verify the calculation efficiency of the proposed method, the calculation efficiency for the proposed CEHMC method and the other methods (in Table 4) are analysed and compared. The unit ramp rate constraints are considered. The computation time growth rate curves with the unit size increasing from 10 to 100 are illustrated in Figure 11. The computation time growth rate is based on the time consumed for the case of 10-unit.
It can be seen that only SHCMIP has a little better calculation efficiency than the proposed CEHMC method as the unit size increases. The time growth rate of CEHMC is significantly lower than the other four methods (MISOCP, BB, HCMIP, MILP) for large-scale UC problems.

6. Conclusions

In this paper, a cross-entropy-based hybrid membrane computing method is proposed to solve the UC problem which is inspired by living cells and their organization in tissues and other higher order structures. In the proposed method, the genetic algorithm-based P system is applied for the unit start-stop plan with embedded generic rules, which can transmit the outer optima into the inner membranes. The biomimetic membrane computing method combined with the cross-entropy is proposed to solve the dynamic economic dispatch problem with strengthened searching ability, which is inspired by the important role of Golgi apparatus in living cells. The 10–100 unit systems for 24 h day-ahead dispatching simulation results showed that the UC problem could be solved by the proposed method with good efficiency and stability. In future research, network constraints will be considered, and the method will be further improved to optimize simulation results.

Author Contributions

Conceptualization, M.X.; methodology, Y.D.; software, Y.D. and P.C.; validation, W.W.; writing—original draft preparation, Y.D.; writing—review and editing, M.X.; supervision, M.X. and M.L.; funding acquisition, M.X.

Funding

The financial supports for this research, from Guangdong Natural Science Foundation Free Application Project (2018A0303130134), are gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pandžić, H.; Dvorkin, Y.; Qiu, T.; Wang, Y.; Kirchen, D. Toward cost-efficient and reliable unit commitment under uncertainty. IEEE Trans. Power Syst. 2016, 31, 970–982. [Google Scholar] [CrossRef]
  2. Mohammad, J.F.; Mitch, C.; Shabbir, A.; Ahmed, S.; Grijalva, S. Large-scale decentralized unit commitment. Int. J. Electr. Power Energy Syst. 2015, 73, 97–106. [Google Scholar]
  3. Bai, Y.; Zhong, H.; Xia, Q.; Kang, C. A decomposition method for network-constrained unit commitment with AC power flow constraints. Energy 2015, 88, 595–603. [Google Scholar] [CrossRef]
  4. Quan, R.; Jian, J.; Yang, L. An improved priority list and neighborhood S. Arch method for unit commitment. Int. J. Electr. Power Energ. Syst. 2015, 67, 278–285. [Google Scholar] [CrossRef]
  5. Yuan, W.; Zhai, Q. Power-based transmission constrained unit commitment formulation with energy-based reserve. IET Gener. Trans. Distrib. 2017, 11, 409–418. [Google Scholar] [CrossRef]
  6. Zheng, H.; Jian, J.; Yang, L.; Quan, R. A deterministic method for the unit commitment problem in power systems. Comput. Oper. Res. 2016, 66, 241–247. [Google Scholar] [CrossRef]
  7. Nasri, A.; Kazempour, S.J.; Conejo, A.J.; Ghandhari, M. Network-constrained AC unit commitment under uncertainty: A benders decomposition approach. IEEE Trans. Power Syst. 2015, 31, 412–422. [Google Scholar] [CrossRef]
  8. GAMS Development Corporation. GAMS, the Solvers’ Manual, 2015. Available online: http://www.gams.com/solvers (accessed on 1 December 2018).
  9. Rouhi, F.; Effatnejad, R. Unit commitment in power system t by combination of dynamic programming (DP), genetic algorithm (GA) and particle swarm optimization (PSO). Indian J. Sci. Technol. 2015, 8, 134. [Google Scholar] [CrossRef]
  10. Magnússon, S.; Weeraddana, P.C.; Rabbat, M.G.; Fischione, C. On the convergence of alternating direction Lagrangian methods for nonconvex structured optimization problems. IEEE Trans. Control Netw. Syst. 2016, 3, 296–309. [Google Scholar] [CrossRef]
  11. Sun, L.; Zhang, Y.; Jiang, C. A matrix real-coded genetic algorithm to the unit commitment problem. Electr. Power Syst. Res. 2006, 76, 716–728. [Google Scholar] [CrossRef]
  12. Shukla, A.; Singh, S.N. Multi-objective unit commitment using search space-based crazy particle swarm optimisation and normal boundary intersection technique. IET Gener. Trans. Distrib. 2016, 10, 1222–1231. [Google Scholar] [CrossRef]
  13. Valenzuela, J.; Smith, A.E. A seeded memetic algorithm for large unit commitment problems. J. Heuristics 2002, 8, 173–195. [Google Scholar] [CrossRef]
  14. Wang, H.; Peng, H.; Shao, J.; Wang, T. A thresholding method based on P systems for image segmentation. ICIC Express Lett. 2012, 6, 221–227. [Google Scholar]
  15. Wang, X.; Zhang, G.; Zhao, J.; Rong, H.; Ipate, F.; Lufticaru, R. A modified membrane-inspired algorithm based on particle swarm optimization for mobile robot path planning. Int. J. Comput. Comm. Control 2015, 10, 732–745. [Google Scholar] [CrossRef]
  16. Xiao, J.H.; Zhang, X.Y.; Xu, J. A membrane evolutionary algorithm for DNA sequence design in DNA computing. Chin. Sci. Bull. 2012, 57, 698–706. [Google Scholar] [CrossRef]
  17. Zhao, J.; Wang, N. A bio-inspired algorithm based on membrane computing and its application to gasoline blending scheduling. Comput. Chem. Eng. 2011, 35, 272–283. [Google Scholar] [CrossRef]
  18. Zhang, G.X.; Cheng, J.X.; Gheorghe, M. A membrane-inspired approximate algorithm for traveling salesman problems. Rom. J. Inf. Sci. Technol. 2011, 14, 3–19. [Google Scholar]
  19. Leporati, A.; Pagani, D. A membrane algorithm for the min storage problem. In International Workshop on Membrane Computing; Springer: Berlin/Heidelberg, Germany, 2006; pp. 443–462. [Google Scholar]
  20. Pereira, S.; Ferreira, P.; Vaz, A.I.F. A simplified optimization model to short-term electricity planning. Energy 2015, 9, 2126–2135. [Google Scholar] [CrossRef]
  21. Păun, G. Computing with membranes. J. Comput. Syst. Sci. 2000, 61, 108–143. [Google Scholar] [CrossRef]
  22. Păun, A. On P Systems with Active Membranes. Unconventional Models of Computation; UMC’2K; Springer: London, UK, 2001; pp. 187–201. [Google Scholar]
  23. Rubinstein, R. The cross-entropy method for combinatorial and continuous optimization. Methodol. Comput. Appl. Probab. 1999, 1, 127–190. [Google Scholar] [CrossRef]
  24. Kroese, D.P.; Porotsky, S.; Rubinstein, R.Y. The cross-entropy method for continuous multi-extremal optimization. Methodol. Comput. Appl. Probab. 2006, 8, 383–407. [Google Scholar] [CrossRef]
  25. Quinn, A.; Kárný, M.; Guy, T.V. Fully probabilistic design of hierarchical Bayesian models. Inf. Sci. 2016, 369, 532–547. [Google Scholar] [CrossRef]
  26. Yang, Z.; Gao, K.; Fan, K.; Lai, Y. Sensational headline identification by normalized cross entropy-based metric. Comput. J. 2015, 58, 644–655. [Google Scholar] [CrossRef]
  27. Senjyu, T.; Shimabukuro, K.; Uezato, K.; Funabashi, T. A fast technique for unit commitment problem by extended priority list. IEEE Trans. Power Syst. 2003, 18, 882–888. [Google Scholar] [CrossRef]
  28. Guo, S.; Guan, X.; Zhai, Q. The necessary and sufficient conditions for determining feasible solutions to unit commitment problems with ramping constraints. IEEE Power Eng. Soc. Gen. Meet. 2005, 1, 344–349. [Google Scholar]
  29. Yuan, X.; Tian, H.; Zhang, S.; Ji, B.; Huo, Y. Second-order cone programming for solving unit commitment strategy of thermal generators. Energy Convers. Manag. 2013, 76, 20–25. [Google Scholar] [CrossRef]
  30. Roy, P.K. Solution of unit commitment problem using gravitational search algorithm. Int. J. Electric. Power Energy Syst. 2013, 53, 85–94. [Google Scholar] [CrossRef]
  31. Juste, K.A.; Kita, H.; Tanaka, E.; Hasegava, J. An evolutionary programming solution to the unit commitment problem. IEEE Trans. Power Syst. 1999, 14, 1452–1459. [Google Scholar] [CrossRef]
  32. Quan, R.; Hua, W.; Jian, J. Solution of large scale unit commitment by second-order cone programming. Proc. CSEE 2010, 30, 101–107. [Google Scholar]
  33. Xie, M.; Zhu, Y.; Wu, Y.; Yan, Y.; Liu, M. Application of ordinal optimization theory to solve large-scale unit commitment problem. Control Theory Appl. 2016, 33, 542–551. [Google Scholar]
  34. Yang, L.; Jian, J.; Han, D.; Zheng, H. A hyper-cube cone relaxation model and solution for unit commitment. Trans. China Electrotech. Soc. 2013, 28, 252–261. [Google Scholar]
  35. Yang, L.; Jian, J.; Zheng, H.; Dan, D. A sub hyper-cube tight mixed integer programming extended cutting plane method for unit commitment. Proc. CSEE 2013, 33, 99–108. [Google Scholar]
  36. Liu, S.; Jian, J. Research on unit commitment considering emission trading. Power Syst. Technol. 2013, 3, 3558–3563. [Google Scholar]
Figure 1. Structure of genetic membrane.
Figure 1. Structure of genetic membrane.
Energies 12 00486 g001
Figure 2. Structure of biomimetic membrane.
Figure 2. Structure of biomimetic membrane.
Energies 12 00486 g002
Figure 3. Flowchart of cross-entropy-based hybrid membrane computing (CEHMC)-based unit commitment (UC) problem.
Figure 3. Flowchart of cross-entropy-based hybrid membrane computing (CEHMC)-based unit commitment (UC) problem.
Energies 12 00486 g003
Figure 4. Generation of the 10 units (not considering unit ramp rate constraint).
Figure 4. Generation of the 10 units (not considering unit ramp rate constraint).
Energies 12 00486 g004
Figure 5. Normalized comparison results of different algorithms (not considering unit ramp rate constraint).
Figure 5. Normalized comparison results of different algorithms (not considering unit ramp rate constraint).
Energies 12 00486 g005
Figure 6. Generation of 10 units (considering unit ramp rate constraint).
Figure 6. Generation of 10 units (considering unit ramp rate constraint).
Energies 12 00486 g006
Figure 7. Normalized comparison results of different algorithms (considering unit ramp rate constraint).
Figure 7. Normalized comparison results of different algorithms (considering unit ramp rate constraint).
Energies 12 00486 g007
Figure 8. Influence of membrane number to convergence (GAPS).
Figure 8. Influence of membrane number to convergence (GAPS).
Energies 12 00486 g008
Figure 9. Influence of membrane number to convergence (BMC).
Figure 9. Influence of membrane number to convergence (BMC).
Energies 12 00486 g009
Figure 10. The sensitivity of Nco′ to convergence (BMC).
Figure 10. The sensitivity of Nco′ to convergence (BMC).
Energies 12 00486 g010
Figure 11. Computation time growth rate curves (considering unit ramp rate constraint).
Figure 11. Computation time growth rate curves (considering unit ramp rate constraint).
Energies 12 00486 g011
Table 1. Initial parameters of the proposed method.
Table 1. Initial parameters of the proposed method.
Unit1020406080100
GAPSN202040506060
No101620203030
Nco224466
Pc0.90.90.90.90.90.9
Pm0.50.50.50.50.50.5
HMCNc102030304050
Nb102020304050
No101010101212
Nco444466
Pc0.950.950.950.950.950.95
Pm0.50.50.50.50.50.5
Pt0.90.90.90.90.90.9
Table 2. Computation characteristics of the UC problem.
Table 2. Computation characteristics of the UC problem.
UnitVariablesConstraints
IntegerContinuousWithout RampWith Ramp
10240240748886
2048048014481724
4096096028483400
601440144042485076
801920192056486752
1002400240070488428
Table 3. Comparison of different algorithms (not considering unit ramp rate constraint), unit: $.
Table 3. Comparison of different algorithms (not considering unit ramp rate constraint), unit: $.
MethodsUnit
1020406080100
SOCP [4]564,5311,124,7132,244,3693,363,7584,484,3575,603,728
IPL [7]563,9771,123,7952,243,5463,360,7644,481,4115,600,108
C&B [30]563,9381,123,7832,243,6873,363,5934,484,4975,603,976
EPSO [31]best563,9381,123,2322,243,4073,365,4804,488,6015,612,742
mean563,9691,124,1272,246,8003,373,8594,501,2545,620,785
worst564,2061,125,8152,250,3643,381,9474,510,9845,633,447
MRCGA [11]best564,2441,125,0352,246,6223,367,3664,489,9645,610,031
mean564,4671,126,1992,249,6093,371,0364,497,3465,616,957
worst565,7561,128,3262,252,0763,375,8154,505,5115,623,248
GA [13]best565,8661,128,8762,252,9093,377,3934,507,6925,626,362
mean567,3291,130,1062,262,5853,394,0444,525,2045,669,362
worst571,3361,131,5652,269,2823,401,8474,552,9825,690,086
MA [13]best565,8271,127,2542,252,9373,388,6764,501,4495,640,543
mean566,4531,128,8242,262,4773,394,8304,527,7795,665,803
worst566,8611,130,9162,270,3163,408,2754,545,3055,698,039
EP [32]best564,5511,125,4942,249,0933,371,6114,498,4795,623,885
mean565,3521,127,2572,252,6123,376,2554,505,5365,633,800
worst566,2311,129,7932,256,0853,381,0124,512,7395,639,148
HMCbest564,5411,127,5942,250,3283,388,0564,522,4915,654,987
mean564,7161,128,1882,251,5043,390,2374,525,3425,659,001
worst564,9491,128,9532,260,9893,392,4454,527,7625,664,080
CEHMCbest563,9301,123,2062,243,3143,360,7794,479,7205,600,004
mean564,0271,123,3092,249,3883,369,9564,480,1225,602,334
worst564,7961,126,7122,260,6843,391,6984,497,3915,609,585
Note: SOCP: second-order cone programming; IPL: improved priority list; C&B: cut and branch; EPSO: elite particle swarm optimization; MRCGA: matrix real-coded genetic algorithm; EP: evolutionary programming; HMC: hybrid membrane computing.
Table 4. Comparison of different algorithms (considering unit ramp rate constraint), unit: $.
Table 4. Comparison of different algorithms (considering unit ramp rate constraint), unit: $.
MethodsUnit
1020406080100
MISOCP [32]565,7771,130,6472,259,2033,382,4704,511,8135,638,456
OO [33]569,7511,139,5042,261,9003,401,8504,570,8085,697,510
BB [33]568,7101,136,6502,260,2143,383,4894,531,8175,658,458
HCMIP [34]566,0841,129,2412,257,2693,379,8524,508,6895,633,984
SHCMIP [35]565,3971,127,4372,251,6173,376,8214,501,4205,625,531
MILP [36]566,1881,127,2182,252,8103,375,9674,501,5325,623,814
PSObest571,7661,141,4302,285,0743,436,2054,590,0275,730,530
mean572,2161,142,6042,307,2583,439,6094,592,2375,732,596
worst572,6231,144,2252,328,4323,441,8074,593,5535,733,525
HMCbest568,6391,134,8352,269,5233,402,5494,542,5995,679,389
mean569,3961,136,1392,271,2923,404,6124,545,6705,683,837
worst570,5131,137,8042,273,0743,406,0024,548,4065,687,362
CEHMCbest565,3981,127,2172,251,6203,375,9604,501,3005,623,716
mean565,4081,127,8822,252,0303,376,3094,501,9765,623,975
worst5,654,8121,130,8612,256,9333,380,6294,503,0355,624,019
Note: MISOCP: mixed integer second-order cone programming; OO: ordinal optimization method; BB: branch and bound method; HCMIP: hyper-cube mixed integer programming; SHCMIP: sub hyper-cube mixed integer programming; MILP: mixed integer linear programming.

Share and Cite

MDPI and ACS Style

Xie, M.; Du, Y.; Cheng, P.; Wei, W.; Liu, M. A Cross-Entropy-Based Hybrid Membrane Computing Method for Power System Unit Commitment Problems. Energies 2019, 12, 486. https://doi.org/10.3390/en12030486

AMA Style

Xie M, Du Y, Cheng P, Wei W, Liu M. A Cross-Entropy-Based Hybrid Membrane Computing Method for Power System Unit Commitment Problems. Energies. 2019; 12(3):486. https://doi.org/10.3390/en12030486

Chicago/Turabian Style

Xie, Min, Yuxin Du, Peijun Cheng, Wei Wei, and Mingbo Liu. 2019. "A Cross-Entropy-Based Hybrid Membrane Computing Method for Power System Unit Commitment Problems" Energies 12, no. 3: 486. https://doi.org/10.3390/en12030486

APA Style

Xie, M., Du, Y., Cheng, P., Wei, W., & Liu, M. (2019). A Cross-Entropy-Based Hybrid Membrane Computing Method for Power System Unit Commitment Problems. Energies, 12(3), 486. https://doi.org/10.3390/en12030486

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

Article Metrics

Back to TopTop