Next Article in Journal
Multilevel Aircraft-Inverter Design Based on Wavelet PWM for More Electric Aircraft
Next Article in Special Issue
Jerusalem Artichoke: Nitrogen Fertilization Strategy and Energy Balance in the Production Technology of Aerial Biomass
Previous Article in Journal
The Impact of Vehicle Technology, Size Class, and Driving Style on the GHG and Pollutant Emissions of Passenger Cars
Previous Article in Special Issue
The Application of Lignocellulosic Biomass Waste in the Iron and Steel Industry in the Context of Challenges Related to the Energy Crisis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Efficient Modeling of Distributed Energy Resources’ Impact on Electric Grid Technical Losses: A Dynamic Regression Approach

by
Alain Aoun
1,*,
Mehdi Adda
1,
Adrian Ilinca
2,*,
Mazen Ghandour
3,
Hussein Ibrahim
2,4 and
Saba Salloum
3
1
Département de Mathématiques, Informatique et Génie, Université du Québec à Rimouski (UQAR), Rimouski, QC G5L 3A1, Canada
2
Ecole de Technologie Supérieure (ETS), Montréal, QC H3C 1K3, Canada
3
Faculty of Engineering, Lebanese University, Beirut 1003, Lebanon
4
Centre National Intégré du Manufacturier Intelligent (CNIMI), Université du Québec à Trois-Rivières (UQTR), Drummondville, QC J2C 0R5, Canada
*
Authors to whom correspondence should be addressed.
Energies 2024, 17(9), 2053; https://doi.org/10.3390/en17092053
Submission received: 11 April 2024 / Revised: 24 April 2024 / Accepted: 25 April 2024 / Published: 26 April 2024
(This article belongs to the Special Issue The Future of Renewable Energy: 2nd Edition)

Abstract

:
Technical losses in electrical grids are inherent inefficiencies induced by the transmission and distribution of electricity, resulting in energy losses that can reach up to 40% of the generated energy. These losses pose significant challenges to grid operators regarding energy sustainability, reliability, and economic viability. Distributed Energy Resources (DERs) offer promising solutions to lower technical losses by decentralizing energy generation and consumption, reducing the need for long-distance transmission and optimizing grid operation. Hence, estimating the impact of DERs on grid technical losses becomes paramount for grid operators and planners. In response, this article proposes the application of regression modeling and nonlinear curve fitting algorithms to provide a more nuanced understanding and better characterize the intricate interplay between DER deployment and technical losses. Through a comprehensive case study based on more than 1080 computer simulations, we demonstrate the effectiveness of our proposed dynamic polynomial varying coefficient regression model in estimating the impact of DERs on technical losses within electrical grids. The proposed model offers a simple and effective methodology that allows grid operators to gain insights into the nonlinear dynamics of DER integration and make quicker and more informed decisions regarding grid management strategies, infrastructure investments, and policy interventions. Also, this research contributes to advancing the field of grid optimization by offering a simple equation that enhances our ability and haste to assess and mitigate technical losses in the context of an evolving energy landscape characterized by increasing DER adoption.

1. Introduction

In today’s constantly changing energy landscape, integrating distributed energy resources (DERs) is gaining popularity as a vital approach to improving the power grid. These DERs, including solar photovoltaic (PV) panels, wind turbines, battery energy storage systems (BESSs), electrical vehicles (EVs), small power generators, and demand side management (DSM), have the potential to transform the conventional concept of energy generation and consumption. In recent years, DERs have emerged as critical components in transforming traditional electric grids into modern, resilient, and sustainable systems. One of the primary roles of DERs lies in enhancing grid resilience. By providing backup power during outages or natural disasters, DERs contribute to the grid’s reliability [1]. These resources can operate independently or in coordination with the main grid, mitigating disruptions and ensuring uninterrupted electricity supply to consumers. Grid stability is also bolstered by DERs, which provide ancillary services such as frequency regulation and voltage support [2]. Battery storage systems, in particular, can respond rapidly to fluctuations in supply and demand, thereby maintaining grid stability and reliability [3].
Furthermore, DERs facilitate demand response programs, allowing consumers to change their electricity usage in response to price signals or grid conditions [4]. This flexibility helps balance supply and demand, reducing the need for peaker plants and expensive infrastructure upgrades. With the rise of EVs, DERs play a pivotal role in integrating EV charging infrastructure into the grid. Smart charging stations and vehicle-to-grid (V2G) technology enable EV batteries to serve as distributed storage units, supporting grid stability and enhancing overall energy management [5]. In addition, DERs are integral to developing microgrids capable of operating equally in off-grid and on-grid modes. Microgrids enhance energy resilience by providing reliable power to critical infrastructure during emergencies or grid outages.
Last, DERs are crucial in managing peak electricity demand and alleviating grid infrastructure strain [6]. Solar panels, wind turbines, and battery storage systems can generate power during periods of high demand, reducing reliance on traditional power plants and minimizing electricity costs for consumers. Integrating renewable energy sources into the grid is another crucial function of DERs. With the increasing adoption of solar and wind energy, distributed generation technologies enable the efficient utilization of these intermittent resources. By generating electricity near where it is consumed, DERs reduce transmission losses and enhance the overall efficiency of the grid. However, integrating DERs into electrical grids presents opportunities and challenges; their full potential can only be realized if they are properly managed and optimized.
Regarding power management, DER integration has matured using advanced control techniques and accurate forecasting of renewable energy generation, load demand, and market prices. Additionally, on the energy level, DER’s contribution is well remunerated via in-place financing mechanisms such as net metering (NM), Feed-in-tariff (FiT), and power purchase agreement (PPA). Yet, other contributions, such as reducing the electric grid’s technical losses, remain less explored.
The operation of electric grids is inherently plagued by technical and nontechnical losses (Table 1), which result in energy dissipation and reduce overall system efficiency [7]. Technical losses in the electric grid refer to the energy lost during transmission and distribution. These losses occur due to resistance in power lines, transformer inefficiencies, and voltage drops. Technical losses lead to a waste of energy and result in increased costs for energy companies and potential disruptions in electricity supply to end-users. By understanding the causes and patterns of technical losses through data analysis, energy companies can take targeted actions to reduce them. Methodologies used today to assess the technical losses rely on simulation software [8,9], measurement and verification (M&V) studies [10], load flow analysis modeling tools [11,12], and Distribution System State Estimation (DSSE) [13,14,15]. However, despite their numerous advantages, these methodologies can be computationally intensive and time-consuming and require specialized knowledge and training to be used effectively. This is where regression models come into play, offering a powerful tool to analyze historical data and identify factors that contribute to technical losses. Regression models enable energy companies to make data-driven decisions and implement strategies to improve grid efficiency.
This study introduces a dynamic polynomial regression equation with varying coefficients to evaluate the impact of Distributed Energy Resources (DERs) integration across different points of the electric grid. The equation factors in DER power, the connected load, and total impedance at the relevant bus, accurately estimating the reduction in technical losses. This dynamic regression model offers a swift and straightforward approach to assessing such reductions as an alternative to traditional methodologies. It boasts quicker results, demands less computational power, and provides comparable precision to simulation models, all without requiring specialized expertise.
This model holds promise for integration into optimization algorithms and resource allocation strategies, facilitating efficient planning while minimizing technical losses. It can also offer valuable insights into the efficacy of different control strategies. Furthermore, it can aid in the development of targeted and cost-effective control mechanisms and help identify inefficiencies within the grid.
By analyzing data on technical losses, energy companies can pinpoint specific areas or factors contributing to these losses, enabling targeted actions to enhance grid resilience, reduce costs, and improve overall operational efficiency.
In the first part of this paper, and within the literature review context, we define the technical losses and the challenges and limitations of DER penetration level in electrical grids and identify the existing methodologies and techniques used to calculate the technical losses. In the Section 3, we introduce the methodology adopted to develop, train, and test our proposed regression model, and in Section 4, we introduce the built model. Finally, in Section 5, the testing results are presented, and an analysis of the accuracy and precision of the proposed model is provided.

2. Literature Review

