• Views 786
• Citations 2
• ePub 44
• PDF 332
`Journal of Applied MathematicsVolume 2014, Article ID 752537, 8 pageshttp://dx.doi.org/10.1155/2014/752537`
Research Article

## Constructing Core Backbone Grid Based on the Index System of Power Grid Survivability and BBO Algorithm

School of Electric Engineering, Wuhan University, Wuhan 430072, China

Received 18 January 2014; Revised 4 May 2014; Accepted 19 May 2014; Published 5 June 2014

Copyright © 2014 Feifei Dong et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Constructing core backbone grid is conducive to strengthen the construction of grid structure and improve the ability of withstanding natural disasters. The survivability index of power grid is made up of four indices, namely, the resistibility, recoverability, security, and connectivity. Based on survivability, a method of constructing core backbone grid with the optimal criteria of the minimal total line length and the largest integrated survivability index is put forward. The biogeography-based optimization (BBO) algorithm is used to search for the core backbone grid. Compared with the traditional algorithms, BBO algorithm shows advantages in fast speed of convergence and high convergence precision. Moreover, the searching results for three kinds of objective functions by BBO algorithm verify the effectiveness of the proposed model of constructing core backbone grid.

#### 1. Introduction

For the past few years, power grid in China is constantly damaged by extreme natural disasters due to the failure of past standards of power facilities unable to resist the increasing natural disasters [1, 2]. Therefore, it is necessary to design the standards of disaster resistance levels differentially, according to the various lines’ geographical landforms and climate conditions where power transmission lines are crossed. The goal of differential planning design is to confirm the core backbone grid, which is made up of important lines that can guarantee continuous power supply of the important load when the major natural disasters attack [3, 4]. Constructing core backbone grid is significant to improve the stability of the power grid’s structure, reduce the secondary investment of repairing, and rebuild the harm of power grid caused by natural disasters, as well as guarantee power grid’s safe and reliable operation under severe natural disasters [5, 6].

The concept of survivability is firstly put forward by Hollway and Neumann  in 1993. Survivability of system refers to the ability that the system can complete its critical services in a timely manner and restore its basic services as soon as possible when it is subjected to an attack, failure, or a sudden incident . Survivability has been widely concerned in the complex network and information system as a new research area  [9, 10], but its application in the related fields of power system is relatively less [11, 12]. Considering survivability is helpful to keep the system alive with supportive network structure. In that case, the constructed core backbone grid has a strong resistance, restorative ability, and connectivity. A key method of line identification based on network survivability evaluation was proposed in . A searching model and a searching method of backbone grid were also presented, which has certain enlightening significance. But the survivability index of this method is relatively simple, and the search algorithm is easy to fall into local optimal solution. The survivability index system that can systematically reflect information system is proposed in . The formalized description and mathematical model were given as well. It has done pioneered work in establishing survivability index of power system.

The search of core backbone grid involves nonlinear, multivariate, and discontinuous optimization problem, which depends on artificial intelligence algorithms [15, 16]. As a new kind of artificial intelligence algorithm, BBO algorithm based on species migration patterns has achieved good results in parameter identification , fault diagnosis , and image classification , as well as transmission network planning  and other fields [21, 22]. It is shown that the algorithm has advantages of less set parameters, light computational work, fast convergence speed, and good stability. Therefore, BBO algorithm is adopted to search the optimal solution of core backbone grid in this paper.

A new method of constructing core backbone grid considering survivability is put forward in this paper. The indices of survivability are built from factors of resistibility, recoverability, security, and connectivity. The largest integrated survivability index and the minimum total length of the backbone grid’s lines are regarded as the optimal solution of objective function, given network connectivity and power grid’s safe operation as constraint conditions. BBO algorithm provided with strong search ability is used to search backbone grid. Empirical simulations show that the proposed method is accurate and effective, with additional advantages of fast convergence speed and high convergence precision over the methods such as particle swarm optimization (PSO), binary ant colony algorithm (BACA), and genetic algorithm (GA).

#### 2. Basic Method of Constructing Core Backbone Grid

