Determining an operational limit of a power transmission line

Information

  • Patent Grant
  • 7107162
  • Patent Number
    7,107,162
  • Date Filed
    Wednesday, December 11, 2002
    21 years ago
  • Date Issued
    Tuesday, September 12, 2006
    18 years ago
Abstract
In a method, computer program and system for determining an operational limit of a power transmission line, time-stamped current phasor information and voltage phasor information for a first end and a second end of the line are determined, an ohmic resistance of the line is computed from the phasor information, and an average line temperature is computed from the ohmic resistance. This allows to determine the average line temperature without dedicated temperature sensors. The average line temperature represents the actual average temperature and is largely independent of assumptions regarding line parameters.
Description

This disclosure is based upon Europe Application No. 01811254.0, filed Dec. 21, 2001, and International Application No. PCT/CH02/00682, filed Dec. 11, 2002, the contents of which are incorporated herein by reference.


FIELD OF THE INVENTION

The invention relates to large-scale electric power transmission networks, and, more particularly, to a method, a computer program and a system for determining an operational limit of a power transmission line according to the preamble of the independent claims.


BACKGROUND OF THE INVENTION

Electric power transmission and distribution systems or networks comprise high-voltage tie lines for connecting geographically separated regions, medium-voltage lines, and substations for transforming voltages and for switching connections between lines. Power generation and load flow in a network with several substations is controlled by an energy management system. For managing the network, it is desirable to determine a state of the network, in particular load flows and stability margins. Measurements of average RMS (root mean square) values of voltages, currents, active power and reactive power flowing in the network and/or measurements of voltage and current phasors are made at various places in the network and collected by substation automation (SA) systems and by a supervisory control and data acquisition (SCADA) system. In recent times, devices and systems for measuring voltage and current phasors at different locations of a network at exactly the same time and for collecting them at a central data processor have become available. Phasors collected from throughout the network in combination provide a snapshot of the electrical state of the network.


In order to maximally utilise a line, it is desirable to maximise the amount of transmitted power. This amount is limited by several factors, in particular a thermal limit and an electrical limit. Currently, these limits are taken into account by making assumptions about the electrical parameters of the line and by calculating a maximum power flow. A temperature of the line, which depends on weather conditions, may be measured by temperature sensors distributed along the line. Such measurements are however not representative for the entire line. Due to the shortcomings of existing systems, power lines are operated with overly conservative limits on transmitted power.


DESCRIPTION OF THE INVENTION

It is therefore an object of the invention to create a method, a computer program and a system for determining an operational limit of a power transmission line of the type mentioned initially, which allows to determine the operational limit with greater precision and to operate the line closer to an actual operational limit.


These objects are achieved by a method, a computer program and a system for determining an operational limit of a power transmission line according to the independent claims.


The inventive method for determining an operational limit of a power transmission line comprises the steps of

    • determining time-stamped current phasor information and voltage phasor information for a first end and a second end of the line,
    • computing an ohmic resistance of the line from the phasor information,
    • computing an average line temperature from the ohmic resistance.


This allows to determine the average line temperature without dedicated temperature sensors. The average line temperature represents the actual average temperature and is largely independent of assumptions regarding line parameters.


An additional advantage is conferred by the fact that phasor measurements incorporate high resolution timing information. It is therefore also possible to obtain a much better temporal resolution for the temperature estimation than with measurements available through conventional SCADA systems. This better temporal resolution allows faster and better quality control of the line and the network.


The computer program for determining an operational limit of a power transmission line according to the invention is loadable into an internal memory of a digital computer, and comprises computer program code means to make, when said computer program code means is loaded in the computer, the computer execute the method according to the invention. In a preferred embodiment of the invention, a computer program product comprises a computer readable medium, having the computer program code means recorded thereon.


The system for determining an operational limit of a power transmission line according to the invention comprises a data processing device configured to

    • determine time-stamped current phasor information and voltage phasor information for a first end and a second end of the line,
    • compute an ohmic resistance R of the line from the phasor information,
    • compute an average line temperature T1 from the ohmic resistance R.


