Technical Loss Calculation in Distribution Grids Using Equivalent Minimum Order Networks and an Iterative Power Factor Correction Procedure

: Energy losses are a fundamental issue in the electricity distribution sector, being an in-evitable consequence of transporting energy from supplying sources to consumers’ installations and are becoming one of the factors to be considered in planning and operation of electrical distribution networks. So, electrical distribution losses must be continuously monitored so that they are kept within acceptable levels to ensure the business proﬁtability as well as the good power quality of supplied energy. In this context, this work introduces a modiﬁed methodology for technical losses calculation with the application of the concept of reduced equivalent networks, via deﬁnition of an Equivalent Operational Impedance, taking as a starting point the electrical network modeling in the Open Distribution System Simulator (OpenDSS). The losses calculation also considers customer’s energy billing measurement data, measurements of injected energy and power factor at the feeder’s coupling bus at the substation, also considering measurement campaigns to characterize the load consumption proﬁles for working days, Saturdays and Sundays. The proposed methodology disag-gregates energy injections in billed energy, technical and non-technical losses parcels, and presents, as the results have demonstrated, a good precision in the proposed calculation procedures.


Introduction
Energy losses are an inherent aspect of the electric energy distribution service, resulting from several factors, mainly due to physical phenomena that originate from the electric current flow or other factors, such as the poor management of electric utilities or unlawful actions by consumers. In this aspect, electrical losses can be classified according to their origin in technical losses and non-technical losses.
Technical losses are associated with the power dissipation in the different elements of the energy distribution system, in the processes of transformation, transport and energy measurement, representing energy losses that occur in delivering energy from the transmission system and substations to final consumers. These losses cannot be eliminated due to the associated physical phenomena, but their minimization is a constant study subject aiming at the optimization of the distribution system [1].
On the other hand, non-technical losses, determined by the difference between total losses and technical losses, correspond to the energy that is distributed and consumed, but does not return any revenue to the electrical utilities due to several factors, such as calculating losses. The same procedure was compared with the methodology for calculating regulatory losses in Brazil, established by the National Electric Energy Agency (ANEEL), through simulations in a test system and in a real distribution feeder. The comparison proved the superior accuracy of the EOI in relation to the ANEEL methodology and the low computational cost in obtaining technical and non-technical losses [21].
Based on the Equivalent Operational Impedance proposal for the calculation of technical and non-technical losses in distribution systems, this work presents a methodology for applying it to a real distribution feeder, taking as a starting point the network modeling in the Open Distribution System Simulator (OpenDSS), considering the customer's billing measurement data, data of feeder's injected energy and power factor measured at the coupling bus with the substation and data of measurement campaigns to characterize the load consumption for working days, Saturdays and Sundays. With all these data, it is possible to calculate technical losses due to supplying regular loads, non-technical losses and also technical losses due to supplying non-technical losses.
Associated with the EOI methodology application, this paper also presents as a contribution an algorithm that aims to adjust the loads power factors in such a way that the calculated power factor at the feeder coupling bus with the substation through the load flow algorithm, is consistent with the measured power factor value at the same bus. In summary, the proposal presented is an iterative process with the purpose of matching the calculated power factor to the measured one, at the feeder coupling bus with the substation, by adjusting the power factor of all loads, with the purpose of gaining precision in technical loss calculation in the electrical network.

Methodology for Calculation of Technical Losses
The calculation methodology proposed in this work is based on three main points: the electric network modeling, with all its constituent elements, in the OpenDSS software and the use of its load flow calculation methodology; the use of Operational Equivalent Impedance for the calculation of total technical losses; and the implementation of a power factor correction algorithm for low voltage loads, a factor that significantly influences the calculation of technical losses.

Network Modeling
The proposed loss calculation methodology is based on a detailed network modeling, using software OpenDSS, covering all distribution network elements, from the substation power transformers to the terminal branches, in which low voltage consumers are connected. In this way, all the necessary parameters for the network's elements modeling were taken from two main sources: the first concerns the BDGD (The Utility Georeferenced Database), which contains geographical, electrical and commercial parameters regarding the network elements; and the second data source comes from measurement campaign data to characterize loads within the scope of the periodic tariff review process.

