Research Article  Open Access
Lei Wang, Ming Cai, Hu Zhang, Fuad Alsaadi, Liu Chen, "Active FaultTolerant Control for Wind Turbine with Simultaneous Actuator and Sensor Faults", Complexity, vol. 2017, Article ID 6164841, 11 pages, 2017. https://doi.org/10.1155/2017/6164841
Active FaultTolerant Control for Wind Turbine with Simultaneous Actuator and Sensor Faults
Abstract
The purpose of this paper is to show a novel faulttolerant tracking control (FTC) strategy with robust fault estimation and compensating for simultaneous actuator sensor faults. Based on the framework of faulttolerant control, developing an FTC design method for wind turbines is a challenge and, thus, they can tolerate simultaneous pitch actuator and pitch sensor faults having bounded first time derivatives. The paper’s key contribution is proposing a descriptor sliding mode method, in which for establishing a novel augmented descriptor system, with which we can estimate the state of system and reconstruct fault by designing descriptor sliding mode observer, the paper introduces an auxiliary descriptor state vector composed by a system state vector, actuator fault vector, and sensor fault vector. By the optimized method of LMI, the conditions for stability that estimated error dynamics are set up to promote the determination of the parameters designed. With this estimation, and designing a faulttolerant controller, the system’s stability can be maintained. The effectiveness of the design strategy is verified by implementing the controller in the National Renewable Energy Laboratory’s 5MW nonlinear, highfidelity wind turbine model (FAST) and simulating it in MATLAB/Simulink.
1. Introduction
With the increasing demand for energy, the importance of new energy development increases with each passing day. Wind energy is an important part of new energy, and actively developing wind energy to improve the energy system structure, ease the energy crisis, and protect the environment is of great significance [1, 2]. For engineering systems, such as wind power systems, their working conditions are poor, and, in the case of longterm operation, a failure of a subsystem or actuator or sensor failure is usually unavoidable. In most cases, we cannot predict the time of system failure. Such a failure may cause the structure and performance of the whole system to change slowly or drastically, which will seriously affect the safety of the wind turbine and its stable power output. Faulttolerant control can guarantee the stable operation of a system under a specified performance index in the event of system component failure, thus opening a new way to improve the reliability of wind turbines in complex and harsh conditions.
In recent years, the application of modelbased faulttolerant control in wind turbines has gradually expanded. Sloth et al. designed an active and passive faulttolerant controller based on the linear variable parameter (LPV) for the pitch fault. The optimization method for the linear matrix inequality was applied to active faulttolerant control, and the bilinear matrix inequality method was applied to passive faulttolerant control [3]. Sami and Patton designed a sliding mode faulttolerant control based on gain adaptive control for a 5 MW wind turbine operating at low wind speeds. A robust observer was used to estimate the state and unknown output (sensor fault and noise), which guaranteed the robustness of the sliding mode to unknown output [4]. Badihi et al. designed two faulttolerant control methods for the model uncertainties and torque actuator faults of wind turbine generator converters. One method was based on fuzzy model reference adaptive control and used fuzzy inference for parameter adaptation without prior information about any faults. The other method involved fuzzy modeling and identification technology to design a method based on the model detection and diagnosis, thus automatically using realtime fault information to correct fault signs and achieving faulttolerant control for the torque actuator [5]. A faulttolerant control strategy based on the unknown input observer array was proposed by Odgaard and Stoustrup. The fault detection and isolation of highspeed and low speed shaft speed sensors can be realized by an unknown input observer [6]. Blesa et al. proposed a method based on the interval observer to detect and isolate the fault of a wind turbine sensor or actuator; according to the fault information obtained, the faulttolerant control method for the corresponding faults can be realized by using a virtual sensor or actuator technology [7].
However, the current faulttolerant control strategies generally aim at only a single fault. Thus when faults appear in both wind turbine actuator and the sensor at the same time, guaranteeing the stable operation of the system in the specified performance index becomes challenging.
To solve the above problems, this paper proposes a faulttolerant control strategy for the pitch system for a wind turbine based on multiple fault reconstruction. The proposed faulttolerant control strategy consists of two parts. The first part adopts the active disturbance rejection control technology to ensure the stable output of the wind turbine power in the case of no fault. In the second part, a descriptor sliding mode observer is designed to realize the continuous estimation of the original system state, pitch actuator failure, and pitch sensor fault, and a faulttolerant controller is designed based on the state estimation and fault estimation information to maintain the stable operation of the system under simultaneous pitch actuator and pitch sensor faults.
The remainder of this paper is organized as follows. Wind turbine modeling is presented in Section 2. A descriptor sliding mode observer is developed in Section 3. An active disturbance rejection controller used for freefaults and a reconfigurable controller used to compensate the pitch system fault effect are developed in Section 4. The proposed FFTC design is then verified in National Renewable Energy Laboratory’s 5MW nonlinear, highfidelity wind turbine model and simulated it in MATLAB/Simulink in Section 5. The conclusions are given in Section 6.
2. Wind Turbine Modeling
A wind turbine system is characterized by nonlinear aerodynamic behavior and the dependence on a stochastic uncontrollable wind force as a driving signal. To conceptualize the system from the analytical and control design requirement standing points, an overall model of the turbines is required. In this section, a generalized nonlinear wind turbine model, its pitch system, and its faults model are presented.
2.1. Nonlinear Wind Turbine Modeling
The model is shown in the form as follows [8] with a generalized nonlinear property:
The state vector , control input , and nonlinear vector are defined aswhere is the speed of low speed shaft; is the generator speed; is the twist angle; is the pitch angle; is the pitch angle control; is the time constant of the pitch actuator; is the generator torque; and are the low speed shaft and generator inertia; is the gear ratio; and and are the drivetrain damping and spring constants, respectively.
The mechanical power captured by the wind turbine is described aswhere is the rotor radius; is the air density; is the wind speed; and is the power conversion coefficient of the wind turbine and is a nonlinear function.
2.2. Pitch System Modeling
The system of hydraulic pitch is modeled as a closedloop transfer function. Essentially, these position servo systems can be modeled quite well via describing a secondorder transfer function as follows [9]:where and are the frequency and damping ratio parameters, respectively; is the pitch angle; and is the pitch angle control.
To facilitate the subsequent controller design, the state space model of the pitch system is as follows:
2.3. Pitch Actuator Fault Modeling
The fault considered for the pitch actuator is hydraulic leakage. The dynamics of the pitch systems will be changed by a drop of oil pressure. The pressure level is modeled as a convex combination of the vertices of the parameters and . Hence, the pitch system can be described according to the socalled fault effectiveness parameter , where corresponds to a faultfree actuator with and relates to a full fault on the actuator with , . Hence, the corresponding damping ratio and damping frequency can be obtained by inverse formula and and can be depicted by the pitch actuator fault [10, 11]:
The corresponding state space model is as follows:where , .
Equation (9) simplifies the design of the observer. However, is the newly defined actuator fault signal and it may not be able to sufficiently reflect the severity of the primordial fault and it is a potential problem. One method to refind the primordial fault signal is to put (10) into use:
2.4. Pitch Sensor Faults Modeling
There are three common faults for the pitch sensor: biased output, fixed output, and no output. The only difference between the front two faults original from the pitch sensor is that detecting a fixed output is necessary [12, 13].
2.4.1. Biased Output on Pitch Sensor
The closedloop pitch system and the pitch angle measurement are affected by a biased pitch sensor measurement. While the bias occurs, the pitch sensor fault model is as follows:where is the pitch angle measurement bias and is the actual pitch angle.
2.4.2. Fixed Output on Pitch Sensor
A fixed output on a pitch sensor is an abrupt fault and can occur at any time leading to the following measurement equation:where is the time when the fault occurs.
2.4.3. No Output on Pitch Sensor
No output on a pitch sensor causes the same changes in the measurement equation as a fixed output. Contrary to a fixed output, the control system is informed when there is no output from the pitch sensor. The fault model is as follows:
3. Descriptor Sliding Mode Observer Design
3.1. Problem Description
Considering the following wind turbine pitch system with four conditions, pitch actuator fault, pitch sensor fault, system uncertainty, and external disturbance [14–17], the corresponding dynamic equations are as follows: where is the system state vector; is the measurement output vector; is the vector of control input; is external and uncertainty disturbance; is the vector of pitch actuator fault; is the vector of pitch sensor fault; and , , , , , and , , are known constant matrices.
Some assumptions are made in this paper for subsequent theoretical analyses and controller design.
Assumption 1. The conditions for the uncertainty and external disturbance are bounded as the following norm, and pitch actuator fault and pitch sensor fault are held:where , , , and are some known constants.
Assumption 2. The pair is stable, and is tested. And a constant exits such that
Assumption 3. The matrices and are full column rank, and the dimension standard of holds.
The purpose of this paper is to get the estimates of pitch system state and fault simultaneously. For the objective, we can define the new augmented variables and matrices as follows:where is a design parameter selected such that Assumption 2 holds.
By the above definitions, we can construct a new descriptor system equivalently as follows:where is the descriptor system state, composed of the initial state, the pitch actuator fault, and the pitch sensor fault, which is imported into equivalently establishing a system for descriptor (18).
3.2. Descriptor Sliding Mode Observer Design
The descriptor observer can be designed by the lemma as the following shows, by the summaries from [18].
Lemma 4. There always exists a matrix designed such that is nonsingular; that is, the matrix exists. Furthermore, it holds that and .
On the basis of Lemma 4, we can design the descriptor sliding mode observer for the descriptor plant (18) as follows:where is a variable in the middle, is the estimator of descriptor system state, , , and are the derivative gain, proportion gain, and sliding mode gain, and is the input of discontinuous sliding mode to get rid of the disturbances of system faults and uncertainties.
We can define the design matrix as follows:where and is a design scalar.
Lemma 5. If Assumption 2 holds, a proportional gain will exist, such that the matrix is Hurwitz, where we can solve via the Lyapunov equation:where satisfies .
The final step is to form the input term of sliding mode to make the input disturbance attenuated. At first, the gain matrix for sliding mode can be designed as
Then, we can define the surface of sliding mode with as follows:where the descriptor error is and the Lyapunov matrix is designed such that the following constraint holds:where is the matrix of design parameter based on matrix .
Remark 6. Constraint (24) is introduced because is unknown in the function (23) of sliding mode; thus, cannot be got. Note that is available, because and can be obtained. Therefore, on the basis of the constraint (24), we can derive it as follows:
It is presented that is available to design the observer. At present, the discontinuous sliding mode input term can be designed in the form as the following shows:where is a small designed value.
3.3. The Descriptor Error and Estimation State Dynamics
First, we can rewrite the observer (19) as follows:
Using Lemma 4, (27) is simplified to
Adding to both sides of the descriptor system (18), it can be rewritten as
Subtracting (29) from (28), the descriptor error dynamics becomes
Then, we can derive the system equation of sliding mode estimation state. First, we can rewrite the observer state Eq. (28) as follows:
which is equivalent to
We can decompose the matrices , , and as follows:
From , it can be easily obtained that
and from (20), we can derive it as
Subsequently, the state Eq. (32) is decomposed to obtain the following estimation state Eq. (36):
Second, an estimation state is constructed based on the integral sliding mode surface as follows:where and are the matrices with parameter designed. is chosen such that is nonsingular, and is designed such that is Hurwitz. From (36) and (37), it yields
Assuming that there exist matrices and such that and , and letting , we can derive the equivalent control scheme as follows:
Substituting (39) into (36), we can obtain the state system equation with the method of sliding mode estimation as the following shows:
Combine the sliding mode estimation state Eq. (40) with the descriptor error system (25) and build a system with the closedloop property as (41) shows:
3.4. The Stability Condition of Observer
Theorem 7. Based on the sliding mode input scheme as (26) presented, the system shown in (41) is stable, if there exist , , and such that the constraint (24) holds and the LMI optimization problem admits a feasible solution [14]: where , , , , , , , , , , .
4. Active FaultTolerant Control Design
The faulttolerant control of a wind turbine in area 3 (above rated wind speed) is mainly studied in this paper. The output power of the generator must be adjusted to the rated power to ensure the output quality and wind turbine safety. There are two general methods for power regulation [19]: With the generator torque fixed and the generator speed adjusted to the rated speed by controlling the pitch angle and the generator torque and pitch angle are adjusted simultaneously. However, regardless of which method is used, the pitch system must be used to regulate the pitch angle. Therefore, the pitch actuator and pitch sensor are very important parts of the power control system. The purpose of this paper to develop a novel active faulttolerant controller to address the case in which the pitch actuator and pitch sensor show faults at the same time, so only a collective pitch method is used. Method will be chosen for the power control in this paper.
An active faulttolerant control for the pitch system (14) comprises a norm controller and a reconfigurable controller, designed aswhere is the norm controller to achieve pitch angle control under the faultfree case and is the reconfigurable controller to compensate the fault effect.
4.1. Norm Controller Design
The active disturbance rejection control (ADRC) technique is used in wind turbine pitch control under the faultfree case in this paper due to its facilitation and robustness. The ADRC controller design procedures are summarized as the following three steps [20, 21]:(1)Firstly, to avoid a large overshoot of the system, the step input signal is transformed into a continuous and smooth signal by designing the tracking differentiator [22].(2)Second, an extended state and disturbance nonlinear observer is designed to estimate and compensate for the unknown time varying nonlinear disturbances in the system online.(3)Finally, the pitch control is realized by a conventional PD controller.
It is assumed that the output of the system is the speed of the low speed shaft and that the input is the pitch angle. From (1) to (3), the dynamic equation of low speed shaft speed is
Its secondorder derivative can be obtained aswhere is the derivative of the wind speed.
Assuming that all nonlinearities of the system represented as and in (45) are unknown, the defined perturbation that contains all nonlinear and time varying dynamics of the system is shown as follows:where is the nominal constant gain.
Then, system (45) can be represented as follows:
4.1.1. Design of Nonlinear State and Perturbation Observer
Defining , , and , the mathematical model of the observer is shown as follows:
Defining , the nonlinear observer is then shown as follows:where is the estimated value of ; is a gain of the observer; is an input error of the nonlinear function; is the precision index; and is the width of the linear area of a nonlinear function.
4.1.2. The Pitch Controller Design
In this controller, the actual perturbation is compensated in realtime by estimating the unknown perturbation of the system with a nonlinear observer. The pitch controller can be expressed aswhere and are two control gains; , , and are the outputs of the observer; and , are the outputs of the tracking differentiator.
4.2. Reconfigurable Controller Design
With the abovementioned stability conditions, we can guarantee the reachability of the sliding surfaces of (25) and of (37). Then presenting the faulttolerant control law as follows:where , is a small design parameter, and is designed such that .
Because it is not as easy to obtain a solution under limit (24) directly in MATLAB, constraint (24) can be transformed into the following inequality condition:where is a small parameter. By using the Schur complement Lemma, we can convert condition (55), equivalently, into
Therefore, we can transform the LMI conditions of (24) and (56) into a problem of LMI optimization to get viable solutions. Recalling the minimization of and , the parameter must be minimized. To this end, we can obtain the LMI optimization problem as follows:where is a fixed parameter to be selected such that . By using the mincx solver in the standard MATLAB LMI Control Toolbox, we can solve the LMI minimization optimization problem.
With the descriptor sliding mode observer (19) proposing, we can now directly estimate the augmented descriptor system state . Note that because is a full column rank, as supposed in Assumption 3, exists. Hence, the estimate of is computed as follows:
Remark 8. To prevent the signals and from becoming chatter, the continuous functions and can be simply used to replace and , where is a small constant [23].
We can summarize the parameters of the descriptor sliding mode method and the faulttolerant control design process into the following four steps:(1)Choose properly such that condition (16) holds, and select the derivative gain matrix in the form (20) such that is nonsingular.(2)Choose such that for . Deal with the Lyapunov Eq. (21) to obtain . Furthermore, compute the proportional gain matrix as .(3)Design the gain matrix of sliding mode, such that . By using the mincx solver in MATLAB, choose a fixed scalar and have the LMI optimization problem (57) solved to get the minimized , , and and the viable solutions , , and . Design as the equation of (26).(4)Design such that . Choose such that is Hurwitz. Let , , and , and design as (43).
5. Simulation Results
In this paper, the faults considered for the pitch system are hydraulic leakage on pitch actuator and a biased output on pitch sensor. Hence, following the design procedures mentioned above, the effective design and verification for the proposed controller can be shown via pitch system (14) with
The simulations are carried out within FASTMATLAB/Simulink combined environment, and the detailed parameters of NREL’s 5MW benchmark wind turbine model are described in Table 1. We can do simulations with the wind speed time signal and its average speed is 16 m/s, as shown in Figure 1.

