Electrical Damping Assessment and Sensitivity Analysis of a Liquefied Natural Gas Plant: Experimental Validation

: Liquefied Natural Gas (LNG) plants are commonly island-operated weak grids where the interaction of high-power Variable Frequency Drives (VFDs) with the Turbine-Generator (TG) units might cause Sub-Synchronous Torsional Interaction (SSTI) phenomena. SSTI phenomena can lead the LNG plant to instability conditions. Each LNG plant configuration is characterized by a risk level, which is considered high when the electrical damping at the TG Torsional Natural Frequencies (TNFs) is negative. Starting from a real case study, a detailed electromechanical model of an LNG plant is presented. The model is comprehensive of the control system of the power conversion stage and of the TG unit. Sensitivity analysis, performed on control system parameters, allows one to detect the parameters that impact the electrical damping and the stability of the overall LNG plant. A complete simulation platform is developed. Experimental results are carried out on a real LNG plant considering four different configurations. The theoretical model and the simulation platform allow one to estimate the electrical damping and the results are confirmed by the experimental validation. It is demonstrated that fine tuning of the power conversion stage control parameters can reduce the risk related to torsional instability.


Introduction
Up-to-date Liquefied Natural Gas (LNG) plants are highly electrified and based on high power electrical motors supplied by Variable Frequency Drives (VFDs). The interaction among the VFDs and the Turbine-Generator (TG) units may cause torsional vibrations known as Sub-Synchronous Torsional Interactions (SSTIs) [1][2][3]. Over the last years some SSTI phenomena have been experienced on site in the LNG plants by the authors. One of the most significant occurred during the commissioning phase and the event is used as a case study for the present paper.
When SSTI phenomena occur, the torsional vibrations measured on the TG shaft-line generate electric disturbances such as voltage fluctuations in the power system. The voltage fluctuations in input to the VFDs imply current fluctuations at the VFD DC link and at the VFD output. Finally, high torsional vibrations can lead to torsional instability.
In the literature the first studies focused on the torsional instability risk assessment are related to the Sub Synchronous Resonance (SSR) phenomenon [4][5][6][7]. In [4], the generator damping, measured at no load condition, is used as an index value to define the risk. In [5] the compensation impact factor is used to detect SSR risks in the case of series compensated transmission lines. In [6] the impact of the transmission expansion on the system damping is evaluated using the frequency-scan approach. In [7] the SSR risk is detected comparing the electrical damping with the modal damping.
The most used risk assessment technique is based on the Unit Interaction Factor (UIF) calculation, originally introduced for HVDC applications [8]. For oil and gas applications, the method implies high variations of the threshold value adopted to evaluate the risk, for this reason the UIF technique is not always feasible.
In this paper, the risk assessment for an LNG plant is based on the electrical damping evaluation. The torsional instability risk can be considered high when the electrical damping at the Torsional Natural Frequencies (TNFs) is negative. Hence the risk assessment relies on the electrical damping assessment.
In [1] a connection between the power electronics analysis and the power systems analysis is established and a complete model of a LNG plant is provided combining the dynamic model of the power conversion stage and the TG electromechanical model. The sole theoretical model allows electrical damping assessment in the case of a basic LNG plant configuration with one TG unit and one VFD. However the analysis provided in [1] lacks a deep investigation about the impact of the synchronization system design on the electrical damping estimation. In the same model the design of the TG controller is not considered.
Some previous studies such as [9,10] have shown that the synchronization system parameters affect the power system stability in the case of weak grids. Besides the impact of the Power System Stabilizer (PSS) on the SSR damping has been previously demonstrated in [11][12][13][14]. Aware of these considerations and starting from the combined electromechanical model developed in [1], the same authors present in this paper a more accurate model providing the electrical damping assessment. Differently from the model presented in [1], both the control systems of the VFD and of the TG unit are included in the overall model. It allows one to detect which control parameters have more impact on the generator units' electrical damping. Moreover a classical tool such as the sensitivity analysis [15][16][17] is adopted for this purpose. In particular the Finite Difference Method (FDM) [18] is employed to achieve a local sensitivity study. The overall advanced model is the first original contribution of the paper compared to [1].
Nevertheless, it has to pointed out that in the case of complex LNG plant configurations, with numerous TG units and several VFDs, the calculation of the electrical damping by the theoretical model can exhibit some practical limitations. Starting from this consideration, a complete simulation platform, which emulates the theoretical model, is developed and the results are shown in the paper. The complete simulation platform allows one to analyze also the complex LNG plant configurations and to extend the results of this study to different plants. The flexibility relies in the possibility to adapt the software to different power ranges and a different kind of loads, including also changes of the operational TG units' number in real time. Finally, an extensive test campaign has been carried out on a real LNG plant in order to take into account also commissioning and potential grid contingencies. These products represent innovative contributions compared to [1].
In conclusion, the theoretical outcomes are confirmed by the simulation and experimental results. The electrical damping assessment allows one to determine the risk of the LNG plant considering different configurations. It is also demonstrated that the proper tuning of the control parameters can modify the electrical damping and, as a consequence, can influence the LNG plant level of risk.
The rest of the paper is organized as follows: in Section 2 the LNG plant configuration is presented; in Section 3 and in Section 4 there are proposed, respectively, the detailed TG unit and VFD models and a preliminary sensitivity analysis is performed; Section 5 treats the LNG plant simulation platform; in Section 6 the simulation and experimental results are shown; in Section 7 a brief discussion about the future insights and possible applications is presented; finally Section 8 is focused on final remarks.

