Numbers, Please: Power- and Voltage-Related Indices in Control of a Turbine-Generator Set

: This paper discusses the proper selection and interpretation of aggregated control performance indices values mirroring the quality of electrical energy generation by a turbine-generator set cooperating with a power system. Typically, a set of basic/classical and individual indices is used in energy engineering to ensure the mirroring feature and is related to voltage, frequency and active or reactive power deviations from their nominal values desired in the power system. In this paper, aggregated indices based on the sum of weighted integral indices are proposed, veriﬁed and built based on the well-known indices originating from control theory. These include an integral of the squared error (ISE) and an integral of the squared error multiplied by time (ITSE), applicable whenever an in-depth analysis and evaluation of various control strategies of the generation system is to be performed. In the reported research, the computer simulation tests veriﬁed their effectiveness in assessing the generated electricity on the example of a turbine-generator set controlled using a predictive control technology as well as applicability, proven by numerous simulation results to take various and different in nature requirements into account efﬁciently, in the form of a single aggregated index.


Preliminaries
The energy policy of the European Union forces progressive modernization changes in the structure of generation sources in power systems related to reducing the impact of energy on the environment, especially for technologies related to the energy combustion of fuels. The generation system in Poland is based 70% on coal-fired power plants. Due to the high emissivity of these sources and the high costs of the fluid gas treatment installation, the Polish Energy Policy by 2040 [1] aims to introduce progressive changes in the energy generation structure. According to the assumptions [1], the transformation of the generation system in Poland is to be based on the replacement of conventional power plants with renewable energy sources and nuclear energy. The first block of the nuclear power plant with a capacity of 1-1.6 GW is to be commissioned in 2033, for which the location was already selected at the end of 2021. The government company, Polish Nuclear Power Plants, selected the coastal location of Lubiatowo-Kopalino in the Choczewo commune in Pomerania for the first nuclear reactor [2]. In the following years, it is planned to launch another five such units at intervals of 2-3 years. As for renewable energy sources, it is expected that by 2040 the installed capacity of renewable sources will increase almost twice, from the current 15% to 28.5% [3]. However, due to rapid climate changes, the transition towards a non-nuclear renewable energy system remains a major political challenge, which can be made even worse by increases in energy supply chain uncertainty given rise by the shift away from fossil fuels [3]. The described development activities are aimed at replacing the shut down of conventional power plants in a way that conditions the fulfillment of the power balance in the system, as well as a significant reduction in the national emission of greenhouse gases and air pollutants. However, from the point of view of introducing a significant share of renewable sources, it is necessary to take the stochasticity of their work and significant fluctuations in generated power into account, which are associated with the need to equip classic coal or nuclear units with modern control systems capable of responding to dynamic changes in the power generated in the system power industry. Therefore, the optimization and analysis of dynamic responses of power unit control systems is a critical issue from the point of view of the planned development and upgrade in the generation structure of the Polish Power System. The latter has been carried out in the reported application by the use of integral-based performance indices.
There is a number of indices for evaluation of electrical energy quality, which differ in their nature and allow simultaneous, and yet independent from one another, evaluation of various characteristics of the electric power generation process. In the case of control system design, a single performance index is usually needed to clearly mirror the control quality. Introduction of the indices is, however, necessary both for the evaluation of the operation of control systems and the synthesis of controllers.
The aim of the research work described in this paper is to find the index which is related to a set of those used hitherto, as per relation to the same phenomena, and which is at the same time useful for implementation and analysis of control systems.
As the quality of electrical energy is a major factor when the development of modern societies is concerned, as well as there is a growing demand on the generation of electric power in accordance with the speed of economic development of societies, one can clearly identify a growing need to increase power plants' efficiency and improve electrical energy quality.
In order to develop solutions improving this quality, it is necessary to have appropriate methods for their assessment. Thanks to the suitably selected indices, it is possible to compare solutions and seek the superior one which satisfies all the requirements. Currently, various indices are used, derived from the control theory and energy engineering. To facilitate the synthesis and analysis of control systems, one aggregated index is needed. This article proposes to use the integral of the squared error (ISE) or integral of the squared error multiplied by time (ITSE) indices as candidates to mirror the quality of operation of the turbine-generator control system. Obviously, the performance of the turbine-generator control loop is a function of a suitable design and tuning of its controller.
In addition, as far specific characteristics of a turbine-generator set are concerned, using the proposed control method (by appropriate tuning of performance indices) clearly leads to improvement in the quality of energy generated by the turbine-generator set. It reduces rotor speed fluctuations in dynamic conditions for arbitrary load changes or for changes in a shaft torque, being typical cases in the system with renewable energy sources present. A proper tuning of controllers is followed by the possibility to fix the properties of transient-and steady-state performance in closed-loop systems [4]. To be able to evaluate the properties of such systems in a qualitative manner, performance measures need to be introduced first. These should be related to the desired properties of specific responses, and thus should be expressed in terms of both transient-and steady-state performance [5].
The obtained closed-loop system, after tuning, is a result of a compromise between the aforementioned factors, as usually steady-state error rejection is the adverse requirement with respect to the stability margins one, and it is impossible to meet contradictory aims simultaneously [6]. A shortlist of these might include: good regulation against disturbances, desirable response to commands, maintaining low amplitudes of critical signals, etc. To override this problem, it is advantageous to introduce integral performance indices to conveniently mirror the performance of the control system [7].
The integral indices assessing the control performance are widely used at the optimization stage in a minimization task of a selected performance-related cost function. To make the results portable across various platforms and applications, asymptotic tracking is usually evaluated for a span of different types of reference signals, such as a unit step, ramp or parabolic inputs, as a basic control task in feedback systems [8]. In such an application, the presented integral performance-related indices, using the information concerning the tracking error in the closed-loop system, are used as quantitative measures of the control quality. In this case, the attained minimum, coupled with selection of controller gains, mirrors the fact that optimal control policy is obtained [9].
Typically, quality indices found in the literature focus either on long-term quality assessment in order to analyze the network or on evaluating electric power suppliers. In [10], the indices are proposed to analyze the quality of electricity every month to give scores to different utility companies, whereas in [11], the authors analyze the situation in the Italian power system in the long term, looking into the continuity of supply and differences in quality between the north and the south of the country. In [12], the usage of wavelet packet transforms to measure non-stationary power quality disturbances, and in [13], wavelets are used to visualize time-varying power quality indices. Some authors use probabilistic methods to assess the power quality using Markov models [14] or the Monte Carlo method [15]. The authors of [16] use neural networks to classify the disturbances occurring in the system offline. Paper [17] contains an overview of commonly used voltage characteristics such as 10 min voltage and total harmonic distortion factor (THD) averages, 10 s frequency averages or 2 h long-term voltage fluctuation analysis. What is more, the authors of [18] use machine learning and artificial neural networks to analyze the power quality in a specific site using common parameters counted as 10 min averages as input. The aforementioned approaches focus on a single power quality index (e.g., total harmonic distortion). On the contrary, the global quality indices (GPQI) [19] proposed in the literature are based on 10/15 min time frames for which quality indices are calculated and then aggregated into a single index showing the global quality.
Despite a very valuable analysis of power quality indices, all of these approaches have, unfortunately, general (site assessment and company assessment) and long-term (10 min, months and years) focus and cannot be used in the case of control system analysis and synthesis where physical processes occur within milliseconds. Therefore, a new approach focused primarily on this particular issue is needed, which constitutes the motivation to undertake the proposed research project. The proposed quality index is expected to be able to properly assess power quality but at the same time must be useful for control systems analysis, synthesis and parameter tuning.
The presented literature review shows some separation of electrical engineering from control engineering. New publications in the field of control engineering analyze new approaches in isolation from the specific case of controlling the turbine-generator set of a nuclear power plant. Publications in the field of electrical engineering, in turn, focus on a detailed analysis of individual energy parameters, or on a global analysis of energy quality but over long periods of time. The presented article connects the achievements of control theory with the achievements of electrical engineering, and thus fills a certain gap between these areas. The conducted research is the starting point for further development of better control systems in the power industry.
The main novelty of the solution presented in the paper is the use of global quality indices with characteristics suitable for use in the synthesis of used control systems. The presented solution takes the analysis of critical parameters from the point of view of the control system of the turbine-generator set and analyzes the data taking very high resolution into account. This is to correctly assess the quality of the control system and to find the optimal controller parameters, which is not possible with the use of widely used individual, classic quality indices.
The main contribution of the work is the confirmation that the integral indices of electric power quality developed in this paper allow to determine the quality of the control system with a single value capturing its viable characteristics. This type of approach significantly facilitates the analysis and synthesis of control systems. The model predictive control technique was applied to control the turbo-generator set, and its synthesis is based on the proposed indices.
The paper is structured as follows: Section 1.2 describes the electrical power quality definition used within this paper. Section 2 discusses the problem including the description of the power quality indices that can be found in the literature. Section 3 describes the turbine-generator set control problem which is the basis for further considerations, while Section 3.2 presents the proposed integral index. Simulation results of computing indices' values for several control systems with different parameters are presented in Section 4. The article ends with a conclusion resulting from the conducted research and lists the directions for further research.

