Modeling of Residential Photovoltaic (PV) System Connected to Low Voltage (LV) Network: Application to Public Distribution Network of Burkina Faso

Abstract

Residential photovoltaic (PV) systems connected to the grid are used for self-consumption. Any surplus production is fed into the grid and contributes to improving the voltage. Several techniques are developed to model their connection. However, studies on methods of injecting energy production into the Low Voltage (LV) network are nowadays a problem. This paper proposes a mathematical model to determine the current to be injected and calculate each node’s voltage. The current equation is a recurrence relation with an initial condition. This initial condition is for the case of a single PV system connected to the LV grid. The equation can also be written in matrix form. Similarly, the voltage solution is a recurrence relation. It also has an initial condition for the first node. Both mathematical formulae with the proposed initial conditions are consistent and can be used for the determination of the current and voltage of the different nodes in the grid.

Share and Cite:

Zongo, A. , Konane, D. and Ouattara, F. (2022) Modeling of Residential Photovoltaic (PV) System Connected to Low Voltage (LV) Network: Application to Public Distribution Network of Burkina Faso. Journal of Applied Mathematics and Physics, 10, 2698-2707. doi: 10.4236/jamp.2022.109180.

1. Introduction

Solar energy captured through photovoltaic (PV) panels represents a viable energy alternative for electricity generation. The latter is a renewable source, both clean, unlimited and with a very low level of accident risk according to B. Boukezata et al., 2014 [1]. In Burkina Faso, current studies seek to propose the optimal option of injecting this renewable energy into the LV network to improve the voltage at the user’s premises. To this end, several methods have been developed.

Guingané, 2018 in his thesis [2] proposes the case of grid-connected PV systems without storage as a solution to enhance local power generation. According to the latter, a generalized integration study of these PV systems on the power grid allows determining with precision the influence of PV production on the power grids of Burkina-Faso. Zongo et al., 2020 in their article [3] carried out theoretical studies that helped to identify the critical lengths of low voltage lines. This allowed them to determine for any LV network the maximum length that should not be exceeded at the risk of not providing a regulated voltage to consumers. Similarly, Zongo et al., 2022 in their paper [4], proved that the peak hours of national electricity consumption are observed between 19:00 and 23:00 Local Time (LT). By analogy M. Du et al., 2017 [5] proposed a voltage drop calculation software by a distribution network node to alert the network operator. Together, these mechanisms prove that power generation from a grid-connected PV system can contribute to substantially improving the voltage quality. Considering that during the high-demand period there is no sunshine in Burkina Faso, it is necessary to use PV systems with storage to optimize the critical lengths of LV lines. In his Ph.D. thesis G.O. Cimuca, 2016 [6] asks whether these PV systems with energy storage, the object of modeling, can increase the efficiency of the existing system of transmission and distribution of electrical energy. Also, with the advent of new technologies, several methods of interconnection between solar PV systems and the power grid have been successfully developed. This is the case for both old and new configurations of grid-connected PV systems. Nowadays, research is ongoing to define the best technique to inject the output of renewable sources into the grid to efficiently improve the voltage of the users. It can be with or without storage. Thus, we will carry out work to study the influence of renewable sources connected to the grid on voltage. The objective of this study is to model a public distribution network with residential PV systems and to analyze the voltage variability.

In the following, we first present the single-phase model of grid-connected residential PV systems, then establish and solve the mathematical equations. Then, we analyze and discuss the results obtained, and finally, we conclude with some perspectives.

2. Single-Phase Model of Grid-Connected Residential PV Systems

The work on critical lengths has made it possible to realize the voltage drops that we have on the public distribution network in Burkina Faso according to Zongo et al., 2022 [4]. However, users of electric current who are connected to isolated networks, as well as those connected to large networks, use the same types of electrical appliances because the requirements of voltage quality are generally the same, G.A. Koucoï, 2017 [7]. To this end, let us ask the question of whether voltage dips can be mitigated by the use of residential grid-connected PV systems. The problem is that there is no sunshine during the periods of high power demand (7 pm to 11 pm) while voltage drops are higher according to Zongo et al., 2022 [4].

The objective of this section is to present the equivalent scheme of grid-connected residential PV systems, Figure 1.

Therefore, we have chosen nodes and connected residential PV systems with storage and performed tests to verify their impacts on the voltage plan As a reminder, these types of tests were carried out among users of the Societé National d’Electricité du Burkina (SONABEL) at substation N˚55 in the Sarfalao district of the city of Bobo-Dioulasso, Burkina Faso according to Zongo et al., 2022 [4]. The data collected was the subject of a published article [4].