LNG Plant Configuration
The LNG plant under analysis is shown in Figure 1. The considered power system is composed of three identical Gas Turbines (GTs) connected to three identical Synchronous Generators (SGs). For the sake of simplicity, Figure 1 represents only the SGs while the complete TG unit is depicted in Figure 2. The SGs are connected to the Point of Common Coupling (PCC) through three step-up transformers, denoted as TTG1, TTG2 and TTG3. Two compression trains operate the natural gas liquefaction. Each compression train is composed of two centrifugal compressors (CC1 and CC2), a GT (the prime mover) and a synchronous Motor (M). M acts as starter and helper motor, which allows one to start-up the entire train and provides additional power when required. The motors are supplied by two power conversion stages denoted as PCS1 and PCS2. The power conversion stages are two Thyristor Variable frequency Drives (TVFDs). Each TVFD is connected to the PCC through a step-down transformer with two secondary windings. The two step-down transformers are indicated with TPCS1 and TPCS2 and adapt the voltage level in order to supply the TVFDs. Each TVFD consists of two Line-Commutated-Converters (LCCs). Each LCC is a double-stage converter since the first stage is a line-commutated-rectifier (LCR), while the second stage is a Line-Commutated-Inverter (LCI). The fundamental frequency of the LCRs is the grid frequency fn. The LCIs supply the synchronous motor M, hence the LCIs fundamental frequency is the motor frequency denoted as fm. Each TVFD is based on 6-pulse H-Bridges.
A Harmonic Filter (HF) is connected at the PCC. The HF consists of resonant circuits connected in parallel. The circuits are designed to cut harmonics at 5 fn, 11 fn and 23 fn.
In order to simplify the analysis, the overall loads connected to the power system are taken into account through a lumped load.  As in [1] the TG unit and the TVFD can be described by the following equations: where i i  is the vector of the current perturbations, PCC v  is the vector of the voltage perturbations,  i x is the state-space vector and the matrixes Ai, Bi, Ci and Di are referred to a generic state-space representation.
Combining the state-space models of the TG and of the TVFD it is possible to obtain the overall model and to estimate the damping associated to each TG. Focusing on the LNG plant shown in Figure 1, the overall damping ξ(fi) related to each TG unit is the sum of the electrical damping ξe(fi) and of the shaft-line inherent mechanical damping ξm(fi) [19].
The mechanical damping is mainly determined by the lube oil bearing actions. For the torsional vibrations, it presents an estimated parameter whose value is normally low. The electrical damping is influenced by all the devices included in the electromechanical system (TG units, power conversion stages, HFs and lumped load) as described in [1,19].
The electrical damping assessment provides information about the torsional instability risk level of an LNG plant. Hence, in this paper, the value of ξe(fi) is chosen as the torsional instability risk index. In particular, if the electrical damping has negative value, high risk can be associated to the LNG plant configuration.

