Time domain distance protection based on polarity comparison for power transmission lines

Abstract

A method can be used for fault detection of a transmission. The method includes obtaining instantaneous measurements of local voltages and local currents of the transmission line, obtaining filtered measurements of the local voltages and local currents by subjecting the measurements of the local voltages and local currents to low pass filtering, obtaining a compensated voltage calculated by subjecting the filtered measurements to a differential equation based algorithm, and performing fault detection by forming a ratio between the compensated voltage and the a reference voltage. An internal fault is determined when the ratio is below a threshold and an external fault is determined when ratio is larger than above the threshold.

Claims

1. A method for fault detection of a transmission line, the method comprising: obtaining instantaneous measurements of local voltages and local currents of the transmission line; obtaining filtered measurements of the local voltages and local currents by subjecting the measurements of the local voltages and local currents to low pass filtering; obtaining a compensated voltage calculated by subjecting the filtered measurements to a differential equation based algorithm; and performing fault detection by forming a ratio between the compensated voltage and a reference voltage, wherein an internal fault is determined when the ratio is below a threshold and an external fault is determined when ratio is above the threshold.

2. The method according to claim 1, wherein the method is performed for each of three voltage phases.

3. The method according to claim 1, wherein the method is performed for both phase-to-phase measurements and phase-to-ground measurements of the local voltages and currents.

4. The method according to claim 1, further comprising providing a result of the fault detection to a distance protection arrangement for making a trip decision.

5. The method according to claim 4, wherein the result causes a trip of a circuit breaker of the distance protection arrangement to be triggered only if the detected fault is an internal fault.

6. The method according to claim 1, wherein the low pass filter used for the low pass filtering has a cutoff frequency between 1 and 5 kHz.

7. The method according to claim 1, wherein the transmission line is represented by a resistor-inductor (RL) circuit and the differential equation based algorithm is an RL based differential equation based algorithm.

8. The method according to claim 1, wherein the transmission line is represented by a resistor-inductor-capacitor in series for series capacitor compensated transmission lines (RLC) circuit and the differential equation based algorithm is an RLC based differential equation based algorithm.

9. The method according to claim 1, wherein the ratio is determined as a normalized convolution between the compensated voltage and the reference voltage.

10. The method according to claim 1, wherein the ratio is formed from measurements collected during a time window.

11. The method according to claim 1, wherein the threshold is set to zero.

12. The method according to claim 1, wherein the threshold is set to −0.2.

13. The method according to claim 1, wherein the transmission line is an ultra high voltage (UHV) transmission line.

14. A method of operating a transmission line, the method comprising: measuring local voltages and local currents of the transmission line to obtain instantaneous measurements of local voltages and local currents; low pass filtering the measurements of the local voltages and local currents to obtain filtered measurements of the local voltages and local currents; subjecting the filtered measurements to a differential equation based algorithm to obtain compensated voltage values; performing fault detection by forming a ratio between the compensated voltage and a reference voltage, wherein an internal fault is determined when the ratio is below a threshold and an external fault is determined when the ratio is greater than the threshold; and providing a result of the fault detection to a distance protection arrangement for making a trip decision.

15. The method according to claim 14, wherein the result causes a trip of a circuit breaker of the distance protection arrangement to be triggered only if the detected fault is an internal fault.

16. A fault detection apparatus for a transmission line, the fault detection apparatus comprising: processing circuitry; a memory coupled to the processing circuitry, the memory storing a computer program comprising computer code that, when run on the processing circuitry, causes the fault detection apparatus to: obtain instantaneous measurements of local voltages and local currents of the transmission line; obtain filtered measurements of the local voltages and currents by subjecting the measurements of the local voltages and local currents to low pass filtering; obtain a compensated voltage measurement by subjecting the filtered measurements to a differential equation based algorithm; and perform fault detection by forming a ratio between the compensated voltage and a reference voltage, wherein an internal fault is determined when the ratio is below a threshold and an external fault is determined when the ratio is above the threshold.

