METHOD AND DEVICE FOR ESTIMATING SOURCE IMPEDANCES ACROSS ONE OR MORE TRANSMISSION LINES

Abstract

The present invention relates to estimating source impedances across one or more transmission lines connecting at least two substations. Each substation is associated with an equivalent source having a respective source impedance. Voltage and current measurements and status signals are obtained. The voltage and current measurements provide terminal or bus voltages and line currents at each terminal, and the status signals are associated with switching events at the one or more transmission lines or at the substations. An event associated with a disturbance or current injection is detected from one or more of the obtained measurements and the obtained status signals. The source impedance of each equivalent source is estimated based on the event, using line parameters and the voltage and current measurements associated with the event.

Claims

1. A method for estimating source impedances across one or more transmission lines connecting at least two substations, each substation associated with an equivalent source having a respective source impedance, the method comprising: obtaining voltage and current measurements and status signals, the voltage and current measurements providing terminal voltages and line currents at each substation, and the status signals being associated with switching events at the one or more transmission lines or at the substations; detecting an event associated with a disturbance in a transmission line or current injection at one of the substations, the event being detected from one or more of the obtained measurements and the obtained status signals; and estimating the source impedance of each equivalent source based on the event, using line parameters and the voltage and current measurements associated with the event, the measurements associated with the event comprising the measurements before and after the event.

2. The method of claim 1, wherein the line parameters comprise one or more of characteristic impedance and propagation constant of a transmission line.

3. The method of claim 2, wherein the line parameters further comprise line length of the transmission line.

4. The method of claim 1, wherein the source impedance being further estimated using a fault current or a current injected at a terminal of the substation.

5. The method of claim 1, wherein the disturbance is a fault on one of the one or more transmission lines, and wherein estimating the source impedance of each equivalent source comprises: estimating a fault location on the transmission line; calculating a fault current using the estimated fault location; calculating two bus impedances, each bus impedance being associated with a respective terminal of the one or more transmission lines, and calculated using a change in the terminal voltages of the corresponding terminal and the calculated fault current; calculating one or more of line and section impedances using the line parameters and the estimated fault location, wherein the section impedances are impedances of line sections located between the terminals and the estimated fault location; and estimating the source impedance of each equivalent source using one or more of the calculated line impedance and the calculated section impedances, and the two bus impedances.

6. The method of claim 5, wherein calculating each bus impedance comprises multiplying the bus impedance with a correction factor, wherein the correction factor is calculated using a ratio of a fault voltage and the calculated fault current, and shunt admittances of the line sections, wherein the fault voltage is calculated using the estimated fault location, and the shunt admittances are calculated using the estimated fault location and the line parameters, optionally wherein estimating the source impedance comprises removing contribution of shunt admittance, wherein the contribution of shunt admittance is calculated using one or more of the shunt admittance of the at least one transmission line, and the shunt admittances of the line sections calculated using on the fault location and the line parameters.

7. The method of claim 1, wherein the estimated source impedance is a complex source impedance having both magnitude and angle.

8. The method of claim 1, wherein the estimated source impedance is used in configuring one of a protection function and a control function in a substation device.

9. A computer program product for estimating source impedances across one or more transmission lines, the computer program product comprising a non-transitory computer readable storage medium comprising instructions that, when executed by a processor, cause the processor to execute the method of claim 1.

10. A device for estimating source impedances across one or more transmission lines connecting at least two substations, each substation associated with an equivalent source having a respective source impedance, the device comprising: a data acquisition unit configured to obtain voltage and current measurements and status signals, the voltage and current measurements providing terminal voltages and line currents at each substation, and the status signals being associated with switching events at the one or more transmission lines or at the substations; a signal processing unit communicatively coupled to the data acquisition unit and configured to: detect an event associated with a disturbance in one of the one or more transmission lines or current injection at one of the substations, the signal processing unit configured to detect the event from one or more of the obtained measurements and the obtained status signals; and estimate the source impedance of the equivalent source based on the event, using line parameters and the voltage and current measurements associated with the event, the measurements associated with the event comprising the measurements before and after the event; and a control unit communicatively coupled to the signal processing unit.

11. The device of claim 10, wherein the line parameters comprise one or more of characteristic impedance and propagation constant of a transmission line.

12. The device of claim 11, wherein the line parameters further comprise line length of the transmission line.

13. The device of claim 10, wherein the source impedance being further estimated using a fault current or a current injected at a terminal of the substation.

14. The device of claim 10, wherein the disturbance is a fault on one of the one or more transmission lines, and wherein the signal processing unit is configured to: estimate a fault location on the transmission line; calculate a fault current using the estimated fault location; calculate two bus impedances, each bus impedance being associated with a respective terminal of the one or more transmission lines, and calculated using a change in the terminal voltages of the corresponding terminal and the calculated fault current; calculate one or more of line and section impedances using the line parameters and the estimated fault location, wherein the section impedances are impedances of line sections located between the terminals and the estimated fault location; and estimate the source impedance of each equivalent source using one or more of the calculated line impedance and the calculated section impedances, and the two bus impedances.

15. The device of claim 14, wherein the signal processing unit is configured to calculate each bus impedance by multiplying the bus impedance with a correction factor, wherein the correction factor is calculated using a ratio of a fault voltage and the calculated fault current, and shunt admittances of the line sections, wherein the fault voltage is calculated using the estimated fault location, and the shunt admittances are calculated using the estimated fault location and the line parameters, optionally wherein the signal processing unit is configured to estimate the source impedance by removing contribution of shunt admittance, wherein the contribution of shunt admittance is calculated using one or more of the shunt admittance of the at least one transmission line, and the shunt admittances of the line sections calculated using on the fault location and the line parameters.