Further preferred embodiments are evident from the dependent patent claims.





BRIEF DESCRIPTION OF THE DRAWING

The subject matter of the invention will be explained in more detail in the following text with reference to the attached drawings.



FIG. 1 illustrates an equivalent circuit for a transmission line in accordance with an exemplary embodiment of the present invention; and



FIG. 2 is a flowchart illustrating a process for determining an operational limit of a power transmission line in accordance with an exemplary embodiment of the present invention.





DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS


FIG. 1 shows an equivalent circuit for a transmission line. The circuit is a standard π-equivalent circuit for a transmission line or line segment connecting a first node n1 of the network with a second node n2. A complex variable representing a voltage phasor v1 corresponds to the voltage at the first node n1, i.e. at a first end of the line, and a complex variable representing a current phasor i1 corresponds to the current flowing into the line from the first end. In analogy, a voltage phasor v2 and current phasor i2 are associated with the second end of the line. The line is represented by electrical line parameters, that is, a line impedance R+jXL and shunt admittances jXC. A power flow s into the second node n2 comprises a real part p and imaginary part jq.


Phasor data is determined with a phasor measurement unit (PMU) residing, for example, at a feeder at the bay level of substations or at branching points along transmission lines. A voltage phasor represents, for example, a voltage of the feeder or line, while a current phasor represents current flowing through the feeder or line. A phasor may also be used to represent electric power.


The phasor data represents a phasor and may be a polar number, the absolute value of which corresponds to either the real magnitude or the RMS value of a quantity, and the phase argument to the phase angle at zero time. Alternatively, the phasor may be a complex number having real and imaginary parts or the phasor may use rectangular or exponential notation. By contrast, conventional sensing devices used in power networks generally measure only scalar, average representations, such as the RMS value of a voltage, current etc.



FIG. 2 is a flowchart illustrating a process for determining an operational limit of a power transmission line in accordance with an exemplary embodiment of the present invention. In Step 10, time stamped current phasor information and voltage phasor information for the first end and a second end of the transmission line is collected from phasor measurement units that are distributed over a large geographic area, i.e. over tens to hundreds of kilometres. Since the phasor data from these disparate sources are analysed in conjunction, they must refer to a common phase reference. Therefore, the different phasor measurement units must have local clocks that are synchronised with each other to within a given precision. Such a synchronisation of the phasor measurement units is preferably achieved with a known time distribution system, for example the global positioning (GPS) system. In a typical implementation, the phasor data 9 is determined at least every 200 or every 100 or preferably every 40 milliseconds, with a temporal resolution of preferably less than 1 millisecond. In a preferred embodiment of the invention, the temporal resolution is less than 10 microseconds, which corresponds to a phase error of 0.2 degrees. Each measurement is associated with a time stamp derived from the synchronised local clock. The phasor data therefore comprises time stamp data.


In another preferred embodiment of the invention, phasor information such as a voltage phasor and a current phasor associated with a node or line is not measured at said node or line, but is inferred from phasor measurements made at a location that is remote from said node or line. For example, with regard to FIG. 1, if the electrical line parameters are known and v1 and i1 at the first node n1 are measured with a PMU, then v2 and i2 at the second node n2 can be computed. In the context of the present invention, this makes sense only if the line whose electrical parameters are assumed to be known and are used to compute phasors at a location that is remote from the PMU is not identical to the line whose electrical parameters and/or temperature need to be estimated.


In step 12, the electric line parameters, R, XL, and XC, or at least the ohmic resistance R of the line, i.e. the real part R of the line impedance Z=R+jXL are determined from measured or computed phasor information representing some or all of the voltage and current phasors at the two ends of line.


In a first preferred variant of the invention, it is assumed that the shunt capacitance jXC remains essentially constant during power line operation and is known from other measurements, design parameters or calculations. This is a valid assumption, since changes in shunt capacitance jXC are relatively small when compared to the ohmic resistance R. Then it is necessary to determine only the two voltage phasors v1 and v2 at either end of the line and one of the current phasors i1 or i2. Let i1 be measured. Then the impedance Z is