The key of differential planning design is to determine the selecting principles and construction standards of affordable power supply, load, and network frame. Then, the quantitative evaluation method is adopted to identify the key components and build core backbone grid, increasing the resistant standards of disasters of element differentially. In this way, the system guarantees the delivery capacity of important load and power supply as well as the ability of the interval exchange capacity, promising the system recovered and reconstructed on the basis of the core backbone grid. The core method lies in that it is the viability and self-healing ability instead of the emergency control and scheduling dispatch during disaster and the repair and reconstruction after disaster that alleviate the damage of major natural disasters to the power grid. Therefore, the core backbone grid needs to satisfy the following conditions: (a) to guarantee the continuous power supply of the important load and the reliable transmission of the important power; (b) to make the rack satisfy the constraints of power grid’s safe operation as well as network’s topology connectivity; (c) to meet the total length of the subrack’s line to be shortest from the perspective of economy; (d) to meet the survivability of the subrack to be stronger from the perspective of the optimal network topology configuration.

Survivability is one of the most important constraints of the core backbone grid, and it has great influence in resistance ability, resilience, and the stability of the overall network topology of the rack against natural disaster. Therefore, establishing an index system of survivability effectively and then determining the comprehensive indicators of network survivability are the key steps to construct core backbone grid.

#### 3. The Index System of Power Grid Survivability

Survivability of power grid is defined as the ability of ensuring continuous power supply to important loads that relied on the high design standards of the core backbone grid, as well as restoring power supply for other loads gradually based on the network frame when major natural disasters attack. The index system of survivability is built from these four aspects as follows: resistibility, recoverability, security, and connectivity, by which the integrated index of survivability is finally obtained.

##### 3.1. The Index of Resistibility

Resistibility of power grid reflects the resistance where the system provides power supply for the important loads with major natural disasters. The three indicators, namely, preserving rate of lines, preserving rate of nodes, and preserving rate of loads, are introduced to evaluate resistibility.

Preserving rate of lines is as follows: where is the number of original rack’s lines and is the number of failure lines of the remaining network frame after natural disasters compared to the original rack.

Preserving rate of nodes is as follows: where is the number of original rack’s nodes and is the number of failure nodes of the remaining network frame after natural disasters compared to the original rack.

Preserving rate of loads is as follows: where is the number of original rack’s lodes and is the number of failure loads of the remaining network frame after natural disasters compared to the original rack.

##### 3.2. The Index of Recoverability

Recoverability of power grid reflects whether the power grid can recover after suffering from natural disasters, and to what extent it can recover. The generator standby indicator and load recovery degree index are introduced to evaluate recoverability.

The generator standby indicator is as follows: where     is the actual output of the generator     in the backbone grid, which is the output power of the generator     obtained by power flow calculation of the core backbone grid. is the largest capacity of the generator in the backbone grid. is the total number of generators.

The load recovery degree index is as follows: where is the active load of node  in the backbone network frame scheme, which is the biggest active load that node    can bear. is the actual active load of node in the backbone network frame scheme, which is the active load that node needs to afford according to the scheme of the backbone grid. They need to meet the safe operation conditions and ensure the rack’s load to be largest.

##### 3.3. The Index of Security

The less important lines of core backbone grid after natural disasters are checked by criteria; the safety margin of core backbone grid is analyzed by the index of security called Sft, mainly including the index of bus voltage fluctuation and index of branch power fluctuation.

The difference between the current bus voltage and upper or lower limits of bus voltage is regarded as the index of bus voltage fluctuation, namely, , which is to measure the current bus voltage with following expression: where , and are the amplitude of current voltage, nominal amplitude, the upper limit, and the lower limit of Bus , respectively.

The bus voltage fluctuations of all nodes in the core backbone grid are taken advantages as the average bus voltage fluctuation, namely, , which can be expressed as

The difference between the current branch power and upper or lower limit of branch power is regarded as the index of branch power fluctuation, namely, , which is used to analyze the current branch power flow. Its expression is as follows: where , and are the amplitude of current power, nominal amplitude, the upper limit, and the lower limit of Branch , respectively.

The average branch power fluctuations, , are defined as the mean of the sum of all lines of the branch power fluctuations in the core backbone grid with the expression of

##### 3.4. The Index of Connectivity

The relative tightness degree and the relative condensation degree of the grid are introduced to evaluate connectivity.

(1) The Relative Tightness Degree of the Grid. If the node has adjacent nodes, there are at most lines among these nodes. Assuming that the actual linking lines are , in that case, the clustering coefficient of the node is defined as follows:

The tightness degree of the grid is obtained by the weighted average of all nodes of   in backbone grid, which is a parameter that shows the connected degree of neighboring nodes, and is expressed as

If the tightness degree of the original grid is , the relative tightness degree of the backbone grid is as follows:

(2) The Relative Condensation Degree of the Grid. The relative condensation degree of the grid is the reciprocal of the product of the number of nodes and the weighted average of the shortest path shown as

The traditional shortest path is the minimum number of edges between two nodes. Considering the actual conditions of power system and assigning the length of the transmission line as the lines’ weight coefficients, the weighted average of the shortest path of backbone grid is as follows: where is the weighted number of edges that the shortest path passes through and is the total number of nodes in the grid.

By substituting formula (14) into formula (13), the condensation degree of the backbone grid is

Given the condensation degree of the original grid as , the relative condensation degree of the backbone grid is as follows:

##### 3.5. Integrated Index of Survivability

There are four indicators included in the index system of power grid survivability. Primary level (Level 1) of index is the integrated index of survivability, while the secondary level (Level 2) of indices is resistibility, recoverability, security, and connectivity index, and the tertiary level (Level 3) of indices is made up by refining indicators of each secondary index. The indexes of Level can be calculated by the indexes of Level according to the following principles.

Vector is set as the individual risk index vector of survivability’s Level indexes. is the total number of Level indicators. is the single index of Level indexes. The survivability index of Level is defined as follows: where and are weight coefficients to satisfy the equation of and and are 1 norm and norm of vector , respectively. Consider

Therefore, the integrated index of survivability is expressed as follows: where and are weighted coefficients that meet the need of ; .

#### 4. The Mathematical Model of Constructing the Core Backbone Grid

The core backbone grid includes the important power source nodes; load nodes; key circuits; the nodes connected to the key circuits that should be guaranteed according to the practical situation; the lines and nodes that need to be preserved considering the constrains of network connectivity, tide of constraints, the shortest length of total lines, and the largest integrated index of survivability. Therefore, mathematical model of constructing core backbone grid is as follows:

In the formula, denotes the on or off state of the lines; 1 means on and 0 means off. is the length of line. is the total length of the backbone grid and is that of the original grid. is the normalized value of the backbone grid’s total length compared with the original grid. The value of is in the (0, 1) and compared with . is the function to test the connectivity. refers to the fact that the network is connected, while indicates that the network is disconnected. is the equality constraint for the power flow. is the inequality constraint for the power flow.

#### 5. Search of the Backbone Grid Based on BBO Algorithm

##### 5.1. The Basic Method of BBO Algorithm

Biogeography-based optimization algorithm is a new type of evolutionary algorithm that uses the biogeography theory to solve optimization problems. Its basic idea is as follows. Firstly, a number of relatively independent habitats are constructed as candidate solutions. Secondly, the species share information through migration between habitats and update the information through mutation. Finally, the most fitted habitat is selected, and the optimal solution of the problem is obtained [23, 24]. For the multiobjective optimization problems of the core backbone grid, the appropriate index HIS is adopted to express the objective function of the mathematical model of constructing core backbone grid, that is, . The variable SIV is used to show variables included in each habitat, namely, the state of lines. There are mainly two operations included in the algorithm, migration and mutation.

(1) Migration Operation. The migration operation is for the information exchange with other habitats, in order to search for the solution in the wide area. Since the linear species migration model cannot accurately simulate the complex process of species migration of the practical biological environment, the cosine migration model which is more similar to the nature process is adopted to simulate migration rate and is shown in Figure 1.

The immigration rate and emigration rate of cosine migration model are shown as follows, respectively:

In this migration model, denotes the largest immigration rate, is the biggest emigration rate, is the number of species, and is the maximum number of species. When there are less or more species in the habitat, immigration rate and emigration rate change smoothly. When habitat has a certain number of species, immigration rate and emigration rate change relatively fast.

(2) Mutation Operation. Mutation operation is to increase the diversity of population. The probability of the habitat, , has a number of species. The variation rate is found out through the following formula: where is the biggest mutation rate and it can be adjusted according to the requirements of different users. is the maximum of .

Mutation operation makes the solution of lower HIS improve by variation and makes the higher HIS obtain the opportunity of improving their solutions.