16. The device of claim 10, wherein the device is a substation server, optionally wherein the substation server is communicatively coupled with one or more measurement equipment and the switching devices; or wherein the device is an intelligent electronic device associated with a terminal of the one or more transmission lines, and communicatively coupled with another intelligent electronic device associated with another terminal of the one or more transmission line.

17. The device of claim 11, wherein the device is a substation server, optionally wherein the substation server is communicatively coupled with one or more measurement equipment and the switching devices; or wherein the device is an intelligent electronic device associated with a terminal of the one or more transmission lines, and communicatively coupled with another intelligent electronic device associated with another terminal of the one or more transmission line.

18. The device of claim 12, wherein the device is a substation server, optionally wherein the substation server is communicatively coupled with one or more measurement equipment and the switching devices; or wherein the device is an intelligent electronic device associated with a terminal of the one or more transmission lines, and communicatively coupled with another intelligent electronic device associated with another terminal of the one or more transmission line.

19. The device of claim 13, wherein the device is a substation server, optionally wherein the substation server is communicatively coupled with one or more measurement equipment and the switching devices; or wherein the device is an intelligent electronic device associated with a terminal of the one or more transmission lines, and communicatively coupled with another intelligent electronic device associated with another terminal of the one or more transmission line.

20. The device of claim 14, wherein the device is a substation server, optionally wherein the substation server is communicatively coupled with one or more measurement equipment and the switching devices; or wherein the device is an intelligent electronic device associated with a terminal of the one or more transmission lines, and communicatively coupled with another intelligent electronic device associated with another terminal of the one or more transmission line.

Description

DRAWINGS

[0056] FIG. 1 is a schematic of a power transmission system with a double-circuit transmission line.

[0057] FIG. 2 shows an equivalent model of the power transmission system of FIG. 1.

[0058] FIG. 3 is a flowchart of a method for estimating source impedances, in accordance with various embodiments of the invention.

[0059] FIG. 4 shows an equivalent two-port model with transmission lines represented by a pi model, in accordance with an embodiment of the invention.

[0060] FIG. 5 is a flowchart of the method for estimating source impedances, in accordance with an embodiment of the invention.

[0061] FIGS. 6 and 7 are example graphical results of bus voltages and line currents following a fault simulated at a mid-point of a line.

[0062] FIG. 8 shows an equivalent two-port model with transmission lines represented by a pi model, in accordance with another embodiment of the invention.

[0063] FIG. 9 is a flowchart of the method for estimating source impedances, in accordance with another embodiment of the invention.

[0064] FIGS. 10 and 11 are example graphical results of bus voltages and line currents following a simulated shunt switching event.

[0065] FIG. 12 is a simplified representation of a power transmission system where a server is provided for estimating source impedances, in accordance with an embodiment of the invention.

[0066] FIG. 13 is a simplified representation of a power transmission system where an Intelligent Electronic Device is provided for estimating source impedances, in accordance with an embodiment of the invention.

[0067] FIG. 14 is a block diagram of a device for estimating source impedances, in accordance with an embodiment of the invention.

[0068] FIG. 15 is a simplified representation of a two-source equivalent source across a transmission line, in accordance with an exemplary scenario.

DETAILED DESCRIPTION

[0069] Various analyses in power transmission lines (e.g. line protection) depend on an equivalent model of a power transmission system. Consider a power transmission system shown in FIG. 1. This power transmission system has a double-circuit transmission line, i.e. circuit-1 and circuit-2 as shown, connecting terminals (or buses) M and N. Further, as shown, the transmission line can be connected to the rest of the network (e.g. through other transmission lines). Also, depending on the operational requirement, there can be shunt elements such as a reactor or capacitor connected at the buses as shown in FIG. 1. The load element may represent the downstream lower voltage level network and/or a direct High Tension (HT) load tapping.

[0070] An equivalent model of the power transmission system of FIG. 1 is shown in FIG. 2, in which the components of the system except the line of concern are all collapsed into two equivalent sources behind corresponding impedances. In FIG. 2, the equivalent sources are shown as Equivalent Source M, and Equivalent Source N. In this embodiment, the equivalent source is a combination of network lines (of FIG. 1). However, the equivalent source can be a combination of other power system equipment or sub-systems. For example, the equivalent source can represent grouping of one or more of a transmission line, power generator, transformer, reactor, and/or load, or other such combinations as would be apparent to those of ordinary skill in the art.

[0071] Typically, an equivalent model similar to the one shown in FIG. 2 is used in various power system devices, for example for setting the operational characteristics of distance relays, for setting the power swing blinders and out of step logic, amongst other system/device configurations etc.

[0072] Source impedances (e.g. impedance of equivalent source M or N), which are the components of a typical equivalent model, are used to define the Source Impedance Ratio (SIR) of the line. This can for instance identify whether it is a short, medium or long line for setting under/overreaching elements and designing pilot protection schemes.

[0073] This invention deals with estimation of impedances Z.sub.SM and Z.sub.SN (see FIG. 2), which are known as the source impedances (i.e. impedances of the equivalent sources). The straightforward approach to solve this problem is to use the system Z.sub.BUS to derive the impedances. However, as a substation level exercise, information of the entire network topology may not always be available. Therefore, in this invention our approach is to find an estimate of the impedances only with the substation level measurements and transmission line data. They are the bus voltage and line current measurements at both ends, currents on the lines incident on the terminal buses, currents through the connected shunt elements and model parameters of the line itself.

[0074] It should be noted that the method and device are described in relation to an equivalent model as shown in FIG. 2. However, it may be possible to extend the invention to other transmission line configurations. For

