Optimizing Power Flow in Northern Cameroon’s Interconnected Grid: Challenges and Solutions

Preprint

Article

Optimizing Power Flow in Northern Cameroon’s Interconnected Grid: Challenges and Solutions

Altmetrics

Downloads

Views

Comments

Jean Ndoumbe^*,Ivan Basile Kabeina

This version is not peer-reviewed

Submitted:

09 August 2024

Posted:

12 August 2024

You are already at the latest version

Alerts

Abstract

This paper presents an analysis of the power flow within the Northern Interconnected Grid of Cameroon. The Newton-Raphson method has been performed, known for its accuracy, under MATLAB software, to model and solve complex power flow equations. This study simulates a series of outage scenarios to evaluate the responsive-ness of the grid. The results obtained underline the crucial importance of reactive power management and highlight the urgent need to consolidate the grid infrastructure of North Cameroon. To increase grid resilience and stability, the paper recommends the strategic integration of renewables and the development of interconnections with other power grids. These measures are presented as viable solutions to meet cur-rent and future energy distribution challenges, ensuring a reliable and sustainable power supply for Cameroon.

Keywords:

Subject: Engineering - Electrical and Electronic Engineering

1. Introduction

Cameroon, like many developing countries, is facing major challenges in managing its electricity system. These challenges are exacerbated by growing energy demand, aging infrastructure, and financial constraints. The situation is particularly critical in the country’s three main electricity grids: the Eastern Interconnected Grid (EIG), the Northern Interconnected Grid (NIG) and the Southern Interconnected Grid (SIG). Each of these networks presents unique problems, influenced by geographical, economic and technical factors.

For example, the NIG suffers from a production deficit and the saturation of its transmission network. Indeed, population growth in northern Cameroon over the past two decades has led to a significant increase in electricity demand, exacerbating access issues, especially in rural areas [1]. Energy production is also affected by unfavourable hydrological conditions, as evidenced by the declining water level in the reservoirs of the Lagdo hydropower plant, leading to the shutdown of turbines due to lack of maintenance [2].

This situation has resulted in the transmission grid being used in ways not originally anticipated, leading to a series of unplanned occurrences. These include power outages, which have been attributed to the escalating demand on the system and alterations in its operational procedures [3]. The increasing load and operational changes have put unexpected stress on the power infrastructure, highlighting the need for adaptive strategies in managing and upgrading the grid to ensure reliability and efficiency [4].

Power flow analysis serves as a critical tool in addressing issues pertaining to power disruptions and ensuring the power grid’s optimal functionality. Its primary objective is to ascertain the electrical status of the grid, a vital aspect for the effective management and operation of power systems [5]. This analysis is instrumental in identifying the voltages and currents across the network segments, confirming that the equipment functions within secure and stable parameters. By pinpointing areas of power losses and potential overloads, it enhances the utilization of generators and transmission lines. Moreover, it furnishes vital information necessary for the strategic planning of new facilities and the grid’s augmentation, thereby guaranteeing an efficient electricity distribution. Furthermore, power flow analysis contributes to the reduction of operational costs by reducing losses and improving efficiency [6]..

A standard power flow analysis yields key parameters such as the voltage magni-tude and phase angle at each bus bar, along with the active and reactive power transfers between buses. Solving the power flow equations determines the network’s steady-state voltage conditions for each bus bar. However, these equations are inherently nonlinear, making mathematical solutions challenging to derive.

Several researchers have tried the linearization approach of power flow equations. Bolognani and Zampieri [7] proposed a linear approximation of the power flow equa-tions, valid for generic line impedances and network topology. However, real distribution networks can have more complex topologies and features than those modelled, making the direct application of theoretical results more difficult. Garces [8] proposed a linear load flow method for three-phase power distribution systems, which is accurate and applicable to both balanced and unbalanced systems, using a complex-plane linear approximation. Limitations of this method include reduced accuracy with high constant power loads and very low voltages, as well as the failure to consider PV nodes and other common controls in power distribution systems. Liu, et al. [9] proposed a data-driven approach to linearize power flow models, using regression algorithms to improve the accuracy of calculations. Although this approach reduces complexity compared to nonlinear methods, it can still be computationally demanding for large power systems.

The Gauss-Seidel and Newton-Raphson techniques remain predominant for the analysis of power flows. The former is preferred for its simplicity and effectiveness in sce-narios where a well-chosen initial assumption can lead to convergence. For example, Chamim, et al. [10] used the Gauss-Seidel approach to examine the energy flow in Bali’s 150 kV radial power grid. However, when it comes to power flow studies, the Newton-Raphson method is often considered more robust and powerful than the Gauss-Seidel method [11,12] This is mainly due to the rapid convergence rate of the Newton-Raphson method, which is particularly advantageous when dealing with very nonlinear and com-plex equations.

