Evaluation of Latvian Power System Static Stability According to a New Development Strategy until 2025

: The current stage of world energy development is characterized by the creation of pow-erful, territorially unified energy systems under free market conditions as well as renewable energy integration. Such transactions challenge networks’ configuration s and their operating modes. Electric power system (EPS) static stability could be considered one of the primary targets in this regard. The aim of this study was to evaluate the undergoing and expected modernization of the Latvian EPS in terms of static stability in a regional scale-based case study. In order to define static stability, a method based on a bus admittance matrix was proposed. A simplified Latvian EPS model was developed and assessed for the following three modes: past (2017), current (2020), and planned (2025), taking into account the Baltic States ’ planned development scenarios including large wind farm integration. The evaluation of computation results provided an opportunity to visually inspect changes in the EPS ’ sensitive elements such as lines and nodes. As a result, positive changes were observed (decrease of several weak points) for the planned mode (2025) as compared to the past mode (2017) under the considered modernization scenario. The detailed analysis and results are presented in the article.


Introduction
Electric power systems are complex objects of dispatching and management, requiring compliance with restrictions related to the reliability, safety, quality, and economy of their operating modes. These requirements must be fulfilled by all participants that make up the technological basis of the electricity market. Such concern is essential for an expansion of existing power systems and the modernization process for future, built-up, large energy systems, moving toward their unification in terms of operation and regulations/legalization within certain regions. Thus, the development of market relations in the power industry drives the necessity of the evaluation of network operation modes, considering their multifactor analysis, and forecasting, leading to market relations becoming a key element in ensuring the smooth functioning of the electricity market as well as the integration of renewable energy resources (RES). Ensuring the reliability of the evaluated modes must be considered one of the priority tasks of operational dispatch management.
Currently, the European Union (EU) is updating its energy policy to facilitate the transition to clean energy and improve energy security. Because Baltic region countries (including the Republic of Latvia) are members of the EU, they must gradually approach their energy and economic policies more similarly to the EU. For example, some reforms are mentioned in [1]. For instance, the EU strategy focuses on the support and protection of investments in the energy sector including RES consideration and energy efficiency patterns and contributes to the secure transit of energy supplies [2]. In this way, the aforementioned energy transition must be evaluated and adapted on the regional level in terms of different aspects including vulnerability (stability studies must be conducted in this regard). The implementation of shared objectives and policies must be accelerated and achieved. For instance, the projects of common interest (PCIs) serve as tools aimed to link the EPSs of the EU countries. A detailed list can be found in [3]. Despite a positive tendency toward achieving set targets, some geopolitical issues still exist. For example, the Nord Stream 2 pipeline project, the twin pipeline of around 1230 km that stretches through the Baltic Sea from Russia to Germany [4], is considered to be having a great impact leading to the reduction of Europe's energy security as well as promoting certain consequences for regional integration in energy security. A study of Southeast Europe's energy transition in terms of energy security, sustainability, resilience, and vulnerability is presented in [5].
In accordance with the agreement between the leaders of the Baltic region, Polish and Baltic transmission system operators (TSOs), a formal procedure for expanding the continental European energy network has been started that includes the Baltic States, which will be managed by ENTSO-E (European Network of Transmission System Operators) [6,7].
Considering the priorities of the EU energy strategy [8] and the impacts of the expected reforms including Latvian energy development guidelines [9], let us highlight the primary targets including:

