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.
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.
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
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
Further preferred embodiments are evident from the dependent patent claims.
The subject matter of the invention will be explained in more detail in the following text with reference to the attached drawing,
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.
The phasor data 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
According to the invention, the temperature of the line is determined in the following way: The electric line parameters, 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
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.
The average line temperature T1 is computed from the ohmic resistance R by modelling 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 resistance 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 the line, differ from the average line temperature T1, in a preferred variant of the invention 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 a further preferred variant of the Invention, a limit to the power delivered through the line is computed from the phasor information v1, v2, i1, i2. This Is done by solving the well-known load flow equation
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 a preferred embodiment of the invention, 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.
Number | Date | Country | Kind |
---|---|---|---|
01811254.0 | Dec 2001 | EP | regional |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/CH02/00682 | 12/11/2002 | WO | 5/3/2005 |