This paper focuses on the analysis of power flow within the NIG, by applying the Newton-Raphson method via the MATLAB software. Various failure scenarios were repli-cated to examine the network’s response to these incidents. The second section describes the northern interconnected network and how it is represented. The third section describes the methods used to model power flow. This same section also discusses the strategy adopted to solve power flow equations, using the Newton-Raphson method. The fourth section is devoted to the presentation of the results obtained, as well as to their analysis and discussion. Finally, the sixth section offers a general conclusion of the study.

2. Grid Overview and Modelling

2.1. Grid Overview

The NIG is instrumental in the distribution of energy generated in the country’s northern region (Figure 1). It facilitates the transfer of power from the Lagdo hydroelectric plant, situated on the Benoue River approximately 66 kilometres from Garoua. The Lagdo facility boasts a considerable storage capacity of 6,300 million cubic meters and spans an area of 700 square kilometres. It consists of four 18 MW units, culminating in a combined capacity of 72 MW and an estimated annual output of 250 GWh. Beyond the Lagdo plant, the NIG also conveys power from the Djamboutou thermal station located near Garoua, which operates on Light Fuel Oil and has a capacity of 20 MW. Additional plants include a 10 MW station in Maroua and a 4.2 MW station in Kousseri, although their output is significantly less. The grid itself is comprised of 400 kilometres of 110 kV transmission lines and 200 kilometres of 90 kV lines, delivering electricity to the urban centres of Ga-roua, Maroua, Ngaoundéré, and Meiganga, thus playing a pivotal role in meeting the en-ergy needs of these cities. Table 1 shows the installed grid connected generation capacities for the NIG.

2.2. Grid Modelling

Modelling large-scale power transmission systems often involves simplification to make calculations more manageable. By decreasing the dimensionality of the grid, the time required for simulations can be significantly reduced while maintaining acceptable accuracy. This simplification is typically achieved by consolidating multiple nodes into one and recalculating the capacities for equivalent transmission lines. This approach al-lows for performance optimization without sacrificing the quality of the results obtained.

The four generators, like the four transformers, share identical characteristics and their parallel configuration allows them to be considered as a single generator or trans-former.

Two lines connect the source node of Lagdo to the node of Garoua. When two power lines run in parallel, their behaviour can effectively be represented by a single line with an equivalent impedance, which is calculated based on the individual impedances of the original lines.

This approach reduces the complexity of the network to make it easier to compute and understand the behaviour of the system. By replacing each component with its theo-retical model, and focusing only on significant loads, such as those at the 110kV and 90kV nodes, a more abstract and manageable representation of the network is obtained. Figure 2 illustrates the simplified topology of the NIG, highlighting the key elements and connec-tions between them, while omitting non-essential details for a clear overview.

2.3. Line Parameters

Calculating the parameters and characteristics of power lines is particularly important for power system analysis.

The calculation of the DC line resistance at a temperature of 20°C is:

R_{C} = ρ \frac{l}{s}

(1)

with:

l: Length of the line (m),

s: Conductor section (m²),

ρ

: Resistivity of metal (aluminium) (Ωm).

The resistivity is a function of the temperature θ and is written by:

ρ (θ) = ρ (θ_{0}) [1 + α (θ - θ_{0})],

(2)

with

α = 4 \cdot 10^{- 3} C^{- 1}

for copper and aluminium,

θ_{0}

being the reference temperature (20°C).

The resistance of a twisted DC conductor exceeds that calculated by the R_C formula. This difference is due to the increased length of the spiral wires, which is longer than the straight conductor. The increase in strength is approximately 1% for conductors composed of three wires and reaches about 2% for those with concentric strands. The AC resistance RA is determined as a function of RC using the Kelvin ratio given by:

q = π d \sqrt{\frac{2 f \cdot 1 0^{- 5}}{ρ}}

(3)

with:

d: Conductor diameter (mm),

f: Frequency of the network (Hz),

ρ: Resistivity of the metal (Ω.mm²/km).

To calculate the reactance of a three-phase power line, it is essential to determine the inductance of the line. In the case of a symmetrical three-phase line, where the conductors are arranged at the vertices of an equilateral triangle, the inductance can be calculated by considering the geometry of the conductors and the uniform distance between them. Considering that the conductors have identical radii, and that the distance between each pair of conductors is equal, then the linear reactance X can be estimated using standardized formulas that consider these parameters.

The inductance L (H/km) of a line is given by:

L = 2.1 0^{- 7} \ln \frac{D}{r e^{- 1 / 4}},

(4)

with

D: the inter-distance between each two conductors,

r: conductor’s radius,

The Inductive reactance X (Ω) is:

X = 2 π f L,

(5)

where

f

is the frequency of the grid.

The capacity C (F/km) of the line can be calculated as:

C = \frac{1}{1.798 \cdot 1 0^{- 7} \ln (\frac{D}{r e^{- 1 / 4}})},