Electrical Power Quality
As per the selected subject of the presented research, it is to be mentioned that the definition of energy generation quality depends on the loads considered and the way in which energy parameter deviations affect their final values. The definition is very customer-oriented, where a quality-related problem is defined as follows [20]: Any power problem manifested in voltage, current, or frequency deviations that results in failure or misoperation of customer equipment.
A few types of such quality problems can be found in the literature, for example, see [17,21], among which one can list amplitude related, waveform disruption related, balance related and frequency related ( Figure 1). As a nearly perfect sine-wave, voltage is generated by the generator [20], the currentand harmonics-related problems are omitted in this paper, since the turbo-generator set control system has no influence on these quality problems, and hence, including these phenomena in the control-oriented indices would not result in any solution improvement. Therefore, only voltage dips, swells, undervoltage, overvoltage, voltage fluctuations and frequency variations are considered in the paper, (highlighted in bold in the previous paragraph).
The turbo-generator set control system should not introduce such problems as they can influence a large group of energy consumers [17], namely: CNC machines, adjustable speed drives, personal computers, programmable logic controllers (PLCs), relays, contactors, motor starters, fax machines, metal-halid and high pressure sodium lighting, telcom switching equipment, electronic ballast fluorescent lighting, etc.
When solving voltage and frequency stabilization control problems resulting with no disruptions in the electrical power system (EPS) given rise by the turbine-generator set, all the considerations presented above must be taken into account. To ensure high quality of stabilization, quality indices must be introduced to properly measure the results of the controller's actions, and mirror the expectations of the designer. From this viewpoint, the power quality index for control purposes is defined in the paper as: The index that can quantify the behavior of the power system assessing the scale of power problems manifested in voltage, current, or frequency and that can be used for control system synthesis and evaluation.
The next Sections describe the details of the introduced control problem and introduce such control-oriented indices in a form of ISE/ITSE integrals.