17. The apparatus according to claim 16, wherein the instantaneous measurements of the local voltages and local currents are performed for each of three voltage phases.

18. The apparatus according to claim 16, wherein the instantaneous measurements of the local voltages and local currents are performed for both phase-to-phase measurements and phase-to-ground measurements of the local voltages and currents.

19. The apparatus according to claim 16, wherein the computer code that, when run on the processing circuitry, causes the fault detection apparatus to provide a result of the fault detection to a distance protection arrangement for making a trip decision.

20. The apparatus according to claim 16, wherein the transmission line is an ultra high voltage (UHV) transmission line.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The inventive concept is now described, by way of example, with reference to the accompanying drawings, in which:

(2) FIGS. 1, 2, 3 are schematic diagrams illustrating transmission systems according to embodiments;

(3) FIG. 4 is a flowchart of methods according to embodiments;

(4) FIGS. 5, 6, 7, 8, 10, 11, 12 illustrate measured electrical signal values and calculated values according to embodiments;

(5) FIG. 13 is a schematic diagram showing functional units of a fault detection arrangement according to an embodiment;

(6) FIG. 14 is a schematic diagram showing functional modules of a fault detection arrangement according to an embodiment; and

(7) FIG. 15 shows one example of a computer program product comprising computer readable storage medium according to an embodiment.

DETAILED DESCRIPTION OF ILLUSTRATIVE EMBODIMENTS

(8) The inventive concept will now be described more fully hereinafter with reference to the accompanying drawings, in which certain embodiments of the inventive concept are shown. This inventive concept may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein; rather, these embodiments are provided by way of example so that this disclosure will be thorough and complete, and will fully convey the scope of the inventive concept to those skilled in the art. Like numbers refer to like elements throughout the description. Any step or feature illustrated by dashed lines should be regarded as optional.

(9) FIG. 1 illustrates an example 500 kV transmission line system. The line length is for illustrative purposes assumed to be 230 km long. The sampling rate is 1 kHz for both current and voltage measurements. The reach setting of the proposed distance protection is 80% of the line (i.e., 184 km). The integration window is 20 ms. Only when more than 3 sample points in a row enter the operate zone, the proposed distance protection will release the trip signal.

(10) In more detail, FIG. 1 schematically illustrates a transmission system 25 of a power distribution system where the herein disclosed embodiments apply. The transmission system 25 comprises at least one fault detection arrangement 10a, 10b for fault detection of a transmission line 20 of the transmission system 25. Two or more fault detection arrangements 10a, 10b may be operatively connected via a communications link 23 for pilot scheme, if needed. In theory, the communication link is not needed for distance protection including the proposed one. But in some cases, two or more distance protections can be connected with communication link or pilot scheme for better performance.

(11) Further, two or more fault detection arrangements 10a, 10b may be part of a common fault detection arrangement 10c for fault detection of the transmission line 20. The fault detection arrangement 10a, 10b, 10c may be part of, or comprise, an intelligent electronic device (IED) operating as a relay. The transmission system 25 further comprises power sources 21a, 21b, current and voltage transformers 22a, 22b, and circuit breakers 23a, 23b. F1 and F2 denote internal and external faults, respective, along the transmission line 20. The transmission line 20 might be an ultrahigh voltage (UHV) transmission line 20. The transmission line 20 might be part of a power distribution system.

(12) FIG. 2 schematically illustrates a simplified example of the transmission system 25 of FIG. 1 with an internal fault occurring at point F1. FIG. 3 schematically illustrates a simplified example of the transmission system 25 of FIG. 1 with an external fault occurring at point F2.

(13) The embodiments disclosed herein relate to mechanisms for fault detection of the transmission line 20 of the transmission system 25. In order to obtain such mechanisms, there is provided a fault detection arrangement 10a, 10b, 10c, methods performed by the fault detection arrangement 10a, 10b, 10c, a computer program product comprising code, for example in the form of a computer program, that when run on a fault detection arrangement 10a, 10b, 10c, causes the fault detection arrangement 10a, 10b, 10c to perform the method.