[0075] The invention utilizes measurements during events such as, but not limited to, line faults, bus faults, short circuits on the line or switching events on the terminal buses for the estimation process. The invention utilizes digital substation measurements at both ends of a transmission line following an event to estimate the source impedances. The measured changes in bus voltages due to the extra current injection (or fault currents) are used to calculate the source impedances of the equivalent sources.

[0076] Referring to FIG. 3, which is a flowchart of a method for estimating source impedances across one or more transmission lines (e.g. for equivalent sources in a power transmission system such as shown in FIG. 1), in accordance with various embodiments of the invention.

[0077] At 302, the method comprises obtaining voltage and current measurements and status signals. The voltages and currents can be obtained with measurement equipment such as voltage or current transformers, Rogowski coils, or other suitable voltage or current sensors.

[0078] The voltage and current measurements provide terminal (e.g. substation terminal or busbar (bus)) voltages and line currents at each substation. In other words, voltage and current measurements at each terminal of the one or more transmission lines are obtained through the corresponding measurement equipment. The line currents include currents incident on, or currents circulated from a terminal. The line currents can also include currents from the shunt elements (if provided at the terminal). The status signals are associated with switching events at the one or more transmission lines or at the substations. These events can include, but need not be limited to, switching (on/off) of a line breaker(s), a breaker at a terminal, a shunt element(s) and so forth, where such switches can be provided on transmission lines or at the substations (e.g. connected at terminals). The status signals (or switching status) can be obtained directly from the switching devices, or via other devices such as an Intelligent Electronic Devices (IEDs) or relays provided for monitoring and/controlling the corresponding switching devices.

[0079] At 304, the method comprises detecting an event associated with a disturbance in one of the one or more transmission lines or associated with current injection at one of the substations (e.g. at a substation terminal or bus). In accordance with various embodiments, the event is detected from one or more of the obtained measurements and the obtained status signals.

[0080] A fault on the line may be detected based on a status output of a protection IED on the line. For example, if a trip output is “1”, then it can be concluded that a fault has occurred on the line, where the fault may have been detected from the obtained measurements. A fault event may also be detected based on status of a line breaker, including bus faults. If the fault is cleared by opening of the line, then the breaker status changes. A flag set by an operator (e.g. substation operator) indicating occurrence of a fault can also be used for detecting a fault.

[0081] Similarly, detection of a switching event of a shunt element can be done based on device status or human flags. For example, if a shunt element is added/removed at the line terminals, then the corresponding breaker status changes. Taking another example, a flag set by an operator can indicate switching (on/off) of a shunt element.

[0082] At 306, the method comprises estimating the source impedance of each equivalent source based on the event. In accordance with different embodiments, the estimated source impedance has a complex value, i.e. has both magnitude and angle. Here, the source impedance is estimated using line parameters and the voltage and current measurements associated with the event, where the measurements associated with the event comprise the measurements before and after the event. The line parameters of a transmission line can comprise one or more of, but not limited to, characteristic impedance, propagation constant and line length of the transmission.

[0083] The following describes embodiments of the method wherein an embodiment, the event is a line fault, and in another embodiment, the event is shunt injection at a substation terminal.

[0084] Let us consider the network shown in FIG. 4, which depicts an equivalent two-port model with the transmission lines represented by a pi-model. Node ‘F’ (Bus F) is a fictitious node which represents the point of occurrence of fault on the line. Note that a double-circuit line is used for the ensuing analysis and discussion as it is very common in High Voltage (HV) transmission. However, the analysis is applicable for a single circuit line also as will be indicated later.

Embodiment: Estimating Source Impedances Using Measurements Related with Line Faults

[0085] Let us consider the network of FIG. 4. A fault at Bus F can be considered as injection of fault current I.sub.flt at the same bus. Due to this extra current injection in the network, the voltages at terminal buses M and N undergo changes. Here, terminals M and N can be seen as substation terminals or busbars (i.e. buses). Let us define the following impedances for each bus as shown in (1) below.

[00001]$\begin{array}{cc}{Z}_{\mathrm{MF}}\text{:=}\frac{V_M^0-V_M^f}{I_{\mathrm{flt}}},Z_{\mathrm{NF}}\text{:=}\frac{V_N^0-V_N^f}{I_{\mathrm{flt}}}& \left(1\right)\end{array}$

[0086] The source impedance of each bus (or bus impedance) is a ratio of the change in the corresponding bus voltage and the fault current. Symbols V and I represent voltage and current phasors respectively, which can be obtained for example by applying a Discrete Fourier Transform (DFT) on the measurements obtained from the measurement equipment (e.g. current and voltage transformers) at the substation. A ‘0’ in superscript denotes steady state and ‘f’ denotes faulted state quantity. Changes in the bus voltages due to the fault current injection can also be expressed in terms of the network parameters as shown in (2) below.

[00002]$\begin{array}{cc}\begin{array}{c}{V}_M^0-V_M^f=I_{\mathrm{flt}}{Z}_{\mathrm{sM}}\times \frac{\frac{Z_{l2}}{\left(\frac{Z_{\mathrm{sum}}}{Z_l}\right)}+Z_{\mathrm{sN}}}{Z_{\mathrm{par}}+\left(Z_{\mathrm{sM}}+Z_{\mathrm{sN}}\right)}\\ V_N^0-V_N^f=I_{\mathrm{flt}}{Z}_{\mathrm{sN}}\times \frac{\frac{Z_{l1}}{\left(\frac{Z_{\mathrm{sum}}}{Z_l}\right)}+Z_{\mathrm{sM}}}{Z_{\mathrm{par}}+\left(Z_{\mathrm{sM}}+Z_{\mathrm{sN}}\right)}\end{array}\}& \left(2\right)\end{array}$