##### 5.2. Search of Core Backbone Grid Based on BBO Algorithm

The multiobjective optimization problem of constructing the core backbone grid is solved by BBO algorithm to obtain the high survivability. Specific steps are as follows.

Step 1. Initiate the original grid’s parameters needed, including the node loads, the output of generators, and the proposed paths. The control parameters of the BBO algorithm are initialized as well. The initial population  that meets the constraint condition is randomly generated.

Step 2. The suitability index of the habitat is calculated and sorted. The individual optimal solution, Xbest, is saved and then judged whether it meets the end condition. If yes, the results are output and then converted to backbone grid scheme and, in that case, the program is over. Otherwise, continue to Step 3.

Step 3. The cosine migration model is established. The habitat’s species number, immigration rate, and emigration rate are calculated.

Step 4. The migration operation is performed to form a new population . The suitability index of the habitat is recalculated, and then the optimal solution Xbest1 is updated.

Step 5. The mutation operation is carried out, and then the optimal solution of the population Xbest2 is updated.

Step 6. Determine whether it meets the maximum number of iterations. If satisfied, the results are output and then converted to core backbone grid scheme. In that case, the program is over. Otherwise, go to Step 2.

#### 6. Case Study

In order to verify the effectiveness of the proposed method, the core backbone grid of IEEE118 node system is constructed based on BBO algorithm, taking the power grid’s survivability into account. The system includes 118 nodes and 179 branches, with 53 generator nodes. Except for the nodes of 17, 22, 23, 57, 58, 84, 102, 108, 109, and 114, other load nodes are all important loads.

The parameters of BBO algorithm are set as follows: the population size is set as 50, the largest number of iterations is 200, the largest mutation rate is 0.01, the biggest immigration rate and emigration rate are both set as 1, and the global mobility is 1. The optimal scheme of core backbone grid considering survivability searched by BBO algorithm is shown in Figure 2. The solid points represent generator nodes. The hollow points denote the load nodes. The reserved lines of the core backbone grid are in the blue solid lines. This scheme of core backbone grid is composed of 71 nodes and 70 lines, including 19 generator nodes and 44 load nodes.

Figure 2: Core backbone grid scheme for IEEE 118-bus power system based on BBO algorithm.

In order to verify the performance of BBO algorithm, PSO algorithm, BACA algorithm, and GA algorithm are used to search the core backbone grid, and their population size and maximum number of iterations are set the same as BBO algorithm. Based on the characteristics of randomness of the algorithm, the statistical results of the four kinds of algorithms are listed in Table 1 with total independence of 50 times. The curves of suitability index with four kinds of algorithms optimal scheme are shown in Figure 3.

Table 1: Comparison of objective function solution by different algorithms.
Figure 3: Contrast of suitability index curve of 4 kinds of algorithm optimal solution.

From Table 1 and Figure 3, the results indicate that the best value and the worst value of BBO search’s objective function are relatively smaller. In that case, the precision of objective function searched by BBO algorithm is higher. Meanwhile, BBO algorithm has got 32 searches of the optimal solution, which shows that this algorithm has a stronger convergence ability compared with other three kinds of algorithms. In addition, it can be concluded from Figure 3 that BBO algorithm has obvious advantages in convergence rate compared with the other three kinds of algorithms.

The optimal core backbone grid schemes of three objective functions searched by BBO algorithm are presented in Table 2.

Table 2: Comparison of the optimal scheme of core backbone grid under different objective functions based on BBO algorithm.

Three kinds of objective functions are set as follows. refers to the fact that the number of lines is the least, where expresses the total number of lines. refers to the shortest length of total lines and the largest integrated survivability index. refers to the minimum length of lines.

From Table 2, under the objective function of , the number of lines of the core backbone grid’s optimal scheme is 67, and the total length of lines is 6.3739. Under the objective function of , there are 70 lines in the lines’ collection and the total length of lines is 5.4490. Meanwhile, 75 lines are in the core backbone grid’s optimal scheme under the objective function with total length of 5.3898. Since in the objective function reflects the content of the objective function and the preserving rate indicators of survivability in the objective function also reflect the content of the objective function , it is reasonable that the searched number of lines and the total length of lines of is between that of and . Thus, the correctness and effectiveness of search algorithm proposed in this paper are verified. In addition, the objective function, , of constructing the core backbone grid takes into account the resistibility, the recoverability, and the configuration of the network frame, compared with the commonly used objective functions, and , which emphasizes economic benefits. This method reflects the characteristic of network frame more comprehensively and gives insights into improving the resistance of power grid against natural disasters and in differentiation design.