(14) Time domain protection based on the polarity information of full component compensated voltage at the reach point is proposed in order to correctly identify internal faults and external faults.

(15) According to the proposed time domain protection scheme, an internal fault is detected by comparing the polarities of compensated voltage at the reach point with a reference voltage (e.g., a local voltage). As shown in FIGS. 2 and 3, the polarities are the opposite for internal faults (F1 fault as in FIG. 2) and the same for external faults (F2 fault as in FIG. 3).

(16) In theory, the principle shown in FIGS. 2 and 3 is valid for every sampling point in the time domain. Furthermore, the principle is valid no matter if the source voltage is a pure fundamental sine wave or any other wave shapes. This is helpful for the operational speed and security of the protection scheme.

(17) It is acknowledged that there might be some measurement errors, parameter errors, calculation errors, etc. or other types of errors in certain practical applications. The herein disclosed methods are designed to consider such errors, making the methods robust and reliable.

(18) Methods for fault detection of a transmission line will now be disclosed with reference to the flowchart of FIG. 4. The methods are performed by a fault detection arrangement 10a, 10b, 10c.

(19) S102: Instantaneous measurements of the local voltages and local currents are obtained. For the voltages this corresponds to measurements of U.sub.A in FIGS. 2 and 3.

(20) S104: Filtered measurements are obtained by subjecting the measurements of the local voltages and local currents to low pass filtering. The cut off frequency of the low pass filter used for the low pass filtering depends on the consideration of different factors, such as the parameters of current transforms (CTs) and capacitor voltage transformers (CVTs), measurement sampling rate of the protection relay, etc. For the voltages this corresponds to measurements of V.sub.ref below. For example, the cutoff frequency can be on the order of several kHz, e.g., 1 to 5 kHz.

(21) S106: Compensated voltage measurements are obtained by subjecting the filtered measurements of the local voltages and local currents to a differential equation based algorithm. For the voltages this corresponds to measurements of U.sub.comp in FIGS. 2 and 3.

(22) In some aspects the transmission line is assumed to be represented by a resistor—inductor circuit (RL circuit) and the differential equation based algorithm is thus an RL based differential equation based algorithm. In other aspects the transmission line is assumed to be represented by a resistor—inductor—capacitor circuit (RLC circuit) in case of series compensation capacitors are included in the transmission line and the differential equation based algorithm is thus an RLC based differential equation based algorithm.

(23) S108: Fault detection is performed by forming a ratio between the compensated voltage measurements and the filtered voltage measurements.

(24) In some aspects the ratio is formed from measurements collected during a time window. For examples, the time window could be 10 ms or 20 ms long.

(25) In some aspects the ratio is determined as a normalized convolution between the compensated measurements and the filtered measurements.

(26) An internal fault is determined if the ratio is below a negative threshold. This corresponds to the situation in FIG. 2. Else, an external fault is determined. This corresponds to the situation in FIG. 3. In some aspects the threshold is set to zero. In other aspects the threshold is set to a negative value close to zero. For examples, the threshold could be set to −0.2.

(27) The above steps are performed for each of three phases.

(28) The above steps are performed for both phase-to-phase measurements and phase-to-ground measurements of the local voltage.

(29) The above steps of obtaining instantaneous measurements and obtaining filtered measurements could also be performed for instantaneous measurements of the local current in each of the three phases and for both phase-to-phase measurements and phase-to-ground measurements of the local current. The compensated measurements (of the local voltage) could then be dependent also on the filtered measurements of the local current in each of the three phases and for both phase-to-phase measurements and phase-to-ground measurements.

(30) S110: In some aspects a result of the fault detection (i.e., if an internal fault was detected or an external fault was detected) is provided to a distance protection arrangement for making a trip decision. A trip is then only triggered if the detected fault is an internal fault. That is, the result causes a trip of a circuit breaker 23a, 23b of the distance protection arrangement to be triggered only if the detected fault is an internal fault.