Utility Georeferenced Database
In the Brazilian electric system, ANEEL, through Module 10 of PRODIST (Electricity Distribution Procedures in the National Electric System), determined the distribution electric utilities to create and maintain a database to standardize the information availability on network elements by type, technical data and locations, to increase reliability and traceability of this information in the periodic tariff review process [22]. Table 1 presents the information taken from the utility's database for the electric network modeling, according to all its constituent elements, according to the segment to which they belong: substation transformers; medium voltage network with capacitor banks and associated medium voltage loads; distribution transformers; low voltage network; low voltage branch circuits; and low voltage loads. For accessing BDGD, a computational tool was developed in order to read the relevant information from network elements to describe them in the OpenDSS programming language and perform the electrical connection, thus creating a concise circuit compatible with the simulation software.

Measurement Campaigns
The second main data source for the electric network modeling are measurement campaigns to obtain load curves to better reproduce daily consumption of distribution customers, according to their usage electricity habits. In Brazil, it is a common practice of electric energy utilities to carry out extensive measurement campaigns to characterize the loads in their area within the scope of the tariff review process, defined by the national regulator. Such measurement campaigns cover all consumers classes according to their tariff group and voltage level and constitute demand data obtained at 15-min intervals (96 points over a measurement day), comprising characteristic consumption curves for a working day, a Saturday and a Sunday for the same consumer type. The load curves obtained in this process were then also used in methodologies for calculating technical losses.
Using the data contained in the measurement campaigns, load curves were obtained, which characterize the energy consumption of a consumer's facility throughout the day, to be used in the simulations. The load curves are contained in the measurement campaign according to the voltage level and the tariff group to which the loads belong. Table 2 shows the number of types of load curves in the measurement campaign, covering medium and low voltage consumers and residential, commercial, industrial, and rural consumers. Still using the measurement campaigns, multiplicative correction factors were implemented for the average load demand, obtained from energy bills, applied to each load, according to its tariff subgroup and load curve typology. For this, the same data obtained from the measurement campaigns were used to characterize the loads. As a first step, billed energy data were obtained for all loads for which load curves were also specified.
The energy data comprise the energy consumed by the load for a working day, a Saturday, and a Sunday: Subsequently, in order to calculate the average monthly energy for loads, the energies obtained in (1) were used, relating them to the number of working days, Saturdays, Sundays and the total days of the respective month for which the calculation of energy loss is performed. e mean = (e load u·n u ) + (e load s·n s ) + (e load d·n d ) n days (2) where: As a result of the used methodology, Figure 1 presents three load curves that describe consumption characteristic of a residential (Figure 1a), commercial ( Figure 1b) and industrial (Figure 1c) consumer. In the figure, the vertical axis variable is the multiplier factor of the average load demand for each instant of the load curve. It can be seen that the curves applied to each low or medium voltage load, depending on the tariff group and voltage level, make the distinction, both in form and in billed energy, for the same load on different days of the week.