To facilitate the study, we have reducepublishedam to a single-phase, Figure 1. To change to a three-phase, we simply multiply the measured voltages and currents by the root of three.

The meaning of the different annotations can be consulted after the bibliographic references.

We notice in Figure 1 that from the Medium Voltage (MV) line, we pass through a transformer that lowers the High Voltage category A to Low Voltage (HTA/BT) to adapt it to the consumers’ needs. According to Zongo et al. [4], most of the voltage drops resulting from the study are higher than the maximum regulatory value of 8% of the nominal values. Thus we make this mathematical study to determine the impact of residential PV systems on critical voltages and lengths.

Figure 1. Model equivalent diagram of grid-connected residential PV systems.

3. Determination and Solution of the Mathematical Equations

The voltage drop depends on the energy flow exchanged between two nodes as well as the impedance of the power line. To reduce the voltage drop, either the power line is replaced by a larger cable section or the exchanged energy flow is reduced according to X.L. DANG, 2014 [8]. In this study, we chose to integrate residential grid-connected PV systems to reduce the energy flow. From Figure 1, we determine and solve the mathematical equations with all the quantities in complex numbers.

3.1. Determination of the Induction Currents of PV Systems

First of all, we calculate the induction differences across each resistor and the current intensity by applying Kirchhoff’s two laws:

From the law of nodes we have:

I S + I G 1 = I C 1 + I 1 (1.1)

From the law of meshes,

If we take mesh 1 we have:

V 1 I 1 Y L 1 V 2 = 0 (1.2)

This gives

I 1 = ( V 1 V 2 ) Y L 1 (1.3)

We also know that

I C 1 = V 1 Y C 1 (1.4)

From relations (1.1), (1.2), (1.3), and (1.4) we find that:

I S + I G 1 = V 1 ( Y L 1 + Y C 1 ) V 2 Y L 1 (1.5)

Then, again applying Kirchhoff’s law for the second node, we have:

I 1 + I G 2 = I C 2 + I 2 (1.6)

V 1 V 2 = I 2 Y L 2 (1.7)

I 2 = ( V 2 V 3 ) Y L 2 (1.8)

We also know that

I C 2 = V 2 Y C 2 (1.9)

and

I G 2 = I C 2 + I 2 I 1 (1.10)

Hence

I G 2 = V 2 Y C 2 + ( V 2 V 3 ) Y L 2 ( V 1 V 2 ) Y L 1 (1.11)

I G 2 = V 1 Y L 1 + V 2 ( Y C 2 + Y L 2 + Y L 1 ) V 3 Y L 2 (1.12)

I C 3 + I 3 = I 2 + I G 3 (1.13)

I G 3 = I 3 + I C 3 I 2 (1.14)

I G 3 = V 3 Y C 3 Y L 2 ( V 2 V 3 ) + Y L 3 ( V 3 V 4 ) (1.15)

I G 3 = V 2 Y L 2 + V 3 ( Y C 3 + Y L 3 + Y L 2 ) V 4 Y L 3 (1.16)

From Equations (1.12) and (1.16), we can generalize by recurrence the equation for injection currents as a function of loads as follows:

For any natural number n > 1, with n the number of nodes in the network, we have Equation (1.17):

I G n = V n 1 Y L n 1 + V n ( Y C n + Y L n + Y L n 1 ) V n + 1 Y L n (1.17)

With all quantities in complex numbers.

When it is a single node the equation is as follows:

I S + I G 1 = V 1 ( Y L 1 + Y C 1 ) V 2 Y L 1 (1.18)

From these equations (1.17) and (1.18) we get the following matrix system;

[ I S + I G 1 I G 2 I G 3 I G n ] = [ Y L 1 + Y C 1 Y L 1 0 0 0 0 0 Y L 1 Y C 2 + Y L 2 + Y L 1 Y L 2 0 0 0 0 0 Y L 2 Y C 3 + Y L 3 + Y L 2 Y L 3 0 0 0 0 0 0 0 0 Y L n 1 Y C n + Y L n + Y L n 1 Y L n ] [ V 1 V 2 V 3 V n ]

To solve this system we can also use numerical methods to determine the voltages V n . Other ways of solving are possible. In this paper, we determine the voltage expression by applying the same Kirchhoff’s laws.