(31) Further details of fault detection of a transmission line 20 will now be disclosed.

(32) As noted above, three actions, might be performed during the fault detection once the instantaneous measurements have been obtained (as in S102): application of low pass filtering (as in S104), running of a differential equation based algorithm (as in S106), and performing voltage comparison (as in S108). Aspects of S104-S108 will now be disclosed in further detail.

(33) First, high frequency noise is removed by low pass filtering of the measurements of the local measurements (voltages as well as currents). In this respect, low pass filtering can be used to decrease the influence of high frequency noise (since the differential equation based algorithm as used thereafter might be sensitive to high frequency noise or error).

(34) The output filtered measurements (voltages as well as currents) of the low pass filtering is used as input to the differential equation based algorithm whereby compensated voltage measurements are obtained.

(35) The compensated voltage in the time domain at the reach point is calculated by means of a differential equation based algorithm. This calculation is valid for any frequency component or any wave shape (in theory). Two examples of a differential equation based algorithm have been mentioned above for a transmission line with typical RL circuits represented as series line parameters and for series compensated transmission line with series capacitors included in the series line parameters as RLC circuits.

(36) Assuming the symbol φ is used to represent any phase A, B, C of a three phase transmission system and the symbol φφ is used to represent any phase to phase of three phase transmission systems, the following definitions could be given:

(37) Uq.sub.φφ is the compensated voltage for the phase-to-phase loop

(38) Uq.sub.φ is the compensated voltage for the phase-to-ground loop

(39) u.sub.φφ is the local phase to phase loop voltage

(40) u.sub.φ is the local phase voltage

(41) i.sub.φφ is the local phase to phase current

(42) i.sub.φ is the local phase current

(43) i.sub.0 is the local zero sequence current

(44) R.sub.0 is the line zero sequence resistance

(45) R.sub.1 is the line positive sequence resistance

(46) L.sub.0 is the zero sequence inductance

(47) L.sub.1 is the line positive sequence inductance

(48) u.sub.SCφφ is the phase to phase voltage across the series compensation capacitor circuits, and

(49) u.sub.SCφ is the one phase voltage across the series compensated capacitor circuit.

(50) For the normal transmission line with series line parameters as RL circuits, the compensated voltage of phase-to-phase loop Uq.sub.φφ and the compensated voltage of phase-ground loop Uq.sub.φ are determined as follows respectively:

(51) $U{q}_{\mathrm{\phi \phi }}=u_{\mathrm{\phi \phi }}-i_{\mathrm{\phi \phi }}.Math.R_1-L_1\frac{{\mathrm{di}}_{\mathrm{\phi \phi }}}{\mathrm{dt}},\mathrm{and}$$U{q}_{\phi }=u_{\phi }-\left(\left(i_{\phi }+k_{R}3i_0\right).Math.R_1+L_1\frac{d}{\mathrm{dt}}\left(i_{\phi }+k_{L}3i_0\right)\right),\mathrm{where}$$k_R=\frac{R_0-R_1}{3R_1},k_L=\frac{L_0-L_1}{3L_1}=\frac{X_0-X_1}{3X_1}.$

(52) For the series compensated line with series line parameters as RLC circuits, the compensated voltage of phase-to-phase loop Uq.sub.φφ and the compensated voltage of phase-ground loop Uq.sub.φ are determined as follows:

(53) $U{q}_{\mathrm{\phi \phi }}=u_{\mathrm{\phi \phi }}-i_{\mathrm{\phi \phi }}.Math.R_1-L_1\frac{{\mathrm{di}}_{\mathrm{\phi \phi }}}{\mathrm{dt}}-u_{\mathrm{SC}\phi \phi },\mathrm{and}$$U{q}_{\phi }=u_{\phi }-\left(\left(i_{\phi }+k_{R}3i_0\right).Math.R_1+L_1\frac{d}{\mathrm{dt}}\left(i_{\phi }+k_{L}3i_0\right)\right)-u_{\mathrm{SC}\phi },\mathrm{where}$$k_R=\frac{R_0-R_1}{3R_1},k_L=\frac{L_0-L_1}{3L_1}=\frac{X_0-X_1}{3X_1}.$