Equivalent Operational Impedance
Equivalent Operational Impedance (EOI) is defined as the complex functional relationship between the total loss of an electrical network divided by the square rms value of the injected electric current at the point of interest for establishing the network equivalent. In radial distribution grids for example, the point of interest can be the feeder's coupling bus at the substation, and in this case, the EOI represents the whole feeder, including the medium and low voltages grids. If the point of interest is a low voltage distribution transformer, then the EOI represents the corresponding low voltage grid. Therefore, Equivalent Operational Impedance is defined as [19,20]: where: In this way, the defined EOI represents a minimal order operational equivalent with the specific purpose of calculating technical losses, the value of which depends on the total calculated losses and the square of the injected rms current. Total losses, in turn, are related to the electrical network topology, the supplied electrical load and the voltage magnitude for an operating point under analysis [18]. Therefore, for each network operational condition, an EOI can be defined, which reflects its operation.
However, taking the definition of the EOI resistive (R EOI ) and reactive (X EOI ) compo- , it is concluded that the load current influences both the numerator, and the denominator, once P T and Q T are given as: P T = R·I 2 inj and Q T = X·I 2 inj ,. Thus, it is reasonable to expect a cancellation effect due to the load current and, consequently, that R EOI and X EOI are not very sensitive to loading variations. Therefore, for variable network loading conditions, but maintaining the same network topology and adequate voltage profile for all loading conditions, for example voltage magnitude variation between ±5%, the EOI is expected to have an approximately constant value for all loading conditions [18][19][20][21].
Also, due to the natural imbalance in distribution networks, in which the loading among phases can vary substantially, a three-phase Equivalent Operational Impedance can also be defined, which is calculated from the relation between the network total loss and sum of the squares of injected currents in each phase [19], according to (5)- (7): where: where: The loss calculation in electric distribution networks has an additional difficulty due to the unbilled energy consumption, which is responsible for non-technical losses. However, non-technical losses generate additional technical losses, and these are difficult to calculate because unbilled energy consumptions are not known with good accuracy. In general, one can write the total active and reactive technical losses of electrical systems as being, where: Losses P T_UC and Q T_UC can be obtained by calculating a load flow to the electric grid considering typical loading curves for costumers' installations as obtained from the customers' energy bills, adjusting the respective power factors, according to the procedure presented in Section 2.3.
Losses P T and Q T , corresponding to total losses, can be obtained for each point of the feeder loading curve, by the corresponding injected current measurement in the substation-feeder coupling bus, using the operational equivalent impedance of the respective electric grid under analysis, that is, P T = R EOI ·I 2 inj and Q T = X EOI ·I 2 inj . Thus, using Equations (8) and (9), non-technical losses are calculated for each point of the feeder loading curve. Therefore, considering the problem presented, of determining the total technical losses, EOI can be used according to the flowchart presented in Figure 2, considering two applications: • First one (Figure 2a): refers to calculation of only one EOL value which is used for all load curve points. This is justified, due to the fact the EOI behavior is little dependent on the load variation along the system load curve, provided that for all operation points voltages are kept adequate, for example, with maximum variations of up to ±5% In these cases, the EOI calculated for a load curve operation point can be used for all other points with good accuracy. This strategy is adopted in references [18][19][20][21]. • Second one (Figure 2b): refers to the new approach of this paper, in which an EOI value is calculated for each operation point following the network loading curve.  (1) Initially, for a given loading condition considering only billed energy consumptions a load flow algorithm is run for the electric grid, to calculate P T_UC and Q T_UC ; (2) For the load flow solution in (1) the EOI is calculated; (3) Measured injected current in second loading condition, with non-technical losses, is obtained. (4) Using the measured injected current for the loading condition under analysis, and the calculated EOI in (2), total active and reactive technical losses are obtained;