3.2. Determination of Node Voltages

To improve the voltage of a given node, it is necessary to inject generation to reduce the energy flow. The solution to the above matrix system allows quantifying the current to be injected at each grid node with the PV system. This injection current influences the node voltage. Thus, to determine the variability of the voltage as a function of the generator currents, we solve the following equations:

From Equation (1.2) we find Equation (2.1)

I C 1 Y C 1 I 1 Y L 1 I C 2 Y C 2 = 0 (2.1)

I 1 Y L 1 = I C 1 Y C 1 I C 2 Y L 2 (2.2)

I 1 = I C 1 Y L 1 Y C 1 I C 2 Y L 1 Y C 2 (2.3)

Replacing I 1 in Equation (1.1) with its expression (2.3), we obtain equation (2.4).

I s + I G 1 = I C 1 + I C 1 Y L 1 Y C 1 I C 2 Y L 1 Y C 2 (2.4)

I s + I G 1 = I C 1 + V 1 Y L 1 Y L 1 Y C 2 I C 2 (2.5)

V 1 Y L 1 = I s + I G 1 + Y L 1 Y C 2 I C 2 I C 1 (2.6)

V 1 = I s Y L 1 + I G 1 Y L 1 + I C 2 Y C 2 I C 1 Y L 1 (2.7)

The relation (2.7) then leads to the equations of V 1 (2.8) and (2.9).

V 1 = 1 Y L 1 ( I s + I G 1 I C 1 ) + I C 2 Y C 2 (2.8)

V 1 = Y L s Y L 1 V L s + I G 1 Y L 1 I C 1 Y L 1 + I C 2 Y C 2 (2.9)

The voltage V 1 can also be expressed as a function of the source voltage V S according to the mesh law:

V 1 = V s I s Y L S (2.10)

Equation (2.10) is then the initial condition.

Applying again the law of nodes and meshes for the second node we find the Equation (2.11).

I 1 + I C 2 = I 2 + I C 2 (2.11)

V L 2 + V 3 V 2 = 0 (2.12)

I 2 Y L 2 + I C 3 Y C 3 = I C 2 Y C 2 (2.13)

I 2 = Y L 2 ( I C 2 Y C 2 I C 3 Y C 3 ) (2.14)

By replacing I 2 in Equation (2.11) with its value in Equation (2.14), we determine equation (2.15).

I 1 + I G 2 = I C 2 + Y L 2 ( I C 2 Y C 2 I C 3 Y C 3 ) (2.15)

I 1 + I G 2 = I C 2 + Y L 2 ( V 2 V 3 ) (2.16)

V 1 V L 1 V 2 = 0 (2.17)

V 1 = V 2 + V L 1 = I C 2 Y C 2 + I 1 Y L 1 (2.18)

I 1 Y L 1 = V 1 I C 2 Y C 2 (2.19)

I 1 = V 1 Y L 1 Y L 1 Y C 2 I C 2 (2.20)

Similarly, by replacing I 1 in Equation (2.15) with its value (2.20), we get the combination (2.21).

I G 2 + V 1 Y L 1 Y L 1 Y C 2 I C 2 = I C 2 + V 2 Y L 2 V 3 Y L 2 (2.21)

V 2 Y L 2 = V 1 Y L 1 I C 2 Y C 2 Y L 1 + I G 2 I C 2 + I C 3 Y C 3 Y L 2 (2.22)

V 2 = I C 2 Y L 2 + I G 2 Y L 2 I C 2 Y C 2 Y L 1 Y L 2 + V 1 Y L 1 Y L 2 + I C 3 Y C 3 (2.23)

Equation (2.24) is for the voltage V 2 .

V 2 = V 1 Y L 1 Y L 2 I C 2 Y L 2 ( Y L 1 Y C 2 + 1 ) + I G 2 Y L 2 + I C 3 Y C 3 (2.24)

Kirchhoff’s law has been reapplied to the third node. Equation (2.25) is found:

I 2 + I G 3 = I 3 + I C 3 (2.25)

V 3 + V L 3 + V 4 = 0 (2.26)

I C 3 Y C 3 = I 3 Y L 3 + I C 4 Y C 4 (2.27)

I 3 = I C 3 Y C 3 Y L 3 I C 4 Y C 4 Y L 3 (2.28)