(54) The output of the differential equation based algorithm is used as input to the fault detection. As disclosed above, in some aspects a comparison of voltage polarities is performed in terms of a normalized convolution. The normalized convolution between two general voltage values v.sub.1(t) and v.sub.2(t) can be expressed as follows:

(55) $s=\frac{\int v_1\left(t\right).Math.v_2\left(t\right)\mathrm{dt}}{\int .Math.v_1\left(t\right).Math.v_2\left(t\right).Math.\mathrm{dt}}$

(56) Here, v.sub.1(t) and v.sub.2(t) are replaced by values of the compensated measurements and local measurement with the filtered measurements for each phase A, B, C. For example, for one phase loop, v.sub.1(t) could be the local voltage u.sub.φ(t) and v.sub.2(t) could be Uq.sub.φ(t); for one of phase to phase loop, v.sub.1(t) could be the local voltage u.sub.φφ(t) and v.sub.2(t) could be Uq.sub.φφ(t).

(57) Although the polarity criterion for internal fault detection in theory is valid for every sampling value, in practice this condition might not be ideal because of measurement errors, calculation errors, etc. Therefore, a normalized convolution might be used to compare the voltage polarities of compensated voltage and reference voltage (e.g., local voltage) for increasing reliability and robustness. The integration window length might depend on system conditions. A short integration window (e.g., 10 ms) could be used, which is beneficial to the sensitivity and speed. In case more serious fault transients (e.g., compensated line) or oscillations occur, a longer integration window (e.g., 20 ms) could be used, which is beneficial to the security and reliability.

(58) In theory, the ratio should be positive for an external fault and negative for an internal fault. In theory, the value zero could therefore be used as the threshold in the criterion for detecting an internal fault. But considering the reliability and errors that might occur in practical applications, a negative setting value (e.g., −0.2) could be used as the threshold to detect the internal fault. Only when the calculated ratio ‘s’ is below the threshold (e.g., −0.2), an internal fault will be detected.

(59) To further enhance the security and reliability, a startup element could be used to open the trip window for some time (e.g., 20 ms). Otherwise, if the startup element does not operate, the final trip signal will be blocked even when the method indicates an internal fault.

(60) The proposed method is based on polarity comparison of voltages, which is not sensitive to voltage amplitude. As a result, it can handle high source impedance ratio (SIR) conditions, while existing time domain protection methods based on amplitude may have weak sensitivity under high SIR conditions.

(61) The proposed method is suitable for both normal power lines and series compensated lines.

(62) The proposed method can work continuously because it is based on full components instead of fault component values.

(63) Simulation results will be disclosed next.

(64) Simulation results for an internal fault will now be disclosed.

(65) An internal phase-to-ground fault (F1 in FIG. 1) for phase A is assumed to occur at 100 ms and be located at the midpoint of the line (i.e., 115 km). The power load angle is 30 degrees and source impedance to line impedance ratio (SIR) is set to SIR=3 in this case. The measurements of this simulation case are shown in FIG. 5.

(66) There are 6 loops (3 phase-to-phase loops and 3 phase-to-ground loops). The calculation of compensated voltage (at 80% of the line) and related polarity comparison of the phase A to ground loop is shown in FIG. 6. The remaining other 5 loops are assumed to remain stable (i.e., no trip) in this simulation case. According to FIG. 6, for the phase A to ground loop the operate value of the ratio ‘s’ enters the “operate” zone, meaning that a signal is issued to the circuit breaker to trip, at 8 ms after the fault occurs and thus the herein disclosed methods detect an internal fault of the phase A to ground loop at 10 ms, after 3 points of confirmation. As mentioned above, the remaining other 5 loops are assumed to remain stable and hence an external fault with very large security margin (s=1) is indicated for these loops and they stay in the “no operate” zone, meaning that they do not issue any signal to the circuit breaker to trip. The “operate” and “no operate” results are correct, fast and secure. The final trip signals of the 6 loops are shown in FIG. 7.

