Online Identification of Distribution Line Parameters by PMUs under Accuracy , Positive Sequence , and Noise Considerations

,is paper proposes an inaccuracy mitigation measure to reduce the error associated with distribution line parameters identification. Additionally, it introduces the concept of positive sequence quantities for determining the line resistance, reactive inductance, and shunt admittance.,e positive sequence-based analysis is required for asymmetrical related studies such as unbalanced fault analysis. ,e paper, also, includes the consideration of noisy distribution networks. It compares the performance of three line parameters identification techniques by using different statistical measures. A total of 12,960 different case studies are simulated and analyzed under six main loading scenarios and four categories with changing line parameters. ,e line parameters are calculated online using voltage and current signals obtained from phasor measurement units (PMUs) placed at the line two terminals. Finally, the study outcomes and the associated recommendations have been summarized for future works considerations.


Introduction
Distribution line (DL) parameters identification forms the basis for distribution power system studies, including dynamic and transient stabilities, state estimate, protection setting, etc. e common practice in the industry, till today, is to determine the parameters using values from design datasheets, manufacture specification sheets, and engineer estimation.
e latter could base the calculation on conductor dimensions, sag, temperature, tower geometries, and other elements.ese elements are used to identify the DL data through different mechanisms such as calculating the geometric mean radius and the geometric mean distance, denoted by GMR and GMD, respectively.Additionally, the official electrical transient analysis program (known as ETAP) model could be utilized to find the DL data, which is an off-line tool.Assumptions and approximations are included in the calculation process steps which reduce the accuracy of results.Basing DL parameter estimation on offline techniques or preidentified information significantly impacts the accuracy level of the power system studies that depend on these values due to the following: (1) Conductor resistance and reactance vary with ambient conditions, conductor situation, and power flow.(2) A number of installed circuits are spliced with other conductors that are different in types and specifications.is represents an inhomogeneity of the line sections.
(3) e overhead conductor arrangement changes due to using different tower configurations and applying the concept of transposition.(4) Cable installation conditions such as grouping, underground, overhead, cable trays, conduits, and submarine, etc., play a major role in line parameter estimation.(5) Cable aging could impact the line parameters due to several factors such as degradation, tension, and life cycle.
e above five factors are sources of conductor impedance and admittance identification errors.With the emergence of PMU technology, it is possible to obtain more accurate data about the system conditions with highfrequency samples along with the corresponding time stamp.Accordingly, it is possible to develop more accurate DL impedance parameters estimation by online measurement techniques using the synchronized PMUs. is online analysis can be used to improve power system operations reliability as detailed below: (a) Power system restoration and reclosing: phasor data is used to bring equipment back into service avoiding the risk of instability or unsuccessful reclosing trials (b) Automated management of voltage and frequency response: the data is used for better system management to frequency and voltage changes (c) Wide-area protection: real-time phasor data allows for improved grid events identification and execution of appropriate system protection measures (d) Planned power system islanding: this is to improve islanding of power system during instability situations (e) Power plant monitoring and integration: real-time data is used for better integration of different power plants that includes intermittent renewables or distributed power sources e majority of research works to estimate the power system line parameters are focused on transmission systems.Numerous techniques have been introduced to calculate the transmission parameters using the synchronized measurement devices.A two-port ABCD parameter identification based technique was introduced in [1]. is method utilized two sets of three samples of sending and receiving terminals' voltage and current signals.is was to find three estimates of ABCD parameters.e ABCD method is referred to in this research work as a "two-port circuit measurement technique".In Reference [2], four methods were discussed to identify short transmission line parameters by synchronized measurements.Reference [3] proposed a novel method to identify transmission line parameters for different cases, including short and long, transposed and untransposed lines with balanced and unbalanced load conditions.e positive sequence line parameters considering the effects of the line shunt capacitance were estimated in [4], employing a two-terminal transmission line model.Likewise, Reference [5] aimed to achieve the same objectives where a new estimation method was presented using synchronized phasor measurements at both line ends.e approach in [6] proposes the use of recursive parameter estimation to find the network branch parameters online and off-line.e least-square technique was leveraged in [7] with the objective of obtaining the line parameters iteratively.
Unlike the abundance of publications on transmission line parameters estimation, the work in distribution is limited.e probability theory, which builds on voltage drop linear equivalent model, was used in [8].
e approach objective was to estimate the DL impedance and get precise parameters.Numerous works discussed the uncertainties of network parameters and inaccuracy of measurements.
In particular, the DL parameters and measurement uncertainties were analyzed in [9].A novel power system uncertainty analysis technique was proposed in [10], where a two-step approach based on static weighted least-squares analysis was used.Reference [11] presents a method to estimate distribution line parameters using only conventional SCADA measurements (voltage magnitude and power measurements).It resulted in a negligible deviation between simulation, experiment, and the actual manufacturer specifications.e key outcomes of the DL parameters estimation studies were that the accuracy of line parameters is crucial for a number of applications including the grid control, stability analysis, and fault location studies.
To the best of the authors' knowledge, the applicability of different methods to identify the sequence DL parameters has not been considered before.is paper proposes the use of PMU to identify the DL parameters under the consideration of accuracy, positive sequence, and noise.e concept of symmetrical components is leveraged to extract the positive sequence of the synchronized phasor voltage and current measurement signals.e online synchronized signals obtained from the PMUs will be used in calculating both the phase and positive sequence DL parameters.In Section 2, three techniques have been developed to measure DL resistance, reactive inductance, and shunt admittance.Section 3 describes the used accuracy statistical measures to evaluate and compare the performance of the three techniques.
e developed case studies along with their results and discussion are presented in Sections 4 and 5, respectively.Finally, the study recommendations and outcomes are stipulated in Section 6. e main data used to support the findings of this study are included within this article.If additional data is required, it could be requested from the corresponding author with proper justification.