Problem Description
Conditions of normal operation usually stipulate nominal values of signals from the networked system. These have also been used here to construct performance indices measuring the degree of deterioration from the expected conditions of work.
In contrast to the energy-generation perspective, the following time-related indices can be listed from a control-engineering-oriented point of view among others from this discipline: • Overshoot-the largest, transient, deviation from the voltage set point; • Stabilization time-the time after which the error is smaller than a certain deviation r ; • Rise time-the time for the controlled variable to rise from 10% to 90% of its final value or ∆Y/∆t rise rate.
The instruction for operation and exploitation of the transmission network [22] gives the value t r ≤ 0.3 s for the generator voltage controller U g for r ≤ 0.5% and a jump of 10% of the value and p ≤ 10% at generator idle state and p ≤ 15% at start-up, and the regulation speed must be greater than 1.5 U n /s. In the literature [23], one can also find a presentation of the criteria for assessing the quality of regulation for electrical energy production, split into the following sub-criteria: • Standard deviation for which (1) is calculated for n measurements every 15 min during one month • Amount and time of deviations greater than 50 mHz from the nominal frequency; • Trumpet characteristics must be met by the control system and are the requirement to contain the frequency waveform inside the trumpet characteristic after a sudden load step change. The curve is described by a pair of Equations (2) and (3), and defines an envelope-like bound on the transients where: t = 900 s is a stabilization time of 15 min; the constant A defines the trumpet characteristic width factor depending on the size of the power disturbance and the characteristics of the power system. The Polish power system belongs to the European UCTE system, for which the abovementioned width of the trumpet characteristic is defined as (4). where: α is a share of the regulatory area in energy production, K fMW/Hz denotes the power frequency equivalent defining the change in power in MW for a change in frequency by 1 Hz determined for the entire system.
On the basis of these requirements, it can be determined that with jumps of several hundred MW (e.g., failure of a generating unit) the permissible temporary frequency change is several dozen mHz, and after 15 min the deviation should be less than 20 mHz [23].
Additionally, a harmonic distortion can be evaluated by calculating the total harmonic distortion (THD) factor [17], but as mentioned, this type of disruption is be analyzed for control purposes in this paper.
In the literature, the aggregated global power quality indices (GPQI) can also be found. In [19,24,25], two such indices are proposed and used in a form of: ADI is a sum of a set of quality parameters calculated in the 10-min time frames: where: i is the number of the quality index; k i are the weights to balance all the components; W i are power quality parameters such as frequency change, voltage level, voltage variation, flicker severity, voltage unbalance, harmonic distortion and voltage change.
The second index, FDI, is used to aggregate all the event-related quality problems, i.e., dips, swells and interruptions. This index determines in percent how many of the 10 min periods were disrupted: where: f is the number of the 10 min period with disruptions, n is the number of all 10 min periods.
ADI and FDI are used in [19,24,25] for the network with distributed generation analysis. The cons of this global power quality indices in this form are that they are calculated over a longer period of time using long time frames and that they also include quality parameters not essential for control purposes.
Therefore, for the control system synthesis, a more precise index is needed. The paper proposes using the integral power quality indices that operate on the shorter time frames and leverage only the control-related parameters. The integral indices assessing the control performance are widely used in the controllers' analysis and synthesis to optimize the performance of the control system and-as a result of that-the whole controlled plant. As mentioned above, optimizing the operation of nuclear power plants is an important issue due to the dynamic development of renewable energy and thus control-oriented indices are needed to properly assess the power quality and performance of the plant. Said indices are presented in the following section.