I 2 + I G 3 = I C 3 + I C 3 Y C 3 Y L 3 I C 4 Y C 4 Y L 3 (2.29)

V 2 = I 2 Y L 2 + I C 3 Y C 3 (2.30)

I 2 = Y L 2 V 2 I C 3 Y C 3 Y L 2 (2.31)

Combining Equations (2.29) and (2.31), we end up with (2.32).

Y L 2 V 2 I C 3 Y C 3 Y L 2 + I G 3 = I C 3 + I C 3 Y C 3 Y L 3 I C 4 Y C 4 Y L 3 (2.32)

V 3 Y L 3 = Y L 2 V 2 I C 3 Y C 3 Y L 2 + I G 3 I C 3 + I C 4 Y C 4 Y L 3 (2.33)

V 3 = Y L 2 Y L 3 V 2 Y L 2 Y L 3 I C 3 Y C 3 + I C 4 Y C 4 + I G 3 Y L 3 I C 3 Y L 3 (2.34)

V 3 = Y L 2 Y L 3 V 2 I C 3 Y L 3 ( Y L 2 Y C 3 + 1 ) + I G 3 Y L 3 + I C 4 Y C 4 (2.35)

V 3 is then the voltage at node 3.

By recurrence, we can generalize the expression for the voltage V n . For n given nodes, we express the corresponding voltages through the general formula (2.36).

V n = Y L n 1 Y L n V n 1 I C n ( Y L n 1 + Y C n Y L n Y C n ) + I G n Y L n + I C n + 1 Y C n + 1 (2.36)

with n > 1.

4. Analysis of Results and Discussion

The relations (1.17) and (2.36) and the matrix system are the three results to which we carry out the analyses and discussions.

On the one hand, the relations (1.17), (1.18), and (2.10), (2.36) can be written as systems of linear equations when we consider the positive integer n as the number of nodes of the residential PV systems connected to the power grid. These are:

I G n = { I S + V n ( Y L n + Y C n ) V n + 1 Y L n , s i n = 1 V n 1 Y L n 1 + V n ( Y C n + Y L n + Y L n 1 ) V n + 1 Y L n , s i n > 1 (3.1)

V n = { V S I s Y L S , s i n = 1 Y L n 1 Y L n V n 1 I C n ( Y L n 1 + Y C n Y L n Y C n ) + I G n Y L n + I C n + 1 Y C n + 1 , s i n > 1 (3.2)

System (3.1) is the equation of the injection currents of the photovoltaic generators. Similarly, system (3.2) represents the equations of the voltages of the different electrical loads. Their recurrence character allows us to deduce that PV systems can be successively connected to several nodes of the LV network; this could be a solution to increase the critical lengths of the lines according to the article of Zongo et al., 2020 [4]. Guinguane, 2018 [2]; Koucoï, 2017 [7] and Dang, 2014 [8] in their respective theses, have proposed methods for connecting PV systems to the grid. All these methods have allowed us to verify the influence of these systems on electricity distribution. In the same logic as these works, the matrix system and the relations (3.1) and (3.2) show that on a power distribution line, it is possible to improve the voltage at the user’s premises by connecting PV generators.

On the other hand, the matrix writing of Equations (1.17) and (1.18) can be put in the form AX = B with:

A = [ Y L 1 + Y C 1 Y L 1 0 0 0 0 0 Y L 1 Y C 2 + Y L 2 + Y L 1 Y L 2 0 0 0 0 0 Y L 2 Y C 3 + Y L 3 + Y L 2 Y L 3 0 0 0 0 0 0 0 0 Y L n 1 Y C n + Y L n + Y L n 1 Y L n ]

X = [ V 1 V 2 V 3 V n ] ; B = [ I S + I G 1 I G 2 I G 3 I G n ]

We can easily notice that A is a square matrix of order n. Also, equation (2.36) explains that matrix A is invertible since it is a solution of the matrix system. From this, we can deduce that the determinant of matrix A is non-zero, which proves that it is possible to solve the linear equation system (3.1) and find a numerical or algebraic solution.

Furthermore, the relations (3.1) and (3.2) are systems of linear equations. They present the current to be injected and the voltage of the grid load. Similar work by Thi Minh Chau Le in his thesis, [9] was used to calculate electrical parameters for the influence and consistency of PV systems connected to the LV grid. He also determined other parameters that should be taken into account to minimize unwarranted disconnections of grid-connected PV systems to optimize their operation.

5. Conclusion

In this paper, the study models a section of the LV grid with five PV system injection nodes. This modeling has allowed me to write mathematical equations that show the proper functioning of the power grid. They allow knowing the currents to be injected into each node of the network to improve the voltage. In addition, the voltages of the nodes were calculated. The different results are recurrence relations with an initial condition. They were analyzed and discussed, which will allow us to use them for further research. Thus, for a better interpretation of the mathematical equations and the obtained solutions, we plan to carry out in our next research work a simulation under Matlab/Simulink, an experiment on an LV network, and a numerical resolution of the linear equation system (3.1).

The Annotations in Figure 1

I S : the secondary current of the transformer;

P V 1 , P V 2 , P V 3 , P V 4 , P V 5 , P V 6 : the PV generators of nodes 1, 2, 3, 4, 5, 6;

Y C 1 , Y C 2 , Y C 3 , Y C 4 , Y C 5 , Y C 6 : the admittances of loads 1, 2, 3, 4, 5, 6;

Y L 1 , Y L 2 , Y L 3 , Y L 4 , Y L 5 , Y L 6 : the admittances of the line sections of nodes 1, 2, 3, 4, 5, 6;

I G 1 , I G 2 , I G 3 , I G 4 , I G 5 , I G 6 : the PV generator currents of nodes 1, 2, 3, 4, 5, 6;

I C 1 , I C 2 , I C 3 , I C 4 , I C 5 , I C 6 : the currents of the loads flowing through nodes 1, 2, 3, 4, 5, 6;

I 1 , I 2 , I 3 , I 4 , I 5 , I 6 : the outgoing currents from nodes 1, 2, 3, 4, 5, 6;

V 1 , V 2 , V 3 , V 4 , V 5 , V 6 : the load voltages at nodes 1, 2, 3, 4, 5, and 6;

V S : the LV voltage of the source;

G S : the synchronous generator.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1] Boukezata, B., Chaoui, A., Gaubert, J.-P. and Mabrouk, H. (2014) Système solaire photovoltaïque connecté au réseau électrique et associé à un filtre actif parallèle. Université de Sétif 1, Algérie. Article présenté au symposium de Génie Éectrique, 9 p.
[2] Guingane, T.T. (2018) Contribution à l’étude de l’influence d’un systeme PV connecte sur le reseau electrique, thèse de doctorat; Université Ouaga I Professeur Joseph Ki-Zerbo.
[3] Zongo, A. and Ouattara, F. (2022) Experimental Study of the Voltage Quality of the Low Voltage Distribution Network: Application to the City of Bobo-Dioulasso. Open Journal of Applied Sciences, 12, 1262-1270.
https://www.scirp.org/journal/ojapps
https://doi.org/10.4236/ojapps.2022.127086
[4] Zongo, A. and Ouattara, F. (2020) Modélisation des Longueurs Critiques des Lignes Basse Tension (BT) des Réseaux de Distribution Publics. European Journal of Scientific Research, 155, 431-439.
http://www.europeanjournalofscientificresearch.com
[5] Du, M.Q., Li, Y., Liu, C.F. and Liang, T.J. (2017) Low-Voltage Analysis of Distribution Network Software Design and Application. Energy and Power Engineering, 9, 183-188.
https://doi.org/10.4236/epe.2017.94B022
[6] Cimuca, G.O. (2016) Systèmes inertiels de stockage d’énergie associé à des générateurs éoliens. Thèse de Doctorat de ‘École Nationale Supérieure d’Arts et de Métiers, Centre de Lille, France, 171 p.
[7] Koucoï, G.A. (2017) Gestion d’énergie dans les systèmes hybrides PV/diesel pour zones isolées et rurales: optimisation et expérimentation. Thèse de Doctorat Co encadrement international, Institut International d’Ingénierie de l’Eau et de l’Environnement de Ouagadougou (2iE), Burkina Faso, Institut National de l’énergie Solaire (INES) en France, 180 p.
[8] Dang, X.L. (2014) Contribution à l’étude des systèmes PV/Stockage distribués: Impact de leur intégration à un réseau fragile. Thèse de Doctorat de ‘École Normale Supérieure de Cachan, France, 191 p.
[9] Le, T.M.C. (2012) Couplage onduleurs photovoltaïques et réseau, aspects contrôle/commande et rejet de perturbations. Autre. Université de Grenoble, 2012. Français.
https://tel.archives-ouvertes.fr/tel-00721980

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.