Let the parameter , which makes sure that condition (16) is available, and select the derivative gain max: ; then, the matrix can be made nonsingular, and and can be calculated.
We select such that for . Via solving the Lyapunov equation (21), we can get the matrix and thus also can calculate the proportional gain matrix as
By selecting the scalar , the LMI optimization problem (57) can be resolved and we can get the practical solutions of and as follows:
We select the matrices and as follows:
The norm controller (53) has the following parameters: , , and .
The simulation results are shown in Figures 2–9. The estimation of the pitch actuator fault signal is redefined according to (9) and shown in Figure 2. However, one problem because of the redefined pitch actuator fault signal is that it is not easy to determine the severity of a fault because it is strongly coupled to the system states. The original pitch actuator fault signal can be constructed by using (10) and the reconstructed original fault is shown in Figure 3. Figures 3–5 show that the proposed descriptor sliding mode observer can obtain a satisfactory estimation of the new augmented pitch system states. With the estimated pitch system augmented states, FTC is carried out with the strategy shown in (43). The FTC results are presented in Figures 6–9. In this part, as for the faulttolerant scene of the faulttolerant controller, the conventional controller is in a divergence state due to the worst case of the fault, so the corresponding effect contrast diagram is not given.
6. Conclusion
In this paper, a new architecture for active FTC controller is proposed for a wind turbine pitch system with simultaneous pitch actuator and sensor faults based on fault estimation and reconstruction to maintain nominal pitching performance. The FTC controller included both a normal controller, which was designed by using the active disturbance rejection control technique for estimating and compensating for the uncertainties, and a reconfigurable controller, which was designed by using the state and fault estimates obtained from an augmented descriptor sliding mode observer. The simulations performed in the FASTMATLAB/Simulink demonstrate that the proposed FTC design can ensure the stability of the generator power output in the event of a fault or in cases without a fault.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The National Natural Science Foundation of China (no. 51205046), Frontier and Applied Basic Research Projects (no. cstc2014jcyjA90005) funded by Chongqing Science and Technology Commission, The 56th Postdoctoral Science Foundation of China (no. 2014M562283), and Special Financial Grant (2015T80956) from the China Postdoctoral Science Foundation and Postdoctoral Special Projects (no. Xm2014003) funded by Chongqing Government are gratefully acknowledged.
References
 T. Burton, N. Jenkins, D. Sharpe, and E. Bossanyi, Wind Energy Handbook, Wiley, New York, NY, USA, 2001.
 T. Jain, J. J. Yame, and D. Sauter, “A novel approach to realtime fault accommodation in NREL's 5MW wind turbine systems,” IEEE Transactions on Sustainable Energy, vol. 4, no. 4, pp. 1082–1090, 2013. View at: Publisher Site  Google Scholar
 C. Sloth, T. Esbensen, and J. Stoustrup, “Robust and faulttolerant linear parametervarying control of wind turbines,” Mechatronics, vol. 21, no. 4, pp. 645–659, 2011. View at: Publisher Site  Google Scholar
 M. Sami and R. J. Patton, “Fault tolerant adaptive sliding mode controller for wind turbine power maximisation,” in Proceedings of the 7th IFAC Symposium on Robust Control Design, ROCOND'12, pp. 499–504, Aalborg, Denmark, June 2012. View at: Publisher Site  Google Scholar
 H. Badihi, Y. Zhang, and H. Hong, “Wind turbine fault diagnosis and faulttolerant torque load control against actuator faults,” IEEE Transactions on Control Systems Technology, vol. 23, no. 4, pp. 1351–1372, 2015. View at: Publisher Site  Google Scholar
 P. F. Odgaard and J. Stoustrup, “Fault tolerant control of wind turbines using unknown input observers,” in Proceedings of the 8th IFAC Symposium on Fault Detection, Supervision and Safety of Technical Processes, SAFEPROCESS 2012, pp. 313–318, Mexico City, Mexico, August 2012. View at: Publisher Site  Google Scholar
 J. Blesa, D. Rotondo, V. Puig, and F. Nejjari, “FDI and FTC of wind turbines using the interval observer approach and virtual actuators/sensors,” Control Engineering Practice, vol. 24, no. 1, pp. 138–155, 2014. View at: Publisher Site  Google Scholar
 Y. Ren, L. Y. Li, J. Brindley, and L. Jiang, “Nonlinear PI control for variable pitch wind turbine,” Control Engineering Practice, vol. 50, pp. 84–94, 2016. View at: Google Scholar
 F. Shi and R. Patton, “An active fault tolerant control approach to an offshore wind turbine model,” Journal of Renewable Energy, vol. 75, pp. 788–798, 2015. View at: Publisher Site  Google Scholar
 J. Lan, R. J. Patton, and X. Zhu, “Faulttolerant wind turbine pitch control using adaptive sliding mode estimation,” Journal of Renewable Energy, pp. 1–13, 2016. View at: Publisher Site  Google Scholar
 L. Chen, F. Shi, and R. Patton, “Active FTC for hydraulic pitch system for an offshore wind turbine,” in Proceedings of the 2nd International Conference on Control and FaultTolerant Systems, SysTol 2013, pp. 510–515, Nice, France, October 2013. View at: Publisher Site  Google Scholar
 M. GálvezCarrillo and M. Kinnaert, “Sensor fault detection and isolation in doublyfed induction generators accounting for parameter variations,” Journal of Renewable Energy, vol. 36, no. 5, pp. 1447–1457, 2011. View at: Publisher Site  Google Scholar
 M. S. Shaker and R. J. Patton, “Active sensor fault tolerant output feedback tracking control for wind turbine systems via TS model,” Engineering Applications of Artificial Intelligence, vol. 34, pp. 1–12, 2014. View at: Publisher Site  Google Scholar
 Y. Song, Y. Lu, and Z. Gan, “Descriptor sliding mode approach for fault/noise reconstruction and faulttolerant control of nonlinear uncertain systems,” Information Sciences, vol. 367368, pp. 194–208, 2016. View at: Publisher Site  Google Scholar
 K.S. Lee and T.G. Park, “An actuator fault reconstruction scheme for linear systems,” Journal of Process Control, vol. 44, pp. 106–119, 2016. View at: Publisher Site  Google Scholar
 L. Chen, X. Huang, and S. Fu, “Observerbased sensor faulttolerant control for semiMarkovian jump systems,” Nonlinear Analysis: Hybrid Systems, vol. 22, pp. 161–177, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 M. Liu and P. Shi, “Sensor fault estimation and tolerant control for Itô stochastic systems with a descriptor sliding mode approach,” Automatica, vol. 49, no. 5, pp. 1242–1250, 2013. View at: Publisher Site  Google Scholar
 T. Wang, H. Gao, and J. Qiu, “A combined adaptive neural network and nonlinear model predictive control for multirate networked industrial process control,” IEEE Transactions on Neural Networks and Learning Systems, vol. 27, no. 2, pp. 416–425, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 O. Bagherieh and R. Nagamune, “Gainscheduling control of a floating offshore wind turbine above rated wind speed,” Control Theory and Technology, vol. 13, no. 2, pp. 160–172, 2015. View at: Publisher Site  Google Scholar
 J. Q. Han, “From PID to active disturbance rejection control,” IEEE Transactions on Industrial Electronics, vol. 56, no. 3, pp. 900–906, 2009. View at: Publisher Site  Google Scholar
 J. Chen, L. Jiang, W. Yao, and Q. H. Wu, “Perturbation estimation based nonlinear adaptive control of a fullrated converter wind turbine for fault ridethrough capability enhancement,” IEEE Transactions on Power Systems, vol. 29, no. 6, pp. 2733–2743, 2014. View at: Publisher Site  Google Scholar
 L. Zhao, Y. Yang, Y. Xia, and Z. Liu, “Active disturbance rejection position control for a magnetic rodless pneumatic cylinder,” IEEE Transactions on Industrial Electronics, vol. 62, no. 9, pp. 5838–5846, 2015. View at: Publisher Site  Google Scholar
 L. Chen, M. Liu, and S. Fu, “Adaptive sliding mode control for stochastic jump systems against sensor and actuator failures,” IET Control Theory & Applications, vol. 10, no. 16, pp. 2000–2009, 2016. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2017 Lei Wang 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.