logoGlobal Energy Interconnection

Contents

Figure(0

    Tables(0

      Global Energy Interconnection

      Volume 6, Issue 3, Jun 2023, Pages 253-259
      Ref.

      Impact of interconnections and renewable energy integration on the Philippine–Sabah Power Grid systems

      Mr. Tristan G Magallones Jr1,2 ,Jai Govind Singh1
      ( 1.Department of Energy, Environment and Climate Change, SERD, Asian Institute of Technology,Pathumthani 12120, Thailand , 2.Department of Electrical Engineering, Central Mindanao University, Bukidnon 8710, Philippines )

      Abstract

      This study presents a comprehensive impact analysis of the rotor angle stability of a proposed international connection between the Philippines and Sabah, Malaysia, as part of the Association of Southeast Asian Nations(ASEAN) Power Grid.This study focuses on modeling and evaluating the dynamic performance of the interconnected system, considering the high penetration of renewable sources.Power flow, small signal stability, and transient stability analyses were conducted to assess the ability of the proposed linked power system models to withstand small and large disturbances, utilizing the Power Systems Analysis Toolbox (PSAT) software in MATLAB.All components used in the model are documented in the PSAT library.Currently, there is a lack of publicly available studies regarding the implementation of this specific system.Additionally, the study investigates the behavior of a system with a high penetration of renewable energy sources.Based on the findings, this study concludes that a system is generally stable when interconnection is realized, given its appropriate location and dynamic component parameters.Furthermore, the critical eigenvalues of the system also exhibited improvement as the renewable energy sources were augmented.

      0 Introduction

      The need for energy security has triggered an initiative for international cooperation to provide a stable and sufficient energy supply [1].As the power industry continues to evolve, a significant number of blackouts worldwide are attributed to instability arising from increased intersystem connections, expansion of power systems,and advancements in power system control technology[2].Reference [3] reported that insufficient damping of oscillations in practical power systems is one of the most common problems in the energy industry.

      One of the anticipated Association of Southeast Asian Nations (ASEAN) international connections is the 500 kV high-voltage direct current (HVDC) link between the Philippine national grid and the Sabah power grid, which is expected to be realized by 2020 [4].The Philippines-Sabah interconnection is a step towards creating the ASEAN Economic Community (AEC), which recognizes the adverse effects of climate change.With a campaign for greener and cleaner energy sources, the community has pledged to use more renewable energy sources alongside fossil fuels.This, in effect, will further complicate the system’s components and controls.The ASEAN Interconnection Master Plan Study (AIMS) II confirmed that the proposed interconnection is technically and economically feasible within the region, even though small-signal stability and the impact of renewable energy sources are absent from its scope of network analyses.

      The archipelagic Philippines has isolated and interdependent power systems.Studying the effects of domestic interconnection aimed at addressing the problems[5] in the southern grid (Mindanao) of the country has been the motivation of this study, as well as investigating the effect of possible international interconnection.Furthermore,the aggressive targets for renewable energy (RE) sources[6] drove the need for this study to examine their immediate effects on existing systems when interconnected, which was not perceived earlier.The authors in [7] have applied the system stability approach to study the impacts of various penetration levels of renewable energy sources on Java-Bali interconnections.Reference [8] reviewed various studies to understand the impact of VRE integration based on different generation types, penetration levels, and power grid characteristics.The authors suggested that balancing issues arise after 20% generation and stability issues occur after more than 50% VRE penetration.The authors in [9]analyzed and quantified system performance in terms of reliability and vulnerability to changing renewable sources.Standard test systems, that is, the IEEE RTS-96 and IEEE RTS-GMLC test networks, were used to analyze their performance.The authors in [10] reviewed the impacts of various technological solutions and challenges in power systems and suggested a matrix with solutions versus challenges for renewable technologies.Reference [11]studied dynamic models of household power consumption demand using solar PV data to demonstrate the impacts of its distribution characteristics on the power grid.

      The authors in [12] revealed that “Australia has struggled previously while maintaining system stability,despite having a high share of natural gas in their domestic electricity production portfolios”.Reference [13] presented a comprehensive evaluation index system to assess the benefits of transnational power grid interconnection projects.The authors in [14] comprehensively reviewed the different operational benefits and challenges of interconnected power grids.The authors in [15] studied and analyzed the impact of Arab counties on the stability of frequency in a Syrian electrical power grid with and without interconnections.Reference [16] presented a load frequency control approach by developing a general state-space model and then a 4-area power system modeled with renewable sources along with electric vehicles.The authors in [17] presented a fuzzy PID controller optimized with an improved ant colony optimization algorithm.The proposed method was tested on a multi-area, connected power system.A decentralized load-frequency control framework was presented in Ref.[18] using dynamic state estimation.This model was tested with the interconnection of four power systems, and it was claimed that the power fluctuations and tie-line deviations were efficiently controlled.

      Most studies have focused on either standard tests or national-level systems.In addition, Philippines-Sabah interconnection data needs to be analyzed for renewable integration in terms of stability studies.

      The contribution of this study is the modeling of the interconnection between the three islands of the Philippines and Sabah, Malaysia.Furthermore, the impacts of the three islands and their interconnections with Sabah were analyzed in terms of dynamic performance.Moreover, the impacts of renewable energy integration on the static and dynamic behavior of the proposed Philippines-Sabah interconnection were analyzed using the Power Systems Analysis Toolbox(PSAT) software.Various case studies were presented to fully understand the stability and performance of the proposed international link.Stability modeling was performed using PSAT.Another objective of this study is to bring the system’s consolidated power grid data to the national/international level and/or the public domain for further analysis by future researchers.However, this work extensively studied most of the operational aspects of the power system, including small signals, transient interconnections, and the integration of renewable energy.Appendix A provides the data file for further analysis.

      1 Proposed methodologies

      1.1 Small signal stability analysis

      The behavior of a dynamic system can be described by a complete set of n first-order nonlinear ordinary differential and algebraic equations (DAE) [18] in the following form:

      The complete set of differential equations with network equations is arranged in state-space form and linearized about an operating point as follows:

      The small-signal perturbation model of the equation can be solved using a Laplace function.The poles of the system are the roots of the characteristic equations, which can be rewritten as:

      The values of λ (λ= λ1, λ2, λ3λn), which satisfy the characteristic equation, are known as the eigenvalues of matrix A.Complex eigenvalues always occur in conjugate pairs as follows:

      In this study, the eigenvalues were calculated along with the damping ratio and frequency of oscillation for different cases, and the participation factor of ith mode.

      1.2 Transient stability analysis

      This study uses numerical integration for the timedomain simulation, specifically the trapezoidal method,which is an implicit method for solving time-domain analysis in PSAT [17].The system was simulated for 20 s in all cases.

      In each generic iteration i and generic time step t + Δt, and q(i) are as follows:

      where, Ac(i) is a matrix depending on the algebraic and state Jacobian matrices of the system, q(i) is a function that depends on the implicit numerical method, Inx is the identity matrix of the same dimension of the dynamic order of the DAE system and all Jacobian matrices and f(i) are computed at the current point (x(i)(t + Δt), y(i)(t + Δt), tt).

      1.3 Solar photovoltaic generator model

      ASEAN region has a very dynamic system; owing to its developing country members who are on the verge of improving their power systems, the inclusion of solar PV is slowly observed in some existing systems.Figure 1 shows a solar PV generator (SPVG) interconnected to the power grid.Solar panels, converters, and controllers, such as MPPT, are some of the major components [18].This model has no moving parts, indicating that the MPPT and controllers are instantaneous in this study.

      The desired active power is the basis of the current setpoint.As such, the current (I) of the terminal voltage has the following dq form:

      Fig.1 SPVG block diagram

      1.4 Wind turbine system model

      A model of the DFIG wind turbines was used in this study [18].In the simulations, all the wind farms were represented by DFIGs, as depicted in Fig.2.

      Fig.2 Basic configuration of a DFIG wind turbine

      The following equation (1) is the main basis for solving the wind turbine’s actual power.It is very useful in the derivation of the capacity of wind power.

      where Cp is the utilization coefficient, ρ is the mass density of air, Ωw is the rotational velocity of the wind turbine, V is the wind velocity, A is the area of the blade.

      It has the boons of fast acting when it comes to flux linkage dynamics, unlike the DFIG rotating masses dynamics.Hence, the following differential algebraic equations express a simplified model of DFIG based wind generation system:

      where Te is the electrical torque, ψdr is the flux linkage due to d-axis rotor current, ψqr is the flux linkage due to q-axis rotor current, Iqr, Iqs are the q-axis rotor and stator current, respectively, Xr, Xm, Xs are mutual, rotor and stator reactances, ψds is the flux linkage due to d-axis stator current, ψqs is the flux linkage due to q-axis stator current,Ids, Idr are the d-axis rotor and stator current, respectively.

      Balanced conditions were assumed for all operating conditions.All values solved in the above equations were considered per unit.A reference frame that was continuously rotating, where the q-axis led the d-axis frame by 90 electrical degrees, was also implemented during the simulation.

      1.5 Power grid model development

      Reference [5] has developed reduced and approximate Philippine models, which are utilized in this work.These models are created separately: Luzon, Visayas and Mindanao and an additional model, Sabah.The PSAT helps to have models for transmission lines, PV generators, PQ load, substation transformers, and reactive sources used for static power system components.The transmission lines in the model were developed by acquiring data such as line lengths, capacities, interconnectivity and conductor impedances.The transformers’impedances are normalized to a common base MVA rating also.All transformer impedances were assumed to follow the values in [16] and converted to a 100 MVA base.The generator’s location or connection point and maximum and minimum real and reactive power were the data needed for the power flow.All the grids were modeled with constant power loads.MW and MVAR data with a 0.9 power factor were considered,and these loads received high power in a high-voltage bus.This was performed to minimize the number of buses used.However, for Mindanao and Sabah, the loading conditions were based on the population of each region.Shunt compensation is rated as acceptable magnitudes needed to maintain acceptable bus voltage levels.Models have MVAr quantities.The voltage levels were maintained in the range of 1.05-0.95 p.u.

      The generation and load data of WESM [19] on September 9, 2012, period 9 were used for the Luzon model, whereas for the Visayas reduced model, the WESM generation schedule of September 10, 2012, period 4, was used.The peak load conditions provided in the NGCP PubCon Presentation were used for the Mindanao Model[5].The Sabah model load was 874 MW on September 23,2013 [17].

      The dynamic components are based on well-documented PSAT dynamic models for synchronous generators,turbine governors, automatic voltage regulators, solar photovoltaic generators, doubly fed induction generators,and wind turbines [20].Simple versions (classical) of synchronous generators were modeled, where the transient voltage in the q-axis was constant.All generators have no mechanical damping, and saturation effects are neglected.In this work, Model 1 TG was used for the conventional thermal generators, while Model 3 was selected for all the hydroelectric power generators, the same as the reference[21].In the latter model, the inelastic penstock, where the inertia of water and the ideal turbine, is considered.However, the parameters of these models are based on previous literature reviewed [22, 23] and that of the IEEE standard data.

      1.6 Philippine power grid

      A 900-ckt-km 500 kV line passing from north to south serves as the backbone of the Luzon power system, with a meshed 230 kV network from the center portion.It has a generation capacity of 7829 MW.

      Visayas predominantly has a 138 kV voltage level for its power transmission system with some 230 kV connections.A mixture of coal, diesel, and geothermal power plants comprises the main generation profile of the Visayas grid,peaking at 1400 MW.Recent developments in the southern grid of the Philippines have added a 230 kV network to its existing 138 kV and 69 kV lines, which originally had a total of 3100 ck-km.Mindanao has experienced lower capacity together with its growing electricity demand.Sabah, Malaysia, has a power grid composed mainly of 275 kV, 132 kV and 66 kV system voltage levels, but 66 kV lines were not included in this study.A total of 2,272 circuit kilometers of both 275 kV and 132 kV were considered in this study.Sabah’s system relies heavily on gas-based generation in terms of both the generation mix and capacity.As shown in Fig.3, the Philippine-Sabah power grid is the realization of an international connection between Malaysia and the Philippines as part of the ASEAN Interconnection Master Plan.This assumes that the connection of the Visaya-Mindanao grid is realized.This includes a 500 MVA 250 kV two-pole HVDC link between the Luzon and Visayas grid (existing) and the isolated Mindanao grid, as planned in the 1995 transmission master plan study of the National Grid of the Philippines[24].The 500 MW, 500 kV link was achieved in PSAT connecting bus 34 (230 kV) and bus 116 (275 kV).In Fig.3,all link buses are denoted by the first letter of each island;for example, Luzon (L), Visayas (V), Mindanao (M), and Sabah (S).

      Fig.3 Single line diagram of the Philippine-Sabah interconnected model

      2 Results and discussion

      2.1 Isolated grid (existing grids)

      A base case was investigated by ensuring that the models without the new interconnection (Luzon-Visayas and isolated grids of Mindanao and Sabah) worked properly.The data in Appendix A were used for this purpose.Although not shown, the load flow results, for example,voltage levels, were observed to be within the limits, as expected.

      More eigenvalues are observed in Fig.4a because this model comprises two groups of islands.The symbol“x”represents the eigenvalues, whereas the dotted lines represent 1% (closest to the y-axis), 5%, and 10% (farthest from the y-axis) damping lines.The eigenvalues nearest the y-axis (0% damping line) were the least damped eigenvalues.The rotor speed and angles were the least damped and prone to instability when small disturbances occurred.This may be caused by the lack of mechanical damping mechanisms (as assumed earlier) in these oscillatory machines.The least damped generator was located in the Visayas area and had the highest rating.

      A fault was also applied to the system to further investigate its behavior when a large disturbance occurred.A three-phase fault, which was assumed to be the worst type of fault, was used in this study.In all models, the fault can be applied anywhere; however, for analysis, the fault is located either at the bus with the highest generator capacity or at the bus with the largest load.Therefore, other configurations should be investigated in future studies.The fault was applied at 0.2 s, followed by clearing after five cycles.

      Since in the Philippines, the frequency is 60 Hz, the cycle is cleared after 0.0833 s.For simplicity, the same clearing time was used for the Sabah model, which is a 50 Hz system.Figures 5a and 5b show the deviations in the voltage profile and rotor angles observed in the Luzon-Visayas model, characterized by the rise in rotor angles after the 8th second.This is also observable in Fig.3a by the unstable fluctuation in voltage, which then increases after the 13th second.This model is the largest of the three existing grids, which is understandable because it is based only on approximate and reduced data.

      Fig.4 Philippine PSAT model eigenvalues

      This instability was resolved by reducing the clearing time to 0.7 s instead of 0.8333 s, which could be an issue if the circuit breaker does not work faster than five cycles.This implies that a practical system should include faster and more efficient protective devices for withstanding large disturbances.This system is shown in Figures 5c and 5d.This demonstrates the sensitivity of the system to the type of fault applied at certain locations.

      2.2 Interconnected models (Philippine Power Grid)

      The Philippine national power grid will be achieved if the proposed link between Visayas and Mindanao is installed.When the proposed link was installed in the model, the system exhibited small signal instability with two positive eigenvalues.Synchronous machine no.10, located in Visayas, contributed to this instability, as indicated by the participation factor.Synchronous machine No.10 accounted for 18.87% of the state of eigenvalue Nos.195 and 196.

      Fig.5 Luzon-Visayas link model time domain simulation results

      Table 1 shows the original critically damped eigenvalues.It can be seen that same, as the values given in Table 1, rotor angle and speed have the least damping ratio compared to other states.It is interesting to note the -1.535% damping ratio of the synchronous machine 10.Examining its mode, 0.565571 Hz, it means that the instability is in inter-area mode.The following sections will discuss the grid’s instability and acquired stability with the help of power system stabilizers (PSS), which is based on the unstable/weak eigenvalue participation factor.

      As shown in Fig.6, the system experienced instability as indicated by the presence of one complex pair of eigenvalues(highlighted in red) located on the right-half of the S-plane.A total of 422 states with 67 complex pairs were initialized and observed.Two positive eigenvalues were observed in the impending interconnection between the Mindanao and main grids.The eigenvalues of the system with machine No.10 became unstable.This is probably because of the approximate and reduced models.For the analysis,participation factors were used to determine the bus that was weakly damped and most exposed to instability.Generators with the highest participation factors are optimal locations where power oscillation dampers (POD), for example,PSS should be placed to improve stability.When PSS was installed in the unstable generator, Syn_10, the eigenvalues(Fig.7b) became stable.The damping ratio improved from a value of -1.12% to a value of 6.4%.With the correct parameters of PSS, the real part of the eigenvalue Nos.195 and 196 became negative (-0.22186).

      For transient stability, unlike the synchronous machines in the area (Fig.7a) where the fault was applied, machine no.10 in Visayas became unstable, as shown in Fig.7b.This is evident in all generators to which this unstable generator belongs.It oscillated after the fault instead of returning to its original value.This is attributed to the unstable states shown in the eigenvalue analysis.The amplitude of the oscillation continues to increase, as shown.The rotor angle has been evident in showing instability as well.

      Fig.6 Philippine PSAT model eigenvalues

      It can be perceived in Fig.7c that the swinging of the faulted generator has minimal fluctuations.It can be seen in the zoomed graph (Fig.7d) that the angle is starting to be unstable, as seen by the fluctuations in the oscillation with increasing amplitude.

      PSS has effects on the transient stability of the system too.The voltage profile of machine No.24, where the fault bus is attached, shows no significant change in its profile after installing the PSS in the system, as shown in Fig.8a.However, the significant improvement was characterized by machine No.10.It can be seen in Fig.8b the improvement in the voltage of the profile.This zoomedin Fig.8c illustrates how the system tried to compensate for the fault by maintaining the voltage levels at the buses.It is noteworthy to examine machine No.10 since this is where the PSS is installed through its AVR to facilitate the voltage regulation of the machine.

      2.3 Interconnected models (Philippine-Sabah Power Grid)

      The stable model of the entire Philippine power grid was utilized when the link to the Sabah power grid was done.The whole system is composed of 126 buses with 531states considered, 80 of which are complex pairs.Figure 9 shows that the system is stable when small disturbances are concerned.

      Table 1 Critical modes of the Philippine power grid

      Eigenvalue No.Most Associated States Real Part Imaginary Part Frequency/Hz Damping Ratio 195 & 196 Rotor Speed & Angle Mach.No.10 0.0546 ±3.5532 0.5656 -0.0154 155 & 156 Rotor Speed & Angle Mach.No.19 -0.1459 ±6.0191 0.9583 0.0242 153 & 154 Rotor Speed & Angle Mach.No.1 -0.2041 ±6.4031 1.0196 0.0319 157 & 158 Rotor Speed & Angle Mach.No.25 -0.2605 ±6.4031 1.0199 0.0406 151 & 152 Rotor Speed & Angle Mach.No.28 -0.3687 ±7.4121 1.1811 0.0497

      Fig.7 Voltage of synchronous machines without PSS

      Fig.8 Voltage of synchronous machine with PSS

      Fig.9 Philippine–Sabah model eigenvalues

      The eigenvalues are listed in Table 2, and the same states, rotor angle, and rotor speed significantly influence the instability of the system.This is also true when studying the Myanmar Grid system [18].This may be why, of all the dynamic parts, it is the one that is closely taken into account.The critical modes are mostly local aside from the eigenvalue No.196 and 197, which is inter-area.This may be due to the tie-lines between these 5 areas.Since this is asynchronous interconnection, a strong HVDC transmission line should be installed (already in the plan [4]) to improve this eigenvalue.The parameters of the HVDC link used in this study are based on the models in IEEE [25].

      The result of this eigenvalue analysis validates that the connection of the Philippines to Sabah is possible and imposes no major problem when the small-signal stability of the system is concerned.

      Fault analysis was also conducted on the system.The voltage behavior of the entire system when a fault is introduced in one of the buses in Mindanao is relatively stable.This location was selected because Mindanao is the most troubled of all the grids owing to its low capacity;therefore, most problems may occur in this area.The system rotor angles did not exhibit any evidence of instability.The fluctuations that occur during the fault smooth with time.

      2.4 Philippine-Sabah Power Grid with High Penetration of Renewable Energy Sources

      The system study also extended subject to high penetration of renewable energy resources, specifically solar photovoltaic and wind turbines in the future.The data are from the Department of Energy as of February 28, 2014 [6], and also available from The National Grid of the Philippines’(NGCP) Transmission Development Plan Volume 1.However, only the wind power plants and solar power plants that are planned to be erected are considered in this work.The assignment of the lumped wind and solar PV farms is based on the geographical location of such power plants, to which they are proximate to a certain bus in the model where the plants are connected.Reference [26]indicated no immediate plans for additional solar PV and wind farms until Sabah in 2023.The load flow results of the entire system complied with the standards set by the Philippines’Energy Regulatory Commission [27].Figure 10 shows the voltage levels of all 133 buses present in the system, with the addition of 7 buses intended for the RE sources injected in the system that has the planned capacity for the year 2023.The total real and reactive loads were 88.4423 and 11.7391 p.u., respectively.Similarly, the total real and reactive power generations were 91.6839 and 24.0990 p.u., respectively.However,total system real and reactive power losses were 3.2317 and 12.3699 p.u., respectively.Small-signal stability analysis was performed by comparing the eigenvalues during the gradual increase in renewable energy sources in the system.The system could withstand up to a 40% increase, as shown in Table 3.Increasing the RE sources by 50% led to malfunctioning of the system components.Table 3 lists the eigenvalues of the critically damped modes.

      Table 2 Critical modes of the Philippine-Sabah power grid

      Eigenvalue No.Most Associated States Real Part Imaginary Part Frequency/Hz Damping Ratio 246 & 247 Rotor Speed & Angle Mach.No.40 -0.0920 ±5.1243 0.8157 0.0180 196 & 197 Rotor Speed & Angle Mach.No.18 -0.0792 ±4.0683 0.6476 0.0195 180 & 181 Rotor Speed & Angle Mach.No.1 -0.2087 ±6.4383 1.0252 0.0324 190 & 191 Rotor Speed & Angle Mach.No.20 -0.2041 ±6.0048 0.9562 0.0340 198 & 199 Rotor Speed & Angle Mach.No.37 -0.2381 ±6.1691 0.9826 0.0386

      Table 3 Damping ratios of the Philippine–Sabah grid with increasing penetration of renewable energy

      Mode/Hz Most Associated States Renewable Penetration Basecase 10% 20% 30% 40%0.6209 Rotor Speed & Angle Mach.No.10 0.0075 0.0166 0.0128 0.0133 0.0204 0.8149 Rotor Speed & Angle Mach.No.40 0.0180 0.0180 0.0180 0.0180 0.0180 0.6471 Rotor Speed & Angle Mach.No.18 0.0194 0.0194 0.0194 0.0194 0.0194 1.0079 Rotor Speed & Angle Mach.No.20 0.0252 0.0299 0.0288 0.0259 0.0301 1.0145 Rotor Speed & Angle Mach.No.1 0.0271 0.0268 0.0265 0.0290 0.0287

      Fig.10 Philippine–Sabah model with high penetration of RE voltage profile and summary of power flow

      The base case is a scenario in where the system is installed with renewable energy resources as planned for 2023, and its damping ratios were considered.The damping ratios of the critically damped modes were compared to visualize which machines have been improving during the growth of RE in the system.

      Table 3 has the damping ratios enlisted.Machines 40 and 18 have no visible changes because these machines are located in the Sabah grid with no RE installed, so even the RE penetration increase will not affect the machines,whereas the other three machines have increased their damping ratios considerably with respect to the base case.

      From Table 3, it can be deduced that the increased RE penetration in the system improved the critically damped eigenvalues, probably because of the reconfiguration of power flows among the corridors/lines after installing distributed RE sources at different locations.

      The time domain simulation results were also evaluated by adding more renewable energy sources to the system.Two subcases were considered in this evaluation.First, the capacity of traditional generators of the buses, where the RE sources are installed is reduced.In the second case, the load increased to the same value as RE generation increased in the region.The load was distributed mostly throughout the area where RE generation was located.The main focus of this study was wind and solar penetration.

      A three-phase fault in the most troubled area (Mindanao)was considered to simulate the behavior of the system when subjected to a large disturbance.In the simulation, the fault was applied at 0.2 s, followed by a clearing time of 5 cycles(0.0833 s).

      When the generation of solar and wind power plants is increased by 20%, the system characteristics, particularly the voltage for Visayas (weak area), are shown in Fig.11.Figure 11a shows the Visayas area was stable when the generation from nearby traditional power sources decreased.However, in the second case, the system became unstable,as shown in Fig.11b.The instability was solved in Fig.11c by acknowledging the role of the wind turbines not being appropriately initialized.Appropriate parameters stabilize the system.The magnetic reactances of DIFG 3 and DFIG 4 located at buses 79 and 77, respectively, increased.This shows the sensitivity of the system parameters when the transient stability is considered.Hence, this signifies the role of the initialization of the dynamic components by power flow analysis.

      Fig.11 Philippine–Sabah model voltage per area with 20%increase in penetration of RE

      3 Conclusions and recommendations

      This study presents a stability study of the Philippines’power grid that is connected to the Sabah power grid as part of the planned ASEAN Power Grid.The results show that all three interconnected islands pose stability threats but can be stabilized using a power system stabilizer.In addition,the Sabah interconnection does not significantly challenge the stability of the grid.However, an interconnected power system exhibits instability issues at moderate penetration of renewable sources, which require a stabilization mechanism for a higher penetration level.

      The author recommends for future researchers:

      ·To include other analyses like voltage stability and line contingency.

      ·To have more practical data and the economic aspect of this interconnection.

      ·To expand the scope of the study among neighboring countries (such as Malaysia, Brunei, and Indonesia) to study international cooperation, such as ASEAN as a whole or the Asian Super Grid.

      Acknowledgements

      The authors would like to thank Engr.Edwin B.Cano of the NYISO for the generous support by providing data for the Philippine power system, PSAT program developer, and all concerned persons/organizations contributed either way.

      Declaration of Competing Interest

      We declare that we have no conflict of interest.

      Appendix A

      Line data

      From To Power/MVA Voltage/kV Freq/Hz Resistance/p.u.Reactance/p.u.Line-charge/p.u.Tap Ratio 1 2 600 500 60 0.01668 0.16668 0.00000 2.173913 2 3 100 230 60 0.01398 0.13876 0.23900 3 4 100 230 60 0.00234 0.02312 0.03983 4 16 100 230 60 0.00370 0.03700 0.06375 2 13 100 230 60 0.01120 0.11100 0.19120 4 14 100 230 60 0.00374 0.03700 0.06373 13 5 100 230 60 0.00233 0.02313 0.03984 13 6 100 230 60 0.00466 0.04626 0.07967 6 9 100 230 60 0.00326 0.03238 0.05577 1 7 100 500 60 0.00124 0.01430 0.71500 7 9 600 500 60 0.00417 0.04167 0.00000 2.173913 7 8 100 500 60 0.00314 0.03640 1.82000 8 10 600 500 60 0.02085 0.20835 0.00000 2.173913 8 41 100 500 60 0.00029 0.03380 1.69000 10 16 100 230 60 0.00280 0.02776 0.04780 16 15 100 230 60 0.00037 0.00370 0.00638 16 14 300 230 60 0.00280 0.02775 0.04780 14 15 600 230 60 0.00280 0.02775 0.04780 16 32 100 230 60 0.00037 0.00370 0.00638 32 31 300 230 60 0.02001 0.20001 0.00000 2.056 32 24 100 230 60 0.00093 0.00925 0.01594 24 25 100 230 60 0.00224 0.02220 0.03824 25 27 100 230 60 0.00122 0.01202 0.02071 27 29 100 230 60 0.00168 0.01666 0.02868 29 30 300 230 60 0.02001 0.20001 0.00000 2.056 27 28 300 230 60 0.01334 0.13340 0.00000 2.056 27 26 100 230 60 0.00210 0.02081 0.03585 27 36 600 500 60 0.01251 0.12501 0.00000 2.173913 27 34 100 230 60 0.00420 0.04162 0.07170 34 33 100 230 60 0.00233 0.02313 0.03984 34 25 100 230 60 0.00256 0.02544 0.04382 25 22 100 230 60 0.00280 0.02776 0.04780 22 21 100 230 60 0.00047 0.00463 0.00797 22 20 100 230 60 0.00196 0.01943 0.03346 21 20 100 230 60 0.00163 0.01619 0.02789

      continue

      From To Power/MVA Voltage/kV Freq/Hz Resistance/p.u.Reactance/p.u.Line-charge/p.u.Tap Ratio 20 17 100 230 60 0.00606 0.06012 0.10357 17 19 100 230 60 0.00303 0.03006 0.05179 17 18 300 230 60 0.02668 0.26668 0.00000 2.056 10 17 100 230 60 0.00420 0.04162 0.71700 8 41 100 500 60 0.00029 0.03380 1.69000 41 40 600 500 60 0.01251 0.12501 0.00000 2.173913 40 39 100 230 60 0.01866 0.18500 0.31867 20 39 100 230 60 0.02332 0.23126 0.39834 41 35 100 500 60 0.00090 0.01040 0.52000 35 36 100 500 60 0.00090 0.01040 0.52000 42 47 150 138 60 0.00147 0.00463 0.00642 43 47 200 138 60 0.03676 0.11578 0.01054 47 48 150 138 60 0.00735 0.02316 0.03211 48 49 150 138 60 0.05881 0.18526 0.25686 49 51 100 138 60 0.08822 0.27788 0.03853 51 52 100 138 60 0.00672 0.02395 0.33118 52 53 100 138 60 0.00000 0.08000 0.00000 49 44 100 138 60 0.02941 0.09263 0.01285 44 45 100 138 60 0.02941 0.09263 0.01285 45 46 100 138 60 0.04411 0.13894 0.19265 45 49 100 138 60 0.00882 0.02779 0.03853 49 50 150 230 60 0.00000 0.24000 0.00000 50 54 175 230 60 0.00172 0.01088 0.10487 54 57 200 230 60 0.00086 0.00544 0.05243 57 58 200 230 60 0.00252 0.02912 0.27687 58 55 220 230 60 0.00020 0.17240 1.66043 55 56 200 230 60 0.00000 0.24000 0.00000 56 59 150 138 60 0.00020 0.09262 0.12843 59 60 175 138 60 0.01470 0.04632 0.06421 60 61 200 138 60 0.00294 0.00926 0.01284 59 62 150 138 60 0.01029 0.03242 0.04495 59 63 150 138 60 0.00010 0.04631 0.06422 62 63 150 138 60 0.00010 0.03474 0.04816 63 64 150 138 60 0.01470 0.04631 0.06422 64 65 100 138 60 0.00882 0.02779 0.03853 65 73 100 138 60 0.01470 0.04632 0.06421 63 68 150 138 60 0.01344 0.04790 0.66235 68 67 150 138 60 0.02940 0.09262 0.12843

      continue

      From To Power/MVA Voltage/kV Freq/Hz Resistance/p.u.Reactance/p.u.Line-charge/p.u.Tap Ratio 68 66 100 138 60 0.03676 0.11578 0.16054 68 70 200 138 60 0.07352 0.23156 0.32107 70 71 175 138 60 0.05882 0.18526 0.25686 71 72 200 138 60 0.10292 0.32420 0.44950 72 69 100 138 60 0.04411 0.13894 0.19265 69 78 100 138 60 0.04411 0.13894 0.19265 72 75 175 138 60 0.00672 0.02395 0.33118 75 79 100 138 60 0.04411 0.13894 0.19265 75 76 100 138 60 0.02940 0.09262 0.12843 75 74 100 138 60 0.04411 0.13894 0.19265 74 77 100 138 60 0.01470 0.04631 0.06422 84 83 150 138 60 0.06616 0.20841 0.03027 83 82 150 138 60 0.01470 0.04631 0.01784 82 81 150 138 60 0.06616 0.20841 0.03027 83 81 150 138 60 0.06616 0.20841 0.03027 81 90 150 138 60 0.02206 0.06948 0.02676 83 85 220 138 60 0.12698 0.28641 0.02788 85 86 220 138 60 0.12698 0.28641 0.02788 86 88 150 138 60 0.07351 0.23157 0.03919 88 111 125 138 60 0.03676 0.11578 0.04460 111 112 125 138 60 0.00735 0.02316 0.00892 112 109 150 138 60 0.04043 0.12736 0.04905 109 108 120 138 60 0.02573 0.08105 0.03122 108 88 100 138 60 0.05514 0.17368 0.01689 108 89 200 230 60 0.00000 0.08000 0.00000 87 89 150 230 60 0.03046 0.19236 0.08527 109 110 120 138 60 0.08086 0.25472 0.09810 110 104 120 138 60 0.08086 0.25472 0.09810 104 105 100 138 60 0.00000 0.08000 0.00000 110 103 120 138 60 0.05881 0.18526 0.02135 103 101 100 138 60 0.12498 0.39366 0.05162 105 102 100 69 60 0.43562 0.75877 0.00783 102 100 150 69 60 0.35875 0.62487 0.00645 102 101 100 138 60 0.00000 0.08000 0.00000 99 100 100 138 60 0.00000 0.08000 0.00000 101 107 100 138 60 0.10292 0.32420 0.02486 99 107 150 138 60 0.09557 0.30104 0.01595

      continue

      From To Power/MVA Voltage/kV Freq/Hz Resistance/p.u.Reactance/p.u.Line-charge/p.u.Tap Ratio 107 106 100 138 60 0.03676 0.11578 0.04459 107 109 120 138 60 0.12498 0.21937 0.05162 106 90 100 138 60 0.17644 0.55576 0.01405 106 92 100 138 60 0.15438 0.48630 0.08729 90 92 100 138 60 0.05881 0.18526 0.02135 90 91 100 138 60 0.06616 0.20841 0.03027 91 92 100 138 60 0.02205 0.06947 0.02676 87 80 150 230 60 0.03104 0.19600 0.08876 80 93 150 230 60 0.03448 0.21776 0.00974 93 92 200 230 60 0.00000 0.08000 0.00000 92 94 100 138 60 0.13232 0.41682 0.06054 94 95 150 138 60 0.05808 0.18294 0.02046 94 96 100 138 60 0.18378 0.57892 0.02297 96 97 100 138 60 0.16174 0.50946 0.09621 97 98 100 138 60 0.04926 0.15516 0.05976 113 114 100 132 50 0.68820 0.00651 0.02761 114 115 150 132 50 0.81400 0.00770 0.03266 115 119 100 132 50 1.12480 0.01064 0.04513 119 118 100 132 50 0.44400 0.00420 0.01781 119 117 100 132 50 0.44400 0.00420 0.01781 118 117 100 132 50 0.29600 0.00280 0.01188 119 120 100 132 50 0.75480 0.00714 0.03028 120 121 100 132 50 0.75480 0.00714 0.03028 121 123 100 132 50 1.05080 0.00994 0.04216 123 122 100 132 50 1.01380 0.00959 0.04068 123 125 100 132 50 0.52540 0.00497 0.02108 125 124 100 132 50 0.59200 0.00560 0.02375 116 119 200 275 50 0.00000 0.08000 0.00000 116 135 150 275 50 0.06350 0.00142 0.00579 135 132 200 275 50 0.00000 0.08000 0.00000 132 133 100 132 50 1.67240 0.01582 0.06710 133 130 100 132 50 0.66600 0.00630 0.02672 130 131 100 132 50 0.66600 0.00630 0.02672 133 129 100 132 50 0.07844 0.00074 0.00315 133 128 100 132 50 1.33200 0.01260 0.05344 128 126 100 132 50 0.48840 0.00462 0.01960 128 127 100 132 50 0.74000 0.00700 0.02969

      Bus data

      Vmin/p.u.Bus 01 0.99805 0.036421 -5.9E-16 1.14E-13 0 0 0 0 PV 0 0 1.1 0.9 Bus 02 1 0 -30.5331 5.798304 0.87 0.14516 0 0 Slack 0.316 -0.316 1.1 0.9 Bus 03 0.991063 0.164911 1.328 0.38014 1.00734 0.1679 0 0 PV 0.229 -0.229 1.1 0.9 Bus 04 0.989955 0.190864 3.664 1.11752 1.00734 0.1679 0 0 PV 0.244 -0.244 1.1 0.9 Bus 05 1 0.251538 3.168 -0.22133 0 0 0 0 PV 0.314 -0.314 1.1 0.9 Bus 06 1 0.191305 5.136 -0.30129 0.88386 0.14731 0 0 PV 0.314 -0.314 1.1 0.9 Bus No.Voltage/p.u.Angle/p.u.Pgen/p.u.Qgen/p.u.Pload/p.u.Qload/p.u.Gcond/p.u.BSuc/p.u.Bus type Qgmax/p.u.Qgmin/p.u.Vmax/p.u.Bus 07 0.99596 0.111748 0 0 0 0 0 0 PV 0 0 1.1 0.9 Bus 08 0.991309 0.234186 1.34E-14 -2.6E-14 0 0 0 0 PV 0 0 1.1 0.9 Bus 09 1 0.166075 4.296 0.861511 0.56541 0.09425 0 0 PV 0.312 -0.312 1.1 0.9 Bus 10 0.98869 0.228132 2.072 0.808598 4.02285 0.67046 0 0 PV 0.3122 -0.3122 1.1 0.9 Bus 100 0.95223 -0.19926 1.11E-16 -2E-15 0.3935 0.03655 0 0 PQ 0 0 1.1 0.9 Bus 101 0.970085 -0.21384 -5E-16 2.75E-15 0 0 0 0 PV 0 0 1.1 0.9 Bus 102 0.954722 -0.25365 4.44E-16 -3.6E-16 0.5935 0.13655 0 0 PQ 0 0 1.1 0.9 Bus 103 1 -0.22028 0.472 -0.13617 0.5935 -0.48855 0 0.0525 PV 0.3 -0.3 1.1 0.9 Bus 104 1 -0.16905 0.928 0.144608 0 0 0 0 PV 0.4 -0.4 1.1 0.9 Bus 105 0.986718 -0.22114 1.11E-16 0 0.5935 0.13655 0 0 PQ 0 0 1.1 0.9 Bus 106 1 0.025675 2.136 -0.14235 0 0 0 0 PV 0.8 -0.8 1.1 0.9 Bus 107 0.986546 -0.05276 3.89E-16 5.16E-15 0.31031 0.07139 0 0 PQ 0 0 1.1 0.9 Bus 108 0.983702 -0.17231 1.44E-15 1.73E-15 0.51622 0.11877 0 0 PQ 0 0 1.1 0.9 Bus 109 0.995639 -0.16994 0.256 0.16 0.51622 -0.50176 0 0.0525 PV 0.5 -0.5 1.1 0.9 Bus 11 0.999806 0.183442 1.05 0.32865 0 0 0 0 PV 0.313 -0.313 1.1 0.9 Bus 110 0.992949 -0.19194 6.66E-16 5.16E-15 0.51622 0.11877 0 0 PQ 0 0 1.1 0.9 Bus 111 0.970186 -0.20759 2.62E-15 -1.1E-15 0 -0.04943 0 0.0525 PQ 0 0 1.1 0.9 Bus 112 0.97081 -0.2086 8.88E-16 -3.6E-15 0.51622 0.11877 0 0 PQ 0 0 1.1 0.9 Bus 113 0.990879 -0.45598 -5.6E-16 1.39E-16 0.233184 0.069951 0 0 PQ 0 0 1.1 0.9 Bus 114 1.009319 -0.42555 2.22E-16 -1.6E-15 0.275911 -0.42659 0 0.5 PQ 0 0 1.1 0.9 Bus 115 0.976037 -0.32098 0.064 0.08 0.254415 0.076325 0 0 PV 1 -0.5 1.1 0.9 Bus 116 0.957927 -0.1406 1.78E-15 -1.1E-15 0 0 0 0 PV 0 0 1.1 0.9 Bus 117 0.972971 -0.06423 3.33E-16 2E-15 0.496317 0.148895 0 0 PQ 0 0 1.1 0.9 Bus 118 1 0 3.934089 0.750254 1.552751 0.465825 0 0 Vθ 0.5 -0.5 1.1 0.9 Bus 119 0.963056 -0.11991 -6.7E-16 2.37E-15 0.360673 0.108202 0 0 PQ 0 0 1.1 0.9 Bus 12 0.991506 0.17767 0.104 0.03939 0.65639 0.1094 0 0 PV 0.303 -0.303 1.1 0.9 Bus 120 0.973656 -0.16032 -2.2E-16 -2.7E-15 0.350671 -0.13938 0 0.258 PQ 0 0 1.1 0.9 Bus 121 0.975829 -0.13596 -9E-17 -7.4E-16 0.05304 0.015912 0 0 PQ 0 0 1.1 0.9 Bus 122 1 0.091767 0.836 -0.08348 0 0 0 0 PV 1 -0.5 1.1 0.9 Bus 123 0.982085 -0.09 -1.1E-15 -2.5E-15 0.278961 0.083688 0 0 PQ 0 0 1.1 0.9 Bus 124 0.955599 -0.20043 0 8.05E-16 0.584098 0.175229 0 0 PQ 0 0 1.1 0.9

      continue

      Vmin/p.u.Bus 125 1 -0.1323 0.528 0.279663 0.242755 -0.18518 0 0.258 PV 0.5 -0.5 1.1 0.9 Bus 126 1 -0.76904 1.1968 0.438445 1.112593 0.33778 0 0 PV 0.5 -0.5 1.1 0.9 Bus 127 0.958433 -0.82476 0.0408 0.0255 0.37643 0.112929 0 0 PV 0.5 -0.5 1.1 0.9 Bus 128 0.986949 -0.77452 0.08 0.1 0 0 0 0 PV 1 -0.5 1.1 0.9 Bus 129 0.962222 -0.71511 3.44E-15 -2.5E-15 0.171606 0.051482 0 0 PQ 0 0 1.1 0.9 Bus 13 1 0.177726 4.176 0.682779 2.11866 0.35311 0 0 PV 0.196 -0.196 1.1 0.9 Bus 130 0.97814 -0.76282 0.2704 0.338 0.564805 0.169442 0 0 PV 1 -0.5 1.1 0.9 Bus No.Voltage/p.u.Angle/p.u.Pgen/p.u.Qgen/p.u.Pload/p.u.Qload/p.u.Gcond/p.u.BSuc/p.u.Bus type Qgmax/p.u.Qgmin/p.u.Vmax/p.u.Bus 131 0.979311 -0.76316 -2.7E-17 1.23E-15 0 0 0 0 PV 0 0 1.1 0.9 Bus 132 0.959769 -0.43669 1.44E-15 -5.1E-15 0.701591 0.210477 0 0 PQ 0 0 1.1 0.9 Bus 133 0.963871 -0.71238 -5.2E-15 4.48E-15 0 0 0 0 PV 0 0 1.1 0.9 Bus 134 1 -0.41494 1.336 0.561667 1.116871 0.335061 0 0 PV 1 -0.5 1.1 0.9 Bus 135 0.954955 -0.38464 2.85E-15 2.51E-15 0 0 0 0 PV 0 0 1.1 0.9 Bus 14 0.98054 0.183946 8.44E-15 -3.5E-14 2.82054 0.4701 0 0 PQ 0 0 1.1 0.9 Bus 15 0.976474 0.186913 -3.2E-14 1.11E-13 3.7174 0.61954 0 0 PQ 0 0 1.1 0.9 Bus 16 0.976585 0.203531 3.456 0.72144 8.68909 1.44818 0 0 PV 0.167 -0.167 1.1 0.9 Bus 17 0.987464 0.264159 5.2 2.002 0 0 0 0 PV 0.308 -0.308 1.1 0.9 Bus 18 0.969554 0.20355 3.55E-15 -5.5E-14 10.61278 1.747431 0 0.0227 PQ 0 0 1.1 0.9 Bus 19 1 0.322765 1.952 0.250768 0 0 0 0 PV 0.292 -0.292 1.1 0.9 Bus 20 1 0.369447 5.472 1.449435 0.39644 0.06607 0 0 PV 0.292 -0.292 1.1 0.9 Bus 21 0.992802 0.365819 7.11E-15 -5.6E-14 1.26729 0.21123 0 0 PQ 0 0 1.1 0.9 Bus 22 0.992259 0.370639 3.584 1.39776 1.35828 0.2264 0 0 PV 0.312 -0.312 1.1 0.9 Bus 23 0.960453 0.290315 7.11E-15 1.99E-15 6.26499 1.040673 0 0.00378 PQ 0 0 1.1 0.9 Bus 24 0.974473 0.339097 5.01E-15 8.99E-14 0 0 0 0 PV 0 0 1.1 0.9 Bus 25 0.978648 0.363372 -4.1E-14 -4.7E-14 2.36561 0.39424 0 0 PQ 0 0 1.1 0.9 Bus 26 0.985586 0.393968 0.352 0.13728 0 0 0 0 PV 0.312 -0.312 1.1 0.9 Bus 27 0.981595 0.386732 9.19E-14 2.18E-13 0 0 0 0 PV 0 0 1.1 0.9 Bus 28 0.964804 0.329641 -1.1E-14 -4.2E-15 4.94569 0.82429 0 0 PQ 0 0 1.1 0.9 Bus 29 0.977283 0.376306 -3E-15 -1.1E-13 0 0 0 0 PV 0 0 1.1 0.9 Bus 30 0.969791 0.357541 -1.8E-15 3.33E-16 2.47177 0.711 0 0 PQ 0 0 1.1 0.9 Bus 31 0.967924 0.219478 -4.9E-15 -6E-15 2.16415 0.36068 0 0 PQ 0 0 1.1 0.9 Bus 32 0.972477 0.236226 2.74E-14 6.34E-14 0 0 0 0 PV 0 0 1.1 0.9 Bus 33 0.980082 0.378493 6.08E-14 1.67E-14 0 0 0 0 PV 0 0 1.1 0.9 Bus 34 1 0.501985 23.328 3.024887 0.91635 0.15275 0 0 PV 0.312 -0.312 1.1 0.9 Bus 35 1 0.442955 10.104 0.640204 0 0 0 0 PV 0.3122 -0.3122 1.1 0.9 Bus 36 0.987113 0.395227 1.08E-14 -2.6E-14 0 0 0 0 PV 0 0 1.1 0.9 Bus 39 1 0.494562 6.04 -0.06906 1.31279 0.21882 0 0 PV 0.165 -0.165 1.1 0.9 Bus 40 1 0.411761 10.04 1.562036 0.45493 0.07579 0 0 PV 0.317 -0.317 1.1 0.9 Bus 41 0.994483 0.385666 1.37E-13 -1.9E-13 0 0 0 0 PV 0 0 1.1 0.9

      continue

      Bus No.Voltage/p.u.Vmin/p.u.Bus 42 1.017322 -0.22579 -3.9E-14 -7.4E-14 0.14481 0.1 0 0 PQ 0 0 1.1 0.9 Bus 43 1.033187 -0.28713 -1.1E-15 -3.1E-15 0.14481 -0.96747 0 0 PQ 0 0 1.1 0.9 Bus 44 0.999643 -0.03967 8.33E-17 -4.4E-15 0.14481 0.1 0 0 PQ 0 0 1.1 0.9 Bus 45 1.003871 -0.03528 1.55E-15 8.69E-15 0.288 0.2 0 0 PQ 0 0 1.1 0.9 Bus 46 1.014534 -0.00277 0.24 -0.0936 0 0 0 0 PV 0.312 -0.312 1.1 0.9 Bus 47 1.017749 -0.22545 5.45E-15 8.04E-14 0 0 0 0 PV 0 0 1.1 0.9 Bus 48 1.015851 -0.20682 9.58E-15 3.33E-15 0.14481 0.1 0 0 PQ 0 0 1.1 0.9 Bus 49 1.007782 -0.03327 -5E-16 -4.5E-15 0.14481 0.1 0 0 PQ 0 0 1.1 0.9 Bus 50 1.000055 -0.00162 1.7E-14 9.65E-14 0 0 0 0 PV 0 0 1.1 0.9 Bus 51 1.034016 -0.08099 -1.9E-15 5.77E-15 0.14481 0.1 0 0 PQ 0 0 1.1 0.9 Bus 52 1.035169 -0.08103 9.85E-16 3.61E-15 0 0 0 0 PV 0 0 1.1 0.9 Bus 53 1.025055 -0.08014 0.256 -0.0224 0.2443 0.10719 0 0 PQ 0.07 -0.07 1.1 0.9 Bus 54 1 0 -9.97384 -0.86165 0 0 0 0 Vθ 0.411 -0.411 1.1 0.9 Bus 55 1.014524 0.273637 5.44E-15 6E-15 0 -1.13218 0 1.1 PQ 0 0 1.1 0.9 Bus 56 1.008171 0.426824 -1.7E-14 1.78E-14 0.50683 -3.08061 0 1.1 PQ 0 0 1.1 0.9 Bus 57 1.001564 0.007402 -2.3E-14 -1.9E-13 0.14481 0.1 0 0 PQ 0 0 1.1 0.9 Bus 58 1.00724 0.047499 4E-15 4.84E-14 0.50683 0.28 0 0 PQ 0 0 1.1 0.9 Bus 59 0.98425 0.618058 -3E-14 1.08E-14 0.50683 0.27354 0 0 PQ 0 0 1.1 0.9 Bus 60 0.993827 0.637115 7.56E-14 1.96E-13 0.50683 0.27354 0 0 PQ 0 0 1.1 0.9 Bus 61 0.995968 0.640925 4 1 0.50683 0.27354 0 0 PV 0.2 -0.2 1.1 0.9 Bus 62 0.998715 0.705833 1.78E-15 -3.3E-15 0.50683 0.27354 0 0 PQ 0 0 1.1 0.9 Bus 63 1 0.807408 4.496 2.430892 0.50683 0.27354 0 0 PV 1.57 -1.57 1.1 0.9 Bus 64 1 0.890302 1.008 -0.20151 0 0 0 0 PV 0.4 -0.4 1.1 0.9 Bus 65 0.998953 0.935699 -5E-15 -2.1E-14 0.14486 0.07819 0 0 PQ 0 0 1.1 0.9 Bus 66 1.007707 0.858189 0.672 -0.168 0 0 0 0 PV 0.2 -0.2 1.1 0.9 Bus 67 1.005382 0.855883 0.936 -0.351 0 0 0 0 PV 0.3 -0.3 1.1 0.9 Bus 68 1.004836 0.84035 1.51E-14 2.48E-14 0 0 0 0 PV 0 0 1.1 0.9 Bus 69 1.028164 1.227105 0.296 -0.0259 0 0 0 0 PV 0.07 -0.07 1.1 0.9 Bus 70 0.999241 0.919991 -2.7E-15 1.27E-14 0.4719 0.25468 0 0 PQ 0 0 1.1 0.9 Bus 71 0.997209 1.008387 -4.4E-16 -1.9E-15 0.4719 0.25468 0 0 PQ 0 0 1.1 0.9 Bus 72 1.002051 1.168356 -2.9E-15 -1.1E-14 0.4719 0.25468 0 0 PQ 0 0 1.1 0.9 Bus 73 1 0.956191 2.144 -0.16957 0.50683 0.27354 0 0 PV 1 -1 1.1 0.9 Bus 74 1.001867 1.299635 0.456 -0.0627 0.34974 0.020899 0 0.3 PV 0.11 -0.11 1.1 0.9 Bus 75 1 1.224193 1.296 -0.13998 0.84974 0.32202 0 0 PV 0.4 -0.4 1.1 0.9 Bus 76 1 1.276678 2.672 -0.34152 0.59658 0.32202 0 0 PV 0.71 -0.71 1.1 0.9 Bus 77 1 1.320512 0.4 -0.19467 0 0 0 0 PV 1 -1 1.1 0.9 Bus 78 1.033283 1.248487 0.16 -0.0624 0 0 0 0 PV 0.312 -0.312 1.1 0.9 Bus 79 1 1.349287 0.832 -0.30422 0 0 0 0 PV 1 -1 1.1 0.9 Angle/p.u.Pgen/p.u.Qgen/p.u.Pload/p.u.Qload/p.u.Gcond/p.u.BSuc/p.u.Bus type Qgmax/p.u.Qgmin/p.u.Vmax/p.u.

      continue

      Bus No.Voltage/p.u.Vmin/p.u.Bus 80 0.99119 -0.09567 1.78E-15 1.04E-14 0.31031 0.07139 0 0 PQ 0 0 1.1 0.9 Bus 81 1 0.134898 1.664 -0.10629 0.2821 0.0649 0 0 PV 0.6 -0.6 1.1 0.9 Bus 82 0.978803 0.087836 -1.9E-15 -1.9E-15 0.25515 0.0587 0 0 PQ 0 0 1.1 0.9 Bus 83 0.977658 0.085129 1.11E-16 6.18E-16 0.25515 0.0587 0 0 PQ 0 0 1.1 0.9 Bus 84 0.979418 0.181018 -4.4E-16 -8.9E-16 0.15515 -0.12437 0 0.17 PQ 0 0 1.1 0.9 Bus 85 0.95636 -0.03316 -1.7E-15 -2.3E-15 0.18067 0.06457 0 0 PQ 0 0 1.1 0.9 Bus 86 0.95831 -0.12858 -9.4E-16 2.22E-15 0.15515 -0.11742 0 0.17 PQ 0 0 1.1 0.9 Bus 87 0.987921 -0.12682 -4.3E-15 -3.2E-15 0.31031 0.07139 0 0 PQ 0 0 1.1 0.9 Bus 88 0.959377 -0.19993 -1.1E-16 2.53E-15 0.51622 0.11877 0 0 PQ 0 0 1.1 0.9 Bus 89 0.985164 -0.14707 4.33E-16 -6.4E-17 0 0 0 0 PV 0 0 1.1 0.9 Bus 90 1 0.126398 1.44 -0.18492 0.31031 0.07139 0 0 PQ 0.6 -0.6 1.1 0.9 Bus 91 0.99092 0.015952 -6.1E-16 3.21E-15 0.31031 0.07139 0 0 PQ 0 0 1.1 0.9 Bus 92 1 0 1.660807 0.585516 0.31031 0.07139 0 0 Vθ 0.2211 -0.2211 1.1 0.9 Bus 93 0.998149 -0.05002 8.81E-17 -7.1E-15 0 0 0 0 PV 0 0 1.1 0.9 Bus 94 0.965712 -0.11379 -6.1E-16 -1.7E-15 0.39362 0.09057 0 0 PQ 0 0 1.1 0.9 Bus 95 0.950056 -0.16671 3.33E-16 5.41E-16 0.39362 0.000309 0 0.1 PQ 0 0 1.1 0.9 Bus 96 0.972405 -0.15582 1.11E-16 1.82E-15 0.39362 0.09057 0 0 PQ 0 0 1.1 0.9 Bus 97 1 -0.14139 0.952 0.245736 0.39362 0.09057 0 0 PQ 0.3 -0.3 1.1 0.9 Bus 98 0.992621 -0.15602 9.44E-16 2.01E-15 0.39362 0.09057 0 0 PQ 0 0 1.1 0.9 Bus 99 0.951333 -0.15671 3.19E-16 3.37E-15 0 0 0 0 PV 0 0 1.1 0.9 Angle/p.u.Pgen/p.u.Qgen/p.u.Pload/p.u.Qload/p.u.Gcond/p.u.BSuc/p.u.Bus type Qgmax/p.u.Qgmin/p.u.Vmax/p.u.

      References

      1. [1]

        Atchatavivan P (2005) ASEAN energy cooperation: an opportunity for regional sustainable energy development.In: Mimeo, Papers on International Environmental Negotiation series, MIT-Harvard Public Disputes Program, Harvard Law School [百度学术]

      2. [2]

        Montaño A M, Ruiz E Y, Morataya J C, et al.(2006) Stability analysis for the electrical integration of Ecuador, Colombia and Panama.Proceedings of IEEE PES Transmission and Distribution Conference and Exposition Latin America.Caracas, Venezuela,2006 [百度学术]

      3. [3]

        Kundur P (1994) Power system stability and control.McGraw-Hill, New York, NY, USA [百度学术]

      4. [4]

        AIMS II Report (2011) ASEAN interconnection master plan study-2, http://www.hapuasecretariat.org/doc2011/Vol_I_Executive%20Summary.pdf.Accessed Nov 2015 [百度学术]

      5. [5]

        Cano E B, Shaikh F A (2013) Reduced and approximate models of philippine major island power grids.Proceedings of International Conference on Power, Energy and Control (ICPEC),2013 [百度学术]

      6. [6]

        Transmission development plan (2013) Final report, national grid corporation of the Philippines https://www.ngcp.ph/savefiles/news/BULLETINS/2013%20 TDP%20Consultation%20Draft_VolI-Major%20Network%20 Development.pdf.Accessed Oct 2015 [百度学术]

      7. [7]

        Marbun M P, Salile A Y, Surya A S (2018) Grid impact study of variable renewable energy integration to Java-Bali system.Proceedings of 2018 8th International Conference on Power and Energy Systems (ICPES), pp 182-189, doi: 10.1109/ICPESYS.2018.8626922 [百度学术]

      8. [8]

        Mararakanye N, Bekker B (2019) Renewable energy integration impacts within the context of generator type, penetration level and grid characteristics, Renewable and Sustainable Energy Reviews, 108: 441-451 [百度学术]

      9. [9]

        Beyza J, Yusta J M (2021) The effects of the high penetration of renewable energies on the reliability and vulnerability of interconnected electric power systems, Reliability Engineering &System Safety, 215: 107881 [百度学术]

      10. [10]

        Erdiwansyah, Mahidin, Husin H, et al.(2021) A critical review of the integration of renewable energy sources with various technologies.Protection and Control of Modern Power Systems,6(1): 1-18 [百度学术]

      11. [11]

        Smith O, Cattell O, Farcot E, et al.(2022) The effect of renewable energy incorporation on power grid stability and resilience.Science Advances, 8(9), doi:10.1126/sciadv.abj6734 [百度学术]

      12. [12]

        International Energy Agency (2021) Countries and regions.https://www.iea.org/.Accessed 25 Dec 2022 [百度学术]

      13. [13]

        Li J, Xu J, Tan X (2018) Dynamic comprehensive benefit evaluation of the transnational tower grid interconnection project based on combination weighting and TOPSIS Grey projection method.Sustainability, 10(12): 4672 [百度学术]

      14. [14]

        Brinkerink M, Gallachóir B Ó, Deane P (2019) A comprehensive review on the benefits and challenges of global power grids and intercontinental interconnectors.Renewable and Sustainable Energy Reviews, 107: 274-287, https://doi.org/10.1016/j.rser.2019.03.003 [百度学术]

      15. [15]

        Alzakkar A, Vladimirovich M V, Samofalov Y, et al.(2022) The impact of electrical interconnection between countries on the stability of electrical power systems.Proceedings of 2022 4th International Youth Conference on Radio Electronics, Electrical and Power Engineering (REEPE), Moscow, Russian Federation,pp.1-6, doi: 10.1109/REEPE53907.2022.9731442 [百度学术]

      16. [16]

        Pham T N, Trinh H, Hien L V (2016) Load frequency control of power systems with electric vehicles and diverse transmission links using distributed functional observers.IEEE Trans Smart Grid, 7(1): 238-252, doi: 10.1109/TSG.2015.2449877 [百度学术]

      17. [17]

        Chen G, Li Z, Zhang Z, et al.(2020) An improved ACO algorithm optimized fuzzy PID controller for load frequency control in multi area interconnected power systems.IEEE Access,8: 6429-6447, doi: 10.1109/ACCESS.2019.2960380 [百度学术]

      18. [18]

        Krishan R, Verma A, Prasad B (2014) Small signal stability analysis of grid connected distributed PV and wind energy system.Proceedings of 6th IEEE Power India International Conference (PIICON), Delhi, India, 5-7 December 2014, p 6 [百度学术]

      19. [19]

        WESM Market Network Model (WESM-MNM).http://www.wesm.ph/inner.php/downloads/market_network_model,Accessed 23 Feb 2016 [百度学术]

      20. [20]

        Milano F (2010) Power system modelling and scripting.Springer,London, Dordrecht Heidelberg, New York [百度学术]

      21. [21]

        Alhelou H H, Golshan M E H, Hatziargyriou N D (2020)Deterministic dynamic state estimation-based optimal LFC for interconnected power systems using unknown input observer.IEEE Transactions on Smart Grid, 11(2): 1582-1592, doi:10.1109/TSG.2019.2940199 [百度学术]

      22. [22]

        Stevenson W D (1984) Elements of power system analysis.4th edition, McGraw-Hill, New York, pp 428 [百度学术]

      23. [23]

        Sabah Electricity Supply Industry Outlook (October 2014)Suruhanjaya Tenaga Energy Commission.http://www.st.gov.my/index.php/component/k2/item/607-sabah-electricity-supplyindustry-outlook-2014.html.Accessed 25 Jan 2016 [百度学术]

      24. [24]

        Lin M K (2014) Small signal stability assessment of MEPE test system in free and open source software.International Journal of Electrical, Computer, Electronics and Communication Engineering, 8(11): 158-166 [百度学术]

      25. [25]

        Chompoobutrgool Y, Li W, Vanfretti L (2012) Development and implementation of a Nordic grid model for power system small-signal and transient stability studies in a free and open source software.Proceedings of IEEE Power and Energy Society General Meeting, San Diego, CA., 22-26 July 2012, p 8 [百度学术]

      26. [26]

        Asian Development Bank (2003) Project completion report on the Leyte-Mindanao interconnection engineering project (Loan 1474-PHI) in the Philippines.http://www.adb.org/sites/default/files/project-document/70576/pcr-phi-30007.pdf.Accessed 15th Jan 2015 [百度学术]

      27. [27]

        Energy Regulatory Commission (2014) Implementing single outage contingency (N-1) operational criterion.www.erc.gov.ph/Files/Render/media/N-1OperationalCriteria.pdf.Acccessed 12 March 2016 [百度学术]

      Fund Information

      The authors would like to thank Engr. Edwin B. Cano of the NYISO for the generous support by providing data for the Philippine power system, PSAT program developer, and all concerned persons/organizations contributed either way.

      Author

      • Mr. Tristan G Magallones Jr

        Mr. Tristan G Magallones Jr received his Master of Engineering degree, Major in Energy at Asian Institute of Technology (AIT), Klong Luang, Pathumthani, Thailand, 2016. He received his Bachelor of Science in Electrical Engineering at Central Mindanao University (CMU), Musuan, Maramag Bukidnon, Philippines last 2011. He is working towards his PhD in Energy Engineering degree at University of the Philippines Diliman, Quezon City, Philippines. Currently, he is also working as a faculty member at Central Mindanao University in Bukidnon, Philippines. His research interests include power system stability studies and renewable energy integration studies.

      • Jai Govind Singh

        Dr. Jai Govind Singh received a Ph.D. degree from IIT Kanpur. Currently, he is working as an Associate Professor at the Department of Energy, Environment and Climate Change, SERD, Asian Institute of Technology (AIT), Bangkok, Thailand. Prior to joining AIT, he worked as a Postdoctoral Research Associate at the KTH Royal Institute of Technology (KTH), Stockholm, Sweden, and as a postdoctoral research fellow at the University of Queensland, Brisbane, Australia. So far, he has supervised seven doctoral and 70 master’s students and has published more than 100 research papers in peer-reviewed international journals and conference proceedings. His teaching and research interests include power system planning, operation, and control, smart grids and microgrids; electric vehicles and battery storage, deregulation, solar and wind integration, and power distribution system planning. He is a fellow of The Institution of Engineers, India.

      Publish Info

      Received:2023-01-07

      Accepted:2023-05-03

      Pubulished:2023-06-24

      Reference: Mr. Tristan G Magallones Jr,Jai Govind Singh,(2023) Impact of interconnections and renewable energy integration on the Philippine–Sabah Power Grid systems.Global Energy Interconnection,6(3):253-259.

      (Editor Yajun Zou)
      Share to WeChat friends or circle of friends

      Use the WeChat “Scan” function to share this article with
      your WeChat friends or circle of friends