[0087] From (2), Z.sub.MF and Z.sub.NF can be calculated as shown in (3) below.

[00003]$\begin{array}{cc}\begin{array}{c}{Z}_{\mathrm{MF}}=Z_{\mathrm{sM}}\times \frac{\frac{Z_{l2}}{\left(\frac{Z_{\mathrm{sum}}}{Z_l}\right)}+Z_{\mathrm{sN}}}{Z_{\mathrm{par}}+\left(Z_{\mathrm{sM}}+Z_{\mathrm{sN}}\right)}\\ Z_{\mathrm{NF}}=Z_{\mathrm{sN}}\times \frac{\frac{Z_{l1}}{\left(\frac{Z_{\mathrm{sum}}}{Z_l}\right)}+Z_{\mathrm{sM}}}{Z_{\mathrm{par}}+\left(Z_{\mathrm{sM}}+Z_{\mathrm{sN}}\right)}\end{array}\}\mathrm{In}\left(3\right)Z_{\mathrm{sum}}=Z_l+Z_{l1}+Z_{l2}\mathrm{and}{Z}_{\mathrm{par}}=\frac{Z_l\left(Z_{l1}+Z_{l2}\right)}{Z_{\mathrm{sum}}}.& \left(3\right)\end{array}$

[0088] In the above relationships of (2), currents in the shunt admittance branches of the line are neglected. The equations when solved to find Z.sub.SM and Z.sub.SN (i.e. source impedances of the equivalent sources M and N respectively) yield the following result as shown in (4) below.

[00004]$\begin{array}{cc}\begin{array}{c}{Z}_{\mathrm{sM}}=\frac{-b+\sqrt{b^2}-4\mathrm{ac}}{2a}\\ Z_{\mathrm{sN}}=\frac{\left(Z_{\mathrm{par}}+Z_{\mathrm{sM}}\right)Z_{\mathrm{MF}}-\frac{Z_{l2}}{\left(\frac{Z_{\mathrm{sum}}}{Z_l}\right)}{Z}_{\mathrm{sM}}}{Z_{\mathrm{sM}}-Z_{\mathrm{MF}}}\end{array}\}& \left(4\right)\end{array}$

[00005]$a=Z_{\mathrm{MF}}-Z_{\mathrm{NF}}-\frac{Z_{l2}}{\left(Z_{\mathrm{sum}}\text{/}{Z}_l\right)},b=Z_{\mathrm{par}}{Z}_{\mathrm{MF}}+\frac{Z_{l1}}{\left(Z_{\mathrm{sum}}\text{/}{Z}_l\right)}\left(Z_{\mathrm{MF}}-Z_{\mathrm{NF}}\right)-\frac{Z_{l1}{Z}_{l2}}{{\left(Z_{\mathrm{sum}}\text{/}{Z}_l\right)}^2},\mathrm{and}c=\frac{Z_{l1}}{\left(Z_{\mathrm{sum}}\text{/}{Z}_l\right)}{Z}_{\mathrm{par}}{Z}_{\mathrm{MF}}.$

[0089] In the above relationships of (4),

[0090] Assuming an equivalent pi circuit model of the transmission line,

Z.sub.l=Z.sub.c sinh(γl),Z.sub.l1=Z.sub.c sinh(γdl),Z.sub.l2=Z.sub.c sinh(γ(1−d)l),

where Z.sub.c is the characteristic impedance and γ is the propagation constant of the line. As evident, Z.sub.l1 and Z.sub.l2 (i.e. impedances of the line sections between the terminals and the fault location) and all related quantities are dependent on the fault location ‘d’, which is a pre-requisite for this embodiment of the method. As an example, this can be found by using the simple principle that the fault point voltage calculated using measurements from either end of the line should be the same. The following can be used for calculating ‘d’.

[00006]$\begin{array}{cc}d=\frac{1}{\gamma l}{\tanh }^{-1}\left(\frac{K_1}{K_2}\right)\mathrm{In}\left(5\right),K_1=V_N^f\cosh \left(\gamma l\right)+\left(I_{N1}^f{Z}_c\right)\sinh \left(\gamma l\right)-V_M^f,K_2=V_N^f\sinh \left(\gamma l\right)-\left(I_{N1}^f{Z}_c\right)\cosh \left(\gamma l\right)-Z_c{I}_{M1}^f.& \left(5\right)\end{array}$

[0091] Quantities Z.sub.MF and Z.sub.NF (i.e. bus impedances) as defined in (1) are calculated by using the steady state and post fault voltage phasors from both end buses and the fault current phasor. Fault current I.sub.flt (considered positive away from the bus ‘F’) can be calculated once the fault location is obtained. To do this, we can first obtain the fault current contribution from each terminal bus M and N as shown below in (6 (6a, 6b)).

[00007]$\begin{array}{cc}{I}_{\mathrm{flt}}^M=-\frac{1}{Z_c}\sinh \left(\gamma \mathrm{dl}\right)V_M^f+\cosh \left(\gamma \mathrm{dl}\right)I_{M1}^f& \left(6a\right)\\ I_{\mathrm{flt}}^N=-\frac{1}{Z_c}\sinh \left(\gamma \left(1-d\right)l\right)V_N^f+\cosh \left(\gamma \left(1-d\right)l\right)I_{N1}^f& \left(6b\right)\end{array}$

[0092] Then the total fault current is obtained by adding up the contributions from bus M side and bus N side as shown in (7).

I.sub.flt=I.sub.flt.sup.M+I.sub.flt.sup.N  (7)