Algorithm for Low Voltage Load Power Factor Adjustment
When calculating technical losses through load flow studies in the low voltage grid, an adequate load modeling procedure is essential for obtaining precise results. Thus, among the necessary parameters for load characterization is power factor that has different ranges of values for different types of consumers in the energy distribution network.
The load power factor determination directly implies the reactive power demand to be supplied by the source or by the reactive compensation elements distributed along the electric grid which impact on the electric current flowing in the distribution lines and, consequently, on the magnitude of the observed technical losses. Thus, there is a sensitive relationship between power factor and technical losses occurrence [23,24].
However, for low voltage consumers the electric utilities commercial databases usually store only information regarding the consumption of active energy for billing purposes, as these customers, in general, still do not follow stricter legislation regarding reactive energy consumption. Therefore, the power factor consideration in calculating losses in the low voltage grid implies a certain level of approximation, since it would be necessary to measure power factor at all low voltage load points, considering the different times of the day.
To solve this problem, several studies have already been carried out in order to arrive at a power factor to be used in calculating load flow that best represents the electric network reality, with an approximation level that returns accurate results in the process of determining losses. The most traditional approach consists on considering an approximate power factor that reflects the characteristics of consumers' equipment in operation for the electric grid under analysis [25] or using power factors considered typical for low voltage consumers [26][27][28][29].
Based on the literature, an algorithm was developed for adjusting power factor of low voltage distribution loads in such a way that the load flow calculated power factor at the substation-feeder coupling bus, be consistent with the power factor measured value at the same point. In summary, the proposal presented deals with an iterative process with the purpose of matching the calculated power factor to the measured one, by adjusting the low voltage loads power factors, with the purpose of gaining precision in technical losses calculation.
In this sense, an algorithm was implemented to estimate the power factor of low voltage consumers, based on the power factor value measured at the feeder-substation coupling point that is usually measured by the electric utility. Figure 3 presents a flowchart to explain the algorithm: Power factor adjustment algorithm: Step 1-The electric network model is loaded in OpenDSS; Step 2-A load flow study is executed in OpenDSS in normal mode, that is, the load flow is run for the system average load condition. Initially, it was considered that all low voltage loads have a power factor equal to 0.92, configured as an initial value. For other loads, consisting of medium voltage loads, the measured value is maintained; Step 3-From the load flow results in step 2, the power factor at the feeder-substation coupling bus is calculated, according to Equations (10) and (11): where: S inj -Total calculated apparent power injected into the feeder [kVA]; P inj -Total calculated active power injected into the feeder [kW]; Q inj -Total calculated reactive power injected into the feeder [kVAr]; PF calc -Calculated power factor at the beginning of the feeder.
Step 4-Having the calculated power factor in step 3, it is compared with the measured power factor provided by the substation's measurement system, thus generating a difference as in (12): where: di f -Difference between calculated and measured power factors; FP calc -Calculated power factor; FP mea -Measured power; Step 5-The calculated power factor difference is checked with the solution tolerance, chosen as 0.000001; Step 6-If the calculated difference is less than the stipulated tolerance, the new power factors calculated by the algorithm will be assigned to the low voltage loads for loss calculation.
Step 7-However, if the difference is greater than the specified tolerance, a new power factor is assigned to the loads, given by subtracting the old power factor by the difference calculated in step 4. The algorithm returns to step 2 again, to run the load flow algorithm with the new power factor for all low voltage loads, repeating the process until the difference between the calculated and measured power factors at the feeder-substation coupling bus is less than the stipulated tolerance. This proposed algorithm seeks to add higher precision to technical losses calculation in electric distribution network, given the sensitive relationship between loss calculation and the choice of the load power factor.
The methodology for defining power factors for low voltage loads for loss calculation can also be implemented considering a substation composed of several connected feeders, provided that each feeder has an active and reactive energy measurement system. For this purpose, a computational routine was implemented in the OpenDSS COM interface that allows tracking all low voltage loads connected to each feeder. Based on this information, the algorithm in Figure 3 was applied to each feeder to define the power factor to be used in the loss calculation procedure. To illustrate this aspect, Figure 4 presents a generic substation with n feeders, for which power factors measurements (PF n) at the coupling points with the substation are available. Each feeder supplies a set of low voltage loads (LV Load n) for which new power factor values (New PF LV Load n) will be adjusted, according to the algorithm proposal presented in the flowchart in Figure 3.

Integrated Loss Calculation Methodology
Using the network modeled in OpenDSS, the defined Equivalent Operational Impedance and the power factor adjustment algorithm for low voltage loads, the integrated methodology for calculating technical and non-technical losses was implemented, as shown in Figure 5.

Flowchart description:
(1) Network modeling in OpenDSS is performed using the BDGD database and load curves of measurement campaigns; (2) The algorithm for low voltage loads power factor adjustment is executed; (3) Using the modeled network, a load flow study is carried out following typical daily load curves for a business day, a Saturday and a Sunday, specified for each individual load, with a simulation step of 15 min between load flow case studies; (4) The Equivalent Operational Impedance (EOI) is calculated using (4); (5) Accessing the BDGD, monthly injected energy is calculated for the feeder or substation; (6) The total inject power curve is calculated using the value of monthly energy, obtained in (5), and the system's load curve, obtained in (3). At this point, the injected power curve considers regular and irregular loads; (7) Total injected currents are obtained at every 15 min intervals from data available in (9) and voltages and power factors obtained in (10), according to relationship (13): where:

I-Total injected current [A]; P-Total three-phase injected active power [kW]; V-Line voltage measurement at the feeder-substation coupling point in [kV];
FP-Power factor measurement at the feeder-substation coupling point.
(1) Total three-phase technical losses calculation is then carried out, which at this point includes also technical losses resulting from the supplied billed and non-billed energy.
where: (2) The total billed power of connected loads is read, at 15 min intervals, which corresponds to the respective consumed energy in this time interval. The active power of each load is given by Equation (15).
where: (3) Total system losses are calculated, from the total injected active power curve, obtained in (6), and from the billed power, obtained in (9), according to Equation (16): (4) Finally, non-technical losses are calculated from the difference between total losses, obtained in (10), and technical losses, obtained in (8).

Results and Discussion
In order to demonstrate the proposed methodology accuracy in calculating losses in electrical network, two electrical systems will be used, being the IEEE 123 bus test system and a real urban distribution substation in Belém City, state of Pará, Brazil, belonging to Equatorial Pará, the local electric utility.

IEEE-123 Test Feeder
The IEEE-123 system shown in Figure 6 is a distribution system with nominal voltage of 4.16 kV and being complex enough for the evaluation of different loss scenarios. It has typical elements of real systems, such as overhead and underground lines, voltage regulators, capacitor banks and switching devices. Considering the original system, non-technical losses were added randomly to the load buses represented in red circles in Figure 6, in a proportion varying between 0 and 40% with respect to the given base loads. Following this procedure, total injected active power curves for a time interval of 24 h, were obtained as shown in Figure 7, for the cases of considering only regular (billed) energy consumption, and irregular consumption, that is, with non-technical losses included. For the simulation scenario presented in Figure 6, the voltage regulator transformers' taps in the IEEE 123 bus test system were maintained fixed, for the two loading conditions presented in Figure 7. This way, due to the load increase when non-technical losses are include, it implies in higher voltage drops along the system loading curve, as can be seen in Figure 8, which shows the voltage profiles for the most loaded condition at 7:45 pm, as seen in Figure 7. Considering that OpenDSS methodology considers each bus phase as a node, Figure 8 shows all IEEE system nodes. Considering the daily mode simulation in the OpenDSS, an EOI curve was calculated representing each loading condition, discretized in 15-min intervals, totaling 96 intervals, as shown in Figure 9. At this point, prior to the results presentation, it is worth highlighting the difference between two loss values that will be presented: real technical loss, that is, the calculated loss from load flow execution in OpenDSS when adding randomly irregular loads to the base system load, according to the location shown in Figure 6; EOI technical loss, that is, the calculated loss using the EOI values shown in Figure 9 obtained for the loading condition considering only the billed demands, after the adjustment of load's power factor, and using the total injected current measurement representing the loading condition with non-technical losses included. Therefore, Figure 10 shows the comparison between the real technical loss and the EOI technical loss curves. It is possible to observe a strong adherence between the two curves at each instant of time, which denotes the high precision of the proposed loss calculation methodology. Taking the power loss values presented in Figure 10, daily energy losses can be obtained using Equation (17): where: The percentual error between the energy values calculated using the EOI procedure and the load flow routine is calculated as in (18). Considering the curves values presented in Figure 10, an error of 0.5827% was obtained. error = e Tec_EOI − e Tec_Real e Tec_Real ·100 (18) where: error-Error between the energy values obtained by the two methodologies [%]; e Tec_EOI -Energy calculated using the EOI methodology [kWh]; e Tec_Real -Energy calculated using the real technical loss curve [kWh]; To better evaluate the EOI performance in loss calculation, 1000 load flow simulations were performed, considering random generation for the fraud percentages to be included in the IEEE 123-bus test system base load condition.
As shown in Figure 11a, the average value of non-technical loss added to the base system was 19.42%, with a variation between 16.81% and 22.03%. On the other hand, the average error obtained between the real loss and the loss calculated by the EOI methodology was 0.7437%, with a variation between −0.4797% and 1.9671%. Still considering a normal distribution for the calculated errors, 95.41% of these values are located in the gray region of Figure 11b, which comprises errors with values in the range specified in (19), which denotes an excellent precision of the proposed methodology. where: error-Percentual error; µ-Error average value; σ-Standard deviation; In order to evaluate more scenarios on the percentual values of non-technical loss to be added randomly to the billed load, these percentual values were varied in the ranges of 0-30%, 0-20%, 0-10% and 0-5%. Thus, the average total non-technical losses, considering these scenarios and after one hundred simulations were equal to 14.76%, 9.76%, 4.87% and 2.44%, respectively. Therefore, for these scenarios, the errors Gaussian curves obtained are shown in Figure 12a-d respectively, for average non-technical losses equal to 14.76%, 9.76%, 4.87% and 2.44%, which resulted in percentual errors equal to 0.6366%, 0.5004%, 0.2863% and 0.1486% respectively for applying the EOI methodology.

A Real Urban Distribution Substation Case Study
The real system used in the simulation studies was the Pedreira substation, in the urban area of Belém city in the state of Pará, Brazil, which belongs to the Equatorial Pará, the local distribution utility. This substation has three 20 MVA, 69 kV/13.8 kV transformers supplying 12 distribution feeders, as shown in Figure 13. In Table 3 are presented data to better characterize the substation feeders as: number of low and medium voltage costumers (loads), number of capacitor banks, and lengths of low and medium voltage grids. In addition, at the substation-feeder coupling points, there is measuring devices to measure voltage, injected current, active power and power factor. Using the substation modeled in OpenDSS and measured electrical data at the substation coupling point, technical and non-technical loss calculation will be performed using the integrated loss calculation methodology, shown in Figure 5.
Following flowchart in Figure 6, in step 3 of the loss calculation methodology, the load power factor is adjusted, according to the algorithm presented in Figure 3. For this purpose, power factor measurements at the beginning of feeders, were used and then the algorithm was executed for each feeder individually. As an example of application, Table 4 presents the monthly average power factor obtained from measurements at the beginning of feeders and the power factor calculated for the corresponding monthly billed demand by load flow simulation, before applying the power factor correction algorithm. The objective of the power factor correction algorithm is to make these two values equal, by adjusting the low voltage loads' power factors. As shown in Table 4, a usual practice in most distribution utilities in Brazil is to assume an initial power factor value equal to 0.92 to all low voltage loads. The algorithm of Figure 3 was then executed individually for each feeder, for the average loading condition for one simulation point. Then, after the algorithm execution, different power factor values were obtained for the low voltage loads for each feeder, to be used in step 4 of Figure 5, as presented in Table 5. Therefore, to perform the loss calculation according to the proposed methodology in Figure 5, the adjusted power factor values as presented in Table 5 are used. To illustrate the power factor correction algorithm applied to a daily load cycle, it is presented in Figure 14 the results obtained for PD01 Feeder's loads for weekday, Saturday and Sunday load curve profiles considering only billed demands.
It is evident from observing Figure 14 that adopting a single power factor value to all low voltage loads, as for example 0.92, may lead to imprecisions when calculating technical and non-technical losses in low voltage distribution grids.
After adjusting the low voltage loads' power factors in step 3, the loss calculation methodology is continued in step 4, by executing load flow cases in OpenDSS, for 96 simulation points for weekday, Saturday, and Sunday. In this step, technical losses due to supplying only billed demands (regular loads) are obtained, as observed in Figure 15a. Then, using this calculated technical loss and the total injected current in the substationfeeder coupling bus, both obtained by the executed load flow cases in OpenDSS, the EOI is obtained in step 5, shown in Figure 15b. Still using the load flow results, the distribution grid load curve is obtained in step 6, shown in Figure 15c. Then, using the feeder's monthly injected energy data, obtained from the commercial database, and having the load curve obtained in step 6, the feeder's total injected power curve is obtained in step 9, shown in Figure 15d, which include regular and irregular loads, that is, billed and non-billed power.  Continuing the methodology, step 11 accomplishes the calculation of total injected current into the feeder's electric network, using the measured injected power curve, given in Figure 15d, measured voltage and power factor curves as shown in Figure 16a,b.
Thus, using the total measured injected current and the EOI calculated in Figure 15b, the total technical loss curve can be obtained, as shown in Figure 17. This calculated technical loss comprises losses resulting from supplying regular and irregular loads as well as other sources of non-technical losses that may have occurred in the electrical network under analysis.
After calculating total technical losses, the methodology is continued to calculate the non-technical losses. For this, the data of monthly energy consumed by each consumer unit are extracted and arranged as demand curves for a weekday, a Saturday and a Sunday, in step 13, as shown in Figure 18a. Then, in step 14, the total loss is calculated by subtracting total injected power curves and power demand curves, shown in Figure 18b.   Finally, in step 15, the electric network non-technical losses are calculated, by subtracting the total losses from the total technical losses, as shown in Figure 19. Aiming to obtain the loss data calculated using the methodology in terms of energy, the daily curves can be integrated to obtain daily, weekly, monthly, and annual energy losses data.