Model and Methods
The classic generator regulation system consists of an excitation controller; its task is to maintain a constant value of the generator voltage and a power system stabilizer (PSS). The additional system of the system stabilizer, through the correction of the set value, minimizes the oscillations of the active power transferred to the power system caused by the operation of the excitation controller. The system stabilizer generates an additional control signal for the generator controller, which dampens potential electromechanical oscillations and improves the dynamic stability of the turbine-generator set connected to the power grid.
In most cases, linear generator controllers are used, which in majority are based on the proportional-integral control (PI) principle. These systems are additionally equipped with a number of modifications and artificial constraints as to ensure safe operation of the turbine set. In addition, these systems receive voltage and power set values from external control systems that ensure appropriate parameters throughout the whole electrical power system (EPS) and appropriate load distribution between generating units.
The steam turbine works with a synchronous generator connected to the power grid. The turbine PI controller regulates the active power of the generator by manipulating the steam flow to the turbine and consequently affecting the torque on the assembly shaft. The active power of the generator is proportional to the torque, so the supplied mechanical energy and the received electrical energy must be in balance with the accuracy of lossesmechanical, thermal and other.
On the other hand, a PI controller of a generator regulates the voltage at the generator terminals by manipulating the excitation voltage (the non-linear relationship between the controller signal and the rectifier voltage is compensated accordingly). As shown in Figure 2:  P g -active power; Q g -reactive power; U g -generator's voltage; I g -generator's current; ω g -generator's angular speed; E f d -excitation voltage; U s -Power system's voltage; ω g -Power system's voltage frequency; p-steam pressure; α s -control valve opening; α u -steam vent valve opening for the heat generation; m-mass flow of the steam for heat generation; α s -control valve opening.
As already mentioned in the study, appropriate models of simplified components were used for the synthesis of the turbine set control systems. The starting point for their development was the structure of the turbine set input-output model (Figure 3), which input quantities were the control valve opening degree (α) and the generator excitation voltage (E f d ), respectively. The outputs can be listed as the power and voltage of the generator (P g and U g ), with the thermal load (Q) being treated as a disturbing input. Such a departure naturally leads to the internal structure of a simplified model built on the basis of four main information processing paths in the following input-output configurations:

Model Predictive Control
In this paper, instead of the typical integral-derivative blocks of the controller, the quadratic dynamic matrix control (QDMC) variant of an MPC controller for the purpose of control of a synchronous generator is suggested [26]. In order to obtain improved quality in the closed-loop system, i.e., the quality in a control loop taking the performance of the turbo-generator, some extent of the exchange in information between the quadratic dynamic matrix controller and the environment form is required. This is performed by introducing an additional signal between the controllers, i.e., rotational speed ω, active power P g or the steam turbine's control valve opening degree α. In addition, in order to also take the continuously changing setpoint of the system, a recursive least-squares (RLS) algorithm is adopted here.
The proposed approach is not standard in the power generation industry or in power plant control, though it offers an easy replacement to the PID control-based algorithms [26]. As per the optimized solution of the MPC problem, not only does it take a complete set of plant-related constraints into consideration [27], but it also offers the optimal solution of the defined control problem. Since this approach can be easily reconfigured (for example by stipulating appropriate horizons of control, prediction, etc), the performance of the control law can easily mirror the expected behavior of the closed-loop system. In this paper, the parameters are adopted as to mimic the requested properties of the process.
The block diagram of the system consisting of a predictive controller (MPC) and the recursive estimation scheme (RLS) is presented in Figure 4. The solutions obtained for this considered control law, are compared among one another for various sets of configuration parameters to distinguish between those offering poor, moderate and superior performance. The linear step-response model of the process is used in classical QDMC methods to estimate its future states [28,29]. Taking both free and forced response components for a multi-variable system with s inputs and r outputs, the following model is considered: with sample number denoted as k and conditional estimate of the future sample as y k+1|k .
The sought sequence of control signal updates calculated in m steps ahead is denoted as ∆u k ∆u k = [(∆u 1(k) , . . . , ∆u s(k) ), . . . , (∆u 1(k+m−1) , . . . , ∆u s(k+m−1) )] T , comprising m values for all s inputs. For the details of the derivation of the model, and for the sake of brevity of this paper, please consult [30][31][32][33]. The optimal control update is found by solving the following problem at every sample k [28,29] min ∆u min ≤ ∆u k ≤ ∆u max , with u k as a vector of control signals calculated on the basis of updates, Γ > 0 as a weight matrix and Λ ≥ 0 as control update penalty term. The notation from Table 1 is used for the described model. The constrained quadratic programming task in (9) is used to calculate optimal control updates, implemented to the plant on the basis of a receding horizon rule. The control updates are found with a sampling period compatible with the documentation of the GTHW-600 generator [34], of which the dominating time constant is estimated atT = 0.0017 s. As per a rule of thumb, to ensure at least 10 sampling periods in a dominating time constant of a model, the sampling period of T = 0.00001 s was selected. In relation to this period, the output prediction horizon should fall within 10 ÷ 20 samples, or 0.0001 ÷ 0.0002 s in the time domain.
The QDMC scheme can be adopted in the structure as in Figure 5. In order to obtain an optimized performance, the time-related response of the model is calculated on the basis of its identified model, following possible changes in operating points. The discussion considering a proper selection of the orders of the model of the system can be found in [30][31][32][33]. At every discrete time instant (i.e., every action taken by the QDMC controller), and for the considered structure of a discrete-time model of the plant, its parameters are estimated via RLS scheme, the step response of the multi-variable model is calculated and adopted into the MPC scheme.
It is to be mentioned here that solving an MPC problem via QDMC framework is not a bottleneck, as the solution is obtained on the basis of either efficient interior-point algorithms or on the basis of an active-set method.
At the moment, the MPC algorithms are implemented in FPGAs or efficient microcontrollers [35] and microcontrollers [36,37] to allow real-time performance, or to offer a distributed approach to MPC control of the nuclear power plant turbine set in real time.

