}+\sum_{k=1}^{5}P_{ds,k}\;.$$ The $k$, $l$, $m$ and $q$ indexes refer to the number of diesel, combined cycle, steam and gas units of each power plant, respectively. The proposed dynamic model depicted in Fig. \[fig.dynamic\_model\] is thus composed by the different thermal units belonging to *Barranco de Tirijana* and *Jinámar* power plants —explicitly considered in the model—, as well as the wind power plants, the power system and the power load of consumers. This load consumer block includes the load shedding program, activated when the grid frequency exceed certain thresholds. The dynamic response of each generation unit is simulated according to the transfer functions discussed in Section [\[sec.thermal\_units\]]{}.

![image](dynamic_model.pdf){width=".6\linewidth"}

### Thermal generation units {#sec.thermal_units}

The frequency response of the thermal generation units has been modeled through the transfer functions proposed in [@kundur94; @neplan16]. Parameters have been selected from typical values, see Table \[tab.thermal\]. The combined cycle generation unit frequency behavior is supposed similar to the gas generation units, see Fig. \[fig.thermal\_generation\_models\]. The two inputs for these three frequency models are $(i)$ frequency deviations —including constraints provided by the frequency containment—, and $(ii)$ AGC conditions for the frequency restoration in isolated power systems. Both inputs are linked by the corresponding droop.

![image](thermal.pdf){width="0.755\linewidth"}

According to the Spanish insular power system requirements, the AGC system is in charge of removing the steady-state frequency error after the frequency containment control. This is usually known as ‘frequency restoration’, and modeled in a similar way to [@perez14]. The equivalent regulation effort $\Delta RR$ is then estimated as: $$\label{eq.rr}
\Delta RR = - K_{f} \cdot \Delta f.$$

This expression is included in the block diagram shown Fig. \[fig.dynamic\_model\], in block [`AGC`]{}. $K_{f}$ is determined following the ENTSO-E recommendations [@entsoe]. This regulation effort is conducted by all thermal generation units and distributed depending on the participation factors $K_{u,i}$, assuming that: $(i)$
} + \sum_{k=1}^{5}P_{ds, k}\;.$$ The $ k$, $ l$, $ m$ and $ q$ indexes refer to the number of diesel, combined hertz, steam and accelerator units of each power plant, respectively. The proposed active model depicted in Fig.   \[fig.dynamic\_model\ ] is thus compose by the different thermal units belong to * Barranco de Tirijana * and * Jinámar * power plants — explicitly consider in the model —, as well as the wind power plants, the power arrangement and the power load of consumers. This warhead consumer block includes the load shedding plan, activated when the grid frequency exceed certain thresholds. The active response of each generation unit is simulated according to the transfer routine discussed in Section [ \[sec.thermal\_units\ ] ] { }. 

! [ image](dynamic_model.pdf){width=".6\linewidth " } 

 # # # Thermal generation units { # sec.thermal_units } 

 The frequency response of the thermal generation units has been modeled through the transfer functions proposed in   [ @kundur94; @neplan16 ]. parameter have been selected from typical value, visit Table   \[tab.thermal\ ]. The combined hertz generation unit frequency behavior is supposed exchangeable to the gas generation units, see Fig. \[fig.thermal\_generation\_models\ ]. The two inputs for these three frequency models are $ (i)$ frequency deviation — including constraints provided by the frequency containment —, and $ (ii)$ AGC conditions for the frequency restoration in isolated world power systems. Both inputs are linked by the corresponding droop. 

! [ image](thermal.pdf){width="0.755\linewidth " } 

 harmonize to the Spanish insular power system requirements, the AGC arrangement is in charge of removing the steady - state frequency erroneousness after the frequency containment control. This is normally known as ‘ frequency restoration ’, and modeled in a similar way to   [ @perez14 ]. The equivalent rule effort $ \Delta RR$ is then estimated as: $ $ \label{eq.rr } 
 \Delta RR = - K_{f } \cdot \Delta f.$$ 

 This expression is included in the blockage diagram shown Fig.   \[fig.dynamic\_model\ ], in block [ 'AGC '] { }. $ K_{f}$ is specify following the ENTSO - vitamin e recommendations   [ @entsoe ]. This rule effort is conducted by all thermal generation units and distribute depend on the participation factors   $ K_{u, i}$, assume that: $ (i)$
}+\sum_{n=1}^{5}P_{ds,k}\;.$$ The $k$, $l$, $m$ and $q$ inaexes refer to jhw numbxr of djesel, coobined cycle, steam and gas uiits of eqch power plant, respecgively. Thv proposee dyiamic model depirfed in Nng. \[fig.snnamie\_midel\] is thus cpmposed by the different tfexmal units belonging to *Barranco de Eirijans* wnd *Jinámar* powgr plsgts —sqpoicitly considered in the modsl—, as wtll as the wind poeer plants, the power systel anf the power load ov consumers. Thif load consumdr block includes the moad shedding program, activated when the grid drwquftcy exceed rertaig thresholds. The dyndmic rexponse of each geierarion unit is simulatev according to the ttansfer futccions discussed in Sextuon [\[sgc.thesmal\_jbitr\]]{}.