TG Units Complete Model
Considering the LNG plant shown in Figure 1, each TG unit includes several inertial masses coupled together via steel shaft sections, special couplings and gears, as discussed in [20]. In [1] a three Degrees of Freedom (DOFs) model is adopted for the TG unit. The first DOF represents the whole gas turbine whose inertia moment is denoted as J1, the second DOF represents the gearbox whose inertia moment is denoted as J2 and the third DOF represents the SG whose inertia moment is denoted as J3. The overall model is characterized by two stiffness coefficients D1 and D2 as shown in Figure 2. In the LNG industry empirical damping assessment is commonly accepted [3]. In the present analysis the mechanical damping for each TG is assumed to be null. The SG model is developed in a d-q reference frame rotating at the grid pulsation ωn as described in [21]. The SG field voltage is assumed constant and the air-gap torque is calculated as function of currents and fluxes. For each TG unit the electromechanical model can be described by [1]: TG  TG  PCC  TG  TG  TG  TG  TG  TG  GT   TG  SG  PCC  TG  TG where δ1, δ2 and δ3 denote the DOFs angular positions and ψfd,, ψkd and ψkq are the rotor fluxes. The model described by Equations (3) and (4) can be modified including the TG unit control loop. The SG excitation circuit is controlled by an Automatic Voltage Regulator (AVR) and a PSS in compliance with [22]. The transfer functions of the AVR and the PSS controllers are shown in Figure  3. The AVR regulates the SG field voltage efd with a PI controller whose proportional gain and time constant are denoted respectively as KAVR and TAVR. Usually, the AVR gain KAVR has high value in order to provide the required air-gap torque when the rotational speed of the rotor ωR deviates from the synchronous value. At very low frequencies (0-2 Hz), the action of the AVR could introduce undamped oscillations, hence the PSS action allows to avoid this instability source. The proportional gain of the PSS is denoted as KPSS. AVR V denotes the signal provided in output by the voltage transducer.  PSS V denotes the signal provided in output by the speed transducer. Differently PSS V denotes the output signal generated by the PSS.  Using the small-signal analysis, the grid voltage magnitude variation can be expressed as: where As a consequence the system defined in Equation (4) can be rearranged as: with The state-space vector and the matrixes, which compose Equation (10), can be expressed as: Independently of the causes, in the presence of SSTI phenomena the generator shaft is led to vibrate. The overall damping of the TG unit coincides with the electrical damping since, as previously declared, the shaft-line inherent mechanical damping is assumed to be null. The TG electrical damping evaluation is based on the matrix ÂTG eigenvalues calculation. Besides the torsional natural frequencies of the shaft-line can be identified. In particular, Equation (10) allows one to detect two modes of torsional vibration for each TG with TNFs equal to 9.2 and 31.5 Hz.
The TG units rated electrical and torsional mechanical parameters are reported respectively in Table 1 and Table 2. 0.14, 0.08, 0.14 pu The sensitivity analysis can be applied to detect the TG control parameters, which impact the torsional stability of the electromechanical system. The Finite Difference Method (FDM) [18] can be adopted to provide the sensitivity analysis of the system shown in Figure 2. FDM implies simple implementation and computational burden proportional to the number of design variables.
Considering the objective function F(P) as a function of a design variable P, its sensitivity coefficient can be approximated from the exact displacement between the initial point P0 and the perturbated point P0 + ΔP, where ΔP is the perturbation of the design variable.
Accuracy can be improved using the central-difference approximation and the sensitivity coefficient dF/dP can be defined as: As discussed in [23] the FDM approximation can lead to an accuracy error, which can be reduced by a proper choice of the perturbation ΔP. In order to reduce the error source, the perturbation ΔP can be assumed high, in particular 50% of the considered design variable.
Considering the data shown in Tables 1 and 2, the sensitivity analysis can be applied to the statespace electromechanical model of the TG units. The analysis is performed assuming rated conditions of the TG unit and starting from the controller parameters provided by the Original Equipment Manufacturers (OEMs) and reported in Table 3. The sensitivity coefficients related to the AVR and PSS control parameters were calculated and the results point out that the AVR and PSS proportional gains (KAVR and KPSS) had more impact than the other design parameters (TAVR, T1PSS and T2PSS). Nevertheless, also the sensitivity coefficients related to the AVR and PSS proportional gains variations were very low as shown in Table 4. Considering the frequency range 5-50 Hz, the torsional stability seemed not influenced by the AVR and the PSS proportional gains variations since the phase of