(6)

To model the average length power line, we use the line impedance

\bar{Z}

and the line admittance

\bar{Y}

such as:

\bar{Z} = (R + j X) l,

(7)

\frac{\bar{Y}}{2} = j 2 π f C

(8)

For the short line model,

\bar{Y}

is negligible.

Table 2 presents the parameters of the simplified network of the NIG. The per unit system is used in this modelling to express the quantities of the system as fractions of a defined base unit quantity. This approach simplifies calculations by ensuring that quantities remain constant when processed on different sides of a transformer, which is especially beneficial in scenarios involving many transformers. It also allows for the comparison of several types of equipment on a common scale, as impedances and other parameters are usually within a narrow range when expressed in unit values, regardless of the actual size of the equipment.

3. Power Flow Modelling

3.1. Buses Classification

An electrical system is defined by datasets about the nodes and lines in the network. Each i node is associated with four essential parameters: the voltage magnitude V_i, the voltage phase angle δ_i, the active power P_i and the reactive power Qi injected. For each node, two of these variables are predetermined while the other two must be calculated. When analysing power flows, buses are divided into three distinct categories as shown in Table 3.

The slack bus is responsible for balancing the active and reactive power in the system by accounting for the difference between the total generated power and the total system load, including losses. Essentially, it acts as a buffer, absorbing or supplying power to maintain equilibrium within the electrical network. This ensures stability and reliable operation of the power grid. The slack bus is commonly considered as the reference bus because both voltage V and angle δ are specified, the powers P and Q are to be determined. Usually identified by the number 1, its voltage magnitude and angle remain constant, while its active and reactive powers require calculation.

The Regulated buses are the rest of generator buses. In this type of node, the real powers Pi and the voltage magnitude V are known and controllable. On the other hand, the reactive power Qi and the phase angle of voltage δ remain to be determined.

The load buses are buses with fixed P and Q powers, the voltage magnitude V and the angle δ are unknown. Most of the buses in practical power systems are load buses [13].

PQ buses are characterized by unknown voltage magnitudes and angles, contrasting with PV buses where only the voltage angle remains undetermined. The Slack bus, however, has predetermined voltage magnitudes and angles, eliminating the need to solve for any variables. Within a network comprising n buses and g generators, the total number of unknowns is calculated as 2(n-1)-(g-1) [13]. The resolution of these unknowns is facilitated through the application of real and reactive power balance equations. These critical equations are formulated based on the transmission network’s representation via the admittance matrix, commonly referred to as Y-bus. This matrix is pivotal in analysing the flow of electrical power through the network and is instrumental in the stability and control of power systems. It serves as the foundation for various algorithms designed to ensure efficient operation and management of the electrical grid. The admittance matrix encapsulates the complex interplay between the network’s buses and lines, providing a comprehensive framework for addressing the intricate challenges inherent in power system analysis [14]. Table 4 presents the bus data of the simplified NIG.

3.2. Admittance Equation

The admittance matrix, or Y-bus, is a fundamental component in power system analysis. It is the complex admittance of the components in a power network, allowing for the analysis of power flow and voltage distribution.

The symmetry along the diagonal is due to the mutual admittance between buses being equal in both directions, reflecting the physical reality of power systems where power can flow bidirectionally between nodes.

Y = [\begin{matrix} Y_{11} & \dots & Y_{1 n} \\ ⋮ & ⋱ & ⋮ \\ Y_{n 1} & \dots & Y_{n n} \end{matrix}] .

(9)

The elements of the diagonal are written as:

Y_{i i} = \sum_{\begin{matrix} j = 0 \\ j \neq i \end{matrix}}^{n} y_{i j}

(10)

The elements outside the diagonal are written as:

Y_{i j} = Y_{j i} = - y_{i j}

(11)

3.3. Power Flow Equation

Equation (12) is the polar form of the power flow equations, formulated for an n-bus system as a function of the intake matrix of the Y bus, with reference to Figure 3.

I_{i} = V_{i} y_{i 0} + (V_{i} - V_{1}) y_{i 1} + (V_{i} - V_{2}) y_{i 2} + \dots + (V_{i} - V_{j}) y_{i j}

(12)

Grouping the elements as a function of the voltage, equation (2) becomes as follows:

I_{i} = V_{i} (y_{i 0} + y_{i 1} + y_{i 2} + . . + y_{i j}) - V_{1} y_{i 1} - V_{2} y_{i 2} - \dots - V_{j} y_{i j}

(13)

Even better,

I_{i} = V_{i} \sum_{\begin{matrix} j = 0 \\ j \neq i \end{matrix}} y_{i j} - \sum_{\begin{matrix} j = 1 \\ j \neq i \end{matrix}} y_{i j} V_{j}

(14)

Either