[0093] It should be noted that the fault location and the fault current phasor can be calculated using other methods, and the method is not limited to the calculations shown above.

[0094] The steps of the above embodiment of the method are illustrated in the flowchart of FIG. 5. As mentioned, the embodiment can be performed after a fault is detected on the line.

[0095] At 502, the fault location ‘d’ is estimated by using two ended voltage and current measurements and line parameters as in (5).

[0096] At 504, the fault current I.sub.flt is calculated by using the fault location d and the two ended voltage and current measurements as in (7).

[0097] At 506, the bus impedances Z.sub.MF and Z.sub.NF are calculated as ratio of the change in respective bus voltage and the fault current as in (1).

[0098] At 508, the line and/or section impedances (impedances Z.sub.l, Z.sub.l1, Z.sub.l2) are calculated using the line parameters and the fault location d.

[0099] At 510, the source impedances of the equivalent sources (Z.sub.SM and Z.sub.SN) are estimated using (4), i.e. using the bus impedances and the line and/or section impedances.

[0100] In equations (1) and (2), the currents in the shunt admittance branches were ignored. In order to take them into account, the following corrections can be done if shunt admittance branches are present in the system.

[0101] Correction 1 (C1):

[0102] To Z.sub.MF and Z.sub.NF as obtained in step 506 above, multiply

[00008]$\xi =\frac{1}{1-Z_{\mathrm{FF}}{Y}_{\mathrm{shd}}}$

as a correction factor where,

[00009]$Z_{\mathrm{FF}}=\frac{V_F^0-V_F^f}{I_{\mathrm{flt}}},\mathrm{and}$$Y_{\mathrm{shd}}=Y_{\mathrm{shl}1}+Y_{\mathrm{shl}2}=\frac{1}{Z_c}\tanh \left(\frac{\gamma \mathrm{dl}}{2}\right)+\frac{1}{Z_c}\tanh \left(\frac{\gamma \left(d-1\right)l}{2}\right).$

[0103] Note that once the fault location d is known, the fault point voltage V.sub.F can be found by using V.sub.M and I.sub.M1 measurements in the pi-equivalent model equations for the line section MF of the faulted circuit.

[0104] Correction 2 (C2):

[0105] From the obtained values of Z.sub.SM and Z.sub.SN as in step 5 above, remove the contribution from the line shunt admittances as:

[00010]$Z_{sM}^{corr}=\frac{1}{\left(\frac{1}{Z_{sM}}-Y_{shl}-Y_{shl1}\right)},Z_{sN}^{corr}=\frac{1}{\left(\frac{1}{Z_{sN}}-Y_{shl}-Y_{shl2}\right)},$

[0106] In the above,

[00011]$Y_{shl}=\frac{1}{Z_c}\tanh \left(\frac{\gamma l}{2}\right),Y_{\mathrm{shl}1}=\frac{1}{Z_c}\tanh \left(\frac{\gamma \mathrm{dl}}{2}\right)\mathrm{and}{Y}_{shl2}=\frac{1}{Z_c}\tanh \left(\frac{\gamma \left(1-d\right)l}{2}\right).$

[0107] A test system of the type shown in FIG. 2 is considered for the numerical illustrations. It is simulated in PSCAD environment. The lines (120 km each) are represented with frequency dependent phase models in the simulation. The characteristic impedance of the line Z.sub.c is 244.16∠−3.4745Ω and the propagation constant γ is 1.0806×10.sup.−06∠86.461°. The line impedance in each circuit is approximately equal to 31.661∠82.922°. The equivalent sources are modelled each as 50 Hz, 220 kV voltage source behind an impedance, with a phase angle difference of 5.5°. The values of source impedances, i.e., Z.sub.SM and Z.sub.SN are kept as 63.32∠84.81° and 158.3∠79° for the purpose of this illustration. A fault is created at the mid-point of Circuit 1 of the transmission line at 0.2 s. Following the fault, the bus voltages and line currents undergo changes as shown in FIGS. 6 and 7.

[0108] For the calculations, we have used only the voltage and current measurements in the window spanning 30 ms to 50 ms (milliseconds) from the time of fault inception. This is done assuming the fault is cleared by operating the circuit breaker at least after two and a half cycles from the fault inception time; after which we will not be having any post fault line current measurements of the faulted circuit. The result of each step of the analysis are presented below. [0109] 1. Fault location estimated using (5): d=0.5002 [0110] 2. Fault current calculated using obtained value of d in (6) and (7): I.sub.flt=2.2835∠14.936° kA [0111] 3. Z.sub.MF calculated as

[00012]$\frac{V_M^0-V_M^f}{I_{flt}}$  and correction factor C1 applied is found to be=46.071∠82.594°Ω, similarly, Z.sub.NF=49.373∠82.193°Ω [0112] 4. Using known values of Z.sub.c, γ and d, we obtain Z.sub.l=31.573∠82.941°Ω, Z.sub.l1=15.823∠82.926°Ω, Z.sub.l2=15.817∠82.926°Ω

[0113] 5. Using calculated values of Z.sub.MF, Z.sub.NF, Z.sub.l, Z.sub.l1, Z.sub.l2 in (4) and then applying correction C2 we obtain Z.sub.sM.sup.corr=63.341∠84.6Ω and Z.sub.sN.sup.corr=158.5∠78.436Ω, which amounts to (total vector) errors of 0.4% and 1% respectively.

Embodiment: Estimating Source Impedances Using Measurements Related with Shunt Injections

[0114] For this embodiment, the equivalent two-port model with transmission lines represented by pi model as shown in FIG. 8 is considered. Let us consider current injection due to a shunt element (e.g. reactor/capacitor/load) switching event on a terminal bus, say bus M. As done for the above embodiment, let us define the following impedances for each bus M and N (i.e. bus impedances).