TVFD Complete Model
The detailed structure of the PCSs shown in Figure 1 is represented in Figure 4. Each TVFD consists of two branches whose DC-links are coupled by the mutual inductance Lm.  In [1] the state-space model of the TVFD first power conversion stage has been presented and based on small-signal linearization. The operation of the power conversion stage can be defined by the following equation.
VFD VFD where ΔT denotes the TVFD commutation period; X5 is the state-space vector when the commutation process is assumed completed, t0 denotes the initial time instant of the first commutation stage; t1 denotes the end of the commutation period; P(θ) is the Park transformation matrix and θ is the Park angle. Further details can be found in [1].
Since the SSTI phenomena are strongly influenced by the TVFD controller parameters [2], a more detailed TVFD model can be developed where the current controllers and the PLLs dynamics can be included. The detailed control scheme of the TVFD and the PLL structure are shown respectively in Figure 5 and Figure 6 where KPI and TPI denote the proportional gain and the time constant of the DC current controllers while KPLL and TPLL denote the proportional gain and the time constant of the PLL PI.
The closed-loop state-space model is described by the following equations system: where the state-space vector ΔXPCS includes the vector of DC currents, the dynamics of the PLLs and the dynamics of the PI current controllers. The inputs of the state-space (24) The PCSs rated parameters and the control parameters provided by the OEMs are reported respectively in Table 5 and Table 6.

LNG Plant Simulation Platform
The theoretical model described in Sections 3 and 4 allows one to obtain electrical damping assessment and stability considerations in case of a basic LNG plant, which consists of a TG unit and a TVFD. In order to manage complex plant configurations, a complete simulation platform has been developed using the software DigSILENT PowerFactory, which allows one to provide detailed simulation results in the time domain. The overall power system shown in Figure 1 has been emulated. In Figure 7 there is represented the scheme of the TG unit integrated model used for the implementation in DigSILENT environment while the TVFD implementation is based on Figure 4.  The TG unit was emulated considering the torsional mechanical model and the control system discussed in Section 3. The SG was modeled through a classic d-q representation provided by the Power Factory software. In Figure 7 the powers related to the TG shaft torque, the air-gap-torque and the driving torque of the TG were indicated respectively with PT, PE and PTUR.
The PCSs and their control system were implemented on the basis of the model discussed in Section 4. The motor M of each compression train was represented as an ideal voltage source connected to the subsynchronous reactance LM and the stator resistance RM.
The TVFDs models and the TG units models were connected in order to obtain a combined electromechanical model. The action of the filter HF was taken into account since it modifies the operative condition of the TG units, changing the amount of active and reactive power supplied by the SGs. The MV transformers were modeled on the basis of their short circuit voltages and the Joule losses. The step-up transformers related to the TG units (TTG1, TTG2 and TTG3) were assumed identical such as the step-down transformers related to the TVFDs (TPCS1 and TPCS2). The transformers parameters are reported in Table 8. The TG shaft-line torque signals provided by the simulation platform are typically affected by time varying components in transient conditions. For this reason the torque signals can be postprocessed through the Prony analysis as already discussed in [1]. In particular, the Prony analysis allows one to identify growing or decaying components of a generic signal f(t) on the basis of the following equation: where Kg is the magnitude, ξg is the damping factor, ωg is the pulsation and g is an integer.
Considering the case study shown in Figure 1, the results provided by the developed LNG plant simulation platform were post-processed through the Prony analysis and successively a proper estimation of the damping factor ξ(fi) was achieved for each TG unit. The stability assessment of the LNG plant and the sensitivity analysis are direct consequences.