#### 7. Conclusions

The index system of survivability is proposed and is made up of resistibility, recoverability, security, and connectivity. Based on the index system of power grid survivability, the model of constructing core backbone grid with the target of the shortest length of total lines and the largest integrated survivability index is built. The results of the case study show that the model can well balance the factors of economy, resistibility, recoverability, and the configuration of the network frame giving insights of differentiated planning with reasonable design schemes.

Compared with the other three kinds of traditional algorithms, BBO algorithm has higher precision and better convergence speed in searching for the optimal solution of core backbone grid.

The system of IEEE118 nodes is studied. The multiple operational results of the four different algorithms and the search results for three kinds of objective functions by BBO algorithm are compared. In that case, the rationality of the proposed model and the superiority of search algorithm are verified.

#### Nomenclature

 Preserving rate of line The number of original rack’s lines The number of failure lines of the remaining grid Preserving rate of node The number of original rack’s nodes The number of failure nodes of the remaining grid Preserving rate of load The number of original rack’s loads The number of failure loads of the remaining grid Generator standby indicator Actual output of generator in backbone grid Largest capacity of generator in backbone grid Load recovery degree index Active load of node in the backbone grid Actual active load of node in the backbone Index of bus voltage fluctuation The amplitude of BUS ’s current voltage Nominal amplitude of Bus The upper limit of Bus The lower limit of Bus Average bus voltage fluctuation Index of branch power fluctuation The amplitude of Branch ’s current power Nominal amplitude of Branch The upper limit of Branch The lower limit of Branch Average branch power fluctuation Clustering coefficient of the node Tightness degree of the backbone grid Tightness degree of the original grid Relative tightness degree of the backbone grid Condensation degree of the backbone grid Condensation degree of the original grid Relative condensation degree of backbone grid The state of line’s input or excision The length of line The total length of the backbone grid The total length of the original grid Judging function of connectivity Immigration rate Emigration rate Probability of the habitat with species Variation rate Biggest mutation rate The maximum of Total number of lines.

#### Conflict of Interests

The authors have declared that no conflict of interests exists.

#### Acknowledgments

This work is funded by the State Grid Corporation of China, Major Projects on Planning and Operation Control of Large Scale Grid (SGCC-MPLG-029-2012), and Natural Science Foundation of China (51207114).

#### References