[00013]$\begin{array}{cc}{Z}_{MF}:=\frac{V_M^0-V_M^f}{I_{inj}},Z_{NF}:=\frac{V_N^0-V_N^f}{I_{inj}}& \left(8\right)\end{array}$

[0115] Here, I.sub.inj represents the shunt current injection at the bus (considered positive when measured in the direction into the bus). Changes in the bus voltages due to the current injection can also be expressed in terms of the network parameters as follows:

[00014]$\begin{array}{cc}{V}_M^0-V_M^f=I_{\mathrm{inj}}{Z}_{\mathrm{sM}}\times \frac{\left(Z_i/2\right)+Z_{\mathrm{sN}}}{\left(Z_i/2\right)+\left(Z_{\mathrm{sM}}+Z_{\mathrm{sN}}\right)},V_N^0-V_N^f=I_{\mathrm{inj}}{Z}_{\mathrm{sN}}\times \frac{Z_{\mathrm{sM}}}{\left(Z_l/2\right)+\left(Z_{\mathrm{sM}}+Z_{\mathrm{sN}}\right)}& \left(9\right)\end{array}$

[0116] From (9), Z.sub.MF and Z.sub.NF can be calculated as shown in (10) below.

[00015]$\begin{array}{cc}{Z}_{\mathrm{MF}}=Z_{\mathrm{sM}}\times \frac{\left(Z_l/2\right)+Z_{\mathrm{sN}}}{\left(Z_l/2\right)+\left(Z_{\mathrm{sM}}+Z_{\mathrm{sN}}\right)}{Z}_{\mathrm{NF}}=Z_{\mathrm{sN}}\times \frac{Z_{\mathrm{sM}}}{\left(Z_l/2\right)+\left(Z_{\mathrm{sM}}+Z_{\mathrm{sN}}\right)}& \left(10\right)\end{array}$

[0117] The equation (10) when solved to find Z.sub.SM and Z.sub.SN (i.e. source impedances of the equivalent sources) yields the result as shown in (11), which can be further simplified to the expressions shown in (12).

[00016]$\begin{array}{cc}{Z}_{sM}=\frac{\left(\frac{Z_l}{2}\right)Z_{MF}}{\frac{Z_l}{2}+Z_{NF}-Z_{MF}},Z_{sN}=\frac{\left(Z_{sM}+\frac{Z_l}{2}\right)Z_{NF}}{Z_{sM}-Z_{NF}}& \left(11\right)\\ Z_{sM}=\frac{Z_{MF}}{1+\frac{I_M^O-I_M^I}{I_{inj}}},Z_{sN}=\frac{Z_{NF}}{\frac{I_N^O-I_N^I}{I_{inj}}}& \left(12\right)\end{array}$

[0118] In (12), I.sub.M=I.sub.M1+I.sub.M2, I.sub.N=I.sub.N1+I.sub.N2, ‘0’ in superscript denotes steady state quantity and ‘I’ denotes value of the entity post the switching event. Note that the expression in (9) can be found by keeping d=0 in equation (2). For shunt injection at bus N, equations similar to (9)-(12) can be derived by keeping d=1 in equations (2).

[0119] The steps of the above embodiment of the method are illustrated in the flowchart of FIG. 9. As mentioned, the embodiment can be performed after a switching event is detected at a terminal bus.

[0120] At 902, the injected current I.sub.inj is obtained directly from the substation measurements.

[0121] At 904, the bus impedances (Z.sub.MF and Z.sub.NF) are calculated as ratio of the change in respective bus voltage and the injected current as in equation (8).

[0122] At 906, the line impedance Z.sub.l is calculated using the line parameters.

[0123] At 908, the source impedances (Z.sub.SM and Z.sub.SN) of the equivalent sources are estimated using the obtained Z.sub.MF and Z.sub.NF in (11). Alternately, the source impedances (Z.sub.SM and Z.sub.SN) are estimated by using Z.sub.MF, Z.sub.NF and measured change in line currents at both the ends as in (12).

[0124] The effect of the line shunt admittances may optionally be considered in the above embodiment. For considering their effect, the following correction needs to be applied.

[0125] Correction:

[0126] From the obtained values of Z.sub.SM and Z.sub.SN as in step 908 above, remove the contribution from the line shunt admittances as

[00017]$Z_{sM}^{corr}=\frac{1}{\left(\frac{1}{Z_{sM}}-2Y_{shl}\right)},Z_{sN}^{corr}=\frac{1}{\left(\frac{1}{Z_{sN}}-2Y_{shl}\right)},\mathrm{where}{Y}_{shl}=\frac{1}{Z_c}\tanh \left(\frac{\gamma l}{2}\right).$

[0127] A test system similar to that for the embodiment of line faults is considered for the numerical illustrations of the above embodiment related to shunt injections. We create a shunt reactor switching on Bus M at 0.2 s and use the obtained measurements to estimate the source impedances. The reactor which is switched onto the bus is rated at 100 MVar at 220 kV nominal line-to-line bus voltage. Following the switching, bus voltages and line currents undergo change as shown in FIGS. 10 and 11. The result of each step of the analysis is presented below. [0128] 1. The injected current due to the reactor switching is measured as I.sub.inj=0.35434∠−91.192° kA [0129] 2. Z.sub.MF calculated as

[00018]$\frac{V_M^O-V_M^I}{-I_{\mathrm{inj}}}=48.699\mathrm{\angle 83}\mathrm{.005}°\Omega ,$  similarly,