Methodology
In order to develop an MPC controller, it is necessary to define a cost function that is used in searching for the optimal control updates which, in the context of this paper, are such that a signal for which a minimal control deviation is attained while stabilizing the voltage and angular speed of the generator while following the set active power trajectory. Each of the criteria presented in Section 2 is of a different nature and are responsible for limiting the impact of disturbances on the controlled quantities (frequency, amplitude and decay time of voltage and power and generator angular speed oscillations) or takes parameters calculated over longer periods of time into consideration. It is necessary to find the appropriate function that allows to aggregate all the requirements resulting from the operation and operating instructions. An additional component of the criterion function, which it can be extended with, are the values of control signal updates (change in the excitation voltage ∆E f d and change in the opening of the control valve ∆α). Such an extension of the criterion allows one to determine the quality of the regulation in relation to its cost (energy expenditure on control effort).
Advanced control systems can use the model of the power plant's turbo generator set as a whole, and thus, the quality indices must take the impact on the quality of both the amplitude and the frequency of the voltage into account. Therefore, ISE-and ITSE-related indices are proposed for assessing the quality of regulation. The ISE criterion takes the following controlled quantities: voltage U g , power P g and frequency ω g combined into a sum of squares-like expressions with a set of weights that determine the share of each component in the sum. Meanwhile, the ITSE index should additionally allow one to assess the accuracy of voltage amplitude and frequency stabilization (minimization of the error in normal operation) and to take the stabilization time into account by increasing the error severity at later times (reduction in disturbances according to trumpet characteristic). Contrary to the ISE index, which is an integral of the squared error (10), it allows not only the amplitude of the disturbances to be taken into account but also the rate of their suppression, as per the weighting factor interpretation of the time product. The proposed ITSE index is used to compare the quality of various analyzed control systems. This criterion can be represented as (11). It is a starting point for further considerations regarding the selection of the best control system selected by means of-first-parameter optimization for a proposed control structure-and next-comparison of a series of different solutions (i.e., selecting the one that provides the lowest value of the index compared with the others).
Criteria (10) and (11) aggregate quality indices (1)-(4) by taking into account deviations from the set points of voltage U g , active power P g and frequency (ω g ). Thanks to this, all deviations and oscillation fluctuations of these quantities are taken into account. The authors have decided to use the above-listed criteria, modified with the use of a linear combination of sub-integral expressions, as per their ability to capture viable properties of the well-tuned control system in power generation tasks. Control performance in this case replaces a number of classical power generation-related criteria, with a single number capturing important characteristics.
In order to assess the deviations of various origin in a single criterion, it is necessary to assign appropriate weights. Therefore, relative values of these quantities related to their nominal values were used. Additionally, the weights (a, b, c, d, e) were used to enable the change in the influence of each component of the sum on the final value of the index in order to increase the impact of one of the component's criteria.
The article analyzes two different sets of weights to illustrate their impact on the results and usability of the proposed indices. As the base case, a set of weights equal to one was selected (a = 1, b = 1, c = 1). Since the sum component responsible for the active power P g (weight a) is significantly dominant, the second case was chosen, in which a lot of emphasis is placed on changes in the ω g speed and changes in the voltage U g , i.e., with the weight a = 1, the weights b = 1000 and c = 1000 were adopted. This is to increase the sensitivity of the index to the appearance of unwanted oscillations in the system. where: The weights a, b and c were selected so that it is possible to:

Results and Discussion
Two types of experiments were performed. The first one was to generate two dimensional surfaces of the ISE/ITSE indices' values for a control system with two parameters (turbine controller's prediction horizon p T and generator controller's prediction horizon p G ). The latter allowed the analysis of the ISE and ITSE indices in a wide range of control cases (672 values-2 indices, 3 sets of weights and 112 different simulations). Based on these results, it is possible to select the best and worst pair of values p G − p T for each of the cases (please refer to Table 2   In the second experiment, the simulations were performed in order to compare the process values from the turbine-generator set for those nine cases: the best and worst values of the different indices. The aim of this experiment was to find the relation between the indices and the performance of the plant (specified as the analysis of the waveforms of active power P g , generator's voltage U g and generator's angular speed ω g ). The results are presented in Figures 9-11.
In the first case (Figure 9), the biggest oscillations are generated by the control system C (the worst ISE/ITSE value). The system B follows the active power trajectory faster but at the cost of the higher angular speed and voltage oscillations. Similarly, in the second case (Figure 10), the system with the worst ISE/ITSE values generates the biggest oscillations. Due to internal interactions between the generator and the turbine, also in this case a slightly better active power control of the system E, worse angular power and voltage stabilization (which is nevertheless better with the system E) is caused. In this case, the differences are not as significant as in the first case. Finally, in the third case (Figure 11), system I introduces the oscillations with the biggest amplitude, while system H follows the active power trajectory best and system G introduces the oscillations with the smallest amplitude. In all the cases (Figures 9 and 10), the oscillations are reflected by the rapid changes of the control signals (control valve opening and excitation voltage), which in real-life application would damage the actuators. Therefore, the systems C, F and I can only serve as an image of the wrong selection of the controller parameters and not systems that may be used in practice. This undesired system behavior is also reflected by the values of the ISE/ITSE indices.   The simulations performed for each group show that the behavior of the system with the worst parameters according to the ISE/ITSE criteria is, as expected, the worst (C, F, I), i.e., there are significant oscillations both in power, amplitude and voltage frequency (angular velocity). The increase in the value of integral indices is clearly reflected in the deterioration of the course of the controlled values. The waveforms for the simulation with the use of the best parameter sets for ISE (sets A, D and G) and ITSE (sets B, E and H) are characterized by much smaller oscillations, lesser overshoot and shorter stabilization time. Improved values of ISE/ITSE indices correspond to better parameters from the point of view of control quality (i.e., classic indices of the control system quality). This also translates into the improvement in the quality of electricity by reducing the amplitude and frequency fluctuations of the voltage.
For the best pairs of ISE/ITSE results for A/B and G/H sets (reduced influence of active power P g on indices in relation to voltage U g and angular velocity ω g ), the power stabilization time is shorter for the best ITSE sets than for ISE sets, but at the cost of increasing oscillations in the amplitude and frequency of the voltage. For the E/F case, i.e., ISE/ITSE, the indices for all weights a, b and c are equal to 1 (equal share of components in the criterion); the waveforms for the ITSE index (set F) are characterized by both shorter regulation times and smaller oscillations of the voltage amplitude and frequency (angular velocity).
The above results are combined together to show the best ISE/ITSE cases for three different sets of weights a, b and c, i.e., comparison for ISE: • The best ISE for weights' set A (a = 1, b = 1000, c = 1000); • The best ISE for weights' set D (a = 1, b = 1, c = 1); • The best ISE for weights' set G (a = 1, b = 100,000, c = 100,000). Comparison for ITSE: • The best ITSE for weights' set B (a = 1, b = 1000, c = 1000); • The best ITSE for weights' set E (a = 1, b = 1, c = 1); • The best ITSE for weights' set H (a = 1, b = 100,000, c = 100,000). The results in this approach are presented in Figure 12. Stabilization time of the active power P g for the set D (weights a = 1, b = 1, c = 1) is shorter as the power takes the dominant part in the index formulation (as mentioned, the active power error has a larger order of magnitude than the other components). Increasing the weights for the generator voltage U g and angular velocity w results in smaller oscillations but also worsens the active power stabilization time. It is caused by aforementioned interactions inside the turbine-generator set.
The ITSE index sets B and E are similar (in both cases the same set of parameters resulted in the smallest ITSE value) and,therefore, the graphs coincide ( Figure 13).
As in the ISE case, higher values of the b and c weights results in smaller oscillations of the amplitude and frequency (angular velocity) of the generator's voltage (set H). For ITSE, smaller differences in weights resulted in equal parameters and outputs. Only after drastic weight changes the expected reduction in oscillations is obtained. As mentioned before, this reduction is the desired behavior of the turbine-generator set's controller and, therefore, this set of parameters results in better quality electrical energy, despite the stabilization time of the active power.
Taking all the results into account, set G mirrors the most desired control system behavior as wave-forms resulting from usage of this parameters' set are characterized by the smallest oscillations of all cases (at the expense of slower power regulation). The change in the turbine's control valve opening and the excitation voltage is also milder for this set of parameters, which has a big impact on the actuators exploitation and their up time.