1. E. Broström, J. Ahlberg, and L. Söder, “Modelling of ice storms and their impact applied to a part of the Swedish transmission network,” in Proceedings of the IEEE Lausanne Power Tech Conference, pp. 1593–1598, Lausanne, Switzerland, July 2007.
2. R. Billinton and G. Singh, “Application of adverse and extreme adverse weather: modelling in transmission and distribution system reliability evaluation,” IEE Proceedings: Generation, Transmission and Distribution, vol. 153, no. 1, pp. 115–120, 2006.
3. R.-A. Hooshmand, R. Hemmati, and M. Parastegari, “Combination of AC transmission expansion planning and reactive power planning in the restructured power system,” Energy Conversion and Management, vol. 55, pp. 26–35, 2012.
4. G. Kjolle, L. Rolfseng, and E. Dahl, “The economic aspect of reliability in distribution system planning,” IEEE Transactions on Power Delivery, vol. 5, no. 2, pp. 1153–1157, 1990.
5. G. Trudel, J.-P. Gingras, and J.-R. Pierre, “Designing a reliable power system: Hydro-Québec's integrated approach,” Proceedings of the IEEE, vol. 93, no. 5, pp. 907–917, 2005.
6. G. Trudel, S. Bernard, and G. Scott, “Hydro-Quebec's defence plan against extreme contingencies,” IEEE Transactions on Power Systems, vol. 14, no. 3, pp. 958–966, 1999.
7. B. A. Hollway and P. G. Neumann, Survivable Computer-Communication Systems: The Problem Working Group Recommendations, US Army Research Laboratory, Washington, DC, USA, 1993.
8. R. J. Ellison, D. A. Fisher, R. C. Linger et al., Survivable Network Systems: An Emerging Discipline, Software Engineering Institute, Carnegie Mellon University, Pittsburgh, Pa, USA, 1997.
9. M. A. Schroeder and K. T. Newport, “Tactical network survivability through connectivity optimization,” in Proceedings of the IEEE Military Communications Conference—Crisis Communications: The Promise and Reality (MILCOM '87), pp. 590–597, Washington, DC, USA, October 1987.
10. A. A. Hagin, “Performability, reliability, and survivability of communication networks: system of methods and models for evaluation,” in Proceedings of the 14th International Conference on Distributed Computing Systems, vol. 11, pp. 562–573, June 1994.
11. A. Dwivedi, X. Yu, and P. Sokolowski, “Identifying vulnerable lines in a power network using complex network theory,” in Proceedings of the IEEE International Symposium on Industrial Electronics (ISIE '09), pp. 18–23, Seoul, Republic of Korea, July 2009.
12. A. Dwivedi, X. Yu, and P. Sokolowski, “Analyzing power network vulnerability with maximum flow based centrality approach,” in Proceedings of the 8th IEEE International Conference on Industrial Informatics (INDIN '10), pp. 336–341, Osaka, Japan, July 2010.
13. W. Yang, T. Bi, S. Huang, A. Xue, and Q. Yang, “An approach for critical lines identification based on the survivability of power grid,” Proceedings of the Chinese Society of Electrical Engineering, vol. 31, no. 7, pp. 29–35, 2011.
14. J. Wang, H. Wang, and G. Zhao, “Index system of quantitative evaluation for information systems survivability,” Computer Engineering, vol. 35, no. 3, pp. 54–56, 2009.
15. A. Soroudi and M. Afrasiab, “Binary PSO-based dynamic multi-objective model for distributed generation planning under uncertainty,” IET Renewable Power Generation, vol. 6, no. 2, pp. 67–78, 2012.
16. R. S. Maciel, M. Rosa, V. Miranda, and A. Padilha-Feltrin, “Multi-objective evolutionary particle swarm optimization in the assessment of the impact of distributed generation,” Electric Power Systems Research, vol. 89, pp. 100–108, 2012.
17. L. Wang and Y. Xu, “An effective hybrid biogeography-based optimization algorithm for parameter estimation of chaotic systems,” Expert Systems with Applications, vol. 38, no. 12, pp. 15103–15109, 2011.
18. M. Ovreiu and D. Simon, “Biogeography-based optimization of neuro-fuzzy system parameters for diagnosis of cardiac disease,” in Proceedings of the 12th Annual Genetic and Evolutionary Computation Conference (GECCO '10), pp. 1235–1242, ACM, Portland, Ore, USA, July 2010.
19. G. Sonakshi, A. Anuja, V. K. Panchal et al., “Extended biogeography based optimization for natural terrain feature classification from satellite remote sensing images,” in Proceedings of the International Conference on Contemporary Computing, pp. 262–269, IC3, Noida, India, 2011.
20. D. Simon, “Biogeography-based optimization,” IEEE Transactions on Evolutionary Computation, vol. 12, no. 6, pp. 702–713, 2008.
21. I. Boussaïd, A. Chatterjee, P. Siarry, and M. Ahmed-Nacer, “Hybridizing biogeography-based optimization with differential evolution for optimal power allocation in wireless sensor networks,” IEEE Transactions on Vehicular Technology, vol. 60, no. 5, pp. 2347–2353, 2011.
22. R. Mukherjee and S. Chakraborty, “Selection of the optimal electrochemical machining process parameters using biogeography-based optimization algorithm,” The International Journal of Advanced Manufacturing Technology, vol. 64, no. 5–8, pp. 781–791, 2013.
23. A. Müller, “On the performance of linear adaptive filters driven by the ergodic chaotic logistic map,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 22, no. 12, Article ID 1250290, 11 pages, 2012.
24. P. K. Roy, S. P. Ghoshal, S. S. Thakur et al., “Optimal reactive power dispatch considering flexible AC transmission system devices using biogeography-based optimization,” Electric Power Components and Systems, vol. 39, no. 8, pp. 733–750, 2011.