Simulation and Experimental Results
A real LNG plant was used as an experimental set-up to carry on an extensive experimental tests campaign. The electrical and mechanical parameters of the experimental set-up were already reported in Tables 1-3 and Tables 5, 6 and 8. The data obtained by the experimental measurements were compared and combined with the simulation results in order to provide a damping assessment of the considered LNG plant. Finally, the simulation platform results were used to perform a sensitivity analysis about the main control systems parameters. The impact of the control system parameters on the TG units overall damping was evaluated and the results could be applied to operate the LNG plant avoiding conditions with high risk of torsional instability.

Experimental Set-Up Configuration
Looking at the LNG plant shown in Figure 1, the three TG units and their control systems were identical. Additionally, the two moto-compressor trains and the related PCSs were assumed equal.
The main components of the experimental set-up are shown in Figures 8 and 9. The speed of the TG units and of the moto-compressor trains was measured by a magnetic pick-up.
In the LNG plant, a Torsional Monitoring System (TMS) was included. It allows one to provide as an output the angular oscillation of the turbine toothed wheel at the point of measurement [24]. The resulting alternating torque was calculated on the basis of the detailed torsional model of the shaft-line.
The SGs stator voltages and currents were measured respectively by means of voltage transformers (VTs) and current transformers (CTs). Probes installed inside the AVR system provide the voltage and the current measurements of the SGs field circuits (Figure 8).
With reference to the PCSs and the moto-compressor trains, the following measurements were available: DC-link's currents, LCR's firing angles, LCI's firing angles, LCR's AC voltages and motor's voltages (Figure 9).  Considering the real LNG plant, the torsional vibration measurements on the TG unit were not available during the first phases of the set-up commissioning. The TMS was installed after the occurrence of the first SSTI phenomena. However, the torsional oscillations of the TG units can be evaluated also by the analysis of the DC-link's currents. Indeed, the currents measured at the DClinks of the TVFDs can be used to estimate the overall damping ξ(TNF1) related to the first TNF. In particular the spectra of DC-link currents can be obtained using the Continuous Wavelet Transform (CWT) [25]. Then, the component related to the first TNF of the TG can be extracted.
The magnitude of the component at the first TNF can be also approximated through the function 1 DC i (TNF ) expressed as: where ATNF1 is the initial value of the current oscillation. Further details can be found in Appendix B.
Hence, the overall damping ξ(TNF1) is estimated comparing 1 DC i (TNF ) (Equation (26)) with the DC-link current component at the first TNF provided by the spectrum.
Finally the derivative-free optimization method [26] can be used to identify the value of ξ(TNF1), which minimizes the difference between the component extracted by the spectrum and the approximation defined by Equation (26).