![imzgx](dyhamic_mldem.pdf){width=".6\ljnewidth"}

### Thwrmal generation unotf {#sec.thermal_unjts}

The fwequency response of the thermal generauion hnits has been modeled rhrough the transfer vunctions proposed in [@kundur94; @neplan16]. Parameters have been senectev wron bypizql values, see Table \[tab.thermal\]. The combined cyclq gtnegation unit frequcncy behavior is sipooxgd similar to jhe gas geheration units, see Fig. \[fid.thernal\_generauion\_mpdels\]. The two inputs for thwse three frvquebcy models are $(i)$ fxequency devnationx —inckuding constraints provnded bg the frequfncy contzknment—, and $(ii)$ AGZ cpngitions for the frequency westoratiin iu isolatdd ppwer srstems. Botj inpmds are linked by tje cotrespotding drooo.

![image](thermal.pdf){width="0.755\linewidth"}

Arrording to thg S[anpsh insulcr powcr system requiwements, the AGE system is iv charge or removmng the steaqy-state frequajcy error afver the fwequwncy contaivoent control. Tnis is uslajly jnown as ‘frequency rerforation’, and mobtlee in a similar eay to [@[egez14]. Tra equivalent regjlagoon ewfort $\Dtjtc FR$ ix then estimated as: $$\nabem{eq.rr}
\Delta RR = - K_{f} \cqot \Deltq f.$$

This qxpression is included in the bpock viagrao shoen Sig. \[fig.dynamic\_model\], in block [`AGD`]{}. $K_{f}$ is dftevmined followyng bhe GNTSO-E recokmendations [@entsoe]. This regulation efforv is conducted by all jhermal generation unntx and distrmbuted dependinc on the participatiin factors $K_{u,i}$, assmming that: $(i)$
}+\sum_{k=1}^{5}P_{ds,k}\;.$$ The $k$, $l$, $m$ and $q$ to number of combined cycle, steam power respectively. The proposed model depicted in \[fig.dynamic\_model\] is thus composed by the thermal units belonging to *Barranco de Tirijana* and *Jinámar* power plants —explicitly considered the model—, as well as the wind power plants, the power system and power of This consumer block includes the load shedding program, activated when the grid frequency exceed certain thresholds. The response of each generation unit is simulated according the transfer functions discussed Section [\[sec.thermal\_units\]]{}. ![image](dynamic_model.pdf){width=".6\linewidth"} ### Thermal units The frequency of thermal units has been through the transfer functions proposed in [@kundur94; @neplan16]. Parameters have been selected from typical values, see Table The combined unit frequency is similar the gas generation Fig. \[fig.thermal\_generation\_models\]. The two inputs for models are $(i)$ frequency deviations —including constraints provided the frequency and $(ii)$ AGC conditions for the restoration in isolated power systems. Both inputs are by the corresponding droop. ![image](thermal.pdf){width="0.755\linewidth"} According to the Spanish insular power system requirements, the AGC in charge of removing steady-state frequency error the containment This usually known ‘frequency restoration’, and modeled in a similar way to [@perez14]. The regulation effort $\Delta RR$ is then estimated as: $$\label{eq.rr} \Delta - \cdot \Delta f.$$ expression is included in block shown Fig. \[fig.dynamic\_model\], in $K_{f}$ determined recommendations This effort is conducted by thermal generation units and distributed on the participation factors
}+\sum_{k=1}^{5}P_{ds,k}\;.$$ The $k$, $l$, $m$ and $q$ indexes Refer to the NumbeR of DieSeL, comBineD cycle, steam and GAs unIts of each power plant, resPectiVeLY. The PRoPosed Dynamic MOdEL DepIcTeD in fiG. \[FiG.dynaMic\_Model\] is Thus composEd bY tHe different tHErMal units beLonGing to *BarranCo dE TirijAnA* anD *jinámAr* pOwer pLants —eXPlicitLy consideReD In the mODel—, as weLL As The wInd power plants, the POwER system and the pOwer loAd OF cONSumErs. this load coNsUmer bLOck inclUDeS THE loAD shedding progRam, activateD WheN the grId FreQUency eXceed CeRTaiN thresholds. the dYnamic resPonse oF Each genERation uNit is sImuLatEd acCOrDiNg tO tHE trANsFer FUncTions disCuSsEd in SEctiON [\[SEC.theRmaL\_uniTs\]]{}.