[00019]$Z_{NF}=\frac{V_N^O-V_N^I}{-I_{\mathrm{inj}}}=44.653\angle 82.735°\Omega$ [0130] 3. Measured change in line currents I.sub.M.sup.o−I.sub.M.sup.f and I.sub.N.sup.o−I.sub.N.sup.f are obtained as 0.082745∠94.675° and 0.099926∠−87.6° kA respectively. Using them with the values of I.sub.inj, Z.sub.MF and Z.sub.NF in (12) and then applying correction C3 we obtain Z.sub.sM.sup.corr=61.344∠84.956°Ω and Z.sub.sN.sup.corr=146.23∠79.983°Ω, amounting to (total vector) error of 1.82% and 1.75% respectively.

[0131] In case of bus faults, a similar approach can be employed, and in place of the current injected due to the switching event, the fault current at the bus can be obtained from the obtained measurements.

[0132] The method described above is not restricted to double-circuit lines and can easily be extended to single circuit lines. For a single circuit transmission line, the following adjustments need to be made. [0133] (a) Line faults: Since there is only one circuit, we can consider impedance of the second circuit as Z.sub.l=∞. This makes the terms

[00020]$\frac{Z_{sum}}{Z_l}=1$  and Z.sub.par=Z.sub.l1+Z.sub.l2. Substitute these values in equation (4) to obtain Z.sub.sM and Z.sub.sN. [0134] (b) Bus injection: Use equation (12) to obtain the source impedances with the only change that I.sub.M=I.sub.M1 and I.sub.N=I.sub.N1 as the second circuit is not present.

[0135] The method can be implemented in a computer program product. Accordingly, in an embodiment, the invention provides a computer program product for estimating source impedances across one or more transmission lines. The computer program product can be, but need not be limited to, an optical disc (e.g. a DVD, CD, or other disc type computer storage product), or a hard-drive or other memory device capable of storing computer program or machine-readable instructions. As an example, the computer program product can be a memory with instructions, that can be part of a server or a substation device.

[0136] The computer program product comprises a non-transitory computer readable storage medium. Examples of such non-transitory computer readable storage medium can be, but need not be limited to, a random access memory (RAM), a read only memory (ROM), an erasable programmable read only memory (EPROM), an external memory device (e.g. USB memory, flash drive etc.) or other similar types of storage medium. The non-transitory computer readable storage medium comprises instructions that, when executed by a processor, cause the processor to execute the method (i.e. part of or the entire method as per the case) as described above. These instructions are program or machine-readable instructions defining the method (i.e. instructions defining the method as described above). These instructions can be executed with a processor which can be one of or a combination of a computer processing unit, a multiprocessor, a microcontroller, and a digital signal processor. When such instructions are executed, one or more steps of the method are performed.

[0137] The method described above can be implemented with a power system or substation device. For example, as shown in FIG. 12, the method can be implemented with a server. This server may be a substation server, or a remotely located server, which receives the obtained voltage and current measurements (e.g. from the measurement equipment, IEDs or relays) and the status signals (e.g. from the switching devices or relays). Taking another example, as shown in FIG. 13, the method can be implemented with an Intelligent Electronic Device (IED), which has local end (terminal) voltage and current measurements and is communicatively coupled with an IED at the remote end to receive measurements/data corresponding to the other end. In this example, the IED also has the status signals.

[0138] FIG. 14 is a block diagram of such a device (e.g. server or IED) for estimating the source impedances in the power transmission system. As shown, the device comprises a data acquisition unit (1402), a signal processing unit (1404), a control unit (1406), a communication unit (1408) and a memory unit (1410) communicatively coupled with each other (e.g. via a communications bus (1412)).

[0139] The various modules of the device are implemented with hardware devices such as a processor, I/O's etc. The processor is configured to invoke each module (individually or simultaneously) by executing respective instructions to perform one or more of the various steps of the method described herein above. Such instructions may be stored in the memory unit. Each module is configured to perform one or more steps of the method. For example, the data acquisition unit can perform step 302, while the signal processing unit can perform steps 304 and 306.

[0140] In accordance with various embodiments, the data acquisition unit is configured to obtain the voltage and current measurements (e.g. as shown by 1414) and the status signals, from the measurement equipment (or IEDs or relays) and the switching devices (or associated IEDs or relays). Further, the signal processing unit is communicatively coupled to the data acquisition unit and configured to detect an event associated with a disturbance in one of the one or more transmission lines, or current injection at one of the terminals.

[0141] The signal processing unit is configured to detect the event from one or more of the obtained measurements and the obtained status signals. Alternately, the event may be detected by another module (e.g. an event detection module (not shown)), and the same can be fed as input to the data acquisition unit and/or the signal processing unit.

[0142] The signal processing unit is also configured to estimate the source impedance of each equivalent source based on the event, using line parameters and the voltage and current measurements associated with the event, the measurements associated with the event comprising the measurements before and after the event.

[0143] The device also comprises a control unit communicatively coupled to the signal processing unit. There may be a need to adapt an equivalent model due to changes in the source impedances. In such situations, the control unit is configured to adapt (e.g. shown by 1416) an equivalent model to perform at least one of a protection function and a control function based on the estimated source impedances. Such adaptations can be pushed to the relevant devices (e.g. relays, when the device is a server) via the communication unit. Alternately, the memory can be updated according to the adaptations. In some cases, as needed, the updated source impedances can be communicated to another device for use at the corresponding devices.

[0144] In an embodiment, the various modules of the device are configured to perform the steps in relation to estimating the source impedances according to the embodiment of FIG. 5. In another embodiment, the various modules of the device are configured to perform the steps in relation to estimating the source impedances according to the embodiment of FIG. 9. In these embodiments, the signal processing unit can also be configured to perform the corrections for shunt admittance branches or line shunt admittances.