The integration of DER into electric grids has emerged as a critical strategy for enhancing grid resilience, sustainability, and efficiency, with publications that address their impact on the grid dating back to the 1970s. The United States Public Utility Regulatory Policies Act (PURPA) [16] encouraged the development of renewable energy and cogeneration facilities by requiring electric utilities to purchase power from qualifying facilities, which included certain types of small-scale, decentralized power generation. This provision effectively incentivized the growth of customer-owned generation, as individuals or businesses could install renewable energy systems and sell excess electricity to utilities at favorable rates.
However, the degree of integration of DERs into the electric grid, defined as the penetration level, remains controversial. A higher DER integration diversifies the energy mix and reduces dependence on fossil fuels. It also helps alleviate stress on the grid during peak periods, reducing the need for costly infrastructure upgrades. However, an excessively high level of DER penetration can pose challenges related to grid stability, voltage regulation, frequency control, and power quality, requiring advanced grid management techniques and infrastructure upgrades [17,18]. Hence, the penetration level of DERs is a key metric used to assess the impact of DERs on the grid and to gauge the extent of their integration into the overall energy system.
Various factors, including technological advancements, regulatory policies, market dynamics, and consumer preferences, influence the penetration level of DERs. The term “DER penetration level” is typically defined as the maximum power capacity of DERs relative to the total installed capacity within a given system. It measures the physical presence of DERs in terms of their maximum potential power output rather than their actual energy production. However, slight variations in how the concept is defined and used may exist in different contexts.
In contrast to capacity penetration, the DER energy penetration’s definition considers the actual energy output or contribution of DERs to the total energy consumed within a specified timeframe (e.g., hourly, daily, annually). It reflects the dynamic nature of DERs’ energy generation or consumption patterns.
Additionally, there are two other definitions: DER market penetration and DER operational penetration. The DER market penetration perspective looks at DER technologies’ market share or adoption rate within a particular market or customer segment. It considers factors such as the number of DER installations, customer participation in DER programs, and the percentage of customers with DER assets. The DER operational penetration definition focuses on the operational impact of DERs on the grid, including their influence on grid stability, reliability, and power quality. It considers factors such as DERs’ ability to provide grid support services, participate in demand response programs, and mitigate grid constraints.
By nature, DERs are located closer to demand centers, hence contributing to reducing electric grid losses by generating electricity closer to where it is consumed. This can mitigate losses incurred during long-distance transmission from centralized power plants. Electric grid technical losses, encompassing both transmission and distribution losses, represent a significant challenge in operating and managing electrical grids worldwide. These losses, arising from conductor resistance, transformer inefficiencies, and voltage drops, contribute to energy wastage, reduced system efficiency, and increased operational costs. Effective modeling of technical losses is essential for optimizing grid performance, enhancing energy efficiency, and ensuring the reliability of the electricity supply [19,20]. This literature review explores various modeling techniques employed to quantify and mitigate electric grid technical losses, encompassing both traditional and advanced approaches.
Traditional methods for modeling electric grid technical losses include empirical formulas, analytical models, and statistical regression approaches. Empirical formulas, such as the I2R formula, estimate losses based on the product of current squared and resistance. Analytical models, such as load flow analysis and network impedance methods [21,22], provide a theoretical framework for calculating losses based on grid topology and operating conditions. Statistical regression approaches, including linear and multiple linear regression, correlate technical losses with load demand, network configuration, and environmental variables. Nevertheless, recent advancements in modeling techniques have expanded the scope and accuracy of electric grid loss estimation. These techniques include machine learning algorithms, such as artificial neural networks (ANNs) [23], support vector machines (SVMs) [24], and decision trees [25]. Machine learning models offer data-driven approaches for predicting technical losses based on historical data and system parameters. These models can capture nonlinear relationships and complex interactions among variables, improving prediction accuracy. Comparably, optimization-based approaches, including genetic algorithms [26], particle swarm optimization [27], and ant colony optimization [28], optimize grid configurations and operational parameters to minimize technical losses by iteratively adjusting system parameters to achieve optimal performance while considering constraints such as voltage limits and equipment ratings. Also, several hybrid models combining multiple techniques exist [29]. A hybrid model may use machine learning for the initial prediction of technical losses and optimization algorithms for fine-tuning system parameters to minimize losses further.
Different methods are also used to allocate the losses to each load connected to the grid. The prorated methodology is a simple empirical method for the electric grid’s loss allocation (pro rata, PR). This methodology aims to equally distribute the grid losses among different components or sections of the distribution system based on their contribution to the overall energy flow, regardless of network configuration and the distance between the generator and the load [30]. The PR approach is straightforward to implement; however, it is unjust for loads located near generating sources. Hence, the distance-adjusted pro rata (DAPR) method may be a more accurate option.
The DAPR considers the distance of the load to a root node and the power demand to define the distance factors. Distance factors are calculated or assigned to each section of the distribution system based on the length of the distribution lines. Longer lines generally have higher losses due to increased resistance, so the distance factor represents the relative contribution of each section to the total line length [31]. The DAPR technique incorporates the distance factor but does not consider the nonlinearity of the power flow.
Another technique that is used is the incremental loss method. By using the linearization in the Newton–Raphson method for solving the power flow, this method calculates the incremental increase in losses caused by the addition of new components, such as transformers, conductors, or other equipment, to the existing system [32,33,34]. By quantifying these incremental losses, utilities and operators can assess the impact of new connections on system efficiency, reliability, and performance. However, high X/R ratio networks may involve complex interactions between resistive and reactive components, making the application of the incremental losses method more challenging. In contrast, network analysis approaches compute losses using the network’s impedance or admittance matrix [35]. Methods under this category include the Z-bus method [36], the modified Y-bus method [37], the branch-current decomposition method [38], and the succinct method [39].
On the other hand, the power tracing method tracks losses through branch power flow and connected nodal injections. Some existing tracing methods use graph theory, as seen in [40,41], while others utilize various algorithms [42,43,44]. Normalization is necessary if the method yields an inaccurate estimate of total losses. A comprehensive overview of loss allocation methods can be found in [45]. Additionally, for radial distribution systems with multiple distributed generation points, article [46] discusses essential considerations such as the characteristics of net generation nodes (power sources) and net load nodes (power sinks), as well as the uncertainty of distributed generator output.
By strategically locating DERs at critical points along the distribution network, utilities can optimize energy delivery, improve system efficiency, and enhance overall grid performance. Optimizing the point of integration of DERs involves selecting the most suitable location within the electrical distribution system to connect DER assets. This optimization aims to maximize the benefits of DERs while minimizing potential negative impacts on system reliability, stability, and performance. Several methods and considerations can be employed to optimize the point of connection (PoC) of DERs. Ehsan and Yang [47] provide an overview of various analytical techniques for optimal integration of DERs in power distribution systems. It reviews a range of analytical methods, including linear programming, nonlinear optimization, heuristic algorithms, metaheuristic algorithms, and artificial-intelligence-based approaches such as genetic algorithms, particle swarm optimization, and ant colony optimization.
Aryani et al. [48] explored the application of genetic algorithms (GAs) to optimize the placement and sizing of DERs in electrical distribution systems, with the main objectives of reducing power losses and improving voltage profiles within the distribution network. Avchat and Mhetre [49] optimized the placement of distributed generation (DG) within distribution networks to achieve loss reduction and voltage improvement objectives using various optimization techniques such as heuristic algorithms, mathematical programming, and evolutionary algorithms to determine the best locations for installing DG units. In the article [50], the authors explore the use of Particle Swarm Optimization (PSO) to determine DG’s optimal sizing and placement in distribution systems.
Thus, numerous case studies and applications demonstrate the effectiveness of modeling techniques for electric grid technical losses across different contexts and scales. These studies encompass urban, rural, and remote grid environments, as well as diverse energy sources and demand profiles. Researchers and practitioners have applied modeling techniques to optimize grid design, improve energy efficiency, integrate renewable energy sources, and improve grid resilience to interruptions and breakdowns. However, despite significant progress, several challenges remain in modeling electric grid technical losses. These challenges include data availability and quality, model complexity, computational requirements, and uncertainty in future scenarios (e.g., demand growth and renewable energy integration).
While simple regression models may not capture all nuances of grid losses’ estimation, their simplicity, interpretability, and efficiency make them valuable tools for initial analysis and as benchmarks for more complex models. Simple regression models often require fewer computational resources, resulting in quicker training and deployment. This can benefit tasks that demand quick predictions in real-time or near-real-time situations. Moreover, simple regression models can be easily scaled to larger data sets and applied to different regions or periods without significant adjustments. This scalability is valuable for analyzing grid losses across diverse geographic areas.
Similarly, basic regression models serve as a benchmark for assessing the effectiveness of more advanced models. By comparing the results of complex models to those of a simple regression model, analysts can assess whether the additional complexity yields significant improvements in accuracy. While dealing with limited data or noisy input variables, complex models are at risk of overfitting. Simple regression models are less susceptible to overfitting, reducing the risk of producing inaccurate estimations. Our proposed model for estimating the technical loss reduction, considering the level of DER penetration at a specific grid node as an independent variable, requires little computational power. Hence, it allows high-frequency sampling, making it perfect for calculating energy losses in real time, considering the fast-changing and intermittent nature of DERs.
Additionally, the proposed model can serve other grid management and optimization algorithms, especially in managing virtual power plants (VPPs), and financially compensates prosumers’ participation in distributed energy generation. Hence, instead of just compensating prosumers for the energy injected back into the grid, our proposed model can calculate the reduction in grid losses resulting from the prosumer’s DER connection to the grid and financially remunerate him for the saved energy. Such a model might be beneficial in increasing the viability of DER deployment in ecosystems where financial incentives, like energy rates and technology prices, are insufficient to incentivize the end users.

3. Methodology