Conclusions
The paper presents the electric energy quality indices typically used to define energy quality and proposes the use of the integral index for the synthesis of control systems. Despite the fact that multiple and various indices of the quality of electricity can be found in the literature, they still most often define quality for the purposes of long-term evaluation of energy suppliers, or in fifteen-minute terms for the purpose of checking compliance with the requirements. Both approaches are not sufficient for the synthesis of control systems.
For these purposes, a preferably single index is needed to take the transient states occurring in the system into account, i.e., changes over a period of milliseconds. This paper proposes ISE/ITSE indices in the synthesis of control systems as they are able to combine a set of different requirements in one value, also taking into consideration fast changing values. The paper consists of a in-depth analysis of the proposed indices taking a number of different approaches into consideration, such as the distinction between ISE and ITSE and three different sets of weights in each of the criteria. The presented simulation results are to illustrate how the proposed solution can be used to compare different power and voltage waveforms and what the relation between the values of the ISE/ITSE indices and the actual quality of electricity is.
The proposed solution was designed for control systems synthesis and their quality assessment and is not designed to replace existing indices used to assess compliance with the regulations. As being control oriented, these indices were used to design a set of controllers of a turbine-generator set in a nuclear power plant: a QDMC model predictive controller [38], a distributed model predictive controller [30], a fuzzy controller [39] and a controller using gain scheduling [40]. They were used to tune the control systems, i.e., to optimize the parameters of controllers, and to compare the quality of several different solutions. The results obtained during the aforementioned studies show the practical usability of the proposed indices in the control system synthesis.
The next step should be to define the range of proper, acceptable and bad values, which would enable an independent assessment of the operation of the control system. Without specific standard values for the index, they are only useful in the comparative analysis of two or more systems and cannot be used as an objective point of reference.
The novelty of the presented paper is the use of knowledge in the field of control theory in practice and the use of known ISE/ITSE indices to solve the problem of unambiguous assessment of electrical quality for the needs of the turbine set control system. Although these indices are widely used in automation, they are not used for this purpose in power engineering. This is to facilitate the synthesis of better control systems and the improvement of the quality of control in the power industry in the face of challenges related to the changing nature of the energy sector, i.e., the growth of renewable energy generation.