Techniques of Distribution Line
Parameters Estimation ree different techniques are discussed in this section with the objectives of identifying the DL parameters.e techniques leverage the PMU voltage and current signals obtained at the two terminals of the line.In order to perform DL parameters estimation, the line is represented in a π-model equivalent circuit as illustrated in Figure 1(a).
e study considers the positive sequence of the voltage and current phasors in addition to the phase values.is aims to explore accuracy enhancement opportunities and compare the results.Additionally, the sequence quantities are required for developing any asymmetrical analysis.e positive sequence equivalent π-model is shown in Figure 1(b).

Ohm's Formula Technique.
e proposed ohm's formula technique (OFT) depends on the ohm's law [12].Under this method, both phase and positive sequence voltage and current phasors are used.is method requires only single set of voltage and current samples of the phasor voltage and current signals produced by PMUs.
e developed OFT equations to calculate the DL parameters are described below: Journal of Electrical and Computer Engineering (1)

Single Measurement Technique.
e proposed single measurement technique (SMT) aims to find DL resistance, reactive inductance and shunt admittance [12].It uses both the phase and positive sequence of the voltage and current signals that are obtained from PMUs at the steady state.e SMT equations are formulated as follows: (2)

Two-Port Circuit Measurement
Technique.e two-port circuit measurement technique (TPCMT) requires two sets of synchronized measurement samples at different loading conditions [12].
e samples are taken from the DL terminals to calculate the two-port circuit parameter known as A, B, C, and D.
e DL impedance and admittance are identified from the ABCD matrix.
e TPCMT is conventionally used to represent transmission lines.Additionally, it provides adequate accuracy for DLs at some cases.Representation of positive sequence TPCMT for DL is shown in Figure 2, where V S1 , V R1 , I R1 , and I S1 are the positive sequence of the sending and receiving ends voltage and current signals, respectively.e following equations form the relation between the sending end and the receiving end quantities: where the parameters A, B, C and D are influenced by the DL resistance, inductance, capacitance, and conductance.e ABDC parameters are complex numbers in which A and D are unit less, B is measured in ohms, and C has a unit of Siemens.e ABCD parameters of the DL equivalent π-model shown in Figure 1 are obtained by the following equations: By applying the Kirchhoff current law (known as KCL) at the sending end, the following equation is obtained: Combining the previous two equations yields Utility  Comparing the last above formula with the ABCD equations yields