I_{i} = V_{i} y_{i i} + \sum_{\begin{matrix} j = 1 \\ j \neq i \end{matrix}} y_{i j} V_{j}

(15)

The power equation at a given bus is expressed in the following manner:

S_{i} = P_{i} + j Q_{i} = V_{i} I_{i}^{*}

(16)

S_{i}^{*} = P_{i} - j Q_{i} = V_{i}^{*} I_{i}

(17)

Thus, we can write, by inserting equation (15) into equation (17):

S_{i}^{*} = P_{i} - j Q_{i} = V_{i}^{*} (V_{i} Y_{i i} + \sum_{\begin{matrix} j = 1 \\ j \neq i \end{matrix}} Y_{i j} V_{j})

(18)

By introducing the polar forms, we obtain:

S_{i}^{*} = P_{i} - j Q_{i} = |V_{i}| ∠ (- δ_{i}) (|V_{i}| |Y_{i i}| ∠ (δ_{i}) + \sum_{\begin{matrix} j = 1 \\ j \neq i \end{matrix}} |Y_{i j}| |V_{j}| ∠ (θ_{i j} + δ_{j}))

(19)

And finally:

S_{i}^{*} = P_{i} - j Q_{i} = \sum_{\begin{matrix} j = 1 \\ j \neq i \end{matrix}} |Y_{i j}| |V_{i}| |V_{j}| ∠ (θ_{i j} - δ_{i} + δ_{j})

(20)

The active and reactive powers are deduced respectively as follows:

P_{i} = \sum_{j = 1}^{n} | V_{i} | | V_{j} | | Y_{i j} | c o s (θ_{i j} - δ_{i} + δ_{j})

(21)

Q_{i} = - \sum_{j = 1}^{n} | V_{i} | | V_{j} | | Y_{i j} | \sin (θ_{i j} - δ_{i} + δ_{j})

(22)

3.4. Newton-Raphson Technique

The Newton-Raphson (N-R) technique, also known as the method of successive approximation, is based on Taylor’s expansion approximation [13]. Applying Taylor’s expansion to the equations (21) and (22), we obtain:

P_{i}^{s c h} = P_{i}^{k} + \sum_{j = 2}^{n} \frac{\partial P_{i}^{k}}{\partial δ_{i}} \partial δ_{i} + \sum_{j = 2}^{n} \frac{\partial P_{i}^{k}}{\partial |V_{i}|} \partial V_{i}

(23)

Q_{i}^{s c h} = Q_{i}^{k} + \sum_{j = 2}^{n} \frac{\partial Q_{i}^{k}}{\partial δ_{j}} \partial δ_{j} + \sum_{j = 2}^{n} \frac{\partial Q_{i}^{k}}{\partial |V_{j}|} \partial |V_{j}|

(24)

In equations (23) and (24), the inconnue variables are both voltage magnitude and angles

|V_{i}| ∠ δ_{i}

at load buses, and the angles

δ_{i}

at regulated buses.

P_{i}^{s c h}

and

Q_{i}^{s c h}

are respectively the scheduled values of active power at the generator bus and the load and reactive power at load buses. We can write:

Δ P_{i}^{k} = P_{i}^{s c h} - P_{i}^{k} = \sum_{j = 2}^{n} \frac{\partial P_{i}^{k}}{\partial δ_{j}} \partial δ_{j} + \sum_{j = 2}^{n} \frac{\partial P_{i}^{k}}{\partial |V_{j}|} \partial V_{j}

(25)

Δ Q_{i}^{k} = Q_{i}^{s c h} - Q_{i}^{k} = \sum_{j = 2}^{n} \frac{\partial Q_{i}^{k}}{\partial δ_{j}} \partial δ_{j} + \sum_{j = 2}^{n} \frac{\partial Q_{i}^{k}}{\partial |V_{j}|} \partial |V_{j}|

(26)

The use of the matrix expression leads us to:

[\begin{matrix} Δ P_{i}^{k} \\ Δ Q_{i}^{k} \end{matrix}] = [\begin{matrix} \frac{\partial P_{i}^{k}}{\partial δ_{j}} & \frac{\partial P_{i}^{k}}{\partial |V_{j}|} \\ \frac{\partial Q_{i}^{k}}{\partial δ_{j}} & \frac{\partial Q_{i}^{k}}{\partial |V_{j}|} \end{matrix}] [\begin{matrix} Δ {δ_{j}}^{k} \\ Δ {V_{j}}^{k} \end{matrix}] = [\begin{matrix} J_{P δ} & J_{P |V|} \\ J_{Q δ} & J_{Q |V|} \end{matrix}] [\begin{matrix} Δ {δ_{j}}^{k} \\ Δ {V_{j}}^{k} \end{matrix}] i, j = 2, . . ., n

(27)

We deduce that:

[\begin{matrix} Δ {δ_{j}}^{k} \\ Δ {V_{j}}^{k} \end{matrix}] = {[\begin{matrix} J_{P δ} & J_{P |V|} \\ J_{Q δ} & J_{Q |V|} \end{matrix}]}^{- 1} [\begin{matrix} Δ P_{i}^{k} \\ Δ Q_{i}^{k} \end{matrix}] i, j = 2, . . ., n

(28)

So, new values are calculated by:

[\begin{matrix} {δ_{j}}^{k} \\ V_{j}^{k} \end{matrix}] = [\begin{matrix} {δ_{j}}^{k - 1} \\ V_{j}^{k - 1} \end{matrix}] + [\begin{matrix} Δ {δ_{j}}^{k} \\ Δ V_{j}^{k} \end{matrix}]

(29)

The iterative process stops when:

[\begin{matrix} Δ {δ_{j}}^{k} \\ Δ V_{j}^{k} \end{matrix}] \leq p r e c i s i o n

(30)

Figure 4 shows the Newton-Raphson algorithm, implemented in the MATLAB/Simulink software.

4. Results and Discussions

4.1. Scenario 1: All Plants Are Operating

Table 5 presents the results of the NIG power flow. All voltage drops do not exceed 10% of their nominal value. So, in this scenario, all busbars meet the standard. To ensure the balance of the system, the slack bus, which is here the Lagdo hydroelectric power plant, must produce a power of 61.02MW or 84.75% of the installed capacity. However, it also must absorb a reactive power of 22.78 MVar, which represents a loss factor of 34.97%. These results show that the Lagdo plant is essential for the balance of the NIG. However, the loss factor of 34.97% indicates that there is significant improvement in terms of reactive power management, including the integration of reactive power compensation devices. The results of the power flow along the lines of the NIG are presented in Table 6. The Lagdo-Ngaoundéré and Garoua-Guider lines are the most overcrowded, but the power transmitted remains acceptable.

4.2. Scenario 2: The Garoua Power Plant is Out of Service

The results of the NIG power flow when the Garoua power plant is out of service are shown in Table 7. There is a drop in voltage in all the buses in the network. The drop in voltage is more significant at the Garoua bus. This generalized voltage drop could indicate a significant dependence of the grid on the Garoua power plant. In addition, it underscores the strategic importance of the Garoua power plant to maintain the stability of the grid voltage. However, the voltage remains acceptable to all buses despite the decrease, which shows that the network has a certain resilience, even if it should be strengthened.

However, this scenario, in which power demand exceeds the production capacity of the Lagdo plant, has several implications and requires specific measures to optimize production while ensuring safety (Table 8). Indeed, the Lagdo plant would be overloaded, which can lead to premature wear and tear of the equipment, increasing the risk of breakdowns and failures. In addition, overloading can cause voltage fluctuations, affecting the quality of the power supply to consumers. To avoid overloading, it may be necessary to set up load shedding, which will certainly lead to power cuts for some users.

4.3. Scenario 3: The Maroua Power Plant is Out of Service

The results of the power flow (Table 9) of the NIG when the Maroua power plant is out of service also show a voltage drop in all the buses in the network, more marked at the Garoua bus. This indicates a significant dependence of the grid on this thermal power plant and highlights its strategic importance for voltage stability.

Table 10 presents the power flow on NIG lines when Maroua Power Plant is out of service. The results show that there is an overload of the Lagdo hydroelectric power plant, thus increasing the risk of outages and voltage fluctuations, which would affect the quality of the power supply. To avoid overloading, it may also be necessary to set up load shedding.

4.4. Scenario 4: The Kousseri Power Plant is Out of Service

Table 11 illustrates the energy flow results for the NIG when the Kousseri power plant is offline. It reveals a minor voltage drop across all network buses. The failure of this low-power plant leads to an increase in reactive power generated by other plants. This indicates that these plants must compensate not only for the loss of active power but also for the increased demand for reactive power, which can impact their performance and efficiency. These findings highlight the electrical network’s resilience. Despite the shutdown of a power plant, the network manages to offset the power loss, albeit with a slight voltage drop and an increase in reactive power.

This outage highlights the vulnerability of the NIG in the face of the loss of large sources of production as illustrated in Table 12. Despite the resilience of the grid, apparent power losses are increasing, which can lead to voltage imbalances and fluctuations. Although the grid is showing some adaptability, the increase in apparent power losses is a cause for concern. This can cause instabilities and voltage variations that could have repercussions on the distribution of electricity. It is therefore crucial to strengthen the robustness of the network to minimize the risk of failure and ensure a stable and reliable power supply for all users.

4.5. Discussions