Conclusions
In this article, two innovative procedures were highlighted in relation to the systematic calculation of technical and non-technical losses in electrical distribution networks. The first refers to using an Equivalent Operational Impedance (EOI) for representing the electrical network in calculating technical losses. It was demonstrated by the results presented that the use of EOI allows the total losses, which include technical and non-technical losses, to be obtained for each operating condition of the distribution network under analysis, simply by multiplying the squared rms value of the current injected in the feeder, by the operational resistances and reactances, resulting in the respective network active and reactive losses, with minimal computational effort. This feature undoubtedly highlights the use of EOI for real-time operation of electrical distribution networks, which will be approached in future works.
For other applications in technical and non-technical losses management, the joint use of EOI with the proposed procedure of adjusting the power factor of customers' loads, so that the power factor calculated in the substation bus meets the measured value at this point, makes the loss calculation procedure more accurate. Also, the proposed loss calculation procedure allows the breakdown of technical and non-technical losses for each operating condition of the electrical grid, which will certainly contribute to a better management of these losses.
The results obtained in calculating technical losses in the IEEE 123-bus system showed the accuracy of the proposed methodology considering different levels of non-technical losses included randomly in network load buses.
The methodology was also applied in a real distribution system, using real operational data obtained from field measurements. The results presented in terms of technical and non-technical losses, voltages, power factors, and demands are adequately detailed to characterize the electrical distribution grid operational conditions, for a daily cycle of operation, involving typical load curves for weekdays, Saturdays, and Sundays. The results were presented by substation, discretized at 15 min intervals, which constitutes relevant information for the distribution grid short term operation planning.
Considering the implementation of the proposed methodology in a real distribution system, some important data must be collected and maintained continuously by utilities, such as: a georeferenced technical database containing updated information on electrical cables, transformers, capacitor banks, loads and other network elements that are part of the electric network topology, and data informing the electrical connections among them for composing mathematical models in OpenDSS; typical load demands curves describing the loads' consumption cycles for weekdays, Saturdays, and Sunday; a commercial database to inform the monthly customers' energy bills; power factor measurement data at feeders coupling points at the substation. If the electric utility does not have available these data or part of it, this may be a limitation to implement the proposed methodology.
As a continuation of this work, there is a proposal to implement a real time application of the proposed methodology to calculate technical losses for every measured real time injected electric current at the distribution feeders. For such implementation, the EOI calculation will be coupled to the output of a real time network configurator program, so that for each topology modification a new EOI is obtained. This proposal is suited to be implemented in a distribution network management system, obtaining technical losses for each feeder at the coupling point with the substation.
In addition, other operating scenarios will be evaluated in test systems, in order to evaluate the loss calculation accuracy in such cases as: system operating with under and over voltage; system operating with and without voltage regulators; insertion of small and medium-sized distributed generation; and uneven distribution of non-technical losses between different feeders or parts of the system, which often occurs in real systems.