Damping Assessment
Considering the experimental set-up shown in Figures 8 and 9 and the LNG plant simulation platform discussed in Section 5, four different Configurations (Cs) are analyzed and denoted as CA, CB, CC and CD. The Cs mainly differ in the number of TG units connected to PCC, the power of the lumped load and the PCSs number. The main data related to the four configurations are reported in Table 9. In the proposed analysis the mechanical damping related to the four Cs was estimated subtracting the electrical damping calculated by the simulations results from the overall damping calculated by the experimental results.
The spectrum of the DC-link current DC i was used to extract the component related to the first TNF as discussed previously. Figures 10-12 show the spectra of DC i in the frequency range (0-100 Hz) related to the configurations CA, CB and CD. As expected, the DC i spectra show that, apart from the direct component, the component at the first TNF exhibited the highest value and it grows with the time.
In Figures 13-15 the components measured at the first TNF were compared with the approximate curves obtained by Equation (26) and represented as dotted lines. A very good match between the two curves was verified in configuration CA ( Figure 13). This provides also information about the proper calculation of the overall damping. The smaller overlap was verified in the case of configurations CB and CD (Figures 14 and 15). This is due to the presence of more TG units, which have the same frequency but operate not exactly in phase. However, in all the examined configurations, Equation (26) guarantees satisfying approximation of the DC i component at the TNF.
As a consequence, also the overall damping related to the TG units was estimated accurately.
The results related to the configuration denoted as CC were not reported since, in this case, the first TNF component of the current DC i exhibited a very low value.       Figures 16-19 show the shaft-line torques related to the first TG unit in configurations CA, CB, CC and CD respectively. The results were obtained by the simulation platform and they are postprocessed through the Prony analysis. The simulations are performed assuming zero mechanical damping and stressing the shaft line of one TG unit with a torque impulse. In this hypothesis the overall damping was calculated using Equation (25) and the overall damping was assumed equal to the electrical component. The fine overlap between the simulation results and Prony curves indicates a good approximation of the electrical damping calculation.    Starting from the overall damping estimated by the experimental measurements and through Equation (26), the mechanical damping related to the four configurations could be evaluated subtracting the electrical damping from the overall damping. The results are summarized in Table  10. Considering the power range of the TG units, [3] suggests that possible values for the mechanical damping are in the range (0.0025-0.01) pu. As a consequence the data reported in Table 10 are in the set range. As expected, the data reported in Table 10 highlight the variability of the mechanical damping related to the same TG unit in the four considered configurations.
The electrical damping is the index of the torsional instability risk. Hence, it has to be pointed out that a high risk level is detected in the four configurations. Just in configuration CC the overall damping is positive and the torsional instability is avoided absolutely.

Sensitivity Analysis
In Sections 3 and 4 the control systems parameters, which impact the SSTI phenomena were detected by means of the sensitivity analysis. In Section 4 it was concluded that the parameters KPI and KPLL are elements with high sensitivity coefficients. In Section 3, it was shown that the torsional stability seems not to be influenced by the AVR and the PSS proportional gains variations.
In this subsection the sensitivity analysis, based on the FDM method, was applied to the results provided by the simulation platform presented in Section 5. The analysis was performed starting from the controller parameters provided by the OEMs and reported in Table 3 and Table 6. Finally, the configuration CC was assessed as the most stiff since the highest electrical damping was verified. Indeed, as summarized in Table 9, three TGs units supply the whole power system in the configurations CC and each TG unit provides around 20% of its rated power.  show the sensitivity coefficients related to the four configurations. In particular, in Figure 24 there is shown how, increasing the value of the parameter KPI, the related sensitivity coefficient tended to be zero. This trend confirms what can be observed in Figure 20 where the electrical damping was almost constant for high values of KPI.
Comparing the sensitivity coefficients of Figures 24 and 25 with the coefficients shown in Figures  26 and 27, it can be noticed that the coefficients related to KPI and KPLL variations were higher than the coefficients associated to KAVR and KPSS variations. This demonstrated that variations of KPI and KPLL affected the LNG plant stability, while the parameters KAVR and KPSS had limited influence. In conclusion, the preliminary analysis developed in Sections 3 and 4 was confirmed and there was demonstrated the good match between the theoretical model and the simulation platform results. Figures 24-27 show also that, despite the different electrical damping values, in the four configurations the sensitivity coefficients were similar and they had a similar trend. This provides hints about the actions to reduce the risk level in all the four configurations. For example, decreasing the values of KPI and KPLL an electrical damping ξe(TNF1) increment was verified.
The torsional instability risk level related to the configurations CB, CC and CD could be reduced through a proper tuning of the parameters KPI and KPLL (Figures 20 and 21). Differently, for the configuration CA, acting on the control parameters did not imply a reduction of the risk level.