The results show that the NIG heavily relies on the Maroua and Garoua thermal power plants to maintain voltage stability. The shutdown of these plants leads to a generalized voltage drop, especially at the Garoua bus, highlighting their strategic importance. The overload of the Lagdo power plant under these conditions is concerning, as it can cause breakdowns, voltage fluctuations, and the need for load shedding. To optimize production and ensure safety, it is crucial to strengthen infrastructure, manage demand efficiently, and integrate additional energy sources [15]. Although the grid shows some resilience, measures are necessary to minimize the risk of failure and ensure a stable and reliable power supply for all users.

Improving the stability of the electricity grid is a major challenge to guarantee a reliable and continuous energy supply. The diversification of energy sources is a fundamental strategy in this regard. By integrating renewable energies such as solar, wind and hydropower, dependence on a single plant is reduced, while contributing to the energy transition to a more sustainable model [16,17]

Smart Grids represent a revolution in energy management. They enable real-time monitoring and optimization of power distribution, using advanced energy management technologies. Although implementation in Cameroon is still subject to several barriers, several researchers have proposed deployment methods [18,19,20].

Finally, interconnections with other power grids provide an additional level of security: NIG is being strengthened and expanded to improve access to electricity and optimize energy exchanges with other networks. First, the interconnection project between the NIG and the SIG aims to create synergy between the two main electricity grids in Cameroon. This will strengthen the stability of the network and optimize the use of resources. Also, an interconnection project between Cameroon and Chad (PIRECT) was launched in November 2023, will allow the construction of 524 km of line between Ntui and Wouro Soua in Ngaoundéré, and 566 km of line from Wouro Soua in Chad via Garoua, Maroua, and Kousseri [21,22].

5. Conclusions

The power flow analysis in the NIG of Cameroon reveals significant challenges related to managing the management of electricity production sources and aging infrastructure. The scenarios studied highlight the crucial importance of the Lagdo, Garoua, Maroua, and Kousseri power plants for grid stability. Diversifying energy sources, integrating renewable energies, and adopting smart grids are essential strategies to improve the resilience and efficiency of the grid. Interconnection projects with other grids, notably with the Southern Interconnected Grid and Chad, offer promising prospects for enhancing stability and optimizing the use of energy resources.

Acknowledgments

Deepest gratitude to Mr. Etienne Gbaweng Belporoh from the Government Technical High School of Ngaoundere for providing the field data that was instrumental in the development of this work.

References