Z
=



v
1

-

v
2




i
1

-



v
1

·
j







X
C









and the desired ohmic resistance R is the real part of Z.


In a second preferred variant of the invention, no assumption on shunt impedances is made, and the two voltage phasors v1 and v2 and the two current phasors i1 or i2 are measured or determined from measurements. Determining the actual electrical line parameters R, XL, XC from these measurements is common knowledge. Since resulting equations for the electrical line parameters are non-linear, numerical methods such as Newton-Raphson approximation are used for determining actual parameter values. The resulting line parameters are actual values in that they are determined online and represent the actual state of the power line, in contrast to average values that are assumed to be constant over all seasons and ambient conditions.


In step 14, the average line temperature T1 is computed from the ohmic resistance R by modeling a relationship between temperature and resistance as linear, i.e.

R=R2+α(T1−T2)


or as quadratic, i.e.

R=R2+α(T1−T2)+β(T1−T2)2


where R2 and T2 are known reference restance and temperature values dependent on the construction of the line and α and β are material constants for the line cables. The linear relationship is typical for common conductor materials such as copper or aluminium. As an example, the parameter values are such that for a line temperature change of ΔT=30° C. the resistance changes by about ΔR=12%. The equation for the chosen relationship is solved for T1, which gives the desired average line temperature.


Since the temperature of the line may, according to location along with the line, differ from the average line temperature T1, in Step 16 the power flow through the line is regulated such that a predetermined maximum average line temperature Tmax is not exceeded. Such a regulation or control regime is implemented with any commonly known control method such as PID control, non-linear control, model predictive control, . . . etc.


In Step 18, the maximum load Smax is determined based on a limit to the power delivered through the line which is computed from the phasor information v1, v2, i1, i2. This is done by solving the well-known load flow equation







p
-

j





q


=


v
2
*

·

(




v
1

-

v
2



R
+

j






X
L




-


v
2


j






X
C




)







where v2* is the complex conjugate of v2: Assuming that v1 is constant and given by a power source or generator at the first node n1, and that consumed power s corresponding to a load at the second node n2 is varied, then there are either two, one or no solutions of the load flow equation for v2. For an increasing load, the power at which there is no solution anymore corresponds to a maximum load smax. This is the maximum amount of power that can be delivered by the line to the second node n2 before the line becomes unstable and the voltage v2 collapses. Since the ohmic resistance R plays an important role in the load flow equation, the maximum load smax depends on actual line conditions, in particular on the line temperature. Determining the maximum load smax according to the invention gives an actual maximum value that allows a less conservative approach than when a maximum load is given and remains constant over a wide range of operating conditions.


Depending on impedance characteristics of the load, power flow is maximised by increasing power flow until the maximum load smax is reached. In Step 20, the power is controlled to increase only until a given safe distance from the maximum load smax is reached. In both cases, according to operator preferences, the load and maximum load are either considered as complex variables, or only the real parts are considered.


If the load has a fixed p/q relationship, then active power p and reactive power q delivered to the load are increased together. If reactive power q of the load is constant, then the active power p is increased alone. As an example, a change of the line resistance ΔR of 10% leads to a loadability change of Δsmax=6.5% for a typical 400 kV line. In the conventional offline worst-case analysis the lower limit has to be considered because R is seen as fixed. With the online approach the loadability can be increased by up to the calculated value, dependent on actual ambient conditions.


In a system according to the invention, the system comprises means for determining phasor information that are configured to receive time-stamped measured phasor data from at least two PMUs located throughout the network, and optionally means for computing phasor information for at least one node from measured phasor data corresponding to other nodes. The system further comprises means for computing an ohmic resistance R of the transmission line from the measured and/or computed phasor information and for computing an average line temperature T1 from the ohmic resistance R.