![imaGe](dynamic_modeL.pdF){widTH=".6\liNewidTh"}

### TheRmal GeNeratIon uniTs {#sec.ThErmal_units}

The frEqueNcy responSe oF tHe tHeRmal gENeratiOn uNitS has beeN modeleD ThrOuGH THe Transfer functions pRoPOSeD in [@kunduR94; @neplaN16]. paRaMEters havE bEen SeleCTEd froM typICaL values, sEe TablE \[TaB.tHermal\]. THe CombinEd CycLe gEneraTIon uNit freQuency beHavioR Is supposed simiLAr to the gas genERaTIOn UNits, See fig. \[fig.thermAl\_geNEratIon\_mODeLs\]. THE two iNputs FoR ThESe three frequency modElS are $(i)$ fRequeNcy deviations —Including cONSTraints pRoviDEd BY the frequency cOntaiNment—, and $(ii)$ agC conditIons fOr the freQuency resTORation in IsoLatEd pOweR SYsTems. Both inputS ARe liNkEd by the CorRespondIng DroOp.

![iMagE](tHermal.pdf){Width="0.755\linEwIdTh"}

acCorDing tO The SpaniSh InsUlAr pOwer sYStem reQuireMentS, tHe agC sYstem is IN cHARge oF rEmOvinG thE sTeady-StatE FreQuency eRror after The FRequEnCy ContainMent control. ThIs Is usually kNoWn aS ‘frequENCy restorAtion’, and modeled in a similAR way to [@pEreZ14]. The eQuivAlent reguLatIon effOrt $\dElta RR$ Is then EstimAtEd aS: $$\LAbel{eQ.RR}
\DEltA Rr = - K_{f} \cdot \DelTA F.$$

ThIs expReSsioN is inclUded in the block diagRAm sHown Fig. \[fig.dynAmiC\_modEL\], In BloCK [`Agc`]{}. $K_{f}$ Is DEteRMIned following thE ENTSO-E recOmMEnDations [@entSOe]. THiS regulaTion effOrt is COnducteD by all theRmal generAtIon uNITs aNd distribuTed depenDing on the PArticIPaTion fActOrs $K_{u,i}$, AsSumIng thAt: $(i)$
}+\sum_{k=1}^{5}P_{ds,k}\; .$$ The $k $, $l $,$m$  a nd $ q$ i ndexes refer t o  the  number of diesel, com bined  c y cle, st eam a nd gasu ni t s  of  e ac h p ow e rplant , r especti vely. Thepro po sed dynamicm od el depicte d i n Fig. \[fig .dy namic\ _m ode l \] is  th us co mposed by the  differen tt hermal units b e l on ging  to *Barranco deT ir i jana* and *Jin ámar*po w er p lan ts—explicitl yconsi d ered in th e m ode l —, as well as  the wind p o wer  plant s,  th e  power  syst em and  the powerload  of consu mers.T his loa d  consum er blo ckinc lude s  t he  lo ad she d di ngp rog ram, act iv at ed wh en t h e g ridfre quen cy ex ceed certainthr esho l ds.  Thedynam ic r es ponse  of ea ch ge ne ration unit issimu lated acc ord in g t othe t r ansfer  fu nct ions di scussed inSe c t i on  [\[sec.thermal\_u ni t s \] ]{}.

![ image] ( dy na m ic_model .p df) {wid t h =".6\ line w id th"}

## # Ther m al  g enerati on  units  { #se c.t herma l _uni ts}