Regression models are statistical tools that analyze the relationship between a dependent variable and one or more independent variables. In the context of energy systems, regression models can be used to predict the performance of DERs based on various factors such as weather conditions, electricity demand, and grid infrastructure.
The main idea is to develop a predictive model capable of calculating the total energy losses in a grid using a DER’s penetration level connected to a specific bus as a predictor.
The DER penetration level X is defined in the following equation:
X = P D E R B u s   L o a d × 100
A significant amount of data are required to build an accurate regression model capable of estimating the impact of DER integration at any grid bus on that grid’s total losses. For this purpose, using ETAP 12.6.0 simulation software, the IEEE-33 bus grid model was used as a reference model to construct a baseline (grid details are provided in Appendix A). The baseline construction methodology was based on the following steps:
Simulate the baseline model without any DER integration, and calculate the total losses of the reference grid model (IEEE-33 bus)
Integrate a certain percentage of DER at a particular bus i, simulate the whole model, and calculate the total grid losses. The total loss output value is then logged.
Increase the DER penetration by a step of 5%, and repeat the simulation and calculation process. Then, this step is repeated until 100% of DER integration at bus i is achieved.
This process is repeated until all busses in our reference model are simulated with DER penetration levels ranging from 0% to 100% with a step of 5%.
Once the learning data are prepared, the regression model can be built. This involves training the model using historical data and optimizing its parameters to minimize the difference between the predicted and actual output. Various algorithms and techniques are available for building regression models, such as ordinary least squares, gradient descent, and machine learning algorithms like random forest and support vector regression. The choice of algorithm depends on the problem’s complexity and data availability. In our case, the least squares method was adopted to develop our regression model. This method aims to minimize the sum of the squares of the differences between the observed and predicted values, also known as residuals. So, the first step is to select a mathematical model that describes the relationship between the independent variables (predictors) and the dependent variable (outcome). For example, in simple linear regression, the model could be represented as Y = β0 + β1X + ε, where Y is the dependent variable, X is the independent variable, β0 and β1 are the coefficients, and ε is the error term. Then, we implement the independent variable data set into the regression equation to generate our predicted data set. Then, we use the gathered data in the observed data set from the constructed baseline to test the precision of the regression model by calculating the difference between the observed and predicted values for each data point. These differences are called residuals. Afterward, we adjust the model coefficients to minimize the sum of the squared residuals; hence, the term “least squares”. To assess how well the model fits the data, we examine various metrics such as the following:
Mean Squared Error (MSE): MSE calculates the average squared difference between predicted and actual values. It penalizes significant errors more heavily than more minor errors.
M S E = 1 n i = 1 n y ^ i y ~ i 2
where y ^ i is the observed value, y ~ i is the predicted value, and n is the number of observations.
Root Mean Squared Error (RMSE): RMSE is the square root of the MSE and represents the average magnitude of the errors in the same units as the dependent variable.
Mean Absolute Error (MAE): MAE calculates the average absolute difference between the predicted and actual values, providing a measure of the average magnitude of the errors.
M A E = 1 n i = 1 n y ^ i y ~ i
The coefficient of determination measures the proportion of the variance in the dependent variable explained by the model’s independent variables. It ranges from 0 to 1, where a higher value indicates a better fit of the model to the data.
R 2 = 1 i = 1 n y ^ i y ~ i 2 i = 1 n y ^ i y ¯ i 2
where y ¯ i is the mean of the observed values.
Mean Percentage Error (MPE): MPE measures the average percentage difference between the predicted and actual values, providing insights into the average directional accuracy of the predictions.
Mean Absolute Percentage Error (MAPE): MAPE calculates the average percentage difference between the predicted and actual values, providing a measure of the overall accuracy of the predictions relative to the observed values.
After the regression model is built, it is important to evaluate its performance to ensure its accuracy and reliability. This can be carried out by comparing the predicted values with the actual values and calculating various metrics such as mean squared error, root mean squared error, and coefficient of determination (R-squared). The performance evaluation provides insights into the effectiveness of the regression model and helps identify any areas for improvement. For this purpose, we decided to use a 14-bus grid model and the IEEE-10 bus model to test the accuracy of the proposed predictive model. Hence, 3 grid models and 1080 simulation results were used to train and test our regression model, as detailed in Figure 1.

4. Model Building

The model was built using an Intel i7-7500U, 2.70 GHz, dual-core CPU with 16 GB RAM, operating on Windows 10. The training and testing data were obtained by simulating three grid models (14-bus grid model, IEEE-10 bus, and IEEE-33 bus models) using the Electrical Transient Analyzer Program version 12.6.0 (ETAP), a software platform used to design, simulate, analyze, control, optimize, and automate electrical power systems. Both models gave practically the same results. The linear regression least square algorithm and nonlinear curve fitting algorithm were programmed using MATLAB R2022a.
After testing various regression equations, we found that a second-degree polynomial equation best fits the data sets from different buses with an acceptable coefficient of determination (R2).
y = a x i 2 + b x i + c
where y is the estimated total technical loss in the electric grid, both variable and fixed as detailed in Table 1, and x i is the DER penetration level at bus i, as defined in Equation (1).
However, it was remarkable that the leading coefficient (a) and the linear coefficient (b) varied from bus to bus. In contrast, the constant coefficient (c) showed minor variations and could be considered constant, as illustrated in Figure 2.
Hence, we needed to identify the independent variables that affect the variabilities of coefficients a and b. It is well known that the load demand (L) directly influences active power losses in the grid. Higher loads increase line currents, leading to higher resistive losses in transmission and distribution lines. Similarly, the impedance of transmission lines contributes to resistive losses as current flows through the conductors. Higher impedance lines experience more significant losses due to increased resistance. Thus, we tried to find a relation between the coefficients a and b from one side and the total impedance Z and load L on the other.
The bus impedance z (in Ohm) is calculated using the resistance r and the reactance x as follows:
z = x 2 + r 2
The total impedance Z of a certain bus constitutes the sum of all the impedances from the utility grid connection to that specific bus. It is defined as the sum of the impedances of all the transmission lines and other components such as transformers, if applicable, that directly connect the power source of the grid to the subject load.
Z = z
To clarify the total impedance formula, let us consider the total impedance for bus 20 of the IEEE-33 bus grid. In this case, the connection between the power source and bus 20 includes bus 0_1, bus 1_18, bus 18_19, and bus 19_20. Hence, Z 20 = z 1 + z 18 + z 19 + z 20 .
Therefore, we moved from a simple polynomial regression model into a dynamic model with varying coefficients. A varying coefficient regression model is considered a type of dynamic regression model. Dynamic regression models are considered when the relationship between the dependent variable and the independent variables changes over time or across different subsets of the data. Varying coefficient regression models fit this definition because they allow the coefficients associated with the independent variables to vary across observations or over time.
In a varying coefficient regression model, the coefficients are not fixed but are allowed to change, which captures the dynamic nature of the relationship between the variables. This flexibility enables the model to capture complex patterns and variations in the data that traditional regression models with fixed coefficients may not adequately capture. Hence, our objective was to find functions f a ( L , Z ) and f b ( L , Z ) that define the variabilities of the leading coefficient (a) and the linear coefficient (b). So, our predictive progression model is presented by Equation (8):
y = f a L i , Z i . x i 2 + f b L i , Z i . x i + c
The values of a i , b i , c i , L i   and   Z i of our training model (IEEE-33 bus) are shown in Table 2.
At first, we tried to find a regression model that predicted the relation between the initial regression model’s coefficients a and b and the Load (L) and total impedance (Z). However, none of the regression models provided an acceptable coefficient of determination. Then, we decided to apply nonlinear curve fitting techniques since, as shown in Figure 3, the relation between the coefficients (a and b) and the load and impedance is nonlinear. A nonlinear curve fitting algorithm is a mathematical technique used to estimate the parameters of a nonlinear model that describes the relationship between variables in the data. Unlike linear regression, which assumes a linear relationship between the independent and dependent variables, nonlinear curve fitting algorithms allow for more complex relationships by fitting a nonlinear function to the data. These algorithms typically involve iteratively adjusting the parameters of the nonlinear model to minimize the difference between the observed and predicted values, often using optimization techniques such as gradient descent, Levenberg–Marquardt algorithm, genetic algorithms, or simulated annealing.
In our case, we used the Levenberg–Marquardt algorithm. It is an optimization technique commonly used for fitting nonlinear least squares curves. The algorithm iteratively adjusts the model’s parameters to minimize the sum of squared differences between the observed and predicted values. The Levenberg–Marquardt algorithm steps are detailed hereafter:
1. Initialization: Start with initial guesses for the model’s parameters.
2. Compute the Jacobian Matrix: Calculate the Jacobian matrix, which contains the partial derivatives of the model equations with respect to each parameter. This matrix provides information about the sensitivity of the model to changes in the parameters.
J i j = r i p j
where p is the parameter vector, and r is the residual vector.
3. Compute the Approximate Hessian Matrix: The Hessian matrix represents the second derivatives of the sum of squared residuals with respect to the parameters. In Levenberg–Marquardt, an approximate Hessian matrix is computed using the Jacobian matrix.
H = J T . J
4. Update the Parameters: Adjust the model parameters using an iterative update rule that considers both the gradient of the objective function (sum of squared residuals) and the curvature of the objective function surface.
5. Control Step Size: The Levenberg–Marquardt algorithm uses a damping parameter, typically denoted as λ (lambda), which controls the step size during each iteration. If the step improves the fit, λ is decreased to allow larger steps. If the step worsens the fit, λ is increased to take smaller steps.
λ k + 1 = λ k × t r a c e ( J T , J ) m × d i a g ( J T . J )
where
λ k is the damping parameter at the k-th iteration;
m is the number of observations (or residuals).
6. Convergence Criterion: Repeat steps 2–5 until a convergence criterion is met, such as reaching a specified tolerance level for the change in the parameters or the objective function.
The Levenberg–Marquardt algorithm was coded in Matlab as a function that tested 269 nonlinear equations (list of equations provided in Appendix B) to find the best nonlinear equation that defines the relation between the initial regression model’s coefficients (a) and (b) and the Load (L) and total impedance (Z). The function calculates the best curve-fitting parameters for each nonlinear equation and the Chi-square (χ2) value. The Chi-square (χ2) is a statistical measure used to evaluate the goodness of fit between observed data and expected values based on a specific model. It is commonly employed in testing the goodness-of-fit tests. The formula for the Chi-square (χ2) is given by Equation (12):
X 2 = i ( O i E i ) 2 E i
where
O i = Observed frequency for category i
E i = Expected frequency for category i
The complete process including the development of the linear regression models using the least square method, as well as using the Levenberg–Marquardt nonlinear curve fitting technique to define the relationship of the leading coefficient (a) and the linear coefficient (b) with respect to the Load (L) and total impedance (Z) is illustrated in Figure 4.
Using the IEEE-33 bus grid model as a learning model, we found the regression equation that allows us to calculate the total grid losses as a function of the DER level of integration at a specific grid node. The equation y = f a L i , Z i . x i 2 + f b L i , Z i . x i + c , is defined as follows:
f a L i , Z i = L i ( 0.64289 × 10 7 0.450559 × 10 6 × Z i )
f b ( L i , Z i ) = 1.472 × 10 4 × Z i × L i
c = I n i t i a l   l o s s e s   o f   t h e   g r i d
If X i is defined as a percentage:
X = P D E R B u s   L o a d
Then, the regression model coefficients will be as follows:
f a L i , Z i = L i ( 642.896 45.0559 × Z i )
f b ( L i , Z i ) = 0.001472 × Z i × L i
The precision of the developed dynamic varying coefficient was tested by comparing the predicted dependent variable outputs (generated by the regression model) to the observed dependent variable values (generated by the simulations using ETAP 12.6.0). The metrics in Section 3 were used to assess the model’s accuracy. Results are presented in Table 3.
Additionally, we compared the performance of our dynamic varying coefficient regression model with a baseline model, which was, in this case, the polynomial regression model developed for each bus separately. The calculated coefficient of determinations for each bus is presented in Figure 5, while the criteria for correlation intensity in regression models according to R2 values is shown in Table 4.
To obtain a more in-depth evaluation of the performance of the proposed dynamic varying coefficient regression model, we decided to calculate the error distribution. The error distribution can help identify any systematic biases in the predictive model. The predictive model error distribution refers to the distribution of errors between the predicted values generated by the model and the actual observed values in the data set. Understanding the error distribution provides insights into the performance and reliability of the predictive model. Additionally, by analyzing the distribution of errors between predicted and actual values, you can understand the spread and variability of prediction errors, thus providing insights into how well the model captures the underlying patterns in the data. A good predictive model error distribution typically exhibits the following characteristics:
  • Symmetry: The error distribution is approximately symmetric around zero, indicating that the model is equally likely to overpredict and underpredict the target variable. A symmetric error distribution suggests that the model is unbiased and does not systematically overestimate or underestimate the outcomes.
  • Normality: The error distribution closely follows a normal (Gaussian) distribution. Normality implies that most prediction errors are minor, with fewer extreme errors. A normal error distribution simplifies interpretation and analysis and is often assumed by many statistical techniques.
  • Constant Variance (Homoscedasticity): The variance of the errors remains relatively constant across different levels of the predictor variables. Homoscedasticity indicates that the model’s predictive performance is consistent across the entire range of the data, and the spread of errors does not systematically change with the magnitude of the predicted values.
  • Zero Mean: The mean of the error distribution is close to zero, indicating that, on average, the model predictions are accurate. A non-zero mean suggests systematic bias in the model predictions, which should be investigated and corrected.
  • No Outliers: The error distribution does not contain extreme outliers or anomalies. Outliers may indicate data points with unusual characteristics or errors in the data collection process. Identifying and addressing outliers is important for improving the overall reliability and performance of the model.
  • Low Dispersion: The dispersion of the error distribution is relatively low, indicating that most prediction errors are concentrated around the mean. Low dispersion suggests that the model provides consistent and precise predictions with little error variability.