Njoh, A.J.; Etta, S.; Ngyah-Etchutambe, I.B.; Enomah, L.E.; Tabrey, H.T.; Essia, U. Opportunities and challenges to rural renewable energy projects in Africa: Lessons from the Esaghem Village, Cameroon solar electrification project. Renewable Energy 2019, 131, 1013–1021. [Google Scholar] [CrossRef]
Djidimbélé, R.; Kidmo, D.K.; Tchaya, G.B.; Djongyang, N. Optimization of the power flow of photovoltaic generators in electrical networks by MPPT algorithm and parallel active filters. Energy Reports 2021, 7, 491–505. [Google Scholar] [CrossRef]
Iweh, C.D.; Gyamfi, S.; Effah-Donyina, E.; Tanyi, E. Analysis of contingency scenarios towards a suitable transmission pathway in the southern interconnected grid (SIG) of Cameroon. e-Prime-Advances in Electrical Engineering, Electronics and Energy 2024, 7, 100486. [Google Scholar] [CrossRef]
Panteli, M.; Trakas, D.N.; Mancarella, P.; Hatziargyriou, N.D. Power systems resilience assessment: Hardening and smart operational enhancement strategies. Proceedings of the IEEE 2017, 105, 1202–1213. [Google Scholar] [CrossRef]
Mouassa, S.; Alateeq, A.; Alassaf, A.; Bayindir, R.; Alsaleh, I.; Jurado, F. Optimal power flow analysis with renewable energy resource uncertainty using Dwarf Mongoose optimizer: Case of ADRAR isolated electrical network. IEEE Access 2024. [Google Scholar] [CrossRef]
Barnawi, A.B. Development and analysis of AC optimal power flow optimization algorithms for minimization of cost and emissions with stochastic renewables. Energy Reports 2024, 11, 2059–2076. [Google Scholar] [CrossRef]
Bolognani, S.; Zampieri, S. On the existence and linear approximation of the power flow solution in power distribution networks. IEEE Transactions on Power Systems 2015, 31, 163–172. [Google Scholar] [CrossRef]
Garces, A. A linear three-phase load flow for power distribution systems. IEEE Transactions on Power Systems 2015, 31, 827–828. [Google Scholar] [CrossRef]
Liu, Y.; Zhang, N.; Wang, Y.; Yang, J.; Kang, C. Data-driven power flow linearization: A regression approach. IEEE Transactions on Smart Grid 2018, 10, 2569–2580. [Google Scholar] [CrossRef]
Chamim, A.N.N.; Putra, K.T.; Al Farisi, M.F. Power Flow Analysis of Electrical Network Systems Using Gauss-Seidel Method and Python. Emerging Information Science and Technology 2023, 4, 28–36. [Google Scholar] [CrossRef]
Mohsin, M.Y.; Khan, M.A.M.; Yousif, M.; Chaudhary, S.T.; Farid, G.; Tahir, W. Comparison of Newton Raphson and Gauss Seidal Methods for Load Flow Analysis. International Journal of Electrical Engineering & Emerging Technology 2022, 5, 01–07. [Google Scholar]
Abaali, H.; Talbi, E.; Skouri, R. Comparison of newton Raphson and gauss seidel methods for power flow analysis. World Acad. Sci. Eng. Technol. Int. J. Energy Power Eng 2018, 12, 627–633. [Google Scholar]
Albadi, M.; Volkov, K. Power flow analysis. Computational Models in Engineering 2020, 67–88. [Google Scholar]
D’orto, M.; Sjöblom, S.; Chien, L.S.; Axner, L.; Gong, J. Comparing different approaches for solving large scale power-flow problems with the Newton-Raphson method. IEEE Access 2021, 9, 56604–56615. [Google Scholar] [CrossRef]
Iweh, C.D.; Gyamfi, S.; Tanyi, E.; Effah-Donyina, E. Distributed Generation and Renewable Energy Integration into the Grid: Prerequisites, Push Factors, Practical Options, Issues and Merits. Energies 2021, 14. [Google Scholar] [CrossRef]
Kidmo, D.K.; Deli, K.; Bogno, B. Status of renewable energy in Cameroon. Renew. Energy Environ. Sustain. 2021, 6. [Google Scholar] [CrossRef]
Fissou Amigue, F.; Ndjakomo Essiane, S.; Perabi Ngoffe, S.; Abessolo Ondoa, G.; Mengata Mengounou, G.; Nna Nna, P.T. Optimal integration of photovoltaic power into the electricity network using Slime mould algorithms: Application to the interconnected grid in North Cameroon. Energy Reports 2021, 7, 6292–6307. [Google Scholar] [CrossRef]
Mbey, C.; Boum, A.T.; Nneme Nneme, L. Roadmap for the Transformation of the South Cameroon Interconnected Network (RIS) into Smart-Grid. Environmental Science, Engineering 2020, 8, 1–8. [Google Scholar]
Dandoussou, A.; Kenfack, P.; Lekunze Toh, G. Distributed Generation and Optimization of Smart Grid Systems: Case Study of Kumba in Cameroon. World Journal of Electrical and Electronic Engineering 2021, 1, 34–43. [Google Scholar] [CrossRef]
Yem Souhe, F.G.; Mbey, C.F.; Boum, A.T.; Ele, P. Forecasting of electrical energy consumption of households in a smart grid. International Journal of Energy Economics and Policy 2021, 11, 221–233. [Google Scholar] [CrossRef]
Morella, E. Disclosable Version of the ISR-Cameroon-Chad Power Interconnection Project-P168185-Sequence No: 03. 2021.
Kidmo, D.K.; Bogno, B.; Ekam, P.-S.N.; Nisso, N.; Aillerie, M. Hydropower generation potential and prospective scenarios for sustainable electricity supply for the period 2022–2042: A case study of the NIN zone of Cameroon. Energy Reports 2022, 8, 123–136. [Google Scholar] [CrossRef]

Figure 1. NIG single-line diagram.

Figure 2. Simplified topology of NIG.

Figure 3. Typical power system bus bar model.

Figure 4. Newton-Raphson algorithm.

Table 1. Installed grid generation capacity.

Plant name	Fuel Type	Installed capacity (MW)
Lagdo	hydro	72
Garoua	HFO	20
Maroua	HFO	10
Kousseri	HFO	4

Table 2. Line parameters.

From bus	To bus	Length (km)	R (pu)	X (pu)
1	5	238	0.1641	0.7586
1	2	50	0.0742	0.1691
2	3	101	0.1670	0.4994
3	4	99	0.1637	0.4896

Table 3. Type of buses in the power flow problem [13]

Bus type	Voltage (\|Vi\|∠δi)		Real power			Reactive power
Bus type	Magnitude	Angle	Generation	Load	Net (Pi)	Generation	Load	Net (Qi)
Slack / Swing	Specified	Specified	Unknown	Specified	Unknown	Unknown	Specified	Unknown
Generator / Regulated/PV	Specified	Unknown	Specified	Specified	Specified	Unknown	Specified	Unknown
Load / PQ	Unknown	Unknown	Specified	Specified	Specified	Specified	Specified	Specified

Table 4. Bus data.

Bus	Type	V	Phase	P_gen	Q_gen	P_load	Q_load
1	Slack	1	0	-	-	0	0
2	PV	1	-	0.2	-	0.29	0.13
3	PQ	-	-	0	0	0.05	0.02
4	PV	1	-	0.142	-	0.42	0.07
5	PQ	-	-	0	0	0.13	0.03