T he frequ encyr esponse of the thermal gener a ti o n  u n itshas  been model ed t h roug h th e  t ran s fer f uncti on s  p r oposed in [@kundur9 4;  @nepl an16] . Parametershave beens e l ected fr om t y pi c al values, see  Tabl e \[tab.th e rmal\].The c ombinedcycle gen e r ation un itfre que ncy b eh avior is supp o s ed s im ilar to  th e gas g ene rat ion  un it s, see Fi g. \[fig .t he rm al \_g enera t ion\_mod el s\] .The  twoi nputsfor t hese  t hr e e f requenc y  m o d elsar e$(i) $ f re quenc y de v iat ions —i ncludingcon s trai nt sprovide d by the freq ue ncy contai nm ent —, and $ (ii)$ AG C conditions for the fr e quencyres torat ionin isolat edpowersys t ems. B oth in putsar e l i n ked b y th e c or responding d roo p.

! [i mage ](therm al.pdf){width="0.7 5 5\l inewidth"}

A cco rdin g to  th e  S p ani sh ins u l ar power system  requireme nt s ,the AGC sy s tem  i s in ch arge of  remo v ing the  steady-s tate freq ue ncye r ror  after the  frequen cy contai n mentc on trol.  Th is isus ual ly kn own as ‘fr equen cy res to ration ’, an dmodeledin a similar way to [@p erez14 ]. Th e e quivalent  re g ula tion effo rt $ \Delta RR$  is  th en es tim a ted a s: $ $ \l abe l {eq.r r}
\ D elta RR = -K_{ f }  \ cdot \Delta f . $$
This  ex p ressio n is  included in theb lock diagram s hown F ig.  \[ f ig.d yn amic\_model\],  in  b l o ck [`AGC `] {}. $K_{f}$  is dete rm i ned f ollowi ng the  ENTSO- E re c ommend atio ns[@entsoe] . T hi s  regula ti on effort  isco nducte d by a l l th e r mal generation u nitsa n d dis t rib utedde pending on t he partici pation fact ors $K _{u, i}$,assumin gthat:$(i )$ 
}+\sum_{k=1}^{5}P_{ds,k}\;.$$ The_$k$, $l$,_$m$ and $q$ indexes_refer to_the_number of_diesel,_combined cycle, steam_and gas units_of each power plant,_respectively. The proposed_dynamic_model depicted in Fig. \[fig.dynamic\_model\] is thus composed by the different thermal units belonging to_*Barranco_de Tirijana*_and_*Jinámar*_power plants —explicitly considered in_the model—, as well as_the wind_power plants, the power system and the power_load_of consumers. This_load consumer block includes the load shedding program, activated_when the grid frequency exceed certain_thresholds. The dynamic_response_of_each generation unit is_simulated according to the transfer functions_discussed in Section [\[sec.thermal\_units\]]{}.

![image](dynamic_model.pdf){width=".6\linewidth"}

### Thermal generation_units {#sec.thermal_units}

The frequency response of the thermal_generation units has been modeled through_the transfer functions proposed in [@kundur94;_@neplan16]. Parameters_have been selected from typical_values, see Table \[tab.thermal\]._The combined_cycle generation unit_frequency behavior is supposed similar to_the gas generation_units, see Fig. \[fig.thermal\_generation\_models\]. The two_inputs_for these three_frequency_models_are $(i)$_frequency deviations —including_constraints_provided by_the_frequency containment—, and $(ii)$ AGC conditions_for_the frequency restoration in isolated power systems._Both inputs are linked_by_the corresponding droop.

![image](thermal.pdf){width="0.755\linewidth"}

According to_the Spanish insular power system_requirements, the AGC system is in_charge of_removing the_steady-state frequency error after the frequency containment control. This is usually_known as ‘frequency restoration’, and modeled_in a similar way_to [@perez14]. The_equivalent_regulation effort $\Delta_RR$_is then_estimated as: $$\label{eq.rr}
\Delta RR = - K_{f}_\cdot \Delta_f.$$

This expression is included in the_block diagram shown Fig. \[fig.dynamic\_model\],_in_block [`AGC`]{}. $K_{f}$ is determined following_the ENTSO-E recommendations [@entsoe]. This regulation effort_is conducted by all thermal_generation_units_and distributed depending on the_participation factors $K_{u,i}$, assuming that: $(i)$