[0145] Thus, the invention utilizes digital substation data during system events for simultaneous estimation of the source impedances across one or more transmission lines. Events such as shunt injections at the buses, bus faults and faults on the line can be utilized for the estimation. Further, the invention provides complex source impedances instead of only providing source impedance magnitudes. The invention utilizes only fault location information as input for the source impedance estimation and does not require information of the external network topology. In other words, it is computationally simple and non-iterative in implementation.

[0146] The invention can help in achieving the following objects. Firstly, an online estimate of the source impedances can facilitate adaptive setting of relays, which can thereby improve its dependability and security. Secondly, a method which depends only substation level measurements can be of immediate use to a substation engineer. Thirdly, a method which can utilize events such as faults on the line and switching of shunt elements at the terminal buses is of more practical importance as compared to staging bus faults, which may not always be possible or acceptable.

[0147] The invention can be used to improve the following protection functions: [0148] It improves the distance relay performance and enables adaptive setting features. [0149] It improves accuracy of single ended fault location algorithms. [0150] It can be used to calculate the system non-homogeneity factor and distance relay reach calculation. [0151] Source impedance information can be used for adaptive switching the different phase selection methods (phasor or super imposed quantities or traveling wave based) to protect the systems with high renewables. [0152] It can be useful in calculation of short circuit ratio (SCR) and it will be used for over current relay co-ordination.

[0153] The above benefits can be understood for example from the description provided below. Consider the two-source equivalent model shown in FIG. 15. It can be shown by circuit analysis that:

[00021]$\begin{array}{cc}{V}_M=I_{M}d{Z}_{line}+\left(\frac{I_{fM}}{D_M}\right)R_f\mathrm{where},D_M=\frac{\left(1-d\right)Z_{\mathrm{line}}+Z_{sN}}{Z_{sM}+Z_{\mathrm{line}}+Z_{sN}}& \left(\mathrm{E1}\right)\end{array}$

represents the current distribution factor for the fault current as seen from the Bus M side. V.sub.M and I.sub.M are the voltage and current measured at the IED during the fault and I.sub.fM is the fault current component in the current I.sub.M or in other words the fault current contribution from the local (with respect to the IED location) Bus M. Equation (E1) can be re-written as:

[00022]$\begin{array}{cc}{Z}_M=\frac{V_M}{I_M}=d{Z}_{line}+\left(\frac{I_{fM}}{I_M}\right)\frac{R_f}{D_M}& \left(\mathrm{E2}\right)\end{array}$

which expresses the impedance measured at the IED. Clearly, estimates of the source impedances can help in analyzing the relay reach for various fault related parameters such as the local/remote fault current contribution, fault location and fault resistance. Based on these analyses, resistive and reactive reach settings of the relay operating characteristics can be adapted.

[0154] The relay impedance trajectory during a power swing is approximated by the expression below:

[00023]$\begin{array}{cc}{Z}_M\left(\delta \right)=\frac{Z_{sM}+Z_{line}+Z_{sN}}{2}\left(1-j\cot \frac{\delta }{2}\right)-Z_{SM}& \left(\mathrm{E3}\right)\end{array}$

in which, the angle δ represents the phase advance of the equivalent source E.sub.M over E.sub.N. The expected swing impedance trajectory during an actual power swing can be analyzed using the model of equation (E3), by plugging in the values of the estimated source impedances and varying angle δ from initial loading angle to 180°. Using this estimated trajectory, power swing blocking parameters such as blinder positions with respect to distance zone elements and delay timer setting can be adaptively determined.

[0155] Equation (E2) also forms the basic principle for the single ended fault location algorithm in which V.sub.M, I.sub.M and I.sub.fM are the measurable inputs and ‘d’ is the unknown to be estimated along with the unknown fault resistance R.sub.f. Source impedance estimates are vital for such algorithms which are used to characterize the term D.sub.M.

[0156] The system non-homogeneity factor represents the variation among the angles associated with the three impedances Z.sub.sM, Z.sub.line and Z.sub.sN. It is a direct result of the source impedance estimates and it can be useful in adapting certain analyses. For example, if the system is homogeneous, i.e., all the three impedances have approximately the same angle, then only the imaginary part of equation (E2) can be solved to directly obtain ‘d’ without the need of estimating R. This is because in such a situation, the imaginary part of the second term on the right-hand side of (E2) tends towards zero.

[0157] The estimated source impedances can be used to decide between a current based or a voltage-based phase selection strategy implemented at the IED bus M. For example, a weak source behind the IED bus (i.e., high source impedance Z.sub.SM) will typically result in lesser incremental (or fault) current contribution from the bus M. In such cases, it is better to switch to an incremental voltage-based phase selection strategy.

[0158] Estimated source impedances can also be used to infer the short circuit MVA at each of the terminal buses M and N as expressed in equations below:

[00024]$\mathrm{SC}\mathrm{MVA}\mathrm{at}\mathrm{Bus}M\left(\mathrm{in}\mathrm{pu}\right)=\frac{Z_{sM}+Z_{line}+Z_{sN}}{Z_{sN}\left(Z_{line}+Z_{sM}\right)}$$\mathrm{SC}\mathrm{MVA}\mathrm{at}\mathrm{Bus}N\left(\mathrm{in}\mathrm{pu}\right)=\frac{Z_{sM}+Z_{line}+Z_{sN}}{Z_{sM}\left(Z_{line}+Z_{sN}\right)}$

[0159] Such information is useful for monitoring of fault current levels and adaptive coordination of overcurrent relay characteristics.