(67) Simulation results for an external fault will now be disclosed.

(68) An external phase-to-phase solid fault between phase A and phases B, C (F2 in FIG. 1) is assumed to occur at 100 ms and be located at the remote busbar (230 km). The power angle is 30 degrees and SIR=3 in this case. The measurements of this simulation case are shown in FIG. 8.

(69) All the 6 loops remain stable in this case. As an example, the calculation of compensated voltage (at 80% of the line) and related polarity comparison results for the phase A to ground loop is shown in FIG. 9. An external fault is detected with large security margin (s=1) and the loop stays in the “no operate” zone.

(70) The final trip signals of the 6 loops are shown in FIG. 10. The “no operate” results are correct and secure.

(71) Further aspects of scenarios of high SIR will now be disclosed.

(72) Most traditional time domain distance protection schemes based on amplitude comparison have an issue with low sensitivity for high SIR conditions. However, since the proposed polarity detection scheme does not need high amplitude to work correctly, the proposed polarity comparison based distance protection scheme has relatively high sensitivity even for high SIR conditions.

(73) One example with a high SIR is shown in FIG. 11. The line length is 230 km and the system model is shown in FIG. 1. In this case, SIR=6 which means the positive source impedance is up to about 350Ω. A phase-to-ground fault for phase A occurs at 70% of the protected line. The reach setting of proposed distance protection is 80% of the line. The voltage and current measurements are shown in FIG. 11 and the characteristics and operate results of proposed distance protection are shown in FIG. 12.

(74) As shown in FIGS. 11 and 12, the proposed distance protection scheme operates correctly. More simulation results with different SIR conditions are shown in Table 1. All the cases in the table are remote faults, where the fault points are located at 70% of the 230 km length line and the proposed distance protection is set to cover 80% of the line.

(75) TABLE-US-00001 TABLE 1 Operate time of remote faults (70%) with different SIR values SIR 0.2 0.5 1 2 3 4 5 6 7 8 9 10 15 20 Operate time (ms) 10 10 11 11 12 12 12 13 13 13 13 13 15 16

(76) As shown in Table 1, the proposed distance protection scheme is able to handle extreme case even when SIR=20. In addition, the operate time of the proposed distance protection is quite stable for different SIR conditions (from 10 ms to 16 ms) besides high sensitivity.

(77) The results have shown that the proposed distance protection has very high sensitivity to handle high SIR scenarios. This is an advantage since existing time domain distance protection schemes based on amplitude may fail to operate under such high SIR scenarios.

(78) FIG. 13 schematically illustrates, in terms of a number of functional units, the components of a fault detection arrangement 10a, 10b, 10c for fault detection of the transmission line 20 of the transmission system 25 according to an embodiment. Processing circuitry 1310 is provided using any combination of one or more of a suitable central processing unit (CPU), multiprocessor, microcontroller, digital signal processor (DSP), etc., capable of executing software instructions stored in a computer program product 1510 (as in FIG. 15), e.g., in the form of a storage medium 1330. The processing circuitry 1310 may further be provided as at least one application specific integrated circuit (ASIC), or field programmable gate array (FPGA).

(79) Particularly, the processing circuitry 1310 is configured to cause the fault detection arrangement 10a, 10b, 10c to perform a set of operations, or steps, as disclosed above. For example, the storage medium 1330 may store the set of operations, and the processing circuitry 1310 may be configured to retrieve the set of operations from the storage medium 1330 to cause the fault detection arrangement 10a, 10b, 10c to perform the set of operations. The set of operations may be provided as a set of executable instructions.