Discussion
The electrical damping is the index of the torsional instability risk. Considering the results reported in Section 6, the high risk level was detected in all the four analyzed configurations. Hence, it was evident that torsional instability could be verified also when the LNG plant was operated with the electrical and control parameters provided by the OEMs. In conclusion, the overall damping was positive and the torsional instability was avoided absolutely just in configuration CC. However, the risk level related to the configurations CB, CC and CD could be reduced through a proper tuning of the parameters KPI and KPLL. While, for the configuration CA, the control system parameters tuning was ineffective. For this configuration, dedicated equipment should be installed in the plant in order to increase the electrical damping of the TG units and to avoid torsional instability. The equipment could be based on a dedicated power converter providing supplementary subsynchronous damping. This solution will be investigated in further publications.
Besides, it has to be pointed out that the theoretical model provided in this paper and the developed simulation platform could be directly applied also to different island-operated power systems such as the combined cycle power plants or to the highly electrified ships since these power systems exhibited similar configurations, power levels and power conversion stages.

Conclusions
Due to the complexity and the numerous devices that compose an LNG plant, torsional instability phenomena can occur also in the case that the plant is operated considering rated electrical, mechanical and control parameters provided by the Original Equipment Manufacturers (OEMs). For this reason a comprehensive approach to the LNG plant stability analysis is required.
In the paper an improved theoretical model was presented to estimate accurately the electrical damping of an LNG plant and sensitivity analysis was performed to assess the impact of the control systems parameters on the risk of torsional instability. On the basis of the theoretical model, a complete simulation platform was developed in order to manage complex LNG plant configurations with numerous TG units and drives. Finally, an extensive tests campaign was carried out on a real LNG plant. Both the theoretical model and the simulation platform provide proper evaluation of the real electrical damping and, as a consequence, of the risk of torsional instability. The results of the sensitivity analysis demonstrated how the tuning of the control systems parameters affected the electrical damping of the LNG plant. Hence fine tuning of the control parameters should be adopted in the LNG plants practice to operate the power systems avoiding conditions with high risk of torsional instability or reducing the risk of torsional instability.
The developed simulation platform represents a valuable and flexible tool to extend the results of this study to different plants. The platform can be arranged to manage other LNG plants with different power ranges, a different kind of loads and a variable number of TG units considering also the occurrence of contingencies and real time variations.
where ψad and ψaq are the flux linkages.
In the small-signal model the variation of rotor circuit fluxes can be expressed as: where ψfd denotes the field winding flux and ψkd and ψkq denote respectively the d-q axes amortisseur winding flux linkages. The fluxes ψad and ψaq are defined as: The SG currents can be defined as: where R  R  R  R  R  R  R  TG  I  I  I  I  I  I  I   I  I  I  I  I  I  I  C  I  I  I  I  I  I R  TG  I  I  I  I   I  I  I  I  D  I  I  I  .
More details about the coefficients used in Equation (A5) can be found in [1]. Considering the Equations (A2-A5) and the torsional model of the TG, the following integrated state-space can be obtained:        0  3  0  3  0  0  0  0  3   0  3  0  3  0  0  0  0  3   3  3  3  0  0  0   3  where q is the modal variable, ωTNF1 is the pulsation relative at the first TNF and KT is the externally applied modal torques. Equation (A8) allows one to identify the damping  related to the first TNF.
The torsional oscillation is transferred to the DC link current, hence the magnitude of the current oscillation at the first TNF can be approximated as 1 DC i (TNF ) : where ATNF1 is the initial magnitude of the oscillation.