Claims
  • 1. A method for determining an operational limit of a power transmission line, comprising the steps of determining time-stamped current phasor information and voltage phasor information for a first end and a second end of the line,computing an ohmic resistance R of the line from the phasor information,computing an average line temperature T1 from the ohmic resistance R.
  • 2. The method according to claim 1, comprising the step of controlling a power flow through the line such that the average line temperature T1 does not exceed a given maximum average line temperature Tmax.
  • 3. The method according to claim 1, comprising the steps of computing actual electrical parameters (R, XL, XC) of the line and computing, from the actual electrical line parameters, a maximum load smax that is deliverable by the line.
  • 4. The method according to claim 3, comprising the step of controlling a power flow through the line such that an actual power s delivered by the line does not exceed the maximum load smax.
  • 5. The method according to claim 4, wherein the step of controlling a power flow through the line ensures that the average line temperature T1 does not exceed a given maximum average line temperature Tmax.
  • 6. A computer program for determining an operational limit of a power transmission line loadable into an internal memory of a digital computer, comprising computer program code means to make, when said program is loaded in the computer, the computer execute the method according to claim 1.
  • 7. A system for determining an operational limit of a power transmission line from time-stamped current phasor information and voltage phasor for a first end and a second end of the line, the system comprising: a data processing device and means forcomputing an ohmic resistance R of the line from the phasor information, andcomputing an average line temperature T1 from the ohmic resistance R.
  • 8. The system according to claim 7, further comprising: means for controlling power flow through the line such that the average line temperature T1 does not exceed a given maximum average line temperature Tmax.
  • 9. The system according to claim 7, further comprising: means for computing actual electrical parameters (R, XL, XC) of the line; andmeans for computing, from the actual electrical line parameters, a maximum load smax that is deliverable by the line.
Priority Claims (1)
Number Date Country Kind
01811254 Dec 2001 EP regional
PCT Information
Filing Document Filing Date Country Kind 371c Date
PCT/CH02/00682 12/11/2002 WO 00 5/3/2005
Publishing Document Publishing Date Country Kind
WO03/055028 7/3/2003 WO A
US Referenced Citations (34)
Number Name Date Kind
4413325 Elfner et al. Nov 1983 A
5181026 Granville Jan 1993 A
6236949 Hart May 2001 B1
6313752 Corrigan et al. Nov 2001 B1
6476521 Lof et al. Nov 2002 B1
6545482 Fedirchuk et al. Apr 2003 B1
6546353 Hallett et al. Apr 2003 B1
6598003 Heino et al. Jul 2003 B1
6618684 Beroset et al. Sep 2003 B1
6662124 Schweitzer et al. Dec 2003 B1
6675100 Hallett et al. Jan 2004 B1
6694270 Hart Feb 2004 B1
6708133 Hallett et al. Mar 2004 B1
6711521 Hallett et al. Mar 2004 B1
6845333 Anderson et al. Jan 2005 B1
6940702 Kojovic et al. Sep 2005 B1
20010021896 Bertsch et al. Sep 2001 A1
20020116139 Przydatek et al. Aug 2002 A1
20030014678 Ozcetin et al. Jan 2003 A1
20030101008 Hart May 2003 A1
20030105608 Hart Jun 2003 A1
20030200038 Schweitzer et al. Oct 2003 A1
20030212512 Hart Nov 2003 A1
20030220752 Hart Nov 2003 A1
20040008461 Kojovic et al. Jan 2004 A1
20040010350 Lof et al. Jan 2004 A1
20040027748 Kojovic et al. Feb 2004 A1
20040059469 Hart Mar 2004 A1
20040133367 Hart Jul 2004 A1
20040172207 Hancock et al. Sep 2004 A1
20040186670 Hart Sep 2004 A1
20050083206 Couch et al. Apr 2005 A1
20050141682 Wells Jun 2005 A1
20050251296 Nelson et al. Nov 2005 A1
Foreign Referenced Citations (2)
Number Date Country
0 795 944 Sep 1997 EP
09 200949 Jul 1997 JP
Related Publications (1)
Number Date Country
20050222808 A1 Oct 2005 US