Patent Application
20210003626
Method for inaccuracy prediction and mitigation of impedance-based fault location in distribution grids were described by Mustafa M. Al Khabbaz and Mohamed A. Abido in “Online Identification of Distribution Line Parameters by PMUs under Accuracy, Positive Sequence, and Noise Consideration”, Journal of Electrical and Computer Engineering Volume 2018, Article ID 8719670, 12 pages, published on Dec. 20, 2018, https://doi.org/10.1155/2018/8719670—incorporated herein by reference in its entirety.
The present disclosure is directed to a method for inaccuracy characterization and prediction, and mitigation of impedance-based fault location in distribution grids.
Fault location accuracy in distribution systems is an important aspect of fault location studies because the fault location accuracy can expedite line repair and system restoration, and subsequently shorten outage and interruption durations.
The majority of phasor measurement unit (PMU)—based fault location techniques were developed for transmission. Fault location applications of PMUs for distribution grids are very limited. Fault locations in distribution networks have the same criticality level, if not more, compared to transmission fault locations for the following four reasons.
Firstly, in numerous cases, distribution lines (DLs) are installed underground, especially in residential and industrial areas for different justifications. The justifications include safety risks associated with aerial lines, such as risks to people and buildings in cities. Underground installations usually cross roads and pass under third-party territories which makes coordination more complicated. This issue is exacerbated further by the fact that underground circuits require manual patrolling and inspection in case of maintenance or abnormal situations. Patrolling to inspect aerial transmission lines is performed by vehicles.
Secondly, live downed-conductors are common problems in distribution networks, which create a safety risk. Therefore, quick and accurate protection and fault location techniques are necessary.
Thirdly, transmission lines are designed with acceptable redundancy to tolerate planned and forced outages. Therefore, faults may cause line outages but do not necessarily cause power interruption to customers. Unlike transmission, most distribution circuits are designed to be radial type. Accordingly, any fault in the radial configuration DL can cause interruptions. These interruptions may cause a tremendous amount of opportunity loss and actual money loss.
Fourth, distribution system faults account for more than 80% of customer interruptions (See D. Chowdhury, Ali; Koval, Power distribution system reliability: practical methods and applications, Second. Wiley-Blackwell, 2009, incorporated herein by reference in its entirety). Such a big percentage confirms the high priority of fault location in distribution grids.
The above four considerations reflect the need for a swift and accurate fault location technique in distribution networks.
Conventional fault location methods used in distribution networks using PMUs can be organized into following four main categories.
In a first category, impedance-based methods consider fundamental voltage and current signals and network parameters to find the fault locations (See J. Ren, S. S. Venkata, and E. Sortomme, “An Accurate Synchrophasor Based Fault Location Method for Emerging Distribution Systems,” IEEE Trans. Power Deliv., vol. 29, no. 1, pp. 297-298, February 2014; Z. Mengsheng, W. Yi, Z. Zhou, and Z. Li, “Research on fault location based on PMU for multi-source distribution network,” in Asia-Pacific Power and Energy Engineering Conference, APPEEC, 2016, vol. 2016-Decem, pp. 1877-1882; A. Rajeev, T. S. Angel, and F. Z. Khan, “Fault location in distribution feeders with optimally placed PMU's,” in Proceedings of IEEE International Conference on Technological Advancements in Power and Energy, TAP Energy 2015, 2015, pp. 438-442; D. Patynowski et al., “Fault Locator approach for high-impedance grounded or ungrounded distribution systems using synchrophasors,” in 2015 68th Annual Conference for Protective Relay Engineers, CPRE 2015, 2015, pp. 302-310; M. M. Ghalei, H. K. Kargar, and M. G. M. M. G. M. Zanjani, “High impedance fault detection of distribution network by phasor measurement units,” in Electrical Power Distribution Networks (EPDC), 2012 Proceedings of 17th Conference on, 2012, no. Proceedings of 17th Conference on, pp. 1-5; J. Mora-Flòrez, J. Meléndez, and G. Carrillo-Caicedo, “Comparison of impedance based fault location methods for power distribution systems,” Electr. Power Syst. Res., vol. 78, no. 4, pp. 657-666, April 2008; K. Ramar and E. E. Ngu, “Generalized Impedance-Based Fault Location for Distribution Systems,” IEEE Trans. Power Deliv., vol. 27, no. 1, pp. 449-451, 2012; R. Dashti and J. Sadeh, “Accuracy improvement of impedance-based fault location method for power distribution network using distributed-parameter line model,” Int. Trans. Electr. Energy Syst., vol. 24, no. 3, 2014; R. Dashti and J. Sadeh, “Applying Dynamic Load Estimation and Distributed-parameter Line Model to Enhance the Accuracy of Impedance-based Fault-location Methods for Power Distribution Networks,” Electr. Power Components Syst., vol. 41, no. 14, pp. 1334-1362, 2013; S. F. Alwash and V. K. Ramachandaramurthy, “New impedance-based fault location method for unbalanced power distribution systems,” Int. Trans. Electr. Energy Syst., vol. 25, no. 6, pp. 1008-1021, 2015; R. H. Salim, K. C. O. Salim, and A. S. Bretas, “Further improvements on impedance-based fault location for power distribution systems,” IET Gener. Transm. Distrib., vol. 5, no. 4, pp. 467-478, 2011, each incorporated herein by reference in their entirety).
In a second category, travelling wave techniques utilize travelling waves of the voltage and current between the network terminals during the fault (See J. Lee, “Automatic Fault Location on Distribution Networks Using Synchronized Voltage Phasor Measurement Units,” in Volume 2: Simple and Combined Cycles; Advanced Energy Systems and Renewables (Wind, Solar and Geothermal); Energy Water Nexus; Thermal Hydraulics and CFD; Nuclear Plant Design, Licensing and Construction; Performance Testing and Performance Test Codes; St, 2014, vol. 2, p. V002T14A008; A. T. Jahromi, P. Wolfs, and S. Islam, “Travelling wave fault location in rural radial distribution networks to reduce wild fire risk,” in 2015 Australasian Universities Power Engineering Conference: Challenges for Future Grids, AUPEC 2015, 2015, pp. 1-6; A. T. Jahromi, “A Travelling Wave Detector Based Fault Location Device and Data Recorder for Medium Voltage Distribution Systems,” in 2016 Australasian Universities Power Engineering Conference (AUPEC), 2016, pp. 3-7; H. Ye, K. Rui, Z. Zhu, X. Zeng, D. Yang, and Y. Cao, “A novel single-phase grounding fault location method with traveling wave for distribution networks,” in 2015 5th International Conference on Electric Utility Deregulation and Restructuring and Power Technologies (DRPT), 2015, pp. 1175-1179; and J. Zhao, T. He, C.-M. Liu, and K. Li, “Travelling Wave Fault Location for Distribution Line Based on Improved Morphological Gradient Algorithm,” in 2016 International Symposium on Computer, Consumer and Control (IS3C), 2016, pp. 156-159, each incorporated herein by reference in their entirety).
In a third category, wavelets, which are mathematical models, are used for digital signal processing with similar principles to Fourier analysis (See E. S. Tag El Din, M. Gilany, M. M. Abdel Aziz, and D. K. Ibrahim, “An PMU double ended fault location scheme for aged power cables,” in IEEE Power Engineering Society General Meeting, 2005, 2005, pp. 423-429; and A. A. P. Bescaro, R. A. F. Pereira, and J. R. S. Mantovani, “Optimal Phasor Measurement Units Placement for fault location on overhead electric power distribution feeders,” in 2010 IEEE/PES Transmission and Distribution Conference and Exposition: Latin America (T&D-LA), 2010, pp. 37-43, each incorporated herein by reference in their entirety).
In a fourth category are other techniques based on different approaches. The other techniques are uncommonly used for fault location such as state estimation (See M. Pignati, L. Zanni, P. Romano, R. Cherkaoui, and M. Paolone, “Fault Detection and Faulted Line Identification in Active Distribution Networks Using Synchrophasors-Based Real-Time State Estimation,” IEEE Trans. Power Deliv., vol. 32, no. 1, pp. 381-392, February 2017; and R. A. F. Pereira, L. G. W. da Silva, and J. R. S. Mantovani, “PMUs optimized allocation using a tabu search algorithm for fault location in electric power distribution system,” in 2004 IEEE/PES Transmision and Distribution Conference and Exposition: Latin America (IEEE Cat. No. 04EX956), 2004, pp. 143-148, each incorporated herein by reference in their entirety).
In general, there is little effort devoted to impedance-based fault location applications of PMUs in distribution grids. Paper J. Ren et al. (See J. Ren, S. S. Venkata, and E. Sortomme, “An Accurate Synchrophasor Based Fault Location Method for Emerging Distribution Systems,” IEEE Trans. Power Deliv., vol. 29, no. 1, pp. 297-298, February 2014) present a method that works for different networks including active and passive, radial, and looped with ranging precision in order of magnitude of 1%. A new PMU-based method was introduced by Z. Mengsheng et al. (see Z. Mengsheng, W. Yi, Z. Zhou, and Z. Li, “Research on fault location based on PMU for multi-source distribution network,” in Asia-Pacific Power and Energy Engineering Conference, APPEEC, 2016, vol. 2016-Decem, pp. 1877-1882) to locate faults in distribution networks. The method is proven via simulations using MATLAB/Simulink for different types of faults. In Reference A. Rajeev (see A. Rajeev, T. S. Angel, and F. Z. Khan, “Fault location in distribution feeders with optimally placed PMU's,” in Proceedings of IEEE International Conference on Technological Advancements in Power and Energy, TAP Energy 2015, 2015, pp. 438-442), an algorithm for optimal PMU placement is validated for a 7-bus test circuit in addition to IEEE 14 and 30-bus networks. The optimal PMU placement is performed to identify the fault location in both a ring and a radial type of distribution system, which is tested by using 11-bus radial and 14-bus ring networks. Patynowski et al. (see D. Patynowski et al., “Fault Locator approach for high-impedance grounded or ungrounded distribution systems using synchrophasors,” in 2015 68th Annual Conference for Protective Relay Engineers, CPRE 2015, 2015, pp. 302-310) disclose a new approach that collects the data from different sources for high impedance grounded or ungrounded lines. A new algorithm for high impedance fault (HIF) was introduced by Kargar and Zanjani (see M. M. Ghalei, H. K. Kargar, and M. G. M. M. G. M. Zanjani, “High impedance fault detection of distribution network by phasor measurement units,” in Electrical Power Distribution Networks (EPDC), 2012 Proceedings of 17th Conference on 2012, no. Proceedings of 17th Conference on, pp. 1-5). This type of fault is difficult to detect by over-current protection relays because of low fault current. The algorithm is sensitive to any change in current phasor. The best accuracy of fault location studies in distribution networks that is achieved is 99.8% by E. S. Tag El Din et al. (See E. S. Tag El Din, M. Gilany, M. M. Abdel Aziz, and D. K. Ibrahim, “An PMU double ended fault location scheme for aged power cables,” in IEEE Power Engineering Society General Meeting, 2005, 2005, pp. 423-429). The method is based on a wavelet technique. The accuracy of J. Ren et al. (see J. Ren, S. S. Venkata, and E. Sortomme, “An Accurate Synchrophasor Based Fault Location Method for Emerging Distribution Systems,” IEEE Trans. Power Deliv., vol. 29, no. 1, pp. 297-298, February 2014) and Z. Mengsheng et al. (see Z. Mengsheng, W. Yi, Z. Zhou, and Z. Li, “Research on fault location based on PMU for multi-source distribution network,” in Asia-Pacific Power and Energy Engineering Conference, APPEEC, 2016, vol. 2016-December, pp. 1877-1882) are 99%, which is the highest for impedance-based fault location techniques in distribution networks.
A number of factors that influence fault location accuracy when impedance-based techniques are used have not been addressed. For example, some main factors can be identified as follows: a) inaccuracy of identifying fault types for fault-locating algorithms, which affects accuracy due to algorithm prejudices are only effective in certain fault types (i.e., a single line-to-ground fault (LG)); b) unaccepted deviation between the line parameters and the actual parameters, which could result from the total length of the line, and could be identified with errors, although the geometry of conductors and parameters are accurate; c) insufficient accuracy due to changes in the line parameters resulting from different conditions, such as weather conditions (i.e., snow, dusts, winter, and summer, etc.), loading conditions, and aging of the network; d) the uncertainty of the zero-sequence impedance when the zero-sequence impedance is impacted by soil resistivity, which varies under different conditions including the weather (e.g. rainy, sunny, and dry, etc.); e) the lack of accurate line models that take into account transposition, inhomogeneity, shunt admittance, mutual effects, and so on.
Accordingly, fixed line parameters can play a major role in increasing the calculation error. The accurate data of DL parameters improves the precision of locating faults.
The present disclosure provides a method for inaccuracy characterization and prediction, and mitigation of impedance-based fault location in DLs by using PMUs. Swift and accurate fault location can advance restoration of faulty lines in an industry in which fault location depends heavily on manual foot patrols with electronic locators and fault finders (See T. W. Stringfield, D. J. Marihart, and R. F. Stevens, “Fault Location Methods for Overhead Lines,” Trans. Am. Inst. Electr. Eng. Part III Power Appar. Syst., vol. 76, no. 3, pp. 518-529, April 1957, incorporated herein by reference in its entirety). In the present disclosure, superposition is introduced. The superposition can be used in developing an asymmetrical fault location method. In addition, the line parameters are identified online through the voltage and current waves obtained from the installed PMUs, which can overcome the deficiencies associated with the fixed predefined line parameters. Further, a percentage error formula is provided, which can be applied along with other statistical measures to evaluate the effectiveness of proposed model in detecting the fault with high accuracy level under different scenarios. Developed case studies along with their results are also presented.
In the present disclosure, novel impedance-based fault location models for homogeneous distribution grids are proposed. The concept of accuracy improvement, namely inaccuracy mitigation measures (IMMs) have also been considered. The IMMs is introduced to reduce the error associated with the fault location estimation and so expedite power restoration. The model associated with homogeneous lines includes asymmetrical and symmetrical fault classification and location algorithms. The algorithms along with the mitigation measures have been evaluated by using different statistical measures. A total of 7,254 different case studies have been simulated and analyzed for different scenarios, cases, fault types and fault distances.
The inputs of the model are calculated online by using voltage and current signals obtained from PMUs placed at the line two terminals.
According to an aspect of the disclosure, a method for predicting a fault location in a distribution system is provided. In the disclosed method, voltage signals and current signals are obtained from phasor measurement units (PMUs) placed at two terminals of a distribution line (DL) in the distribution system. The two terminals include a sending terminal and a receiving terminal. The voltage signals and current signals of the DL are converted into phasors. A fault type of a faulty line of the DL is subsequently classified in the distribution system based on the converted phasors. The fault type includes a symmetric type and an asymmetric type. A fault location of the faulty line of the DL is predicted based on the fault type through an impedance-based fault location model. Further, inaccuracy mitigation measures on the predicted fault location are applied to improve prediction accuracy.
In some embodiments, when the fault type is the symmetric type, the predicting the fault location of the faulty line of the DL further includes calculating parameters of the faulty lines, and applying an symmetric method of the impedance-based fault location model on the calculated parameters to predict the fault location of the faulty line.
In some embodiments, when the fault type is the asymmetric type, the predicting the fault location of the faulty line of the DL further includes calculating parameters of the faulty lines, determining symmetrical components of the voltage signals and current signals of the DL, identifying superposed quantities of the voltage signals and current signals of the DL, determining an equivalent impedance of a source and an equivalent impedance of a load in the faulty line of the DL, and applying an asymmetric method of the impedance-based fault location model onto the calculated parameters, the identified symmetrical components, the determined superposed quantities, and the determined equivalent impedances to predict the fault location of the faulty line.
In some embodiments, the inaccuracy mitigation measures includes ÿ={dot over (y)}(1+ε({dot over (y)})). {dot over (y)} is an originally predicted value of the fault location based on the impedance-based fault location model. ÿ is an enhanced prediction, and ε({dot over (y)}) is taken from pre-developed inaccuracy mitigation measures in which the originally predicted value {dot over (y)} corresponds to a respective value of ε.
According to another aspect of the disclosure, an apparatus is provided. The apparatus has processing circuitry. The processing circuitry is configured to perform the disclosed method for predicting a fault location in a distribution system.
Aspects of the disclosure also provide a non-transitory computer-readable medium storing instructions which when executed by a computer cause the computer to perform the method mentioned above.
Aspects of the present disclosure are best understood from the following detailed description when read with the accompanying figures. It is noted that, in accordance with the standard practice in the industry, various features are not drawn to scale. In fact, the dimensions of the various features may be arbitrarily increased or reduced for clarity of discussion.
Linear systems obey a superposition principle, which separates the faulty system (post-fault) into two circuits: steady state (pre-fault) and during fault circuits (See J. J. Grainer and W. Stevenson, Power System Analysis. McGraw-Hill, 1994, incorporated herein by reference in its entirety). The superposition principle results in a new circuit called a superposed network (e.g., a superposed distribution system for electrical distribution). The superposed distribution system is converted into a sequence electrical distribution system that is illustrated in
As shown in
Consider dVSk and dISk to be a superposed phase voltage and a phase current at a sending terminal of a distribution line (DL) respectively, and dVRk and dIRk to be a superposed phase voltage and a phase current at a receiving terminal of the DL respectively, where k reflects the phase a, b or c. Applying concepts of symmetrical components transformation method, the following yields:
ΔVSi=M−1×dVSk (1)
ΔVRi=M−1×dVRk (2)
ΔISi=M−1×dISk (3)
ΔIRi=M−1×dIRk (4)
Where, M−1 is a reciprocal of M. M is known as a symmetrical component transformation matrix (SCTM), which converts the phasor quantities into their symmetrical components. The M−1 can be described as follows:
ZSSi and ZSRi are equivalent impedances at the source (sending) and load (receiving) sides (terminals), respectively. Values of the ZSSi and ZSRi can be directly related to a change of generation mode and a load of the distribution system. In the present disclosure, the equivalent impedances ZSSi and ZSRi can be obtained from an administration agency, such as a Power Bureau, and the equivalent impedances ZSSi and ZSRi can change with time. However, it is beneficial to calculate the equivalent impedances ZSSi and ZSRi online in order to apply the equivalent impedances ZSSi and ZSRi in a fault location model to reflect more practical and synchronized values. In the present disclosure, the equivalent impedances ZSSi and ZSRi can be calculated online at the two terminals of the DL by using the superposed sequence voltage and current phasors. From
The equivalent impedances ZSSi and ZSRi can reflect the generation and load modes of the distribution grids during faults.
By utilizing a superposition system configuration such as that shown in
In the present disclosure, to overcome line parameters related influencing factors on the accuracy of the determination of the fault location, the line parameters can be identified online through the distribution line (DL). Identification of the line parameters is based on voltage and current signals obtained from PMUs that are placed at the two terminals of the DL. A ±20% inaccuracy of line parameters could result in substantial errors in determining fault location. The resulted error can reach up to 6.6% as per Joe-Air Jiang et al. (See Joe-Air Jiang, Jun-Zhe Yang, Ying-Hong Lin, Chih-Wen Liu, and Jih-Chen Ma, “An adaptive PMU based fault detection/location technique for transmission lines. I. Theory and algorithms,” IEEE Trans. Power Deliv., vol. 15, no. 2, pp. 486-493, April 2000, incorporated herein by reference in its entirety).
In the present disclosure, a single measurement technique (SMT) can be used to identify a DL resistance, a reactive inductance and a shunt admittance. A DL capacitance can also be considered despite the fact that the DL capacitance does not play the same role as in electrical power studies in electrical power transmission. A consideration of the DL capacitance can ensure that a high accuracy is maintained. The SMT can apply symmetrical components of the voltage and current signals that are obtained from PMUs at the steady state. The SMT equations are formulated as follows:
As shown in equations (9) and (10), ZDLi can be a ith sequence of distribution line impedance. VSi can be a ith sequence of voltage collected from the PMUs at the sending terminal. VRi can be a ith sequence of voltage collected from the PMUs at the receiving terminal. ISi can be a ith sequence of current collected from the PMUs at the sending terminal. IRi can be a ith sequence of current collected from the PMUs at the receiving terminal. YDLi can be a ith sequence of distribution line admittance.
The effectiveness of the SMT method in finding the symmetrical components of DL parameters has been investigated by M. M. Al-khabbaz (See, M. M. Al-khabbaz, “Fault Location In Power Distribution Grids Using Phasor Measurement Units,” King Fahd University Of Petroleum & Minerals, 2018, incorporated herein by reference in its entirety).
In some embodiments, as show in equations (9) and (10), the pre-fault symmetrical components of voltage and current at the sending and receiving ends can be used in the determination of the line parameters. In some embodiments, in case of symmetrical faults, either phase or positive-sequence voltage and current quantities can be used to find the phase line parameters, which is recommended by M. M. Al-khabbaz (See, M. M. Al-khabbaz, “Fault Location In Power Distribution Grids Using Phasor Measurement Units,” King Fahd University Of Petroleum & Minerals, 2018, incorporated herein by reference in its entirety).
The impedance-based fault location method for distribution grids of the invention which uses PMUs is demonstrated in
At step S204, the obtained voltage and current signals can be converted into phasors. For example, a four-sample method (see D. M. Abido, “EE541 Power System Protection II (Course Handouts).” Dhahran, p. 150, 2010) can be used to convert the obtained voltage and current signals into phasors. In some embodiments, the four-sample method is similar to Fourier analysis in term of obtaining the fundamental (60 Hz) component from noisy signals. The advantage of using the four-sample method is its simplicity and no need for using trigonometric functions.
At step S206, a fault type of a faulty line of the DL in the distribution system can be classified based on the converted phasors. The fault type can include a symmetric type and an asymmetric type. In some embodiments, for the symmetrical fault, all phases are short-circuited to each other and often to earth. The symmetric fault is balanced in a sense that the systems remain symmetrical, or the lines are displaced by an equal angle (i.e. 1200 in three phase line). The symmetrical fault is a most severe type of fault involving largest current, but symmetrical faults rarely occur.
An asymmetrical fault can involve one or two phases. In the asymmetrical fault, the three phase lines become unbalanced. The asymmetrical fault can occur between line-to-ground or between lines. The asymmetric fault can include a series fault and/or a shunt fault. The asymmetrical series fault can occur between phases or between phase-to-ground. The asymmetrical shunt fault can be unbalanced in the line impedances. The asymmetric fault can include a single line-to-ground fault (LG), a line-to-line fault (LL), a double line-to-ground fault (LLG), and the symmetric fault can include a three-phase short circuit fault (LLL), and a three-phase-to-ground fault (LLLG).
In the single line-to-ground fault (LG), one conductor comes in contact with the ground or a neutral conductor. The line-to-line fault (LL) occurs when two conductors are short circuited. The double line-to-ground fault (LLG) occurs when two conductors make electrical contact with the ground or come in contact with the neutral conductor.
The process 200 then proceeds to step S208. When the fault is classified as the asymmetric fault, the process 200 proceeds to step S210 where parameters of the faulty line can be calculated. The parameters of the faulty line can include the VSi, VRi, ISi, and IRi, which are described at equations (9) and (10). The parameters of the faulty line can also include dVSk, dISk, dVRk and dIRk, which are described at equations (1)-(4).
The process 200 then proceeds to S212 where the symmetrical components of the DL parameters, ZDLi and YDLi, can be determined based on equations (9) and (10).
At step S214, the superposed quantities ΔVSi, ΔVRi, ΔISi, and ΔIRi of the voltage signals and current signals of the DL can be identified based on equations (1)-(4). Subsequently, the equivalent source impedance ZSSi and the equivalent load impedance ZSRi are determined based on equations (7) and (8) respectively.
The process 200 proceeds to step S218 where the asymmetric impedance-based fault location method can be applied. The asymmetrical impedance-based fault location method leverages a sudden change of the system conditions resulted from the fault. The asymmetrical impedance-based fault location method is used to solve for the symmetrical components of voltages at the sending and receiving ends. The sudden change of voltage and current signals of the faulty DL shown in
The sudden voltage change of the VSFi at the sending terminal and the sudden voltage change of the VRFi at the receiving terminal must be equal at the fault location when calculated from the sending and receiving terminals, respectively, when data of the two terminals is considered to be absolutely synchronous. The absolutely synchronous data can be obtained through synchronous measurement units, such as PMUs, which shows an advantage of applying PMUs for the fault location. Therefore, the asymmetric impedance-based fault location method can be described as follows:
The above fault location method described by (13) can be used to solve for D, which is the distance between the fault point and the sending bus. The solution can, also, be obtained graphically by plotting VSFi and VRFi for the full line length. An intersection point is the solution for the above fault location model and can provide the estimated D. Alternatively, the solution can be found through the following objective function:
D=min(|VSFi|−VRFi|) (14)
The process 200 then proceeds to step S220 where inaccuracy prediction measures can be applied, and subsequently proceeds to step S228 to print out the fault location of the fault.
Still referring to
The process 200 subsequently proceeds to step S224 where the symmetrical impedance-based fault location method can be applied. In general, the symmetrical impedance-based fault location method analysis is less involved compared to the asymmetrical impedance-based fault location method. The symmetrical impedance-based fault location method can be derived based on the principle of voltage drop. First, the method identifies a voltage rise that is calculated from the receiving end by using the receiving end voltage and current and moving toward the sending terminal. Secondly, the method calculates a voltage reduction that moves from the sending end toward the receiving end by using the sending terminal quantities.
The symmetrical impedance-based fault location method performs fault location steps by using post-fault phase quantities, and applying the fault location steps to the pre-fault circuit arrangement (i.e., assuming there is no fault in the DL). The equations developed for symmetrical fault location can be described as follows:
V
SFj
=V
Spj−(ISpj−DVSpjYDL)DZDL (15)
V
RFj
=V
Rpj+(IRpj−(L−D)VRpjYDL)(L−D)ZDL (16)
The calculated fault distance D can be obtained by one of the following equations:
|VSFj|=|VRFj|, solving for D (17)
D=first min(|VSFj|−|VRFj|) (18)
The process 200 then proceeds to step S226 where the inaccuracy prediction measures can be applied, and further to step S228 to print out a fault location of the fault.
In the present disclosure, the steps S220 and S226 of the inaccuracy prediction measures can be demonstrated by
ÿ={dot over (y)}(1+ε({dot over (y)})) (19)
Where {dot over (y)} is an originally calculated value and ÿ is an enhanced estimation. The symbol ε is taken from the pre-developed inaccuracy mitigation measures demonstrated in
The inaccuracy mitigation measures (IMMs) can be organized into two categories. The first category can be a global-based inaccuracy mitigation measure (GIMM), which considers the average of all fault type errors for a specific case and under a specific algorithm. The second category can be a local-based inaccuracy mitigation measure (LIMM), which develops a specific measure for each fault type. The two IMMs can be developed for each of distribution lines as appropriate by the design consultant firms during the design phase of the lines, and part of the protection and coordination studies.
The exemplary electronic device 2300 of
The memory 2350 includes but is not limited to Read Only Memory (ROM), Random Access Memory (RAM), or a memory array including a combination of volatile and non-volatile memory units. The memory 2350 can be utilized as working memory by the controller 2310 while executing the processes and algorithms of the present disclosure. Additionally, the memory 2350 can be used for long-term storage, e.g., of image data and information related thereto.
The electronic device 2300 includes a control line CL and data line DL as internal communication bus lines. Control data to/from the controller 2310 can be transmitted through the control line CL. The data line DL can be used for transmission of voice data, display data, etc.
The antenna 2301 transmits/receives electromagnetic wave signals between base stations for performing radio-based communication, such as the various forms of cellular telephone communication. The wireless communication processor 2302 controls the communication performed between the electronic device 2300 and other external devices via the antenna 2301. For example, the wireless communication processor 2302 can control communication between base stations for cellular phone communication.
The speaker 2304 emits an audio signal corresponding to audio data supplied from the voice processor 2303. The microphone 2305 detects surrounding audio and converts the detected audio into an audio signal. The audio signal can then be output to the voice processor 2303 for further processing. The voice processor 2303 demodulates and/or decodes the audio data read from the memory 2350 or audio data received by the wireless communication processor 2302 and/or a short-distance wireless communication processor 2307. Additionally, the voice processor 2303 can decode audio signals obtained by the microphone 2305.
The exemplary electronic device 2300 can also include a display 2320, a touch panel 2330, an operations key 2340, and an antenna 2306 connected to the short-distance communication processor 2307. The display 2320 can be a Liquid Crystal Display (LCD), an organic electroluminescence display panel, or another display screen technology. In addition to displaying still and moving image data, the display 2320 can display operational inputs, such as numbers or icons which can be used for control of the electronic device 2300. The display 2320 can additionally display a GUI for a user to control aspects of the electronic device 2300 and/or other devices. Further, the display 2320 can display characters and images received by the electronic device 2300 and/or stored in the memory 2350 or accessed from an external device on a network. For example, the electronic device 2300 can access a network such as the Internet and display text and/or images transmitted from a Web server.
The touch panel 2330 can include a physical touch panel display screen and a touch panel driver. The touch panel 2330 can include one or more touch sensors for detecting an input operation on an operation surface of the touch panel display screen. The touch panel 2330 also detects a touch shape and a touch area. Used herein, the phrase “touch operation” refers to an input operation performed by touching an operation surface of the touch panel display with an instruction object, such as a finger, thumb, or stylus-type instrument. In the case where a stylus or the like is used in a touch operation, the stylus can include a conductive material at least at the tip of the stylus such that the sensors included in the touch panel 2330 can detect when the stylus approaches/contacts the operation surface of the touch panel display (similar to the case in which a finger is used for the touch operation).
According to aspects of the present disclosure, the touch panel 2330 can be disposed adjacent to the display 2320 (e.g., laminated) or can be formed integrally with the display 2320. For simplicity, the present disclosure assumes the touch panel 2330 is formed integrally with the display 2320 and therefore, examples discussed herein can describe touch operations being performed on the surface of the display 2320 rather than the touch panel 2330. However, the skilled artisan will appreciate that this is not limiting.
For simplicity, the present disclosure assumes the touch panel 2330 is a capacitance-type touch panel technology. However, it should be appreciated that aspects of the present disclosure can easily be applied to other touch panel types (e.g., resistance-type touch panels) with alternate structures. According to aspects of the present disclosure, the touch panel 2330 can include transparent electrode touch sensors arranged in the X-Y direction on the surface of transparent sensor glass.
The touch panel driver can be included in the touch panel 2330 for control processing related to the touch panel 2330, such as scanning control. For example, the touch panel driver can scan each sensor in an electrostatic capacitance transparent electrode pattern in the X-direction and Y-direction and detect the electrostatic capacitance value of each sensor to determine when a touch operation is performed. The touch panel driver can output a coordinate and corresponding electrostatic capacitance value for each sensor. The touch panel driver can also output a sensor identifier that can be mapped to a coordinate on the touch panel display screen. Additionally, the touch panel driver and touch panel sensors can detect when an instruction object, such as a finger is within a predetermined distance from an operation surface of the touch panel display screen. That is, the instruction object does not necessarily need to directly contact the operation surface of the touch panel display screen for touch sensors to detect the instruction object and perform processing described herein. Signals can be transmitted by the touch panel driver, e.g. in response to a detection of a touch operation, in response to a query from another element based on timed data exchange, etc.
The touch panel 2330 and the display 2320 can be surrounded by a protective casing, which can also enclose the other elements included in the electronic device 2300. According to aspects of the disclosure, a position of the user's fingers on the protective casing (but not directly on the surface of the display 2320) can be detected by the touch panel 2330 sensors. Accordingly, the controller 2310 can perform display control processing described herein based on the detected position of the user's fingers gripping the casing. For example, an element in an interface can be moved to a new location within the interface (e.g., closer to one or more of the fingers) based on the detected finger position.
Further, according to aspects of the disclosure, the controller 2310 can be configured to detect which hand is holding the electronic device 2300, based on the detected finger position. For example, the touch panel 2330 sensors can detect a plurality of fingers on the left side of the electronic device 2300 (e.g., on an edge of the display 2320 or on the protective casing), and detect a single finger on the right side of the electronic device 2300. In this exemplary scenario, the controller 2310 can determine that the user is holding the electronic device 2300 with his/her right hand because the detected grip pattern corresponds to an expected pattern when the electronic device 2300 is held only with the right hand.
The operation key 2340 can include one or more buttons or similar external control elements, which can generate an operation signal based on a detected input by the user. In addition to outputs from the touch panel 2330, these operation signals can be supplied to the controller 2310 for performing related processing and control. According to aspects of the disclosure, the processing and/or functions associated with external buttons and the like can be performed by the controller 2310 in response to an input operation on the touch panel 2330 display screen rather than the external button, key, etc. In this way, external buttons on the electronic device 2300 can be eliminated in lieu of performing inputs via touch operations, thereby improving water-tightness.
The antenna 2306 can transmit/receive electromagnetic wave signals to/from other external apparatuses, and the short-distance wireless communication processor 2307 can control the wireless communication performed between the other external apparatuses. Bluetooth, IEEE 802.11, and near-field communication (NFC) are non-limiting examples of wireless communication protocols that can be used for inter-device communication via the short-distance wireless communication processor 2307.
The electronic device 2300 can include sensors 2308. The motion sensors 2308 can include PMUs to obtain the voltage and current signals of the DL.
Electronic device 2300 can include a data processor 2309, which is configured to receive the inputs from the operator via the touch panel 2330 and the signals from the sensors 2308, and operate a method, such as the method described by the process 200, to determine the fault location based on the signals and inputs.
As shown in
Further, the claimed advancements may be provided as a utility application, background daemon, or component of an operating system, or combination thereof, executing in conjunction with CPU 2401 and an operating system such as Microsoft® Windows®, UNIX®, Oracle® Solaris, LINUX®, Apple macOS® and other systems known to those skilled in the art.
In order to achieve the computer 2400, the hardware elements may be realized by various circuitry elements, known to those skilled in the art. For example, CPU 2401 may be a Xenon® or Core® processor from Intel Corporation of America or an Opteron® processor from AMD of America, or may be other processor types that would be recognized by one of ordinary skill in the art. Alternatively, the CPU 2401 may be implemented on an FPGA, ASIC, PLD or using discrete logic circuits, as one of ordinary skill in the art would recognize. Further, CPU 2401 may be implemented as multiple processors cooperatively working in parallel to perform the instructions of the inventive processes described above.
The computer 2400 in
The computer 2400 further includes a display controller 2408, such as a NVIDIA® GeForce® GTX or Quadro® graphics adaptor from NVIDIA Corporation of America for interfacing with display 2410, such as a Hewlett Packard® HPL2445w LCD monitor. A general purpose I/O interface 2412 interfaces with a keyboard and/or mouse 2414 as well as an optional touch screen panel 2416 on or separate from display 2410. General purpose I/O interface 2412 also connects to a variety of peripherals 2418 including printers and scanners, such as an OfficeJet® or DeskJet® from Hewlett Packard.
The general purpose storage controller 2420 connects the storage medium disk 2404 with communication bus 2422, which may be an ISA, EISA, VESA, PCI, or similar, for interconnecting all of the components of the computer 2400. A description of the general features and functionality of the display 2410, keyboard and/or mouse 2414, as well as the display controller 2408, storage controller 2420, network controller 2406, and general purpose I/O interface 2412 is omitted herein for brevity as these features are known.
In
Referring again to
The PCI devices can include, for example, Ethernet adapters, add-in cards, and PC cards for notebook computers. The Hard disk drive 2560 and CD-ROM 2566 can use, for example, an integrated drive electronics (IDE) or serial advanced technology attachment (SATA) interface. In one implementation the I/O bus can include a super I/O (SIO) device.
Further, the hard disk drive (HDD) 2560 and optical drive 2566 can also be coupled to the SB/ICH 2520 through a system bus. In one implementation, a keyboard 2570, a mouse 2572, a parallel port 2578, and a serial port 2576 can be connected to the system bus through the I/O bus. Other peripherals and devices can be connected to the SB/ICH 2520 using a mass storage controller such as SATA or PATA, an Ethernet port, an ISA bus, a LPC bridge, SMBus, a DMA controller, and an Audio Codec.
The mobile device terminals can include a cell phone 2710, a tablet computer 2712, and a smartphone 2714, for example. The mobile device terminals can connect to a mobile network service 2720 through a wireless channel such as a base station 2756 (e.g., an Edge, 3G, 4G, or LTE Network), an access point 2754 (e.g., a femto cell or WiFi network), or a satellite connection 2752. In one implementation, signals from the wireless interface to the mobile device terminals (e.g., the base station 2756, the access point 2754, and the satellite connection 2752) are transmitted to a mobile network service 2720, such as an EnodeB and radio network controller, UMTS, or HSDPA/HSUPA. Mobile users' requests and information are transmitted to central processors 2722 that are connected to servers 2724 to provide mobile network services, for example. Further, mobile network operators can provide service to mobile users for authentication, authorization, and accounting based on home agent and subscribers' data stored in databases 2726, for example. The subscribers' requests are subsequently delivered to a cloud 2730 through the Internet.
A user can also access the cloud through a fixed terminal 2716, such as a desktop or laptop computer or workstation that is connected to the Internet via a wired network connection or a wireless network connection. The mobile network service 2720 can be a public or a private network such as an LAN or WAN network. The mobile network service 2720 can be wireless such as a cellular network including EDGE, 3G and 4G wireless cellular systems. The wireless mobile network service 2720 can also be Wi-Fi, Bluetooth, or any other wireless form of communication that is known.
The user's terminal, such as a mobile user terminal and a fixed user terminal, provides a mechanism to connect via the Internet to the cloud 2730 and to receive output from the cloud 2730, which is communicated and displayed at the user's terminal. In the cloud 2730, a cloud controller 2736 processes the request to provide users with the corresponding cloud services. These services are provided using the concepts of utility computing, virtualization, and service-oriented architecture.
In one implementation, the cloud 2730 is accessed via a user interface such as a secure gateway 2732. The secure gateway 2732 can for example, provide security policy enforcement points placed between cloud service consumers and cloud service providers to interject enterprise security policies as the cloud-based resources are accessed. Further, the secure gateway 2732 can consolidate multiple types of security policy enforcement, including for example, authentication, single sign-on, authorization, security token mapping, encryption, tokenization, logging, alerting, and API control. The cloud 2730 can provide to users, computational resources using a system of virtualization, wherein processing and memory requirements can be dynamically allocated and dispersed among a combination of processors and memories to create a virtual machine that is more efficient at utilizing available resources. Virtualization creates an appearance of using a single seamless computer, even though multiple computational resources and memories can be utilized according to increases or decreases in demand. In one implementation, virtualization is achieved using a provisioning tool 2740 that prepares and equips the cloud resources, such as the processing center 2734 and data storage 2738 to provide services to the users of the cloud 2730. The processing center 2734 can be a computer cluster, a data center, a main frame computer, or a server farm. In one implementation, the processing center 2734 and data storage 2738 are collocated.
Embodiments described herein can be implemented in conjunction with one or more of the devices described above with reference to
In the present disclosure, the accuracy of the proposed methods is evaluated using different statistical measures. The different statistical measures are applied to ensure that the measures can converge for various case studies disclosed in the present disclosure. That is, in case one statistical measure fails to perform in one of the case studies, the evaluation can be achieved by other measures.
In the present disclosure, a first step toward accepting or rejecting the disclosed method is assessing the accuracy of the method by using a percentage error given by the following equation:
Secondly, the following statistical measures can be used to evaluate the performance of the impedance-based fault location method: a) the coefficient of determination (CoD), denoted by R2. The range of coefficient of determination varies between 0 and 1. The higher the number used in the proposed formula the more descriptive and reflective it is of actual values. b) 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). c) 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. e) root mean square error (RMSE), obtained by applying the square root to the RMSE. f) mean absolute percentage error (MAPE), which is the average of absolute errors over the actual records.
In the present disclosure, the effectiveness of the impedance-based fault location method can be validated by using 25 kV distribution systems shown in
The large number of case studies is developed to test the robustness and accuracy of the proposed fault location model. The 7,254 simulations/case studies differ in the line lengths, fault locations, loading conditions, line parameters and applied IMMs.
The selected DL can be described as a structure of three-phase DL with a n-type. The structure includes of one set of resistance and inductance elements that are in series connected between sending and receiving terminals. Two sets of shunt capacitances that are lumped are also included at both ends as illustrated in
In MATLAB, two sets of simulated PMUs are placed at both terminals of the selected DL to measure the voltages and currents waveforms simultaneously. The recorded waveforms are in the shape of sinusoidal signals and then converted into phasor equivalents.
In the present disclosure, the simulation results of the 7,254 case studies can be organized into two categories: asymmetrical and symmetrical. Under each category, three simulation scenarios, which are shown in Table 2, are considered. In addition, three inaccuracy mitigation measures are applied: 1) case 1: base case, which includes the results obtained directly from the proposed fault location methods shown in
In the present disclosure, one inaccuracy mitigation measure is applied to the symmetrical impedance-based fault location method where the three phases and three phases to ground faults use the same balanced faults analysis.
The calculation is performed based on the voltage and current signals obtained from PMUs that are installed at both ends of the line.
Table 2 provides three scenarios that can be applied to three types of unbalanced/asymmetric faults (LG, LL and LLG). As shown in Table 2, the faults can be placed at different distances starting from 0.50 km to approximately the total line length with a step size of 0.5 km. In each fault distance, the fault resistance has been varied to be 0.001 Ω, 0.01Ω, 0.1Ω, 1Ω, 10Ω and 100Ω. The fault location at each distance has been averaged for the aforementioned fault resistance values.
Simulated results can be further applied with the three inaccuracy mitigation measures as mentioned above: 1) case 1: without IMM; 2) case 2: GIMM; and 3) case 3: LIMM.
The scenarios of Table 2 can be simulated by the proposed asymmetrical method for different fault distances and resistances. The average errors of fault resistances for fault distances and unbalanced fault types are demonstrated in
The results in
Table 3 provides statistical measures results for all asymmetrical-based simulations (in %) that are categorized into fault types (LG, LL, and LLG).
Table 4 provides statistical measures results for all symmetrical-based simulations (in %) that are categorized into three scenarios.
The maximum and average errors for the three scenarios considering the variation of both the distance and resistance of fault are displayed in Table 5. It is observed that the line length decreases as the maximum and average errors drops. The best numbers are obtained for the LG and LL and then for LLG. The base case overall average error for all scenarios is 0.62%, which reflects the accuracy of proposed asymmetrical fault location method by using PMUs.
The GIMM considers the average of all unbalanced fault types errors for a specific scenario. This is under different loading conditions and the change of both distance and resistance of the fault. The results are averaged for fault resistances at specific distance and shown in
Global inaccuracy prediction and mitigation has tremendously enhanced fault location determination accuracy. The inaccuracy reduction could reach up to 96% of the base case. This is observed in the maximum error of LL fault under Scenario 2 as shown in Table 5 where the maximum error has been reduced from 1.12% to 0.0500 (above 22 times of reduction). In simple terminologies, the error for the base case is equivalent to 224 meters. By introducing the global inaccuracy prediction and mitigation, the uncertainty can drop to only 10 meters.
The overall average error obtained by the GIMM for all scenarios is 0.030%, which is over 20 times less than the base case. This proves the robustness and effectiveness of the disclosed GIMM.
The LIMM is almost similar to the global IMM except that the LIMM uses a specific measure for each type of faults. The LIMM is developed based on the change of both distance and resistance of the fault. The results are averaged for fault resistances at specific distance and presented in
The concept of local inaccuracy prediction and mitigation has resulted in a similar enhancement in the accuracy of fault location determination as is the case for global inaccuracy prediction and mitigation. The LIMM can achieve same level of reduction in inaccuracy as the GIMM, for example more than 22 times comparing to the case 1 (base case). The overall average error obtained by the LIMM measure for all scenarios is 0.027%. This number is less than both the base and global cases by 22.38 and 0.13 times, respectively. This proves the robustness and effectiveness of the LIMM. The symmetrical algorithms described in
The results are categorized into two cases: case 1 (base case) and case 2 (inaccuracy mitigation measures, IMM). The global IMM will be the same as the local IMM as the three phases and three phases to ground faults use the same balanced faults analysis.
The scenarios of Table 2 can be simulated by the symmetrical impedance-based fault location method for different fault distances and resistances. The average errors of fault resistances for fault distances and unbalanced fault types are demonstrated in
The results shown in
The maximum and average errors for the three scenarios considering the variation of both the distance and resistance of fault can be displayed in Table 6. It is observed that the line length decreases as the maximum and average errors drops. The base case overall average error for all scenarios is 0.067%. This very low number reflects the high accuracy of proposed symmetrical fault location method by using PMUs.
When the inaccuracy mitigation measures are applied, the resulted calculation errors of varying fault resistances for the three scenarios by using the proposed inaccuracy mitigation measure can be displayed in
The application of inaccuracy prediction and mitigation to the determination of fault location has provided an enhancement of the fault location accuracy. The inaccuracy reduction could reach up to 43% of the base case maximum errors and 29% of the average of the base case maximum errors. This confirms the robustness and effectiveness of the proposed inaccuracy mitigation measure.
Table 6 demonstrates the maximum and average errors for the three scenarios considering the variation of both distance and resistance of the faults.
In order to perform the model evaluation, the MAD, MSE, RMSE, MAPE and CoD can be applied to the three scenarios, three cases (i.e., case 1, case 2 and case3) and different resistance values. The results for the latter are averaged into one value for each scenario and each case. The results are tabulated in Table 3 and Table 6 to evaluate the robustness of the proposals. It is noticed from the Tables 3 and 6 that generally the values under the proposed IMMs categories are improved compared to those in the base case (Case1). This shows the strength of the proposed inaccuracy mitigation concept. Case 3 is superior to all cases and can provide very accurate results.
With the consideration of maximum and average errors presented in Table 5, the proposed IMMs concept is between 99.81% and 99.97%. The accuracy could reach more than four 9s if the minimum errors are considered.
Considering maximum and average errors shown in Table 6, the proposed IMMs concept is between 99.74% and 99.95%. The accuracy could reach more than four 9s if the minimum errors are considered.
The present disclosure provides a method and system to advance DL repair and system restoration which ensures prompt and accurate fault location determination. An impedance-based fault location method was developed and evaluated for different simulations by statistical measures. The method included an asymmetric impedance-based fault location method and a symmetric impedance-based fault location method. Additionally, global and local IMMs were proposed to improve the fault location. In the present disclosure, the proposed IMMs can significantly improve the fault location estimate up to 23 times compared to the base case without IMMs. The accuracies of the IMMs that are characterized through the maximum errors are 99.81% and 99.74% for the asymmetrical and symmetrical respectively. By applying the IMMs, fault location can be identified within a few meters of the exact location if not at the actual location. These accuracies demonstrate the effectiveness and the superiority of the proposed method comparing to related examples (see, J. Ren, S. S. Venkata, and E. Sortomme, “An Accurate Synchrophasor Based Fault Location Method for Emerging Distribution Systems,” IEEE Trans. Power Deliv., vol. 29, no. 1, pp. 297-298, February 2014; Z. Mengsheng, W. Yi, Z. Zhou, and Z. Li, “Research on fault location based on PMU for multi-source distribution network,” in Asia-Pacific Power and Energy Engineering Conference, APPEEC, 2016, vol. 2016-Decem, pp. 1877-1882). In the related examples mentioned above, a 99% of accuracy is disclosed, which is a highest accuracy for impedance-based fault location techniques in distribution networks by using PMUs among the related examples.
In the present disclosure, both the asymmetrical and symmetrical methods result in high-level of accuracy with overall average error of 0.62% and 0.067%, respectively. It is observed that the line length decreases as the maximum and average errors drops for both the symmetrical and asymmetrical methods. The two IMMs (i.e., GIMM and LIMM) can be developed for each line (DL) by design consultant firms during the design phase of the lines, and part of the protection studies. In addition, by calculating the line parameters online (in situ), the disclosed method does not suffer from the line impedance inaccuracy issues under the idea situation (homogeneous, radial, and uncompensated, etc.).