Table 5. Results of power flow of NIG

No_Bus	V	Phase	P_gen	Q_gen	P_load	Q_load
1	1	0	61.0175	-22.7774	0	0
2	1	-5.4040	20.0000	29.3023	29	13
3	0.9796	-16.6583	0	0	5	2
4	1	-26.2229	14.2000	22.1575	42	7
5	0.9803	-5.8643	0	0	13	3
Total			95.2175	28.6824	89	25

Table 6. Power flow on NIG lines

From bus	To bus j	Pij	Qij	P_loss	Q_loss	S_loss
1	5	13.2899	-3.6989	0.2899	-6.6989	6.7051
1	2	47.7276	-19.0785	1.9326	2.7937	3.3970
2	3	36.7949	-5.5699	2.294	4.6653	5.1987
3	4	29.5009	-12.2352	1.7009	2.9223	3.3812
Total				6.2175	3.6824	7.2262

Table 7. Power flow behaviour of NIG when Garoua Power Plant is out of service.

No_Bus	V	Phase	P_gen	Q_gen	P_load	Q_load
1	1	0	84.0371	6.3937	0	0
2	0.9359	-6.8421	0	0	29	13
3	0.9423	-20.2635	0	0	5	2
4	1	-31.0465	14.2000	30.4009	42	7
5	0.9803	-5.8643	0	0	13	3
Total			98.2371	36.7946	89	25

Table 8. Power flow on NIG lines when Garoua Power Plant is out of service.

From bus	To bus j	Pij	Qij	P_loss	Q_loss	S_loss
1	5	13.2899	-3.6989	0.2899	-6.6989	6.7051
1	2	70.7472	10.0926	3.7918	7.1401	8.0844
2	3	37.9554	-10.0475	2.9035	6.7072	7.3087
3	4	30.0519	-18.7547	2.2519	4.6462	5.1632
Total				9.2371	11.7946	14.9812

Table 9. Power flow behaviour of NIG when Maroua Power Plant is out of service.

No_Bus	V	Phase	P_gen	Q_gen	P_load	Q_load
1	1	0	76.0840	-27.5014	0	0
2	1	-7.0731	20.0000	38.9393	29	13
3	0.9639	-22.5126	0	0	5	2
4	1	-36.1358	4.2000	31.4820	42	7
5	0.9803	-5.8643	0	0	13	3
Total			100.2840	42.9199	89	25

Table 10. Power flow on NIG lines when Maroua Power Plant is out of service.

From bus	To bus j	Pij	Qij	P_loss	Q_loss	S_loss
1	5	13.2899	0.2899	0,2899	-6.6989	6.7051
1	2	62.7941	3.3091	3,3091	5.9373	6.7971
2	3	50.485	4.2684	4,2684	10.6037	11.4305
3	4	41.2166	3.4166	3,4166	8.0778	8.7706
Total				11.284	17.9199	21.1766

Table 11. Power flow behaviour of NIG when Kousseri Power Plant is out of service.

No_Bus	V	Phase	P_gen	Q_gen	P_load	Q_load
1	1	0	67.0801	-24.7363	0	0
2	1	-6.0775	20.0000	32.9034	29	13
3	0.9739	-19.0197	0	0	5	2
4	1	-30.2233	10.0000	25.7458	42	7
5	0.9803	-5.8643	0	0	13	3
Total			97.0801	33.9129	89	25

Table 12. Power flow on NIG lines when Kousseri Power Plant is out of service.

From bus	To bus j	Pij	Qij	P_loss	Q_loss	S_loss
1	5	13.2899	-3.6989	0.2899	-6.6989	6.7051
1	2	53.7902	-21.0374	2.4439	3.9613	4.6545
2	3	42.3463	-5.0953	3.0211	6.852	7.4884
3	4	34.3253	-13.9473	2.3253	4.7986	5.3323
Total				8.0801	8.9129	12.0303

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer

Optimizing Power Flow in Northern Cameroon’s Interconnected Grid: Challenges and Solutions

Abstract

1. Introduction

2. Grid Overview and Modelling

2.1. Grid Overview

2.2. Grid Modelling

2.3. Line Parameters

3. Power Flow Modelling

3.1. Buses Classification

3.2. Admittance Equation

3.3. Power Flow Equation

3.4. Newton-Raphson Technique

4. Results and Discussions

4.1. Scenario 1: All Plants Are Operating

4.2. Scenario 2: The Garoua Power Plant is Out of Service

4.3. Scenario 3: The Maroua Power Plant is Out of Service

4.4. Scenario 4: The Kousseri Power Plant is Out of Service

4.5. Discussions

5. Conclusions

Acknowledgments

References

MDPI Initiatives

Important Links

Subscribe