The error distribution of our predictive model is shown in Figure 6. The curvature of the error distribution meets the criteria of a normal error distribution curve, as detailed here.

5. Testing and Results Analysis

Following the learning phase, it is important to test the predicting model to evaluate its performance and accuracy using a separate data set not used during the training process. This process helps assess how well the model generalizes to new, unseen data and provides insights into its predictive capabilities. Hence, to obtain a separate data set that was not used during the regression model training, we used two grid models: IEEE-10 bus model and a 14-bus grid model (grid models data are provided in Appendix C and Appendix D). The testing data were generated and processed the same way as the training data, using the same simulation software, ETAP 12.6.0. Then, we used the trained regression model to predict the dependent variable using the same independent variables as the testing data set. And finally, we applied the same comparison criteria to assess the accuracy and precision of our models. Both grid model results are, respectively, shown in Table 5 and Table 6.
The evaluation metrics indicate that the predictive model performs well, with low error values and a high coefficient of determination (R2). Compared to a baseline model using simple polynomial regression, the developed predictive model demonstrates comparable performance, with slightly higher error metrics and lower R2 value (Appendix CTable A5).
While the developed predictive model demonstrates strong performance in predicting the impact of a DER integration on the total technical losses of an electric grid, several limitations should be acknowledged to provide a comprehensive understanding of its capabilities and constraints. Firstly, the model’s predictive accuracy may be influenced by the quality and representativeness of the training data. Hence, a larger data set for training might help improve the predictivity of the model.
Additionally, the analysis of the predicted values across the three tested models showed that the model tends to divert from the baseline model, thus resulting in higher error when the load connected to the bus is high and the total impedance is high. This divergence was less noticed once just one of these factors (high load or high total impedance) was applicable. Lastly, the developed dynamic varying coefficient regression model was not tested for radial topology. Additionally, the model was tested on a symmetrical balanced load, considering a pure sine wave for voltage and current. Hence, further analysis is required to test the validity of the proposed dynamic regression model for non-symmetrical unbalanced loads. Therefore, while the predictive model offers valuable insights into the impact of DER integration on the grid’s technical losses, it should be used judiciously to make informed decisions in real-world scenarios. Ongoing refinement and validation of the model are essential to address these limitations and enhance its robustness and applicability in diverse contexts.

6. Conclusions

As the energy landscape continues to evolve, optimizing DERs and minimizing technical losses in the electric grid will play an increasingly crucial role. Regression models offer a powerful tool for achieving these goals by analyzing vast amounts of data and predicting the performance of DERs. By leveraging the power of predictive modeling, grid operators, energy providers, and end consumers can unlock the true potential of DERs. This will lead to a more sustainable and efficient energy system and pave the way for a greener and more resilient future.
This article introduces a novel regression model for calculating technical losses in electric power distribution systems. Developing accurate methods for estimating technical losses is crucial for utilities to optimize their operations, enhance grid efficiency, and reduce financial losses. The regression model proposed in this study offers a data-driven approach to estimate technical losses based on key system parameters and operational characteristics such as DER integration level, line impedance, and active load connected to a certain bus. One of the strengths of the proposed dynamic varying coefficient regression model is its flexibility and simplicity, making it suitable for easy implementation in optimization algorithms, control loops, or calculation models. The main advantage of such a model is that it does not require considerable computational power or time and can be used without unique technical expertise. Moreover, we have demonstrated the model’s ability to accurately predict technical losses across various scenarios and operating conditions through comprehensive data analysis and regression modeling techniques.
While the regression model shows promising results, it is essential to acknowledge its limitations and areas for further improvement. Future research could focus on enhancing the model’s predictive capabilities by training the model with a more extensive data set and incorporating additional factors such as voltage levels, equipment conditions, and grid modernization initiatives. Furthermore, validation studies in real-world utility settings are valuable for assessing the model’s performance and practical utility.
Overall, the regression model presented in this article represents a significant advancement in technical loss calculation in electric power distribution systems as a function of the DER integration level. By providing utilities with a reliable tool for assessing the impact of DER integration level on technical losses, the model contributes to the optimization of grid operations, better selection of DER integration points for higher grid resilience and efficiency, improved energy management, and, ultimately, the delivery of reliable and cost-effective electricity to consumers. Continued research and collaboration in this area will further advance our understanding of technical losses and support the development of innovative solutions for a sustainable energy future.

Author Contributions

Conceptualization, A.A.; methodology, A.A.; software, S.S.; validation, A.I., M.A. and M.G.; formal analysis, A.A.; investigation, A.A.; resources, H.I.; data curation, A.A.; writing—original draft preparation, A.A.; writing—review and editing, A.I. and M.A.; visualization, A.A.; supervision, M.A.; project administration, H.I.; funding acquisition, A.I. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. IEEE-33 Bus Grid Model Data

Energies 17 02053 i001
Table A1. Buses resistance and reactance data.
Table A1. Buses resistance and reactance data.
Impedance ReferenceConnectionx (Ohm)r (Ohm)Impedance z (Ohm)
Z1Bus 0_10.09220.4770.485829023
Z2Bus 1_20.4930.25110.553263238
Z3Bus 2_30.3660.18640.410732224
Z4Bus 3_40.38110.19410.427682148
Z5Bus 4_50.8190.7071.081947318
Z6Bus 5_60.18720.61880.646496156
Z7Bus 6_71.71141.23512.110535944
Z8Bus 7_81.030.741.268266534
Z9Bus 8_91.0440.741.279662455
Z10Bus 9_100.19660.0650.207066559
Z11Bus 10_110.37440.12380.394337165
Z12Bus 11_121.4681.1551.867899623
Z13Bus 12_130.54160.71290.895297141
Z14Bus 13_140.5910.5260.791174443
Z15Bus 14_150.74630.5450.924115085
Z16Bus 15_161.2891.7212.150200456
Z17Bus 16_170.7320.5740.930215029
Z18Bus 1_180.1640.15650.226689766
Z19Bus 18_191.50421.35542.02477821
Z20Bus 19_200.40950.47840.629727568
Z21Bus 20_210.70890.93731.175189559
Z22Bus 2_220.45120.30830.546470795
Z23Bus 22_230.8980.70911.144214495
Z24Bus 23_240.8960.70111.137698207
Z25Bus 5_250.2030.10340.227816944
Z26Bus 25_260.28420.14490.319007288
Z27Bus 26_271.0590.93371.411834512
Z28Bus 27_280.80420.70061.066573017
Z29Bus 28_290.50750.25850.56954236
Z30Bus 29_300.97440.9361.351129661
Z31Bus 30_310.31050.36190.47684574
Z32Bus 31_320.3410.53020.63039118
Table A2. Buses load and total impedance values.
Table A2. Buses load and total impedance values.
Bus NumberLoad (KW)Load (KVA)Total Impedance Z (Ohm)
11001170.4858
29098.51.0391
31201441.4498
46067.11.8775
56063.22.9595
62002243.6060
72002245.7165
86063.26.9848
96063.28.2644
104554.18.4715
116069.58.8658
126069.510.7337
1312014411.6290
146060.812.4202
156063.213.3443
166063.215.4945
179098.516.4247
189098.50.7125
199098.52.7373
209098.53.3670
219098.54.5422
22901031.5856
234204652.7298
244204653.8675
2560653.1873
2660653.5063
276063.24.9181
281201395.9847
292006326.5542
301501667.9054
312102338.3822
326072.19.0126