(80) Thus the processing circuitry 1310 is thereby arranged to execute methods as herein disclosed. The storage medium 1330 may also comprise persistent storage, which, for example, can be any single one or combination of magnetic memory, optical memory, solid state memory or even remotely mounted memory. The fault detection arrangement 10a, 10b, 10c may further comprise a communications interface 1320 at least configured for obtaining current values from the transmission system 25, to provide current values to the distance protection fault detection arrangement 1210, and for communications with another fault detection arrangement 10a, 10b, 10c. As such the communications interface 1320 may comprise one or more transmitters and receivers, comprising analogue and digital components. The processing circuitry 1310 controls the general operation of the fault detection arrangement 10a, 10b, 10c, e.g., by sending data and control signals to the communications interface 1320 and the storage medium 1330, by receiving data and reports from the communications interface 1320, and by retrieving data and instructions from the storage medium 1330. Other components, as well as the related functionality, of the fault detection arrangement 10a, 10b, 10c are omitted in order not to obscure the concepts presented herein.

(81) FIG. 14 schematically illustrates, in terms of a number of functional modules, the components of a fault detection arrangement 10a, 10b, 10c for distance protection of the transmission line 20 of the transmission system 25 according to an embodiment. The fault detection arrangement 10a, 10b, 10c of FIG. 14 comprises a number of functional modules; an obtain module 1410a configured to perform step S102, an obtain module 1410b configured to perform step S104, an obtain module 1410c configured to perform step S106, and a fault detection module 1410d configured to perform step S108.

(82) The fault detection arrangement 10a, 10b, 10c of FIG. 14 may further comprise a number of optional functional modules, such as a provide module 1410e configured to perform step S110. In general terms, each functional module 1410a-1410e may in one embodiment be implemented only in hardware and in another embodiment with the help of software, i.e., the latter embodiment having computer program instructions stored on the storage medium 1330 which when run on the processing circuitry makes the fault detection arrangement 10a, 10b, 10c perform the corresponding steps mentioned above in conjunction with FIG. 4. It should also be mentioned that even though the modules correspond to parts of a computer program, they do not need to be separate modules therein, but the way in which they are implemented in software is dependent on the programming language used. Preferably, one or more or all functional modules 1410a-1410e may be implemented by the processing circuitry 1310, possibly in cooperation with the communications interface 1320 and/or the storage medium 1330. The processing circuitry 1310 may thus be configured to from the storage medium 1330 fetch instructions as provided by a functional module 1410a-1410e and to execute these instructions, thereby performing any steps as disclosed herein.

(83) FIG. 15 shows one example of a computer program product 1510 comprising computer readable storage medium 1530. On this computer readable storage medium 1530, a computer program 1520 can be stored, which computer program 1520 can cause the processing circuitry 1310 and thereto operatively coupled entities and devices, such as the communications interface 1320 and the storage medium 1330, to execute methods according to embodiments described herein. The computer program 1520 and/or computer program product 1510 may thus provide means for performing any steps as herein disclosed.

(84) In the example of FIG. 15, the computer program product 1510 is illustrated as an optical disc, such as a CD (compact disc) or a DVD (digital versatile disc) or a Blu-Ray disc. The computer program product 1510 could also be embodied as a memory, such as a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM), or an electrically erasable programmable read-only memory (EEPROM) and more particularly as a non-volatile storage medium of a device in an external memory such as a USB (Universal Serial Bus) memory or a Flash memory, such as a compact Flash memory. Thus, while the computer program 1520 is here schematically shown as a track on the depicted optical disk, the computer program 1520 can be stored in any way which is suitable for the computer program product 1510.

(85) The inventive concept has mainly been described above with reference to a few embodiments. However, as is readily appreciated by a person skilled in the art, other embodiments than the ones disclosed above are equally possible within the scope of the inventive concept, as defined by the appended patent claims.