Two-port circuit
From the simple DL (only series impedance representation) analysis and derivation B is obtained to be e above A, B, C and D equations are solved to find Z DL1 and Y DL1 which will be as follows: is method could be extended to accommodate two sets of PMU measurements.e two sets could be obtained from two different redundant PMUs or from two readings recorded at different timing or loading conditions.e ABCD equations for the two sets are as follows: e samples of the voltages and currents for the receiving and sending ends are as the following: (i) V S1 ′ , V R1 ′ , I S1 ′ , and I R1 ′ are for the first set (ii) V ″ S1 , V ″ R1 , I ″ S1 , and I ″ R1 are for the second set e ABCD parameters are calculated to account for the two sets to be as follows:

Accuracy Statistical Measures
e accuracy of the proposed methods is evaluated using different statistical measures.is is to ensure that the measures will converge for all case studies analyzed in this paper.at is, in case one statistical measure fails to perform in one of the cases, the evaluation will be achieved by the other measures.

Percentage Error.
e first step toward accepting or rejecting the proposed methods is assessing its accuracy using the percentage error given by the following equation:

Coefficient of Determination.
e coefficient of determination (CoD), denoted by R 2 , is used to indicate the difference of the obtained values by a proposed formula compared to the actual ones.It measures the strength of the proposed formula and benchmarks it with the ideal situation which will result in a coefficient of determination of 100%.It is, also, called the squared error which is the error between the curve obtained by the proposed formula and the actual curve.
e range of coefficient of determination varies between 0 and 1. e higher the number means the proposed formula is more descriptive and reflective to the actual values.Figure 3 is an explanatory sketch for calculating the CoD.
e coefficient of determination equation is formulated as follows: (13) e parameters are described as follows: CoD, coefficient of determination; SE, total square error between the calculated points and the actual values; TV, total variation between the calculated points and the actual values; y act , mean of the actual values; y calc , calculated value; and y act , actual value.

Other Accuracy Statistical Measures.
Other accuracy statistical measures are required to be integrated with the percentage error and CoD. is is due to the fact that the percentage error does not represent the correlation and the CoD has certain shortfalls, especially for small scientific numbers.
e following additional statistical measures are used to evaluate the proposals presented in this paper:

Case Studies
A 25-kV distribution system (refer to Figure 4) is modeled in MATLAB/Simulink to verify the effectiveness of the three line parameters identification techniques.A total of 12,960 different case studies have been performed under six main loading scenarios (stated in Table 1) and four categories (presented in Table 2) with changing of the line parameters.e line parameters have been varied in 60 steps.e loading scenarios and categories considered under this study are tabulated in below tables.e large number of case studies has been developed to test the robustness and accuracy of this paper proposals.e 12,960 simulations differ in the loading conditions, line lengths, noise, and inaccuracy mitigations.
e selected DL is modeled as three-phase DL with a π-type.
e model consists of one set of resistance and inductance elements in series connected between sending and receiving terminals.Two sets of shunt capacitances lumped are, also, included at both ends as illustrated in Figure 1. e initial DL parameters are stated in Table 3. e total series resistance, reactive inductance, and shunt admittance are given by the following formulas, respectively: where R, L, and C are the total DL resistance, inductance, and capacitance, and ℓ is the total length of the line.
In MATLAB, two sets of simulated PMUs are placed at both terminals of the selected DL to measure the voltages and currents waveforms simultaneously.
e recorded waveforms are in the shape of sinusoidal signals and then converted into phasor equivalents.

Results and Discussions
e simulation results of the 12,960 cases are summarized in this section and organized into four categories.Under each category, the resistance, reactive inductance, and shunt admittance are calculated using the three methods for different loading conditions and parameter values.e calculation is based on the voltage and current signals obtained from PMUs that are installed at both ends of the line.Figure 5 shows the voltage and current signals obtained from PMU devices considering noise-free system.