Appendix B. Non Linear Fitting Equations

Ref.EquationRef.Equation
1Y = A*(X1^B)*X2^C136Y = (A+X2)/(B+C*X1^2)+D*X1
2Y = X1/(A+B*X2)137Y = A*X1^(B*X2^C)+D*X1
3Y = X2/(A+B*X1)138Y = A*X2^(B*X1^C)+D*X1
4Y = A*X1^(B*X2)139Y = A*X1+B*X2+C*X1^4
5Y = A*X2^(B*X1)140Y = A*X1+B*X1^3+C*X2+D*X2^3
6Y = A*X1^(B/X2)141Y = A*(X1^B)*LOG(X2+C)+D*X2
7Y = A*X2^(B/X1)142Y = A*(X2^B)*LOG(X1+C)+D*X2
8Y = A*X1+B*X2^2143Y = A/X1+B*EXP(C/X2)+D*X2
9Y = A*X2+B*X1^2144Y = A/X2+B*EXP(C/X1)+D*X2
10Y = X1/(A+B*X2^2)145Y = A/X1+B*EXP(C*X2)+D*X2
11Y = X2/(A+B*X1^2)146Y = A/X2+B*EXP(C*X1)+D*X2
12Y = A*(B^X1)*X2^C147Y = A*X1^(B+C*X2)+D*X2
13Y = A*(B^X2)*X1^C148Y = A*X2^(B+C*X1)+D*X2
14Y = A*(X1*X2)^B149Y = A*X1^(B+C/X2)+D*X2
15Y = A*(X1/X2)^B150Y = A*X2^(B+C/X1)+D*X2
16Y = A*(X1/X2)^B+C151Y = A*X1^(B+C*LOG(X2))+D*X2
17Y = A*(B^(1/X1))*X2^C152Y = A*X2^(B+C*LOG(X1))+D*X2
18Y = A*(B^(1/X2))*X1^C153Y = A*X2^(B+C/LOG(X1))+D*X2
19Y = A+B/X1+C/X2^2154Y = A*X1^(B+C/LOG(X2))+D*X2
20Y = A+B/X2+C/X1^2155Y = A*EXP(B*X1+C*X2^2)+D*X2
21Y = A+B*X1+C/X2156Y = A*EXP(B*X2+C*X1^2)+D*X2
22Y = A+B*X2+C/X1157Y = A*EXP(B/X1+C*X2)+D*X2
23Y = A*((X1/B)^C)*EXP(X2/B)158Y = A*EXP(B/X2+C*X1)+D*X2
24Y = A*((X2/B)^C)*EXP(X1/B)159Y = (A+X1)/(B+C*X2)+D*X2
25Y = A+B*X1+C*X2^2160Y = (A+X2)/(B+C*X1)+D*X2
26Y = A+B*X2+C*X1^2161Y = (A+X1)/(B+C*X2^2)+D*X2
27Y = A*(X1^B)*X2^C+D162Y = (A+X2)/(B+C*X1^2)+D*X2
28Y = X1/(A+B*X2)+C163Y = A*X1^(B*X2^C)+D*X2
29Y = X2/(A+B*X1)+C164Y = A*X2^(B*X1^C)+D*X2
30Y = A*X1^(B*X2)+C165Y = A*X1+B*X2+C*X2^4
31Y = A*X2^(B*X1)+C166Y = A*X1+B/X1^3+C*X2+D/X2^3
32Y = A*X1^(B/X2)+C167Y = A*(X1^B)*LOG(X2+C)+D/X2
33Y = A*X2^(B/X1)+C168Y = A*(X2^B)*LOG(X1+C)+D/X2
34Y = A*(B^(1/X1))*X2^C+D169Y = A/X1+B*EXP(C/X2)+D/X2
35Y = A*(B^(1/X2))*X1^C+D170Y = A/X2+B*EXP(C/X1)+D/X2^2
36Y = X1/(A+B*X2^2)+C171Y = A/X1+B*EXP(C*X2)+D/X2
37Y = X2/(A+B*X1^2)+C172Y = A/X2+B*EXP(C*X1)+D/X2^2
38Y = A*(B^X1)*X2^C+D173Y = A*X1^(B+C*X2)+D/X2
39Y = A*(B^X2)*X1^C+D174Y = A*X2^(B+C*X1)+D/X2
40Y = A*((X1/B)^C)*EXP(X2/B)+D175Y = A*X1^(B+C/X2)+D/X2
41Y = A*((X2/B)^C)*EXP(X1/B)+D176Y = A*X2^(B+C/X1)+D/X2
42Y = 1/(A+B*X1+C/X2)177Y = A*X1^(B+C*LOG(X2))+D/X2
43Y = 1/(A+B*X2+C/X1)178Y = A*X2^(B+C*LOG(X1))+D/X2
44Y = A+B*X1+C/X2^2179Y = A*X2^(B+C/LOG(X1))+D/X2
45Y = A+B*X2+C/X1^2180Y = A*X1^(B+C/LOG(X2))+D/X2
46Y = A*X1^(B+C*X2)181Y = A*EXP(B*X1+C*X2^2)+D/X2
47Y = A*X2^(B+C*X1)182Y = A*EXP(B*X2+C*X1^2)+D/X2
48Y = A*X1^(B+C/X2)183Y = A*EXP(B/X1+C*X2)+D/X2
49Y = A*X2^(B+C/X1)184Y = A*EXP(B/X2+C*X1)+D/X2
50Y = A*X1^(B+C*LnX2)185Y = (A+X1)/(B+C*X2)+D/X2
51Y = A*X2^(B+C*LnX1)186Y = (A+X2)/(B+C*X1)+D/X2
52Y = A*X2^(B+C/LnX1)187Y = (A+X1)/(B+C*X2^2)+D/X2
53Y = A*X1^(B+C/LnX2)188Y = (A+X2)/(B+C*X1^2)+D/X2
54Y = A*EXP(B*X1+C*X2^2)189Y = A*X1^(B*X2^C)+D/X2
55Y = A*EXP(B*X2+C*X1^2)190Y = A*X2^(B*X1^C)+D/X2
56Y = A*EXP(B/X1+C*X2)191Y = A*X1+B*X2+C/X2^4
57Y = A*EXP(B/X2+C*X1)192Y = A*X1+B/X1^2+C*X2+D/X2^2
58Y = (A+X1)/(B+C*X2)193Y = A*(X1^B)*LOG(X2+C)+D/X1
59Y = (A+X2)/(B+C*X1)194Y = A*(X2^B)*LOG(X1+C)+D/X1
60Y = (A+X1)/(B+C*X2^2)195Y = A/X1+B*EXP(C/X2)+D/X1^2
61Y = (A+X2)/(B+C*X1^2)196Y = A/X2+B*EXP(C/X1)+D/X1
62Y = A*EXP(B*X1)+C*EXP(D*X2)197Y = A/X1+B*EXP(C*X2)+D/X1^2
63Y = A*(EXP(B*X1)-EXP(C*X2))198Y = A/X2+B*EXP(C*X1)+D/X1
64Y = A*X1^B+C*X2^D199Y = A*X1^(B+C*X2)+D/X1
65Y = A*X1^B+C*EXP(D*X2)200Y = A*X2^(B+C*X1)+D/X1
66Y = A*X2^B+C*EXP(D*X1)201Y = A*X1^(B+C/X2)+D/X1
67Y = A*(X1^B)*(C-X2)^D202Y = A*X2^(B+C/X1)+D/X1
68Y = A*(X2^B)*(C-X1)^D203Y = A*X1^(B+C*LOG(X2))+D/X1
69Y = (A+B*X1^C)/(D+X2^C)204Y = A*X2^(B+C*LOG(X1))+D/X1
70Y = (A+B*X2^C)/(D+X1^C)205Y = A*X2^(B+C/LOG(X1))+D/X1
71Y = (A+B*X1)/(1+C*X2+D*X2^2)206Y = A*X1^(B+C/LOG(X2))+D/X1
72Y = (A+B*X2)/(1+C*X1+D*X1^2)207Y = A*EXP(B*X1+C*X2^2)+D/X1
73Y = A*X1^(B*X2^C)208Y = A*EXP(B*X2+C*X1^2)+D/X1
74Y = A*X2^(B*X1^C)209Y = A*EXP(B/X1+C*X2)+D/X1
75Y = X1/(A+B*X2+C*SQ210Y = A*EXP(B/X2+C*X1)+D/X1
76Y = X2/(A+B*X1+C*SQ211Y = (A+X1)/(B+C*X2)+D/X1
77Y = A*X1^B+C*EXP(D/X2)212Y = (A+X2)/(B+C*X1)+D/X1
78Y = A*X2^B+C*EXP(D/X1)213Y = (A+X1)/(B+C*X2^2)+D/X1
79Y = A*X2^2+B*X2+C*X1+D214Y = (A+X2)/(B+C*X1^2)+D/X1
80Y = A*X1^2+B*X1+C*X2+D215Y = A*X1^(B*X2^C)+D/X1
81Y = A*X1^3+B*X1^2+C*X1+D*X2216Y = A*X2^(B*X1^C)+D/X1
82Y = A*X2^3+B*X2^2+C*X2+D*X1217Y = A*X1+B*X2+C/X1^4
83Y = EXP(A+B/X1+C*LOG(X2))218Y = A*X1+B/X1^3+C*X2+D*Ln(X2)^3
84Y = EXP(A+B/X2+C*LOG(X1))219Y = A*(X1^B)*LOG(X2+C)+D*Ln(X2)
85Y = EXP(A+B/X1+C*LOG(X2))+D220Y = A*(X2^B)*LOG(X1+C)+D*Ln(X2)
86Y = EXP(A+B/X2+C*LOG(X1))+D221Y = A/X1+B*EXP(C/X2)+D*Ln(X2)
87Y = A*(X1^B)*LOG(X2+C)222Y = A/X2+B*EXP(C/X1)+D*Ln(X2)^2
88Y = A*(X2^B)*LOG(X1+C)223Y = A/X1+B*EXP(C*X2)+D*Ln(X2)
89Y = A*(X1^B)*LOG(X2+C)+D224Y = A/X2+B*EXP(C*X1)+D*Ln(X2)^2
90Y = A*(X2^B)*LOG(X1+C)+D225Y = A*X1^(B+C*X2)+D*Ln(X2)
91Y = A/X1+B*EXP(C/X2)+D226Y = A*X2^(B+C*X1)+D*Ln(X2)
92Y = A/X2+B*EXP(C/X1)+D227Y = A*X1^(B+C/X2)+D*Ln(X2)
93Y = A/X1+B*EXP(C*X2)+D228Y = A*X2^(B+C/X1)+D*Ln(X2)
94Y = A/X2+B*EXP(C*X1)+D229Y = A*X1^(B+C*LOG(X2))+D*Ln(X2)
95Y = A*X1^(B+C*X2)+D230Y = A*X2^(B+C*LOG(X1))+D*Ln(X2)
96Y = A*X2^(B+C*X1)+D231Y = A*X2^(B+C/LOG(X1))+D*Ln(X2)
97Y = A*X1^(B+C/X2)+D232Y = A*X1^(B+C/LOG(X2))+D*Ln(X2)
98Y = A*X2^(B+C/X1)+D233Y = A*EXP(B*X1+C*X2^2)+D*Ln(X2)
99Y = A*X1^(B+C*LOG(X2))+D234Y = A*EXP(B*X2+C*X1^2)+D*Ln(X2)
100Y = A*X2^(B+C*LOG(X1))+D235Y = A*EXP(B/X1+C*X2)+D*Ln(X2)
101Y = A*X2^(B+C/LOG(X1))+D236Y = A*EXP(B/X2+C*X1)+D*Ln(X2)
102Y = A*X1^(B+C/LOG(X2))+D237Y = (A+X1)/(B+C*X2)+D*Ln(X2)
103Y = A*EXP(B*X1+C*X2^2)+D238Y = (A+X2)/(B+C*X1)+D*Ln(X2)
104Y = A*EXP(B*X2+C*X1^2)+D239Y = (A+X1)/(B+C*X2^2)+D*Ln(X2)
105Y = A*EXP(B/X1+C*X2)+D240Y = (A+X2)/(B+C*X1^2)+D*Ln(X2)
106Y = A*EXP(B/X2+C*X1)+D241Y = A*X1^(B*X2^C)+D*Ln(X2)
107Y = (A+X1)/(B+C*X2)+D242Y = A*X2^(B*X1^C)+D*Ln(X2)
108Y = (A+X2)/(B+C*X1)+D243Y = A*X1+B*X2+C*Ln(X2)^4
109Y = (A+X1)/(B+C*X2^2)+D244Y = A*X1+B/X1^2+C*X2+D*Ln(X2)^2
110Y = (A+X2)/(B+C*X1^2)+D245Y = A*(X1^B)*LOG(X2+C)+D*Ln(X1)
111Y = A*X1^(B*X2^C)+D246Y = A*(X2^B)*LOG(X1+C)+D*Ln(X1)
112Y = A*X2^(B*X1^C)+D247Y = A/X1+B*EXP(C/X2)+D*Ln(X1)^2
113Y = A*X1+B*X2+C248Y = A/X2+B*EXP(C/X1)+D*Ln(X1)
114Y = A*X1+B*X1^2+C*X2+D*X2^2249Y = A/X1+B*EXP(C*X2)+D*Ln(X1)^2
115Y = A*(X1^B)*LOG(X2+C)+D*X1250Y = A/X2+B*EXP(C*X1)+D*Ln(X1)
116Y = A*(X2^B)*LOG(X1+C)+D*X1251Y = A*X1^(B+C*X2)+D*Ln(X1)
117Y = A/X1+B*EXP(C/X2)+D*X1252Y = A*X2^(B+C*X1)+D*Ln(X1)
118Y = A/X2+B*EXP(C/X1)+D*X1253Y = A*X1^(B+C/X2)+D*Ln(X1)
119Y = A/X1+B*EXP(C*X2)+D*X1254Y = A*X2^(B+C/X1)+D*Ln(X1)
120Y = A/X2+B*EXP(C*X1)+D*X1255Y = A*X1^(B+C*LOG(X2))+D*Ln(X1)
121Y = A*X1^(B+C*X2)+D*X1256Y = A*X2^(B+C*LOG(X1))+D*Ln(X1)
122Y = A*X2^(B+C*X1)+D*X1257Y = A*X2^(B+C/LOG(X1))+D*Ln(X1)
123Y = A*X1^(B+C/X2)+D*X1258Y = A*X1^(B+C/LOG(X2))+D*Ln(X1)
124Y = A*X2^(B+C/X1)+D*X1259Y = A*EXP(B*X1+C*X2^2)+D*Ln(X1)
125Y = A*X1^(B+C*LOG(X2))+D*X1260Y = A*EXP(B*X2+C*X1^2)+D*Ln(X1)
126Y = A*X2^(B+C*LOG(X1))+D*X1261Y = A*EXP(B/X1+C*X2)+D*Ln(X1)
127Y = A*X2^(B+C/LOG(X1))+D*X1262Y = A*EXP(B/X2+C*X1)+D*Ln(X1)
128Y = A*X1^(B+C/LOG(X2))+D*X1263Y = (A+X1)/(B+C*X2)+D*Ln(X1)
129Y = A*EXP(B*X1+C*X2^2)+D*X1264Y = (A+X2)/(B+C*X1)+D*Ln(X1)
130Y = A*EXP(B*X2+C*X1^2)+D*X1265Y = (A+X1)/(B+C*X2^2)+D*Ln(X1)
131Y = A*EXP(B/X1+C*X2)+D*X1266Y = (A+X2)/(B+C*X1^2)+D*Ln(X1)
132Y = A*EXP(B/X2+C*X1)+D*X1267Y = A*X1^(B*X2^C)+D*Ln(X1)
133Y = (A+X1)/(B+C*X2)+D*X1268Y = A*X2^(B*X1^C)+D*Ln(X1)
134Y = (A+X2)/(B+C*X1)+D*X1269Y = A*X1+B*X2+C*Ln(X1)^4
135Y = (A+X1)/(B+C*X2^2)+D*X1

Appendix C. The 14-Bus Grid Model Data

Energies 17 02053 i002
Table A3. Buses resistance and reactance data.
Table A3. Buses resistance and reactance data.
BusImpedance ReferenceR (Ohm)X (Ohm)Z (Ohm)
Bus 1_2Z10.12330.41270.430725179
Bus 2_3Z20.0180.06090.063504409
Bus 3_4Z30.84631.11.38788461
Bus 4_5Z40.69840.60840.926235996
Bus 5_6Z51.5311.6222.230436056
Bus 6_7Z60.90.781.190965994
Bus 7_8Z72.21.052.437724349
Bus 8_9Z84.522.14.984014446
Bus 9_10Z95.536.264982043
Bus 4_11Z100.540.921.066770828
Bus 11_12Z111.291.061.66964068
Bus 9_13Z120.930.81.226743657
Bus 13_14Z132.112.32594067
Table A4. Buses load and total impedance values.
Table A4. Buses load and total impedance values.
Bus ReferenceLoad (KVA)Load (KW)Total Impedance Z (Ohm)
2107710000.430725179
39018500.494229588
4145614001.882114198
5156212002.808350194
6164315305.03878625
77887806.229752245
8105210508.667476593
971270013.65149104
101434142019.91647308
116365502.948885027
12106310204.618525706
1387380014.8782347
1476163517.20417537
Table A5. Regression equation coefficients and coefficient of determination (CoD).
Table A5. Regression equation coefficients and coefficient of determination (CoD).
R2abc
Bus 20.9997731822.35 × 10−5−0.065991116.619
Bus 30.9998126441.91 × 10−5−0.065241116.616
Bus 40.9999985790.000353−0.69581116.585
Bus 50.9999991990.000479−0.933651116.58
Bus 60.9999996940.001858−2.072431116.599
Bus 70.9999531470.000752−1.271431116.459
Bus 80.9999997240.002145−2.373691116.556
Bus 90.9999911090.001346−2.275621117.953
Bus 100.9999931470.012075−5.530861115.824
Bus 110.9986440717.86 × 10−5−0.292561116.588
Bus 120.9986440710.001013−0.677681117.429
Bus 130.9999994550.002947−2.774891116.533
Bus 140.9999994230.002144−2.250151116.548

Appendix D. IEEE-10 Bus Grid Model Data

Energies 17 02053 i003
Table A6. Buses resistance and reactance data.
Table A6. Buses resistance and reactance data.
BusImpedance ReferenceX (Ohm)R (Ohm)Z (Ohm)
Bus 1_2Z10.12330.41270.430725
Bus 2_3Z20.0140.60570.605862
Bus 3_4Z30.74631.2051.417388
Bus 4_5Z40.69840.60840.926236
Bus 5_6Z51.98311.72762.630074
Bus 6_7Z60.90530.78861.200607
Bus 7_8Z72.05521.1642.361936
Bus 8_9Z84.79432.7165.51017
Bus 9_10Z95.34343.02646.14093
Table A7. Buses load and total impedance values.
Table A7. Buses load and total impedance values.
Bus ReferenceLoad (KW)Load (KVA)Total Impedance Z (Ohm)
2184018970.4307
398010370.6059
4179018451.4174
5159824370.9262
6161017182.6301
77507881.2006
8115011522.3619
99809895.5102
10164016526.1409

References

  1. Navidi, T.; El Gamal, A.; Rajagopal, R. Coordinating distributed energy resources for reliability can significantly reduce future distribution grid upgrades and peak load. Joule 2023, 7, 1769–1792. [Google Scholar] [CrossRef]
  2. Hu, C.; Zhang, X.; Wu, Q. Collaborative Active and Reactive Power Control of DERs for Voltage Regulation and Frequency Support by Distributed Event-Triggered Heavy Ball Method. IEEE Trans. Smart Grid 2023, 14, 3804–3815. [Google Scholar] [CrossRef]
  3. Eid, C.; Codani, P.; Perez, Y.; Reneses, J.; Hakvoort, R. Managing electric flexibility from Distributed Energy Resources: A review of incentives for market design. Renew. Sustain. Energy Rev. 2016, 64, 237–247. [Google Scholar] [CrossRef]
  4. Patnam, B.S.K.; Pindoriya, N.M. Demand response in consumer-centric electricity market: Mathematical models and optimization problems. Electr. Power Syst. Res. 2021, 193, 106923. [Google Scholar] [CrossRef]
  5. Kempton, W.; Tomić, J. Vehicle-to-grid power implementation: From stabilizing the grid to supporting large-scale renewable energy. J. Power Sources 2005, 144, 280–294. [Google Scholar] [CrossRef]
  6. Padullaparti, H.; Pratt, A.; Mendoza, I.; Tiwari, S.; Baggu, M.; Bilby, C.; Ngo, Y. Peak Load Management in Distribution Systems Using Legacy Utility Equipment and Distributed Energy Resources. In Proceedings of the 2021 IEEE Green Technologies Conference (GreenTech), Denver, CO, USA, 7–9 April 2021; pp. 435–441. [Google Scholar] [CrossRef]
  7. Agüero, J.R. Improving the efficiency of power distribution systems through technical and nontechnical losses reduction. In Proceedings of the PES T&D 2012, Orlando, FL, USA, 7–10 May 2012; pp. 1–8. [Google Scholar] [CrossRef]
  8. Recalde, D.; Trpovski, A.; Troitzsch, S.; Zhang, K.; Hanif, S.; Hamacher, T. A Review of Operation Methods and Simulation Requirements for Future Smart Distribution Grids. In Proceedings of the 2018 IEEE Innovative Smart Grid Technologies—Asia (ISGT Asia), Singapore, 22–25 May 2018; pp. 475–480. [Google Scholar] [CrossRef]
  9. Sokolova, E.S.; Martynyuk, M.V.; Dmitriev, D.V.; Tyurin, A.I. Optimization of the Parameters of the Distribution Network Computer Model to Reduce Losses. In Proceedings of the 2018 International Multi-Conference on Industrial Engineering and Modern Technologies (FarEastCon), Vladivostok, Russia, 3–4 October 2018; pp. 1–5. [Google Scholar] [CrossRef]
  10. Chen, Z.; Amani, A.M.; Yu, X.; Jalili, M. Control and Optimisation of Power Grids Using Smart Meter Data: A Review. Sensors 2023, 23, 2118. [Google Scholar] [CrossRef] [PubMed]
  11. Farag, H.E.; El-Saadany, E.F.; El Shatshat, R.; Zidan, A. A generalized power flow analysis for distribution systems with high penetration of distributed generation. Electr. Power Syst. Res. 2011, 81, 1499–1506. [Google Scholar] [CrossRef]
  12. Ibrahim, K.A.; Au, M.T.; Gan, C.K.; Tang, J.H. System wide MV distribution network technical losses estimation based on reference feeder and energy flow model. Int. J. Electr. Power Energy Syst. 2017, 93, 440–450. [Google Scholar] [CrossRef]
  13. Kapoor, S.; Blackhall, L.; Sturnaberg, B.; Shaw, M. Distribution System State Estimation With Losses. In Proceedings of the 2021 IEEE PES Innovative Smart Grid Technologies Europe (ISGT Europe), Espoo, Finland, 18–21 October 2021; pp. 1–6. [Google Scholar] [CrossRef]
  14. Nainar, K.; Iov, F. Three-Phase State Estimation for Distribution-Grid Analytics. Clean Technol. 2021, 3, 395–408. [Google Scholar] [CrossRef]
  15. Majdoub, M.; Boukherouaa, J.; Cheddadi, B.; Belfqih, A.; Sabri, O.; Haidi, T. A Review on Distribution System State Estimation Techniques. In Proceedings of the 2018 6th International Renewable and Sustainable Energy Conference (IRSEC), Rabat, Morocco, 5–8 December 2018; pp. 1–6. [Google Scholar] [CrossRef]
  16. Kebir, N.; Maaroufi, M. Technical losses computation for short-term predictive management enhancement of grid-connected distributed generations. Renew. Sustain. Energy Rev. 2017, 76, 1011–1021. [Google Scholar] [CrossRef]
  17. Ashish, V.; Shivya, S. Solar PV Performance Parameter and Recommendation for Optimization of Performance in Large Scale Grid Connected Solar PV Plant—Case Study. J. Energy Power Source 2015, 2, 40–53. [Google Scholar]
  18. Bezerra, U.H.; Soares, T.M.; Vieira, J.P.A.; Tostes, M.E.L.; Manito, A.R.R.; Paye, J.C.H. Equivalent operational impedance: A new approach to calculate technical and nontechnical losses in electric distribution systems. In Proceedings of the 2018 Simposio Brasileiro de Sistemas Eletricos (SBSE), Niteroi, Brazil, 12–16 May 2018; pp. 1–6. [Google Scholar] [CrossRef]
  19. Oureilidis, K.O.; Demoulias, C.S. A decentralized impedance-based adaptive droop method for power loss reduction in a converter-dominated islanded microgrid, Sustainable Energy. Grids Netw. 2016, 5, 39–49. [Google Scholar] [CrossRef]
  20. Monteiro, R.V.A.; Guimarães, G.C.; Silva, F.B.; da Silva Teixeira, R.F.; Carvalho, B.C.; Finazzi, A.D.P.; de Vasconcellos, A.B. A medium-term analysis of the reduction in technical losses on distribution systems with variable demand using artificial neural networks: An Electrical Energy Storage approach. Energy 2018, 164, 1216–1228. [Google Scholar] [CrossRef]
  21. Toma, R.N.; Hasan, M.N.; Nahid, A.-A.; Li, B. Electricity Theft Detection to Reduce Nontechnical Loss using Support Vector Machine in Smart Grid. In Proceedings of the 2019 1st International Conference on Advances in Science, Engineering and Robotics Technology (ICASERT), Dhaka, Bangladesh, 3–5 May 2019; pp. 1–6. [Google Scholar] [CrossRef]
  22. Yao, M.; Zhu, Y.; Li, J.; Wei, H.; He, P. Research on Predicting Line Loss Rate in Low Voltage Distribution Network Based on Gradient Boosting Decision Tree. Energies 2019, 12, 2522. [Google Scholar] [CrossRef]
  23. Queiroz, L.M.O.; Lyra, C. Adaptive Hybrid Genetic Algorithm for Technical Loss Reduction in Distribution Networks Under Variable Demands. IEEE Trans. Power Syst. 2009, 24, 445–453. [Google Scholar] [CrossRef]
  24. Tuzikova, V.; Tlusty, J.; Muller, Z. A Novel Power Losses Reduction Method Based on a Particle Swarm Optimization Algorithm Using STATCOM. Energies 2018, 11, 2851. [Google Scholar] [CrossRef]
  25. Ahuja, A.; Pahwa, A. Using ant colony optimization for loss minimization in distribution networks. In Proceedings of the 37th Annual North American Power Symposium, 2015, Ames, IA, USA, 25 October 2005; pp. 470–474. [Google Scholar] [CrossRef]
  26. Nagi, J.; Yap, K.S.; Tiong, S.K.; Ahmed, S.K.; Mohammad, A.M. Detection of abnormalities and electricity theft using genetic Support Vector Machines. In Proceedings of the TENCON 2008—2008 IEEE Region 10 Conference, Hyderabad, India, 19–21 November 2008; pp. 1–6. [Google Scholar] [CrossRef]
  27. Public Utility Regulatory Policies Act of 1978, 16 USC §2601. Available online: https://www.ferc.gov/media/public-utility-regulatory-policies-act-1978 (accessed on 25 February 2024).
  28. Alam, M.S.; Al-Ismail, F.S.; Salem, A.; Abido, M.A. High-Level Penetration of Renewable Energy Sources Into Grid Utility: Challenges and Solutions. IEEE Access 2020, 8, 190277–190299. [Google Scholar] [CrossRef]
  29. Cole, W.J.; Greer, D.; Denholm, P.; Frazier, A.W.; Machen, S.; Mai, T.; Vincent, N.; Baldwin, S.F. Quantifying the challenge of reaching a 100% renewable energy power system for the United States. Joule 2021, 5, 1732–1748. [Google Scholar] [CrossRef]
  30. Conejo, A.J.; Arroyo, J.M.; Alguacil, N.; Guijarro, A.L. Transmission loss allocation: A comparison of different practical algorithms. IEEE Trans. Power Syst. 2002, 17, 571–576. [Google Scholar] [CrossRef]
  31. Happ, H.H. Cost of wheeling methodologies. IEEE Trans. Power Syst. 1994, 9, 147–156. [Google Scholar] [CrossRef]
  32. Costa, P.M.; Matos, M.A. Loss allocation in distribution networks with embedded generation. IEEE Trans. Power Syst. 2004, 19, 384–389. [Google Scholar] [CrossRef]
  33. Galiana, F.D.; Conejo, A.J.; Kockar, I. Incremental transmission loss allocation under pool dispatch. IEEE Trans. Power Syst. 2002, 17, 26–33. [Google Scholar] [CrossRef]
  34. Carpaneto, E.; Chicco, G.; Sumaili Akilimali, J. Loss partitioning and loss allocation in three-phase radial distribution systems with distributed generation. IEEE Trans. Power Syst. 2008, 23, 1039–1049. [Google Scholar] [CrossRef]
  35. Atanasovski, M.; Taleski, R. Energy Summation Method for Loss Allocation in Radial Distribution Networks With DG. IEEE Trans. Power Syst. 2012, 27, 1433–1440. [Google Scholar] [CrossRef]
  36. Conejo, A.J.; Galiana, F.D.; Kockar, I. Z-bus loss allocation. IEEE Trans. Power Syst. 2001, 16, 105–110. [Google Scholar] [CrossRef]
  37. Parastar, A.; Pirayesh, A.; Mozafari, B.; Khaki, B.; Sirjani, R.; Mehrtash, A. A new method for power loss allocation by modified Y-Bus matrix. In Proceedings of the 2008 IEEE International Conference on Sustainable Energy Technologies, Singapore, 24–27 November 2008; pp. 1184–1188. [Google Scholar] [CrossRef]
  38. Carpaneto, E.; Chicco, G.; Akilimali, J.S. Branch current decomposition method for loss allocation in radial distribution systems with distributed generation. IEEE Trans. Power Syst. 2006, 21, 1170–1179. [Google Scholar] [CrossRef]
  39. Fang, W.L.; Ngan, H.W. Succinct method for allocation of network losses. Gener. Transm. Distrib. IEEE Proc. 2002, 149, 171–174. [Google Scholar] [CrossRef]
  40. Strbac, G.; Kirschen, D.; Ahmed, S. Allocating transmission system usage on the basis of traceable contributions of generators and loads to flows. IEEE Trans. Power Syst. 1998, 13, 527–534. [Google Scholar] [CrossRef]
  41. Lim, V.S.C.; McDonald, J.D.F.; Saha, T.K. Development of a new loss allocation method for a hybrid electricity market using graph theory. Electr. Power Syst. Res. 2009, 79, 301–310. [Google Scholar] [CrossRef]
  42. Rao, M.S.S.; Soman, S.A.; Chitkara, P.; Gajbhiye, R.K.; Hemachandra, N.; Menezes, B.L. Min-max fair power flow tracing for transmission system usage cost allocation: A large system perspective. IEEE Trans. Power Syst. 2010, 25, 1457–1468. [Google Scholar] [CrossRef]
  43. Savier, J.S.; Das, D. An exact method for loss allocation in radial distribution systems. Int. J. Electr. Power Energy Syst. 2012, 36, 100–106. [Google Scholar] [CrossRef]
  44. Jagtap, K.M.; Khatod, D.K. Loss allocation in radial distribution networks with various distributed generation and load models. Int. J. Electr. Power Energy Syst. 2016, 75, 173–186. [Google Scholar] [CrossRef]
  45. Kalambe, S.; Agnihotri, G. Loss minimization techniques used in distribution network: Bibliographical survey. Renew. Sustain. Energy Rev. 2014, 29, 184–200. [Google Scholar] [CrossRef]
  46. Exposito, A.G.; Santos, J.M.R.; Garcia, T.G.; Velasco, E.A.R. Fair allocation of transmission power losses. IEEE Trans. Power Syst. 2000, 15, 184–188. [Google Scholar] [CrossRef]
  47. Ehsan, A.; Yang, Q. Optimal integration and planning of renewable distributed generation in the power distribution networks: A review of analytical techniques. Appl. Energy 2018, 210, 44–59. [Google Scholar] [CrossRef]
  48. Aryani, N.K.; Abdillah, M.; Negara, I.M.Y.; Soeprijanto, A. Optimal placement and sizing of Distributed Generation using Quantum Genetic Algorithm for reducing losses and improving voltage profile. In Proceedings of the TENCON 2011—2011 IEEE Region 10 Conference, Bali, Indonesia, 21–24 November 2011; pp. 108–112. [Google Scholar] [CrossRef]
  49. Avchat, H.S.; Mhetre, S. Optimal Placement of Distributed Generation in Distribution Network Using particle Swarm Optimization. In Proceedings of the 2020 International Conference for Emerging Technology (INCET), Belgaum, India, 5–7 June 2020; pp. 1–5. [Google Scholar] [CrossRef]
  50. Bhumkittipich, K.; Phuangpornpitak, W. Optimal Placement and Sizing of Distributed Generation for Power Loss Reduction Using Particle Swarm Optimization. Energy Procedia 2013, 34, 307–317. [Google Scholar] [CrossRef]
Figure 1. Regression model training and testing process.
Figure 1. Regression model training and testing process.
Energies 17 02053 g001
Figure 2. Regression equations of buses 2 to 5.
Figure 2. Regression equations of buses 2 to 5.
Energies 17 02053 g002
Figure 3. Plot of the leading coefficient (a) and the linear coefficient (b) as a function of L and Z.
Figure 3. Plot of the leading coefficient (a) and the linear coefficient (b) as a function of L and Z.
Energies 17 02053 g003
Figure 4. The full process logic diagram.
Figure 4. The full process logic diagram.
Energies 17 02053 g004
Figure 5. The CoD of the predictive regression model tested on the IEEE-33 bus.
Figure 5. The CoD of the predictive regression model tested on the IEEE-33 bus.
Energies 17 02053 g005
Figure 6. The predictive model error distribution.
Figure 6. The predictive model error distribution.
Energies 17 02053 g006
Table 1. Classification of electric grid losses.
Table 1. Classification of electric grid losses.
Level 1Level 2Level 3Components of Level 3
LossesTechnical LossesFixed LossesHysteresis losses
Eddy current losses
Variable LossesOhmic losses
Nontechnical LossesNetwork Equipment issuesTheft and fraud
Measurement errors
Network information issuesMissing or unregistered connection points
Incorrect location or energization status of connection points
Incorrect information on measurement equipment
Energy data processing issuesEstimation of unmetered consumptions
Estimation of consumptions between meter readings and calculations
Estimation of technical losses
Estimation of detected issues
Other energy data processing issues
Table 2. Values of a, b, c, L, and Z of the buses in the IEEE-33 bus model.
Table 2. Values of a, b, c, L, and Z of the buses in the IEEE-33 bus model.
BusabcZ (Ohm)L (kW)BusabcZ (Ohm)L (kW)
20.0000040.01072110.4858100180.000080.134221116.424790
30.0000050.02312111.039190190.0000090.00532110.712590
40.000010.04332111.4498120200.000010.00922112.737390
50.0000040.02822111.877560210.000010.00982113.367090
60.0000020.04242112.959560220.000010.01072114.542290
70.000070.15132113.6060200230.0000080.02792111.585690
80.00010.19482115.7165200240.00020.16682112.7298420
90.0000080.06392116.984860250.00030.18642113.8675420
100.000020.07072118.264460260.000030.04772113.187360
110.0000080.05482118.471545270.000040.04832113.506360
120.000020.07272118.865860280.0000080.05392114.918160
130.000020.079321110.733760290.000050.11892115.9847120
140.000090.163121111.6290120300.00020.21182116.5542200
150.000030.083521112.420260310.00010.16782117.9054150
160.000030.085521113.344360320.00020.23692118.3822210
170.000030.086321115.494560330.000030.06782119.012660
Table 3. Regression model testing results using IEEE-33 bus grid.
Table 3. Regression model testing results using IEEE-33 bus grid.
IEEE-33 Bus ModelMSERMSEMAEMAPER2
Value3.281.811.120.550.89
Recommendation <10~1
Table 4. Characterization of correlation intensity in regression models according to R2 values.
Table 4. Characterization of correlation intensity in regression models according to R2 values.
R2 ValueCorrelation Intensity
0.00(N) null
(0.00–0.09)(L) low
(0.09–0.36)(M) moderate
(0.36–0.81)(H) high
(0.81–0.98)(VH) very high
1.00(P) perfect
Table 5. Regression model testing results using IEEE-10 bus grid.
Table 5. Regression model testing results using IEEE-10 bus grid.
IEEE-10 Bus Model MSERMSEMAEMAPER2
Value618.2424.8613.332.320.95
Recommendation <10~1
Table 6. Regression model testing results using 14-bus grid model.
Table 6. Regression model testing results using 14-bus grid model.
14-Bus Grid Model MSERMSEMAEMAPER2
Value970.5531.1521.892.240.95
Recommendation <10~1
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.

Share and Cite

MDPI and ACS Style

Aoun, A.; Adda, M.; Ilinca, A.; Ghandour, M.; Ibrahim, H.; Salloum, S. Efficient Modeling of Distributed Energy Resources’ Impact on Electric Grid Technical Losses: A Dynamic Regression Approach. Energies 2024, 17, 2053. https://doi.org/10.3390/en17092053

AMA Style

Aoun A, Adda M, Ilinca A, Ghandour M, Ibrahim H, Salloum S. Efficient Modeling of Distributed Energy Resources’ Impact on Electric Grid Technical Losses: A Dynamic Regression Approach. Energies. 2024; 17(9):2053. https://doi.org/10.3390/en17092053

Chicago/Turabian Style

Aoun, Alain, Mehdi Adda, Adrian Ilinca, Mazen Ghandour, Hussein Ibrahim, and Saba Salloum. 2024. "Efficient Modeling of Distributed Energy Resources’ Impact on Electric Grid Technical Losses: A Dynamic Regression Approach" Energies 17, no. 9: 2053. https://doi.org/10.3390/en17092053

APA Style

Aoun, A., Adda, M., Ilinca, A., Ghandour, M., Ibrahim, H., & Salloum, S. (2024). Efficient Modeling of Distributed Energy Resources’ Impact on Electric Grid Technical Losses: A Dynamic Regression Approach. Energies, 17(9), 2053. https://doi.org/10.3390/en17092053

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop