8/9/2019 PSSE Model User Manual wind Vestas http://slidepdf.com/reader/full/psse-model-user-manual-wind-vestas 1/53 Vestas Wind Systems A/S Alsvej 21 DK-8940 Randers SVWWW VEST S COM USER MANUAL PSS/E Model for Vestas GridStreamer TM Wind Turbines Version 8.0 : s o cu me nt c on ta ns vaua e con ent a n ormat on o e st as n y st em s / . t s protecte y c op yr g t aw as an unpu s e wor . e st as r es er ve s a p at en t, c op yr g t, trade secret, and other proprietary rights to it. The information in this document may not be used, reproduced, or disclosed except if and to the extent rights are expressly granted by Vestas in writing and subject to applicable conditions. Vestas disclaims all warranties except as expressly granted by written agreement and is not responsible for unauthorized uses, for which it may pursue legal remedies against responsible parties.
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
USER MANUALPSS/E Model for Vestas GridStreamer TM Wind Turbines
Version 8.0
: s ocument conta ns vaua e con ent a n ormat on o estas n ystems / . t s protecte y copyr g t aw as an unpu s e wor . estas reserves a patent, copyr g t ,trade secret, and other proprietary rights to it. The information in this document may not be used, reproduced, or disclosed except if and to the extent rights are expressly granted by Vestas in writing and subject to applicableconditions. Vestas disclaims all warranties except as expressly granted by written agreement and is not responsible for unauthorized uses, for which it may pursue legal remedies against responsible parties.
VESTAS MAKES NO WARRANTY OR REPRESENTATION EITHER EXPRESS OR IMPLIED IN RESPECTOF THE PSS/E MODEL, INCLUDING WITHOUT LIMITATION AS TO ACCURACY, COMPLETENESS,FUNCTIONALILTY, PRECISION, USEFULNESS, FITNESS FOR A PARTICULAR PURPOSE OF THE PSS/EMODEL OR OTHERWISE. THE PSS/E MODEL IS PROVIDED “AS IS” AND VESTAS SHALL HAVE NORESPONSIBILITY OR LIABILITY WHATSOEVER FOR RESULTS OF USE OR PERFORMANCE OF THEPSS/E MODEL. TO THE MAXIMUM EXTENT PERMITTED BY LAW, IN NO EVENT VESTAS SHALL BELIABLE FOR ANY CONSEQUENTIAL DAMAGES, DIRECT, INDIRECT, SPECIAL, PUNITIVE OR OTHERDAMAGES WHATSOEVER ARISING OUT OF OR IN ANY WAY RELATED TO THE USE OF OR INABLITYTO USE THE PSS/E MODEL, WHETHER BASED IN CONTRACT, TORT, NEGLIGENCE, STRICT LIABILITY
OR OTHERWISE.
For the avoidance of doubt, Vestas makes no warranty or representation either express or implied as tothe performance of the wind turbine model in terms of it being in accordance with the performance ofthe actual wind turbine generator, as other circumstances, including, but not limited to deviations in themarkets and optional features might have influence on the performance of the actual wind turbine gen-erator. The performance of the wind turbine model is expected only to be indicative to the performanceof the actual wind turbine generator.
Copyright Notice The documents are created by Vestas Wind Systems A/S and contain copyrighted material,trademarks, and other proprietary information. All rights reserved. No part of the documentsmay be reproduced or copied in any form or by any means—such as graphic, electronic, or
mechanical, including photocopying, taping, or information storage and retrieval systems—without the prior written permission of Vestas Wind Systems A/S, and its respective parentcompanies, subsidiaries, affiliates, successors, assigns, licensees, representatives andagents (together "Vestas"). The use of these documents by you, or anyone else authorizedby you, is prohibited unless specifically permitted by Vestas. You may not alter or removeany trademark, copyright or other notice from the documents.
00 2013-04-01 This document is based on a previous document with docu-ment number 0026-5729_V03.Updated document to reflect changes of Pre-sales model ofGridStreamer TM 3.3MW Turbines.Only parameters of GridStreamer TM 3.3MW Turbines areavailable in this version.
01 2013-04-05 Parameters of V164 8.0MW Mk0 concept model are included
in this version.02 2013-07-04 Parameters of V80/V90 2.0MW Mk8 model are included in this
version
03 2013-09-06 Updated document to reflect changes of As-built model ofGridStreamer TM Turbines.Only parameters of 3 MW Turbines are available in this ver-sion.
The description given in this document covers the design of the Vestas PSS/E wind turbine
user written model. The model is applicable to all Vestas GridStreamer TM wind turbines.
This document is a user manual that explains how to implement the Vestas GridStreamer TM wind turbines in an existing grid and states a number of steps that must be performed to geta successful simulation.
The document is a general overview of the PSS/E user written model version 8.0 for VestasGridStreamer TM wind turbines. The model is developed for the purpose of dynamic stabilityanalysis and therefore generally omits effects with time constants shorter than one AC cycle. A PSS/E time step (DELT) in the range 1ms to ½ cycle is recommended.
Vestas PSS/E Model Package
This user manual is part of a model package including the following files:
PSS/E Model for Vestas GridStreamer TM Wind Turbines Version 8.0 User Manual
A library file (.LIB) that includes all the sub models in the wind turbine is compiled for PSS/Eversion 29, 30, 31 and 32.
VestasGS_8_0_n_PSSE<m>.lib – Model library file for PSS/E ver. <m>
A Dynamic data file (.DYR).
V<xx>.GS.<yy>MW.<zz>Hz. Mk<k>.Df.v8.0.n.dyr – Dynamic Data file (.dyr)
This PSS/E user written model has been developed to support the Vestas GridStreamer TM wind turbines topology. This is explained – in short - in the following sections.
2.1 Vestas GridStreamer TM ConceptIn the full-scale converter system the generator is decoupled from the grid-side convertercompletely by the DC-link. In this case the converter will be in total control of the powerand the GridStreamer ™ WTG is coupled via a step-up transformer to the connection point.
2.2 Model Structure
The PSS/E user written model developed for the Vestas GridStreamer ™ wind turbine is asimplified ‘dynamic’ model of the WTG and omits some detail of mechanical componentsand the generator-side converter. However, this model was shown to accurately representthe complete GridStreamer™ turbine, because the mechanics and generator -side con-verter have very little impact on the grid-side output due to decoupling through the DC-link.
Figure 1 shows the generic model structure for Vestas GridStreamer TM wind turbines.Here the various signals and blocks may be observed as well as input and output of themodel. Note that the grey blocks are not part of the PSS/E model, and have directfeedthrough signals.
Figure 1 Illustration of the generic PSS/E model structure. Blocks marked as grey are not a part ofthe model
As mentioned in the Introduction, this model is intended for grid dynamic stability analysis.It may also be designated as a performance model, since it receives grid voltage and fre-quency and delivers active and reactive power as response.
2.3 Model FeaturesIt is important to decide the types of analysis and studies needed to be carried out beforedeciding which approach to use in modelling the physical system. The modelling approachshould reflect the intended types of analysis and studies.
These questions led to shape the modelling requirements and to give a better overview ofthe user’s needs in order to develop an analysis-dedicated model. The modelling of thewind turbine needs to be made according to the performance of the simulation tool, whichis the PSS/E.
The pre-modelling analysis has been done for using PSS/E simulation tool and the out-come is listed in Table 1.
Table 1 Outcome of the pre-modelling analysis
Model purpose Dynamic stability studies – 1s to 5 min.
Implementation in large bus system.
Aggregated Wind-power plant performance.
Model operational range Full power range of the WTG.
PF control mode. Reactive power control mode.
Full voltage and frequency range of WTG.
Studies Dynamic stability analysis.
Protection Voltage settings.
Frequency settings.
LVRT settings.
Model requirements Positive sequence based. Model bandwidth 0-10 Hz.
This section explains how to add Vestas Wind Turbines Model and modify the user’s exist-ing system simulation files. With the package files provided by Vestas Wind Systems A/S,the end users should modify their load flow data (RAW or SAV file) as well as their dynam-ic data (DYR file). The suggested modifications are stated below.
3.1 Load Flow DataThe PSS/E user should begin by preparing a solved load flow case which includes a gen-erator and unit transformer for each WTG or WTG aggregate.
-The WTG’s Generator :
The wind power plant units can be included as separate units or as a single aggregatedgenerator, with appropriate steady-state real power (P) and reactive power (Q) for eachwind power plant (turbine). The relevant WTG generator should be implemented by theuser according to the parameters listed in Appendices and the set up instructions in Chap-ter 6.
NOTE : In newer versions of PSS/E a generator is added to the load flow case as either a‘classical machine’ or a ‘wind machine’. The Vestas wind turbines should be added as a
classical machine not a wind machine, as the model is designed to be back compatible toPSS/E version 29.
- The WTG’s Transformer :
In addition to the substation transformer, a step-up unit transformer should also be includ-ed in the load flow case by the user. The delivered model package does not include theWTG transformer and therefore this transformer should be implemented by the user ac-cording to the parameters of the transformer as given in Appendices and the set up in-structions provided in Chapter 6.
The user may refer to the raw file template in Template 1.
/*add WTG transformer data here */****************************************************************
0 / END OF TRANSFORMER DATA, BEGIN AREA DATA
0 / END OF AREA DATA, BEGIN TWO-TERMINAL DC DATA
0 / END OF TWO-TERMINAL DC DATA, BEGIN VSC DC LINE DATA
0 / END OF VSC DC LINE DATA, BEGIN SWITCHED SHUNT DATA
0 / END OF SWITCHED SHUNT DATA, BEGIN IMPEDANCE CORRECTION DATA
0 / END OF IMPEDANCE CORRECTION DATA, BEGIN MULTI-TERMINAL DC DATA
0 / END OF MULTI-TERMINAL DC DATA, BEGIN MULTI-SECTION LINE DATA
0 / END OF MULTI-SECTION LINE DATA, BEGIN ZONE DATA
0 / END OF ZONE DATA, BEGIN INTER-AREA TRANSFER DATA
0 / END OF INTER-AREA TRANSFER DATA, BEGIN OWNER DATA
0 / END OF OWNER DATA, BEGIN FACTS DEVICE DATA
0 / END OF FACTS DEVICE DATA Template 1: The PSS/E raw file template
3.2 Short Circuit DataUnder fault conditions the behaviour of the WTG is determined by the converter control.During a fault the converter controls will limit the magnitude of the fault current contributionto 1.05 p.u. (Ik”) and 1.45 p.u. (peak) on WTG rated current.
The WTG should be modelled as a current source for short circuit studies. The currentsource should contribute a maximum of 1.05 p.u. (Ik”) fault current.
In PSS/E the value of Zsource is obtained from the maximum short circuit current contri-bution (Ik”=1.05 p.u.). The recommended values are listed in Appendices and the set upinstructions provided in Chapter 6.
3.3 Dynamic DataThe user should prepare a DYR file (.dyr) containing the dynamics data. This will generallyinclude dynamic data records for other power system plant(s), as well as the dynamic datafor Vestas WTG model. In this package the user will find a parameter template(V<xx>.GS.<yy>MW.<zz>Hz.Mk<k>.Df.v8.0.n.dyr) as shown in Template 2. The template
The dynamic data file (DYR file) Template 2 supplied requires user modification, to insertthe details specific to each WTG. If there is more than one WTG in the PSS/E simulation,the user should duplicate the DYR file contents for each WTG added, and to modify each
of the user inputs denoted by the “<>” to reflect the particulars of each WTG. User inputparameters are listed in Table 2. The details that must be entered by the user for eachWTG are:
The PSS/E bus number (<bus>) and machine ID (<mach>) corresponding to thesingle or aggregate WTG generator.
The WTG reactive power control mode (<WTGmode>): 0 for direct reactive power(Q) control, 1 for power factor (PF) control.
Table 2 User input parameters to the dynamic data file (dyr file)
Parameter Description
<bus> The generator bus number
<mach> The machine ID
<WTGmode> WTG control mode: 0 for Q control, 1 for PF control
In case several WTGs are to be aggregated into one scaled WTG the set of data shouldnot be duplicated, this type is handled through the WTG load flow inputs as explained inChapter 6, where the Pgen, Pmax, Qgen, Qmax, Qmin, & Mbase, will be scaled to repre-sent the required amount of aggregated WTGs as per the instructions provided.
The model is developed for a maximum simulation time step of half a cycle. The model isdeveloped for PSS/E versions from 29 upwards.
The wind turbine(s) impacts on the system (grid) will be dependent on the short circuit ratio(SCR) and X/R at the bus to which the wind turbine(s) is (are) connected. In this sense it isvery important to assess the model performance taking into consideration the grid strengthseen from the connection point.
Section 4.1 gives the recommended settings for the dynamic simulation. In Section 4.2,the simulation procedure is given (dependent on the version of PSS/E). The list should befollowed exactly in order to get a successful simulation.
1. Create a working directory – where files related to this simulation set-up shallbe placed – including all relevant files listed.
2. Start up a command prompt in the working directory, ensuring the PATH isset correctly to be able to run PSS/E programs.
3. Start-up PSS/E dynamics (PSSDS4) in the working directory.4. Enter PSS/E load flow (LOFL).5. In case the user’s preferred method of representing the load flow case is a
raw file then read the load flow raw file (READ) to which the Vestas wind tur-bines have been added. If the user prefers to use an already created savedcase then open (CASE) the load flow case to which the Vestas wind turbineswill be added.
6. In case the user has chosen to use a Saved file (*.sav) then the following maybe necessary:
i. Modify this load flow case by adding the desired number of new buses,feeders and transformers necessary to model the WPP collector sys-tem up to the point of common coupling with the existing power sys-tem.
ii. Add the WPP aggregated generator(s), specifying P and Q output foreach.
iii. Add the wind turbine transformers as per the implemented structure:1. Aggregated transformer for aggregated WTGs (if applicable).
2. Separate transformer for separate WTGs (if applicable). 7. Solve the load flow case and save the solved case file.8. Convert the load flow case for dynamic simulations by using CONG and
CONL, followed by ORDR, FACT and TYSL.9. Make a copy of the relevant DYRE parameter template, and fill in the details of
each WPP. Add the DYRE records to any records for other system plant tocreate a system wide ‘.dyr’ file.
10. Enter PSS/E dynamics (RTRN).11. Read (DYRE) the new ‘.dyr’ file, and specify CONEC.FLX and CONET.FLX
files. Specify the name of the compile-file.12. Save the snapshot.
13. Exit the PSS/E program (STOP).14. Ensure the command prompt is pointing to the correct working directory.15. Compile the CONEC.FLX and CONET.FLX subroutines – by typing the name
of the compile-file.16. Link the object codes using a command line of the following form:
CLOAD4 VestasGS_8_0_n_PSSE29.lib <other user models>17. You are now ready to run simulations: enter PSS/E dynamics (PSSDS4), open
1. Create a working directory – where files related to this simulation set-up shall
be placed – including all relevant files listed.2. Ensure the system PATH is set correctly to be able to run PSS/E programs
(using the PSS/E Environment Manager if necessary).3. Start-up PSS/E in the working directory.4. In case the user’s preferred method of representing the load flow case is a
raw file then read the load flow raw file (READ) to which the Vestas wind tur-bines have been added. If the user prefers to use an already created savedcase then open (CASE) the load flow case to which the Vestas wind turbineswill be added.
5. In case the user has chosen to use a Saved file (.sav) then the following maybe necessary:
i. Modify this load flow case by adding the desired number of new buses,feeders and transformers necessary to model the WPP collector sys-tem up to the point of common coupling with the existing power sys-tem.
ii. Add the WPP aggregated generator(s), specifying P and Q output foreach. Ensure the g enerator(s) spec ify ‘not a wind machine’.
iii. Add the wind turbine transformers:1. Aggregated transformer for aggregated WTGs (if applicable).2. Separate transformer for separate WTGs (if applicable).
6. Solve the load flow case and save the solved case file.7. Convert the load flow case for dynamic simulations by using CONG and
CONL, followed by ORDR, FACT and TYSL.8. Make a copy of the relevant DYRE parameter template, and fill in the details of
each WPP. Add the DYRE records to any records for other system plant tocreate a system wide ‘.dyr’ file.
9. Read (DYRE) the new ‘.dyr’ file, and specify CONEC.FLX and CONET.FLXfiles, and a name for the compile-file as appropriate.
10. Save the snapshot.11. Use the ‘Create User DLL’ function in the Dynamics menu to compile the
CONEC and CONET subroutines and link the model code into PSS/E. Specifythe library ‘VestasGS_8_0_n _PSSE32.lib’ when r equired.
12. You are now ready to run simulations: enter PSS/E dynamics (PSSDS4), openthe snapshot and load flow case.
Note : It is not necessary to specify the WPP real or reactive power output in the *.dyr fileas these will be taken from the load flow case at initialization (i.e. when using the activitySTRT).
Once the snapshot is taken and the model compiled and linked, the same dynamics snap-shot can be used with multiple load flow cases, provided the network topology and con-nected models did not change from one case to another. It is therefore possible to assem-ble a number of load flow cases, each with a different level of initial WTG power outputand/or voltage, without reassembling the snapshot or recompiling before using each differ-
5.1 System ModelVestas has built a simple example to show how to setup and simulate with the WTG mod-els. The test grid is modelled with four buses and two branches, with the Vestas WT modelat one end, and a slack bus generator at the other. The middle bus is where the fault isapplied. The system is shown in Figure 3.
Figure 3 Network Diagram of Vestas test system
The files/folders listed in Table 4 are delivered for setting up the test cases for both, thereactive power control and the power factor control examples (for WTG level control).
Note that the example files provided represent a simple test case for Vestas GridStream-er TM wind turbines and do not reflect any particular site-specific configuration.
Table 4 Files included in the delivered package
Simulation Files DescriptionV<xx>.GS.<yy>MW.<zz>Hz.Mk<k>.Df.v8.0.n_Base_29.raw RAW file for Base case in Ver29
V<xx>.GS.<yy>MW.<zz>Hz. Mk<k>.Df.v8.0.n_Base_32.sav SAV file for Base case in Ver32
Simulate_QControl.psa PSA file for QContorl Simulation
Simulate_PFControl.psa PSA file for PFContorl Simulation
V<xx>.GS.<yy>MW.<zz>Hz. Mk<k>.Df.v8.0.n.dyr DYR file for Dynamic model
TG Z 20 to 30 Z 30 to 40
Wind Generator Bus Bus # 10 Bus Volt = 0.65 kV Bus Type ”2"
Wind turbine Transformer Bus Bus # 20 Bus Volt = 20.0 kV
5.2 Simulation StepsTo perform a simulation – following steps has to be carried out:
In PSS/E versions 29 to 31
1. Create a Working directory – where files related to this simulation set-up shallbe placed – including all relevant files listed.
2. Start up a command prompt in the working directory, ensuring the PATH isset correctly to be able to run PSS/E programs.
3. Start up PSS/E dynamics (PSSDS4) in the working directory.4. Read the load flow raw file “*.raw”, convert the load flow case for dynamic
simulations by using CONG and CONL, followed by ORDR, FACT and TYSL.5. Save the converted case file.6. Read (DYRE) the ‘*.dyr’ file, and specify CONEC.FLX and CONET.FLX files.
Specify the name of the compile-file.7. Save the snapshot.8. Exit the PSS/E program (STOP).9. Ensure the command prompt is pointing to the correct working directory.10. Compile the CONEC.FLX and CONET.FLX subroutines – by typing the name
of the compile-file.11. Link the object codes using a command line of the form:
CLOAD4 VestasGS_8_0_n_PSSE29.lib <other user models>12. Start up PSS/E dynamics (PSSDS4) and run “Simulate_QControl.psa” or
“Simulate_PFControl.psa”
13. Exit PSS/E14. Start plot program (PSSPLT)
In PSS/E version 32 or above
1. Create a Working directory – where files related to this simulation set-up shallbe placed – including all relevant files listed.
2. Ensure the system PATH is set correctly to be able to run PSS/E programs, byrunning the PSS/E Environment Manager if necessary.
3. Start up PSS/E in the working directory.4. Read the case file “*.sav”, convert the load flow case for dynamic simulations
by using CONG and CONL, followed by ORDR, FACT and TYSL.5. Save the converted case file.6. Read (DYRE) the ‘*.dyr’ file, and specify CONEC.FLX and CONET.FLX files.
Specify the name of the compile-file.7. Save the snapshot.8. Use the ‘Create User DLL’ function in the Dynamics menu to compile the
CONEC and CONET subroutines and link the model code. Specify the codelibrary “VestasGS_8_0_n _PSSE32.lib” when prompted.
9. Run “Simulate.idv” 10. Exit PSS/E11. Start plot program (PSSPLT)
Figure 4 shows the PSS/E plant and machine editor. The WTG is to be set on a bus hav-ing a voltage equal to the generator rated voltage (at the low voltage side of the transform-er).
The machine parameters in PSS/E load flow should be entered as indicated in the tablebelow. Note that for the purpose of the PSS/E WTG model, the MBASE value for a singleturbine is taken as the maximum MW output of the turbine.
Table 5 Variables as in plant & machine editor of Figure 4
# Name Units Description RecommendedValue for a Single
Turbine
RecommendedValue for an Ag-gregated Turbine
1 Pgen MW Is the generated active power User supplied value(0≤ P≤ Prate)
N*P
2 Pmax MW Is the maximum allowed activepower
Prate N*Prate
3 Pmin MW Is the minimum allowed activepower
Always ZERO 0
4 Qgen MVAR Is the generated / absorbed reac-tive power
User supplied value*(Qind≤ Q≤ Qcap)
N*Q
5 Qmax MVAR Is the absolute maximum allowedreactive power
Equals Qgen* N*Q
6 Qmin MVAR Is the absolute minimum allowedreactive power
Equals Qgen* N*Q
7 Mbase MVA Is the base apparent power for n Prate N*Prate
8 ZSource pu Is machine positive sequence im-pedance
Refer to Appendices Refer to Appendices
9 Rtran pu Is the Resistance of the internal
transformer
Always ZERO 0
10 Xtran pu Is the Reactance of the internaltransformer
Always ZERO 0
Where Prate is the rated MW value; Qind and Qcap are the turbine reactive power inductive and capacitive lim-its.
*Vestas WTGs are not configured for terminal voltage control. Therefore the user should select the appropri-ate reactive power condition and configure the load flow case accordingly.
In PSS/E the value of ‘Zsource’ is calculated from the maximum short circuit current con-tribution (Ik”=1.05 p.u.). Recommended values of source resistance and reactance aregiven in Appendices.
The bus to which the wind turbine(s) is (are) connected should have the same level ofvoltage as the generator terminals, in fact this bus is the low voltage side of the 2-windingstur bine’s transformer. The turbine’s transformer is not included in the model and is left tothe user to implement in the simulation as a PSS/E standard 2 winding transformer formore details see section 6.2.
For aggregated WTGs, parameters 1, 2, 4, 5, 6 & 7 needs to be scaled to represent thedesired number N of aggregated units.
6.2 WTG 2-Winding Transformer Load Flow InputsFigure 5 shows the PSS/E load flow 2-winding transformer editor, all data are to remain asdefault except the 9 encircled fields that are to be changed as specified in Table 6.
Figure 5 The PSS/E load flow – 2-winding transformer editor
Table 6 Input values for the 2-winding transformer editor of Figure 5
# Name Units Description
1 Name The transformer name (Optional)2 Impedance data I/O
codeMVA Base to use for the p.u. val-ues
3 Line R pu Winding resistance
4 Line X pu Leakage reactance
5 Winding MVA base MVA Winding MVA base
6 Rmax Maximum Tap
7 Rmin Minimum Tap
8 Tap Positions No of Taps
9 Connection Code Vector Grouping
** Vector grouping assumes that the LV side is declared as the primary winding.
It is important to be aware that aggregating WTGs also means aggregating WTG unittransformers, which means that voltage ratings remain the same but the MVA rating will bescaled.
The different transformer parameters are as per the data sheet in Appendices.
6.3 Output Channels for PlottingOne of the possible ways to call a channel in PSS/E is through the CHAN activity in theuser interface (Dynamics). Table 7 shows the possible useful channels that can be calledto represent the relevant inputs/outputs of the Vestas wind turbine model in PSS/E.
This section describes different sub models in the dynamic data file.
7.1 Reactive Power ControlThe turbine can operate within certain limits for active and reactive power. The limits areprimarily determined by the converter and the generator. There are two possible controlmodes:
Q-Control.
Power Factor Control.
The PSS/E user written model will independent of which of the modes that is selected, pri-oritize the active power higher than the reactive power. This means that the highest possi-
ble active power – which is the user-defined (requested) active power – will be supplied tothe grid even though it might compromise the control mode performance.
Active Power PriorityThe model implemented always gives highest priority to the active power. That means thatif both P and Q are outside the valid operating area, Q is limited in order to bring the tur-bine inside the safe operating area. This is illustrated in Figure 7 (A typical PQ chart forVestas turbines) where the references (P* and Q*) are both outside the chart. In this caseQ is limited in order to bring the turbines operating point inside the operating area again.
Q
P
(P*,Q*)
Q*Lim
P*
P over Q
Q over P
Figure 7 Example of limiting function where P is given priority over Q
The value of P* and Q* are set in the load flow interface.
Q-Control A certain reactive power reference can be set. The reference must be within its lower andupper limits – vertical limits for each P according to Figure 11.
When activated, the required Q will be calculated based on the P and the PF set point.Depending on the working point, whether it is a full or partial load operating point and alsowhether it is lead or lag, the PF set point is subject to possible minimum or maximum as
shown in Figure 11.
Changing Q and PF set point (Qset)Reactive power set point (Qset), defined in ‘GSVARS’ sub-routine. This set point is as-signed to VAR L+1 can be changed to desired value using ‘ALTR’ activity.
Table 8 Description of VAR set points P and Q or PF
VARs Description
L Available real power, Pset
L+1 Q or PF set point, Qset
7.2 LVRT Control SettingsDuring low-voltage-ride-through (LVRT), Vestas turbines switch from power control to cur-rent control mode. The following settings can be changed by the user in order to carry ondifferent system studies.
Active current Priority As default, Vestas turbines are designed to maintain and control reactive current duringLVRT. However the user could also select the active current priority, If supplying activepower is required during LVRT.
The Active current priority setting can be enabled by changing CON J+7 of ’GSLVRT’ to1.0, and be disabled by changing CON J+7 to 0.0.
Q-offset SettingIn case the turbine reactive output is not zero at pre-fault condition, the Q-offset can beused to remove the step change when Vestas turbines switches between power controland current control.
The Q-offset setting can be enabled by changing CON J+31 of ’GSLVRT’ to 1.0, and be
disabled by changing CON J+31 to 0.0.
Asymmetrical fault
If asymmetrical fault study is required, asymmetrical fault settings should be enabled andthe user should provide negative sequence magnitude (U-) input in VAR L+12 of’GSLVRT’.
The asymmetrical fault setting can be enabled by changing CON J+32 of ’GSLVRT’ to 1.0,and be disabled by changing CON J+32 to 0.0.
Figure 8 Voltage protection limits defined in the protection model VWVPR6
The WTG will be disconnected if the voltage limits are exceeded and the protection timesout. Table 10 shows the protection settings for a Vestas GridStreamer TM wind turbine with-out the LVRT option.
Table 10 Voltage protection settings (LVRT not enabled)
Voltage limit Setting Timeout Setting
UV1 0.90 tUV1 60 s
UV2 0.85 tUV2 11 s
UV3 0.85 tUV3 11 s
OV1 1.10 tOV1 3600 s
OV2 1.21 tOV2 2 s
OV3 1.36 tOV3 150 ms
OV4 1.36 tOV4 150 ms
OV5 1.36 tOV5 150 ms
Note : the voltage protection settings may vary with different turbine models. The user issuggested to check the General Specification of the specified turbine for the voltage pro-
tection settings.
7.3.2 LVRT Logic
As default, the Vestas GridStreamer TM wind turbines is a LVRT enabled turbine, and theVestas PSS/E model will represent this. The turbine is designed to stay connected on thegrid in case of a number of severe faults.
The LVRT option is affected by two user defined parameters (LVRT_flag and AGO_enable) as explained in Table 11.
I+2 AGO_enable GSLVRT AGO enableLVRT mode can be enabled (set both “LVRT_flag” and “AGO_enable” to 1) or disabled(set both “LVRT_flag” and “AGO_enable” to 0) by changing these two parameters in thedyr file template.
The LVRT under voltage protection curve is represented by the parameters shown inTable 12.
Table 12 LVRT settings (LVRT enabled)
Voltage limit Setting Timeout Setting
ULVRT1 0.00 tLVRT1 0.55 sULVRT2 0.70 tLVRT2 2.6 s
ULVRT3 0.70 tLVRT3 10 s
ULVRT4 0.70 tLVRT4 10 s
Note : the LVRT protection settings may vary with different turbine models. The user issuggested to check the General Specification of the specified turbine for the LVRT protec-tion settings.
The default low voltage ride-through protection settings for a connected turbine is illustrat-ed in Figure 9.
Figure 9 LVRT protection limits for Vestas GridStreamer TM
wind turbine defined in the protection sub-routine VWVPR6
7.3.3 Frequency Protection
The frequency protection is implemented having three protection levels for high frequen-
cies and three protection levels for low frequencies.
Each protection level can be accepted within a certain time-window. If this window is ex-ceeded, the turbine disconnects. Setting frequency and time limit equal to zero disables
the protection level.
The frequency protection setting may be seen below in Figure 10.
Frequency variation
Time
Hz50 of 1 of 2 of 3
tof1
tof2
tof3
uf 1uf 2uf 3
tuf1
tuf2
tuf3
Figure 10 Frequency protection limits defined in the protection model VWFPR6
Table 13 Frequency protection settings
Frequency limit Setting for 50Hz Setting for 60Hz Timeout Setting
OF1 53 Hz 63.6 Hz tOF1 200 ms
OF2 53 Hz 63.6 Hz tOF2 200 ms OF3 53 Hz 63.6 Hz tOF3 200 ms UF1 47 Hz 56.4 Hz tUF1 200 ms UF2 47 Hz 56.4 Hz tUF2 200 ms UF3 47 Hz 56.4 Hz tUF3 200 ms
Note : the frequency protection settings may vary with different turbine models. The user issuggested to check the General Specification of the specified turbine for the frequency protection settings.
VESTAS R&D Nonstandard Model Data SheetModel Version 8.0.1 GSCORE
GSCORE
VESTAS WIND TURBINE GENERATOR (BEHAVIOURAL MODEL)
This model is located at system bus # ________ (IBUS)Machine ID # ________ (IMACH)This model uses 2 ICONs starting with # ________ (I)and 45 CONs starting with # ________ (J)and 23 STATEs starting with # ________ (K)and 104 VARs starting with # ________ (L)
ICON # Value Description
I Not in use
I+1 Flag for internal loop mode,FI 0 = not in internal loop mode
1 = in internal loop
CON # Vaule Description
J Generator kVA rating, S
J+1 Gen voltage rating, V
J+2 Rated rotor current, IR
J+3 Time constant for 10msmoving average filter
J+4 Disp.
J+5 Disp. J+6 Disp. J+7 P proportional gain
J+8 P integral gain
J+9 Q proportional gain
J+10 Q integral gain
J+11 Lower limit for P control loop
J+12 Upper limit for P control loop
J+13 Lower limit for Qref
J+14 Upper limit for QrefJ+15 Limit for P (LVRT) J+16 Limit for P at zero voltage
(LVRT) J+17 Slope for the power derate
curve (LVRT)
J+18 P positive slope J+19 P negative slope J+20 Disp. J+21 Disp. J+22 Disp. J+23 Model interface, MI
J+24 Stator volt. filter time const
(current injection)
J+25 Stator volt. threshold (currentinjection)
J+26 Current injection mode.
J+27 Reserved for Future Version J+28 Reserved for Future Version J+29 Reserved for Future Version J+30 Reserved for Future Version J+31 Reserved for Future Version J+32 Reserved for Future Version J+33 Reserved for Future Version J+34 Reserved for Future Version J+35 Reserved for Future Version J+36 Reserved for Future Version J+37 Reserved for Future Version
J+38 Reserved for Future Version
J+39 Reserved for Future Version
J+40 Reserved for Future Version
J+41 Reserved for Future Version
J+42 Reserved for Future Version
J+43 Reserved for Future Version
J+44 Reserved for Future Version
STATE # Description
K Integ for active powerK+1 Integ for reactive power
K+2 Integ for measured voltage (10msmoving average)
K+3 Transient volt. real part
K+4 Transient volt. imag part
K+5 Disp. K+6 Disp. K+7 Stator volt. filtered real part (current
injection)
K+8 Stator volt. filtered imag part (currentinjection)
K+11 Reserved for Future Version K+12 Reserved for Future Version K+13 Reserved for Future Version K+14 Reserved for Future Version
K+15 Reserved for Future Version K+16 Reserved for Future Version K+17 Reserved for Future Version K+18 Reserved for Future Version K+19 Reserved for Future Version K+20 Reserved for Future Version K+21 Reserved for Future Version K+22 Reserved for Future Version
VARs # Description
L Disp.
L+1 Active current
L+2 Storage for raw real powerL+3 Storage for raw reactive power
L+4 Saved angle of stator voltage.
L+5 Reactive current
L+6 Voltage in power flow, real part
L+7 Voltage in power flow, imag prt
L+8 P from P_loop after limiter (pu)
L+9 Qref after limiter (pu)
L+10 Q from Q_loop after limiter (pu)
L+11 Active current from P_loop (pu)
L+12 Rective current from Q_loop (pu)
L+13 Memory for delay.
L+14 Memory for delay.
L+15 Memory for delay.
L+16 Memory for delay.
L+17 Memory for delay.
L+18 Memory for delay.
L+19 Memory for delay.
L+20 Memory for delay.
L+21 Memory for delay.
L+22 Memory for delay.
L+23 Memory pos. for delay. .
L+24 Required no. of pos. for delay.
L+25 Delta time for delay.
L+26 Next time step for update delay.
L+27 Delayed angle of voltage.L+28 Limited Pref in normal control L+29 Pref after slope limiter L+30 Saved time for slope limiter L+31 P before LVRTslope limiter L+32 Limited Pref in LVRT L+33 Pref before slope limiter L+34 Max value for active current in LVRT L+35 Reserved for Future Version L+36 Reserved for Future Version L+37 Reserved for Future Version L+38 Save internal loop store(K) for
GSCORE
L+39 Save internal loop store(K+1) for
GSCORE
L+40 Save internal loop store(K+2) forGSCORE
L+41 Save internal loop store(K+3) for
GSCOREL+42 Save internal loop store(K+4) forGSCORE
L+43 Save internal loop store(K+5) forGSCORE
L+44 Save internal loop store(K+6) forGSCORE
L+45 Save internal loop store(K+7) forGSCORE
L+46 Save internal loop store(K+8) forGSCORE
L+47 Save internal loop store(K+9) forGSCORE
L+48 Save internal loop store(K+10) for
GSCOREL+49 Save internal loop store(K+11) forGSCORE
L+50 Save internal loop store(K+12) forGSCORE
L+51 Save internal loop store(K+13) forGSCORE
L+52 Save internal loop store(K+14) forGSCORE
L+53 Save internal loop store(K+15) forGSCORE
L+54 Save internal loop store(K+16) forGSCORE
L+55 Save internal loop store(K+17) for
GSCOREL+56 Save internal loop store(K+18) forGSCORE
L+57 Save internal loop store(K+19) forGSCORE
L+58 Save internal loop store(K+20) forGSCORE
L+59 Save internal loop store(K+21) forGSCORE
L+60 Save internal loop store(K+22) forGSCORE
L+61 Save internal loop store(K) forGSLVRT
L+62 Save internal loop store(K+1) for
GSLVRTL+63 Save internal loop store(K+2) forGSLVRT
GSPWRCL+73 Save internal loop store(K+2) forGSPWRC
L+74 Save internal loop store(K+3) forGSPWRC
L+75 Save internal loop store(K+4) forGSPWRC
L+76 Save internal loop store(K+5) forGSPWRC
L+77 Save internal loop store(K+6) forGSPWRC
L+78 Save internal loop store(K) forGSMEAS
L+79 Save internal loop store(K+1) for
GSMEASL+80 Save internal loop store(K+2) forGSMEAS
L+81 Save internal loop store(K+3) forGSMEAS
L+82 Save internal loop store(K+4) forGSMEAS
L+83 Save internal loop store(K+5) for
GSMEAS
L+84 Save internal loop store(K+6) forGSMEAS
L+85 Save internal loop store(K+7) for
GSMEASL+86 Reserved for Future Version L+87 Reserved for Future Version L+88 Reserved for Future Version L+89 Reserved for Future Version L+90 Reserved for Future Version L+91 Reserved for Future Version L+92 Reserved for Future Version L+93 Reserved for Future Version L+94 Reserved for Future Version L+95 Reserved for Future Version L+96 Reserved for Future Version L+97 Reserved for Future Version L+98 Reserved for Future Version L+99 Reserved for Future Version
L+100 Reserved for Future Version L+101 Reserved for Future Version L+102 Reserved for Future Version L+103 Reserved for Future Version
Note: The synchronous generator arrays (SPEED, PMECH, XADIFD and ECOMP) are used for identifyingthe CONEC models placement in the ICON, CON, STATE and VAR arrays.Note: The synchronous generator arrays (EFD) is used for identifying models connected to the same bus.
VESTAS R&D Nonstandard Model Data SheetModel Version 8.0.1 GSVARS
GSVARS
VESTAS WIND TURBINE VARIABLES INTERFACE
This model is located at system bus # ________ (IBUS)Machine ID # ________ (IMACH)This model uses 2 ICONs starting with # ________ (I)and no CONsand no STATEsand 30 VARs starting with # ________ (L)
ICON # Value Description
I No. of wind farm bus, IBUS
I+1 Wind farm machine ID, IMACH
VARs # Description
L Available real power, Pset
L+1 Q or PF setpoint, Qset
L+2 Real power after PQ chart, Ptrim
L+3 React pwr after PQ chart, Qtrim
L+4 Real power after slope lim, Plim
L+5 React pwr after slope lim, Qlim
L+6 Real power request, Pref L+7 Reactive power request, Qref
L+8 Real power output, Pactual
L+9 Reactive power output, Qactual
L+10 Real power measurement, Pmeas
L+11 React pwr measurement, Qmeas
L+12 Disp.
L+13 Disp. L+14 Disp. L+15 AGO (LVRT) status, S AGO
L+16 Active current for LVRT, IP
L+17 Reactive current for LVRT, IQ
L+18 Reserved for Future Versions L+19 Reserved for Future Versions L+20 Reserved for Future Versions L+21 Reserved for Future Versions L+22 Reserved for Future Versions L+23 Reserved for Future Versions
L+24 Reserved for Future Versions L+25 Reserved for Future Versions L+26 Reserved for Future Versions L+27 Reserved for Future Versions
L+28 Reserved for Future Versions
L+29 Reserved for Future Versions
Note: All power quantities are in per unit on the machine base.
VESTAS R&D Nonstandard Model Data SheetModel Version 8.0.1 GSLVRT
GSLVRT
VESTAS WIND TURBINE GENERATOR LOW VOLTAGE RIDE THROUGH
This model is located at system bus # ________ (IBUS)Machine ID # ________ (IMACH)This model uses 3 ICONs at # ________ (I)and 65 CONs starting with # ________ (J)and 10 STATEs starting with # ________ (K)
and 35 VARs starting with # ________ (L)
ICON # Value Description
I WTG bus number, IBUS
I+1 WTG machine ID, IMACH
I+2 AGO enabler (1=enable)
CONs # Value Description
J AGO threshold, V AGO
J+1 Disp.
J+2 RegainPQ delay
J+3 LVRT IP pos slope, RIP+, Ifequal zero then disabled
J+4 LVRT IP neg slope, RIP-, Ifequal zero then disabled
J+5 LVRT IQ pos slope, RIQ+, Ifequal zero then disabled
J+6 LVRT IQ neg slope, RIQ-, Ifequal zero then disabled
J+7 Active Current Priority1- Active current priority0- Reactive current priority
J+8 Current Overload Factor (Ippriority)
J+9 Offset (Ip priority)
J+10 Gain (Ip priority)J+11 Short term current overload
J+63 HVRT Leaving ThresholdJ+64 LVRT instantaneous Voltage
Threshold
STATE # Description
K State for 10ms FRT Filter
K+1 State for 1 min Average Filter
K+2 Not in use K+3 Not in use K+4 Not in use K+5 Not in use K+6 Not in use K+7 Not in use K+8 Not in use K+9 Not in use
VARs # Description
L Saved time for slope limiters
L+1 Slope limit value for IP
L+2 Slope limit value for IQ
L+3 Disp.
L+4 Input voltage for Ip calculation (Ippriority)
L+5 Active current (Ip priority)
L+6 Upper limit of Iq (Ip priority)L+7 Lower limit of Iq (Ip priority)
L+8 Voltage signal used to compare withthe threshold
L+9 Limited value of One minute movingaverage
L+10 Disp.
L+11 Memory for old value of realVoltStat
L+12 negative sequence magnitude (U-)
L+13 Recorded Time for FRT entry L+14 Reserved for Future Version L+15 Reserved for Future Version L+16 Reserved for Future Version
L+17 Reserved for Future Version L+18 Reserved for Future Version L+19 Reserved for Future Version L+20 Reserved for Future Version L+21 Reserved for Future Version L+22 Reserved for Future Version L+23 Reserved for Future Version L+24 Reserved for Future Version L+25 Reserved for Future Version L+26 Reserved for Future Version L+27 Reserved for Future Version L+28 Reserved for Future Version
L+29 Reserved for Future Version L+30 Reserved for Future Version L+31 Reserved for Future Version L+32 Reserved for Future Version L+33 Reserved for Future Version L+34 Reserved for Future Version
Note: All power quantities are in per unit on the machine base. Slope limits are in p.u. per second.Note: For the parameters that are not used for GridStreamer
VESTAS R&D Nonstandard Model Data SheetModel Version 8.0.1 GSPWRC
GSPWRC
VESTAS WIND TURBINE GENERATOR POWER CONTROL
This model is located at system bus # ________ (IBUS)Machine ID # ________ (IMACH)This model uses 3 ICONs at # ________ (I)and 30 CONs starting with # ________ (J)and 7 STATEs starting with # ________ (K)
Note: All power quantities are in per unit on the machine base. Slope limits are in p.u. per second.DYRE input line:0 'USRMDL' 0 'GSPWRC' 8 0 3 30 7 10 ICON(I) ICON(I+1) ICON(I+2) CON(J)…CON(J+29)/
VESTAS R&D Nonstandard Model Data SheetModel Version 8.0.1 GSMEAS
GSMEAS
VESTAS WIND TURBINE MEASUREMENTS MODEL
This model is located at system bus # ________ (IBUS)Machine ID # ________ (IMACH)This model uses 2 ICONsand 10 CONs starting with # ________ (J)and 8 STATEs starting with # ________ (K)
and 5 VARs # ________ (L)
ICON # Value DescriptionI No. of wind farm bus, IBUS
I+1 Wind farm machine ID, IMACH
CON # Value DescriptionJ Real pwr time const, TP
J+1 React pwr time const, TQ J+2 Reserved for Future Version J+3 Reserved for Future Version J+4 Reserved for Future Version J+5 Reserved for Future Version
J+6 Reserved for Future VersionJ+7 Reserved for Future VersionJ+8 Reserved for Future Version
J+9 Reserved for Future Version
STATE # DescriptionK Active power filter
K+1 Reactive power filter
K+2 Reserved for Future Version K+3 Reserved for Future Version K+4 Reserved for Future Version K+5 Reserved for Future Version K+6 Reserved for Future Version K+7 Reserved for Future Version
VARs # Description
L Saved AGO stateL+1 Reserved for Future Version L+2 Reserved for Future Version L+3 Reserved for Future Version L+4 Reserved for Future Version
VESTAS R&D Nonstandard Model Data SheetModel Version 8.0.1 GSVPRT
GSVPRT
VOLTAGE PROTECTION FOR VESTAS WIND TURBINE GENERATORS
This model is located at system bus # ________ (IBUS)Machine ID # ________ (IMACH)This model uses 7 ICONs starting with # ________ (I)and 30 CONs starting with # ________ (J)and no STATEs
and 18 VARs starting with # ________ (L)
ICON # Value Description
I No. of wind farm bus, IBUS
I+1 Wind farm machine ID, IMACH
I+2 Protection enable flag,
SPROT
I+3 LVRT protection flag, SLVRT
I+4 Lockout state (memory)
I+5 Relay state (memory)
I+6 Relay state (memory)
CONs # Value Description
J Extreme UV limit
J+1 Extreme UV timeout
J+2 Short-term UV limit
J+3 Short-term UV timeout
J+4 Continuous UV limit
J+5 Continuous UV timeout
J+6 Continuous OV limit
J+7 Continuous OV timeout
J+8 Short-term OV limit
J+9 Short-term OV timeout
J+10 Extreme OV1 limitJ+11 Extreme OV1 timeout
J+12 Extreme OV2 limit
J+13 Extreme OV2 timeout
J+14 Extreme OV3 limit
J+15 Extreme OV3 timeout
J+16 Reserved for Future Version J+17 Reserved for Future Version J+18 Reserved for Future Version J+19 Reserved for Future Version J+20 Reserved for Future Version
J+21 Reserved for Future Version
J+22 LVRT extreme UV limit
J+23 LVRT extreme UV limittimeout
J+24 LVRT short-term UV
J+25 LVRT Short-term UV timeout
J+26 LVRT mid-term UV
J+27 LVRT mid-term UV timeout
J+28 LVRT continuous UV limit
J+29 LVRT continuous UV
timeout
VARs # Description
L Saved time for extreme UV
L+1 Saved time for short-term UV
L+2 Saved time for mid-term UV
L+3 Not used, but holds last UV recoverytime
L+4 Saved time for mid-term OV
L+5 Saved time for short-term OV
L+6 Saved time for extreme OV1
L+7 Saved time for extreme OV2
L+8Saved time for extreme OV3
L+9 Reserved for Future Version L+10 Reserved for Future Version L+11 Reserved for Future Version
L+12 Reserved for Future Version
L+13 Reserved for Future Version
L+14 Saved time for LV extreme UV
L+15 Saved time for LV mid-t UV
L+16 Saved time for LV short-t UV
L+17 Not used, but holds last LVRT recoverytime
Note: All voltage quantities are in pu, time quantities in seconds. ‘UV’ = undervoltage, ‘OV’ = overvoltage. Iffewer than three distinct undervoltage and overvoltage levels are required, set two consecutive limit values
VESTAS R&D Nonstandard Model Data SheetModel Version 8.0.1 GSFPRT
GSFPRT
FREQUENCY PROTECTION FOR VESTAS WIND TURBINE GENERATORS
This model is located at system bus # ________ (IBUS)Machine ID # ________ (IMACH)This model uses 3 ICONs starting with # ________ (I)and 12 CONs starting with # ________ (J)and no STATEs
and 7 VAR at # ________ (L)
ICON # Value Description
I No. of wind farm bus, IBUS
I+1 Wind farm machine ID, IMACH
I+2 Relay state (memory)
CONs # Value Description
J Extreme Under-Freq limit
J+1 Extreme Under-Freqtimeout
J+2 Short-term Under-Freq limit
J+3 Short-term Under-Freqtimeout
J+4 Continuous Under-freq limit
J+5 Continuous Under-freqtimeout
J+6 Continuous Over-freq limit
J+7 Continuous Over-freq
timeout
J+8 Short-term Over-Freq limit
J+9 Short-term Over-Freqtimeout
J+10 Extreme Over-Freq limit
J+11 Extreme Over-Freq timeout
VARs # Description
L Saved time for Extreme Under-Freq
L+1 Saved time for Short-term Under-Freq
L+2 Saved time for Continuous Under-Freq
L+3 Not_used, but can hold last Under-Freqrecovery time
L+4 Saved time for Continuous Over-Freq
L+5 Saved time for Short-term Over-Freq
L+6 Saved time for Extreme Over-Freq
Note: All frequency quantities are in pu, time quantities in seconds. If fewer than three distinct under frequencyand over frequency levels are required, set two consecutive limit values equal, and ensure both have the sametimeout figure. Also, the following relationships must hold for declaring the frequency limits if they are utilised:
Figure 11 Vestas V112 GS 3.075MW wind turbine PQ Capability Chart with active power (P) on X-axisand reactive power (Q) on Y-axis. The above plot refers to LV (0.65kV) side of the WTG transformer and
with the LV side voltage between 0.9 p.u. and 1.1 p.u.
Table 15 Parameters of the PQ capability for V112 GS 3.075MW
Appendix B. Data for V112/V117/V126 GS 3.3MW 50/60Hz
B.1 PQ Capability Chart
Figure 12 Vestas V112/V117/V126 GS 3.3MW wind turbine PQ Capability Chart with active power (P) onX-axis and reactive power (Q) on Y-axis. The above plot refers to LV (0.65kV) side of the WTG trans-
former and with the LV side voltage between 0.9 p.u. and 1.1 p.u.
Table 18 Parameters of the PQ capability for V112/V117/V126 GS 3.3MW Offshore