•
Increasing energy efficiency; • Developing the electricity market; • Increasing the reliability of the power supply; • Using and integrating RES. For instance, there is a possibility of increasing wind energy share from the existing 8% to 20% when considering the technical and economic factors of the Republic of Latvia, thus increasing the competitiveness of energy prices; • Preserve the environment while producing, transmitting, and distributing different types of energy.
To address the above-highlighted targets, one important considerations must be devoted to power system stability issues. Stability analysis of any power system leads to defining its performance and weakness indicators in terms of current and planned operation modes. In that turn, it justifies following up on the feasibility of established strategies, regulations, and planned implementation measures for adaptation to the examined power network.
As it is known, the stable operation of a power network depends on its capacity for continuous provision of the generated and consumed power balance as well as the required quality level of electric power. In the past, the stability of proposed systems has been estimated by extrapolating the experience obtained from existing systems and using the rule-of-thumb method. In the future, however, more precise computation methods for predicting and evaluating stability will be required [10] that consider the prospects of growth of the share of renewables.
Usually, a list of the most "dangerous" systems is obtained by simulation of single or joint failures, using a complete or simplified mathematical model of the EPS based on expert judgements and heuristics. For instance, the vulnerability of the EPS may be estimated on the basis of complex systems' general theory [11]. Using this approach, in compliance with the overload level (by comparison with other lines), only one weak point (or branch) of the network, the losing of which brings about a relative decrease in the operation and general efficiency of the whole network for a 34-node diagram, will be determined. Moreover, different optimization techniques have been applied to access power system vulnerability; for instance, a genetic algorithm presented in [12] defined critical system components, a dynamic vulnerability assessment described in [13] used the coherency concept and a fuzzy clustering algorithm considering large power networks, and a nonsequential Monte Carlo approach was applied for power system analysis with protection failure, shown in [14].
Another approach [15] provided a formalized technique for analysis of cascade-developing emergency processes, which allowed obtaining a fairly adequate depiction of the actual events, states, and processes in the EPS that appeared in the cascade development of the accident but did not allow analyzing emergency situations before they occurred. A lot of other studies have been devoted to cascading failure and blackouts analysis; for example, using fault chain theory in power transmission networks [16], a cascading failure graph-based model taking into account continuous transmission line temperature impact [17], a PageRank-based rapid screening method for large-scale power grids [18], a stochastic model under uncertain events (renewable generation and load variation) [19], using new metrics such as the critical moment for a DC power flow model and transient stability analysis [20], and a cascading failure model presented in [21], which included both the structural vulnerability and the state vulnerability of a power network.
A further study [22] aimed to identify weak points and allowed finding the critical (weak) sections (by static stability), based on simplified models. The disadvantage of this approach was the difficulty of analysis and systematization of the results of numerous calculations. Some other research works are related to the same issue; for instance, the geodesic vulnerability approach in [23] explains how to identify power network critical buses, and an adjacent graph is used in [24] to indicate branch importance and at the same time, define fault adjacent relationships between branches. Moreover, to assess the stable operation of a power system, a whole range of different software packages have been designed, both market available and research-focused, which can be used to model a variety of operating modes of power systems.
All the above-mentioned factors highlight the importance of stability analysis for power system operation modes.
This paper aimed to evaluate the Latvian power network in terms of static stability issues for three operation modes by considering their models as follows: past (before modernization, according to system infrastructure in 2017), current (after some modernization was implemented and in particular, the "Kurzeme Ring" project's second and third stages, wherein newly constructed 110 kV and 330 kV transmission lines expanded the existing transmission network, laid out in 2020), and planned (the expected modernization includes new interconnections of transmission lines and the construction of new large wind parks according to the established plan by 2030 [25], when the expected break of continental Europe from the energy ring of BRELL (Belarus, Russia, Estonia, Latvia, and Lithuania) will be fully realized.
An overview of Latvian energy system specifications is presented in Section 2. Next follows the theoretical background of the proposed approach in Section 3. The obtained results and analysis are given in Section 4. A brief discussion of highlights is in Section 5. The final section, Section 6, draws the conclusions.

An Overview of Latvian Energy System Specifics
Latvia (officially the Republic of Latvia) is a sovereign state in the Baltic region of Northern Europe. Since the restoration of independence (1990), Latvia has been referred to as one of the Baltic States. It is bordered by Estonia in the northern region, Lithuania in the southern; to the east is Russia, Belarus is to the southeast, and it shares a maritime border with Sweden to the west. Latvia has about 1,900,000 inhabitants and a territory of 64,589 km 2 . The country has a temperate seasonal climate.
The Latvian transmission power system historically developed as a deficit power system and then developed the transmission system to a greater extent, paying less attention to the construction of new generating capacities because there was a constant opportunity to buy necessary volumes of electricity from neighboring power systems. Currently, the Latvian transmission power system operates in a single energy ring with the energy systems of Belarus, Lithuania, Russia, and Estonia (BRELL). However, in view of the existing economic and political situation, all three Baltic States (BALTSO) are expected to exit BRELL, joining the Union for the Coordination of Transmission of Electricity (UCTE) (Figure 1) [26].

Figure 1.
Preliminary operation scheme for the Baltic countries' system (BALTSO) in an isolated model [26].
Latvia depends on the import of primary energy resources. Without fossil resources, Latvia is highly dependent on oil and gas imports. Hydropower and gas provide almost all the domestic supply of electricity and in recent years have been supplemented by wind energy and energy derived from biomass. Self-sufficiency in energy supply has reached 40% and therefore, both the security of supply and the liberalization of the energy market are vital. Historically, Latvia has been strong in the renewable energy sector due to its high percentage of hydropower use (Table 1). More than a third of total energy consumption in Latvia comes from RESs, and Latvia ranks second in terms of renewable energy consumption among EU countries, but mainly, its strength is related to the significant role of water energy resources (60% of total electricity production; 92% of electricity derived from renewable sources is produced in hydroelectric power stations) and a large amount of biomass ( Figure 2) [27].  The most significant local energy resources used are fuelwood and hydro energy (the Daugava Hydroelectric Power Plant (HPP) cascade). Solid fuel, oil products, and electricity are imported from several countries and supply regions, but there is only one supplier for natural gas-Russia.
Overview of the Baltic States' electricity generation, consumption, and average Nord Pool electricity exchange prices in 2020 are given in Figure 3 and Table 2. As shown, in Latvia, variable hydro resources influence the amount of electricity produced by the Daugava HPP cascade. Because Latvian power plants are not capable of producing the required amount of electricity throughout the whole year, a shortage in electricity supply (or sale of surplus produced electricity in full-water periods of the year) is realized using interconnections. According to Nord Pool information in 2020, Sweden (via NordBalt cable) exported to Lithuania 4,559,426 MWh and to Poland (through the Lit-Pol Link), 458,046 MWh, while Finland exported to Estonia 6,438,480 MWh. Additionally, 3,949,672 MWh were supplied from third-party countries (Belarus and Kaliningrad (Russia)) through Nord Pool Lithuanian and Latvian trade areas. As result, the total consumption of electricity was 7,135,520 MWh in 2020, 2.2% lower than consumption in 2019 [29]. Table 2. Electricity generation and import/export, consumption indicators in Latvia, 2020 [29]. Using local generation, Latvia covered its electricity consumption by 77%. In 2020, primary energy consumption was 5.4 Mtoe, but final energy consumption was 4.5 Mtoe (Table 3) [33]. Table 3. Absolute level of energy consumption in 2020 (Mtoe) in EU countries [33]. * Adjusted with an estimate for international aviation energy consumption of 6.5 Mtoe.

Theoretical Background of the Proposed Approach
Computation methods used for practical purposes require experimental verification. To be able to assess rapidly the performability, efficiency, and reliability of new equipment and new devices; test new theories and methods; and evaluate the assumptions and premises of already existing methods, it is necessary to have the ability to test experiments in the system, allowing for simulating many modes and reproducing any accidents as well as observing and identifying their consequences. In the real power system, the possibilities of conducting experiments are very limited. All observations of "naturally occurring" accidents or in rare cases, accidents "specially produced" in systems, must be processed and presented in criterial dependencies with the application of data from a single (in essence, unique) experiment to a number of (class) similar systems and conditions. However, the difficulty and fundamentally limited possibilities of obtaining such results lead to the need to create an experimental system. It is obvious that only a model system can be used as such, in connection with the various analogue and computing devices finding extensive use in the electric power industry.
The construction of reliable mathematical models of power systems provides the possibility of controlling the electric power system (EPS) regime based on a single approach of structural analysis, simultaneously associated with the widest possible set of regime indicators and electrical system characteristics.
Application of this approach allows solving the actual management tasks related to: • Determination of the admissibility of adding new loads to the existing power system and the need to modernize the energy system to meet current electricity demand; • Identification of power system weak connections with regimes closest to the stability limit; • Visibility of marked nodes where installation of active mode control devices is rational, leading to increased reliability and economy in the transmission of electricity.
The technological processes of electricity production and transmission are connected with many other processes in industry, transport, and other life sectors. Therefore, the implementation of economic development plans in Latvia should provide the possibility for construction of new electric power facilities, including high-voltage interconnections. Currently, JSC Augstspriegumu Tikls is implementing several projects: for instance, the third 330 kV interconnection between Estonia and Latvia, the 330 kV transmission line between Riga TEC-2 and Riga HES, and synchronization of the Baltic States with continental Europe, as shown in Figure 4.
Furthermore, when addressing issues related to the safety of the energy supply, safety aspects should be carefully evaluated, both nationally and regionally, to ensure a flexible and secure energy supply on a national scale.
Analysis of the regime characteristics of existing overhead power lines with a nominal voltage of 330 kV requires specification of the limiting regimes. Thus, the carrying capacity is not determined by all indicators characterizing the stable operation of power transmission lines. Reserve coefficients of static stability in some cross sections are lower than normatively permissible. Regulatory properties of the receiving power system are not fully considered. In these conditions, let us consider the influence of various factors affecting the ultimate power transmission modes and the possibility of improving carrying capacity.
For certain values of both the regime and system parameters, the limiting regimes in the cross sections can be estimated based on a criterion analysis. In the first approximation, the limiting regimes, obviously, are determined by the carrying capacity under the stability conditions of the normal preset operating modes. When calculating modes, steady-state calculation methods are used. To do this, one can use the node potential method or other iterative calculation methods.
The uniqueness of the solution of the steady-state electric network with successive weighting of the regime, as it is known, characterizes the steady-state regime. The deterioration of convergence characterizes the presence of "heavy" regimes. The violation of convergence determines the approach to the limiting regimes for both aperiodic stability and other constraints. The limit of regime existence is close to the limit in aperiodic stability.
Consecutive weighting of the regime, as a rule, must begin with the "weak" ties in this cross section. Next, it is necessary to determine the limit regime for other ties in this cross section. The following situations are possible: First, a solution can exist, but the regime is aperiodically unstable, or secondly, the state of the power system cannot be limited in stability.
The EPS sensitivity issue in steady-state regimes was first studied in relation with the analysis of errors' influence in initial data. From the EPS modeling point of view, errors' influence in initial data on the results of an accurate model is quite relevant to the effects of the real disturbances in the values of the regime parameters. Thus, the errors' study in results is interesting when considering the problem of weak points. For instance, in [34], attention was drawn to role of the Jacobian matrix condition number in estimating the calculation results error, showing that the main reason for poor conditionality was the branches' resistance difference in the electric network. The disadvantage of most approaches is vast numerous calculations. At the same time, the main difficulty is related not so much to the modeling of situations and the execution of calculations, but rather to the analysis and systematization of the obtained results. In addition, it should be noted that the calculation of the roots is very simple for the characteristic equations of the first and second degrees. There are general expressions for the roots of equations of the third and fourth degrees, but these expressions are cumbersome and useless. As for the equations of higher orders, it is generally not possible to write a general expression for them. Therefore, a direct solution of the characteristic equation has not been widely used thus far in stability calculations.
The calculation based on the bus admittance matrix and Jacobian matrix provides a more accurate calculation of small eigenvalues in comparison with competing algorithms.
To determine sensitivity based on the Jacobian matrix and the bus admittance matrix, a smaller number of operations are required, while sensitivity analysis under normal modes requires less resources, providing the accuracy of the results of computations. The advantage of using the Jacobian inverse matrix is that it is not necessary to solve a system of linear equations in this case. However, the use of the Jacobian matrix has advantages, which will be considered further.
Sensitive elements are those elements of the electrical network in which a given (single) disturbance of the primary parameter causes the greatest change (reaction) of the secondary parameter. At the same time, the definition of sensitivity (sensority) is the degree of regime parameters' reactions to a single disturbance, which can be determined taking into account numerical experiments or indirect indicators, specifically by the following proposed indicators connected with the singular values and eigenvalues of sensitivity matrices. As a result, operational parameters and elements of the EPS that are related to them, with sensitivity that is significantly higher than the sensitivity of others, are called sensory [36].
The most interesting question in the research of EPS vulnerability is the localization of the diagram's branches, changes in the parameters of which influence sensitivity.
In mathematics, in particular functional analysis, the singular value decomposition of the matrix A allows one to calculate the singular numbers of a given matrix as well as  Therefore, to assess the impact of the diagram parameters and operational parameters on the sensitivity of its elements, the matrix derived from the Jacobian matrix by conductance of branches yij of the analyzed diagram was studied.
The presence and placement of sensors and weak points was determined using numerical and analytical methods for studying the scheme of the EPS and its parameters of singular analysis with the Jacobian matrix formulated in Equation (2) and the bus admittance matrix given in Equation (3):  [36,37].
Based on singular analysis, the following indicators   J and   Y were used as generalized indicators to identify weak branches [36,37]: -were derivatives of Jacobian matrix by the conductance of branches, The algorithm for the calculation of steady-state parameters of the EPS is shown in the flowchart in Figure 5. Appendix A shows detailed description of the computational procedure. It should be noted that the considered approach to the analysis of the power system was very convenient for using. It made it easy to identify the weak points of the EPS and to predict the balance/imbalance of the power system in conditions of shortage and excesses of electricity as well as to take measures to achieve a normal balanced regime. To compare sensitivity estimation obtained using the singular analysis of the Jacobian matrix and the singular analysis of the nodal conductance matrix, we concluded that almost the same nodes and branches were sensitive to the deviation of the voltage modules. The distinction in results may have been explained by the fact that the sensitivity of some elements of EPS were more dependent on system parameters and others on regime parameters [38].
The greatest interest from the point of view of the research of the vulnerability of the EPS was the localization of the branches of the network, the change in the parameters of which maximized the sensitivity. The indicator of increasing sensitivity was the decrease of the minimum singular value σ1 of the Jacobian matrix. Therefore, to evaluate the influence of the parameters of the network scheme and the parameters of the regime on the sensitivity of its elements, the derivative matrix of the Jacobian matrix, with respect to the conductivities of the branches of the analyzed network, was studied.

The Obtained Results and Analysis
It should be noted that many scientific articles and monographs are devoted to the problem of calculating the stationary modes of an EPS [15,[39][40][41][42][43][44][45]. IEEE test schemes [46][47][48][49][50][51] are often used for testing, comparing, and demonstrating methods for solving problems of analysis and modeling the static modes of various power systems. In this paper, the authors use the EPS model, which approximates real conditions, partially reflecting the situation in the EPS of the Baltic States [11,[52][53][54][55][56][57].
Research of the Latvian power system was carried out based on three model schemes obtained for the past mode (2017) presented in Figure 6, the current mode (2020), and the planned mode (2025) including development scenarios (Figures 7 and 8) considering the share of the growth prospects of renewable power generating capacity and the disconnection of continental Europe from the BRELL. The Latvian power system has interconnection lines with Estonian, Lithuanian, and Russian networks and is a part of the Nord Pool Spot market. The transmission power system of Latvia operates in two voltage classes: 110 kV and 330 kV.   The normal balanced mode of operation of the studied EPS model was considered the initial one. The mathematical model included the parameters of the stationary mode of the EPS (voltage class, loads, and generation) and network parameters (the combination of internal and external communication for power transmission, branch resistances).
The Nord Pool statistical data from 2017 have been taken for initial computations that are given in Table 4. At the same time, a 2% increase in consumption and generation was taken into account, considering forecasts until 2025. In addition, the expected redistribution between energy types over the examined period, until 2025, was considered [25].
Difference estimates corresponding to the nodes at the ends of the branches of the components of the first right singular vector σmin of the matrix made it possible to identify the branches sensory to the loss of voltage: 2, 4, 6, 10, 12, 18, 23 (Figure 9b).
Evaluation of the computation results of the sensitivities obtained using the singular analysis of the Jacobian matrix ( Figure 9) and the singular analysis of the bus admittance matrix ( Figure 10) showed that practically the same nodes and branches had sensitivity in terms of deviation of voltage modules. Analysis of the obtained results allowed us to establish that branches 2' and 18', whose components were maximal, were weak.
The difference in results was explained by the fact that the sensitivity of some elements of the EPS depended only on the topological parameters of the network, while others on the regime parameters as well.
Using the methodologies "Singular analysis of the Jacobian matrix" and "Singular analysis of the bus admittance matrix", the calculations for the schemes of 2020 ( Figure 7) and 2025 (Figure 8) were provided. Obtained results are presented in Table 5.
As followed by the computation results (Figures 9-14) as well as the results of Table  5, the received results of weak point localization for different development and modernization stages of the Latvian EPS model as well as vulnerability assessment in the design phase allowed us to provide activities associated with the construction of full-strength networks, thereby pre-emptively eliminating dangerous or unwanted situations, and in the operational phase, to develop protective measures minimizing damage to the power system and the country's economy through possible emergencies.           Summarizing the obtained results of the examined network in terms of the sensors and weak points, the following was observed: • For past mode-2017 ( Figure 6), it could be stated that nodes 3' (CHP-1 (Riga)), 4' (CHP-2 (Riga)), 5' (Bishiciems (Riga)), and 6' (HPP (Riga), and branches 2 and 18 were sensory in terms of voltage ( Figure 15); • For current mode-2020 (Figure 7), sensory nodes were 1' and 10' and 20', and branch 2 was sensory ( Figure 16).

Discussion
The formation of focused directions for power industry development can be carried out based on an investigation of both the stationary and transient states of the EPS, including a study of the static stability. The complicity, structure, impacting factors, and possible state of the EPS must be considered. The initial task was set up for the EPS' requirements as well as the equipment and control systems.
The most important phase in the EPS vulnerability analysis was the examination of the generation impact scenarios leading to the development of an EPS crash. When a power system has a relatively small number of elements, impact scenarios can be created manually by listing different variants. However, the elements' number increases rapidly due to EPS expansion, leading to network configuration changes.
As a result, the proposed method provided the possibility of investigating and discovering the most vulnerable parts of the power system that caused EPS crashes by searching so-called sensitive points (nodes, lines) without performing statistical tests and counting a huge number of external exposure scenarios.
The computations considered the following parameters: • EPS stationary mode parameters such as voltage magnitude and angle, power, load, and generation at nodes; • Network parameters such as the internal and external communication sets for power transmission and branch resistance.
In order to ensure convergence (balance between generation and consumption) of the EPS operational mode, both the balance node and reactive power were set arbitrarily under the examined model and assumed to be close to the real operating conditions of the EPS' models. The computations of the respective modes were approximate as well, and thus the estimation of active and reactive power flows could not be considered with high reliability.
The task of analyzing EPS heterogeneities included use of the matrices based on generalized values of EPS elements and data about the applied points of perturbations for searching and localizing examined EPS weak points. The generalized values were expressed via parameters of the nodal conductance matrix and the Jacobi matrix of equations in a steady-state mode serving as the main data sources about the EPS' sensitivity and heterogeneity. Sensitive points and their locations were defined using numerical and analytical methods considering the power network's scheme and parameters.
Due to the limited scope of this work, computations were performed for only one mode, aiming to demonstrate the effectiveness of the proposed approach rather than a detailed study of the behavioral properties of the system or namely, parameters of the EPS' stationary mode and network.
This study included evaluation of three modes of the Latvian EPS model, in particular the past (2017), current (2020), and planned (2025), by analyzing and comparing obtained computation results based on the proposed method. The performance of the Latvian EPS revealed positive results under realized and planned modernization scenarios, leading to an increase in the stability of the power system by eliminating some of the existing weak points. The examined EPS models were close to real-world conditions and partially reflected the situation of the Baltic States' EPSs (taking into account development strategy) without claiming detailed similarity to the latter.
The obtained results of EPS vulnerability localization allowed us to provide measures aimed at constructing an equally stable power network and thus preventing dangerous or undesirable situations as well as developing protective measures to minimize EPS damage during its operational phase from possible external exposures.
Currently, the authors are working on data library creation devoted to collecting information about operating modes that are the most typical for a particular time of year, considering climatic conditions (including long-term forecast of their changes), time of day, market requirements, etc. in order to increase the reliability of the performed computations.

Conclusions
Latvian EPS configuration has been experiencing great changes because of modernization processes (new interconnections of transmission lines, construction of large new wind parks), leading to a need for evaluation of the Latvian power system's static stability for both current and planned modes. For this purpose, the detection methods of the weakest network parameters in relation to external impacts must be applied, highlighting the identification of such nodes as a priority task. In this way, the nodes' sensitivity will increase, contributing to the occurrence of their significant violations in the deviation and content of higher harmonic voltage components.
This article demonstrated an application of analytical techniques for node conductance matrix Y and Jacobian matrix J of power system equations of the steady-state mode. Such an application resulted in the development of an original instrument for identification of the sensitivities and vulnerabilities of a heterogenous power system.
Scientific novelty lies in finding the "hidden values" in aforementioned matrixes and enhancing the potential of their usage for the purposes of power system research. The possibility and efficiency of these methods' applications, leading to the ability to solve analysis problems and control modes by considering voltage limits at buses, are what justifies a wide range of applications: The proposed approach allows detecting the most vulnerable parameters of the power system such as voltage drop and external impacts (sensors) as well as their relationships with EPS parameters. The availability of data of sensor locations provides the possibility of determining and controlling nodes, in which the biggest oscillations of operational parameters are observed due to disturbances in the system. As result, EPS behavioral properties are improved.
The model of the Latvian EPS for the periods of 2017, 2020, and 2025 has been analyzed in terms of an effect on the nodes' sensitivity degree (sensory "leaf" nodes) based on a weighting scheme of the reactive load power. By identifying weaknesses, predicting balance/imbalance of the power system under deficits and surpluses of electricity will force taking further measures to achieve a balanced system regime.
Simulations of the developed Latvian EPS model have been performed considering new interconnection lines and an increased share of WPS energy; thus, to keep the system balanced, a part of the reacting power has been increased.
Comparison of the past mode (2017) and the planned mode (2025) showed positive changes in the revealed weak places/nodes of the Latvian power network. To conclude, the computation results of the planned and ongoing modernization of the Latvian EPS predict an obvious decrease of several weak spots.
It is worth mentioning that the existence of weak spots is inherent in any power supply system, but their sensitivity can be reduced. As the obtained computation results showed, weak point localization as well as vulnerability evaluation of the Latvian EPS in the design stage allows for defining further measures to increase system static stability.

Data Availability Statement:
The datasets used and/or analyzed in the current study are available from the corresponding author upon reasonable request.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A
This appendix presents the computation example that clarifies how the proposed mathematical model was applied to define EPS sensitivity elements, evaluation of the EPS model based on the current mode-2020. The method (with nodal conductance matrices and Jacobi matrices) and the process of identifying these sensitive points are presented below.
The following steps were performed (for details, kindly see Figures A1-A6 below): 1. Set up the initial data of the considered power network in a steady state; 2. Calculate the main parameters' values of the examined EPS' sensitive elements; 3. Build the matrix of own and mutual node conductance; 4. Localize the EPS' sensitive points using singular decomposition of the Jacobian matrix and obtain results; 5. Localize the EPS' sensitive points using singular decomposition of the node admittance matrix and obtain results; 6. Analyze the obtained results and evaluate the EPS' stability.