Phase Quantities.
In this category, the phase quantities of voltage and current are used to perform the analysis.is type of analysis is required for asymmetrical related studies such as unbalanced fault analysis.e parameters identification errors of the six loading scenarios are averaged for the three methods.e voltage and current waveforms are assumed to be noise-free.Results of the average errors for the resistance, reactive inductance, and shunt admittance are shown in Figure 6. e maximum errors for each method are stated in Table 4. e results reveal that SMT is more effective in calculating the DL parameters.
e TPCMT shows weakness in calculating the shunt admittance for short lines.is is expected as the method was developed specifically for medium transmission lines.However, it performs very well when the DL length is ranging between 10 and 30 km which is a common sort of DLs.

Positive Sequence Quantities.
Both OFT and SMT have excellent performance in identifying the DL parameters using positive sequence quantities.e average and maximum errors recorded in the simulated studies are presented in Figure 7 and Table 4, separately.It is observed form the results that TPCMT fails to calculate the line parameters using positive sequence voltage and current quantities.
erefore, the results were excluded from Figure 7. e results demonstrate that SMT is superior to OFT in calculating the line parameters using positive sequence quantities.

Phase Quantities with
Noise.Actual voltage and current signals of any distribution system are not pure sinusoidal.Noise is always impeded in the signals due to several factors, e.g., harmonics produced from electronic based devices.e electronic devices could be at residential areas such as televisions, computers, laptops, electronic games, and so on.ere are a number of applications that produce harmonics at the industrial sector, for example capacitor bank, variable frequency drives, and other electronic based equipment.
Accordingly, all input signals to PMUs will be associated with additional harmonics beside the fundamental frequency (60 Hz) as in the Kingdome power system.Although PMU measurements showed an improved accuracy compared to other devices, this performance is not fully materialized in the actual field due to errors from other channels such as instrumentation, CT, and potential transformer (shortly PT) and etc. Figure 8 shows the voltage and current signals obtained from PMU devices considering noisy system.
e OFT and SMT have extraordinary performance when applying the phase values to noisy system (Figure 9).TPCMT still shows weakness in estimating the line parameters, especially for short lines capacitance.As the line length increases as TPCMT converges for identifying the X C .e maximum error recorded in the simulated studies is shown in Table 4. From the calculated average and maximum errors of the six loading scenarios for the three methods considering phasor quantities, it is concluded that SMT is superior to the other techniques for noisy system.

Phase Quantities with Inaccuracy Mitigation for Noisy
Systems.It is observed form the simulated case studies that the error follows specific trend under different line parameters, irrespective of the loading conditions.Knowing the error trend will ease predicting the error magnitude and hence mitigating it.is category proposes to apply inaccuracy mitigation measures to improve the line parameter calculation errors.e measures are developed based on line characteristics and possible loadings.
e proposed inaccuracy mitigation measure concept is illustrated in Figure 10 and given by the following formula: where _ y is the originally calculated value and € y is the enhanced measurement.
e symbol ε is taken from the    10. e inaccuracy mitigation curve could take different shapes based on line loading and characteristics.e proposed concept has been applied to Category 3, and the simulation results are illustrated in Figure 11.e results reveal significant improvements of Category 4 compared to Category 3 in Figure 9.
e inaccuracy mitigation measures will result in accuracy improvement up to 98% of the maximum error of Category 3. e maximum errors for the four categories and six loading scenarios considering the variation of the line parameters are tabulated in Table 4.
e MAD, MSE, RMSE, MAPE, and CoD have been applied to the four categories and six loading scenarios.e results for the latter are averaged into one value for each category and parameter.e results are tabulated in Table 5 to evaluate the robustness of this paper proposals.It is noticed from the table that generally the values under the proposed inaccuracy mitigation measures category (Category 4) are improved compared to those in Category 3. is shows the strength of the proposed inaccuracy mitigation concept which could be applied for ideal and noisy systems.e use of positive sequence quantities will perform very well when using OFTand SMT.However, the phase quantities will result in more accurate line parameters estimation.Unlike OFT and SMT, TPCMT does not function when using the positive sequence values.erefore, ABCD should not be used for any asymmetrical related studies in DLs.
MAPE is found to be the only method applicable for calculating the line shunt admittance since the values of the capacitances are very small scientific numbers.

Conclusions
To carry out any asymmetrical related analysis at DLs such as asymmetrical fault studies, the symmetrical components should be leveraged to identify the positive, negative, and zero sequences.erefore, robust and accurate line parameters calculation techniques are required.Based on that, three line parameters identification techniques have been applied to different case studies and evaluated using different statistical measures.e outcomes of this analysis along with the associated recommendations are as follows: (1) e proposed inaccuracy mitigation concept will result in accuracy improvement up to 98% of the maximum error.erefore, it is recommended to use this concept for any online impedance and admittance calculations using PMUs.(2)   -0.50% 0.00% 0.50%     (5) SMT is superior to OFT in calculating the line parameters using positive sequence quantities.(6) TPCMT does not perform when the line resistance is small (short line) and using phase quantities. is is expected as the method was developed for medium transmission lines.As the line impedance or length increase, TPCMT will boost up its resistance calculation accuracy.(7) It is expected the TPCMT will not perform very well for capacitance identification of short DLs.erefore, it is unrecommended to use this method for short DLs.(8) TPCMT fails to produce result using positive sequence voltage and current signals.erefore, it should not be applied for any asymmetrical studies at the distribution level.(9) Some statistical measures do not function under certain conditions, such as in case of small scientific figures.erefore, there is a need for a wide range of statistical measures to ensure covering all study cases.(10) SMT is ranked to be the most robust technique for identifying all DL parameters under different conditions and OFT comes the second.erefore, it is recommended to use SMT for any distribution related case studies.

Figure 1 :
Figure1: Distribution line equivalent model (π-Type).(a) One-line diagram using the phasor quantities (b) One-line diagram using the positive sequence quantities.e parameters of the above circuits are described as follows: Z S , equivalent impedance at the source side; Z S1 , positive sequence equivalent impedance at the source side; Z R , equivalent impedance at the receiving end; Z R1 , ppositive sequence equivalent impedance at the receiving end; Z DL , distribution line impedance; Z DL1 , positive sequence distribution line impedance; Y DL , distribution line admittance; Y DL1 , positive sequence distribution line admittance; V S , phase voltages at sending end; V S1 , positive sequence phase voltages at sending end; V R , phase voltages at receiving end; V R1 , positive sequence phase voltages at receiving end; I S , Phase current at sending end; I S1 , positive sequence phase current at sending end; I R , phase current at receiving end; I R1 , positive sequence phase current at receiving end.

Figure 2 :
Figure 2: Representation of positive sequence two-port circuit for distribution line.

( 1 )
Mean absolute deviation (MAD), which is the summation of the absolute deviation between the actual and calculated values over the number of records (or the length of the range) (2) Mean square error (MSE), which is considered as the most common error metric.It is mainly the summation of the squared errors over the number of records 4 Journal of Electrical and Computer Engineering (3) Root mean square error (RMSE) is obtained by applying the square root to the MSE (4) Mean absolute percentage error (MAPE) is the average of absolute errors over the actual records

Figure 5 :
Figure 5: PMU current and voltage signals with noise-free.

Figure 6 :
Figure 6: Calculation average errors of the six loading scenarios for the three methods under Category 1.(a) Resistance.(b) Reactive inductance.(c) Shunt admittance.

Figure 7 :Figure 8 :
Figure 7: Calculation average errors of the six loading scenarios for the OFT and SMT techniques under Category 2. (a) Resistance.(b) Reactive inductance.(c) Shunt admittance.

Figure 9 :
Figure 9: Calculation average errors of the six loading scenarios for the three methods under Category 3. (a) Resistance.(b) Reactive inductance.(c) Shunt admittance.

Figure 10 :
Figure 10: Indicative graphs of the proposed inaccuracy mitigation measure concept.(a) Actual vs. calculated line parameters.(b) Proposed inaccuracy mitigation measure.

Figure 11 :
Figure 11: Calculation average errors of the six loading scenarios for the three methods under Category 4. (a) Resistance.(b) Reactive inductance.(c) Shunt admittance.

Table 1 :
e six loading scenarios for simulations.

Table 2 :
e four categories for this study.

Table 3 :
Initial parameters of the distribution line test circuits.

Table 4 :
Maximum errors of the six loading scenarios considering the variation of the line parameters.