Login| Sign Up| Help| Contact|

Patent Searching and Data


Title:
ELECTRICAL INVERTER AND METHOD OF OPERATION
Document Type and Number:
WIPO Patent Application WO/2014/060065
Kind Code:
A2
Abstract:
The present invention relates to means for adapting a power rotation angle of an electrical inverter to match the unknown grid impedance at the point of connection. The invention further relates to an electrical inverter comprising such means. The present invention also relates to a method of adapting power rotation angle of an electrical inverter to match the unknown grid impedance at the point of connection.

Inventors:
JOSTOCK MARKUS (DE)
Application Number:
PCT/EP2013/002908
Publication Date:
April 24, 2014
Filing Date:
September 27, 2013
Export Citation:
Click for automatic bibliography generation   Help
Assignee:
UNIV LUXEMBOURG (LU)
International Classes:
H02J3/38; H02J3/46
Foreign References:
GB2483879A2012-03-28
Other References:
See references of EP 2907214A2
Attorney, Agent or Firm:
LAMPE, SIGMAR (Université du Luxembourg162a avenue de la Faïencerie, Luxembourg, LU)
Download PDF:
Claims:
CLAIMS

Means for adapting a power rotation angle (3) of an electrical inverter to an unknown grid impedance angle (Θ) of the electric grid or busbar at a point of connection of the electrical grid or busbar with the inverter output, the output power angle adapting means comprising

means for detecting a momentary active power (P(t)) of the inverter; means for detecting a momentary reactive power (Q(t)) of the inverter; means for generating a rotated momentary active power ( rot( and a rotated momentary reactive power (Qmt(t)) by carrying out a rotation operation of the momentary active power (P(t)) and the momentary reactive power (Q(t)) by a power rotation angle (#); and

means for providing a reference frequency ( (t)),

wherein the reference frequency providing means further comprise means for providing a nominal frequency (/¾), wherein the nominal frequency corresponds to the nominal grid or busbar frequency, and

means for applying a momentary jitter signal (jitter(t)) to the nominal frequency ( <j) for generating a reference frequency (/¾ ))> so that the reference frequency ( (t)) varies in time; and

the inverter further comprises means for adapting the power rotation angle (#) in response to the rotated momentary active power (PTOt(t) and the rotated momentary reactive power (Qrot(t))-

The means of claim 1 , further comprising means for detecting fluctuations in time in one or both of the rotated momentary active power and the rotated momentary reactive power.

The means of claim 1 or 6, further comprising means for adapting the power rotation angle ($) until it approaches, or optionally matches, the grid impedance angle (Θ).

The means of claim 3, further comprising means for judging the quality of the match between the power rotation angle (&) and the grid impedance angle (Θ).

The means of claim 4, wherein the means for judging the quality of the match further comprise means for observing the rotated momentary active and reactive power which are altered through a jittering inverter frequency. The means of any one of the preceding claims, further comprising means for controlling a momentary inverter output frequency ( vsi(t)) in response to the rotated momentary active power (PTOt(t) and the reference frequency ( (t)).

The means of claim6, wherein the means for controlling the momentary inverter output frequency ( vsi( comprise frequency droop control means.

The means of any one of the preceding claims, further comprising means for controlling a momentary inverter output voltage (wvsi( in response to the rotated momentary reactive power (Qxoi(t);

The means of claim 8, wherein the means for controlling a momentary inverter output voltage («vsi( comprise voltage droop control means.

The means of any preceding claim, wherein the means for adapting the power rotation angle (9) comprise a minimum correlation finding means in response to only the rotated momentary reactive power (Qrot(t))-

The means of claim 10, wherein the minimum correlation finding means further comprise means for calculating at least one harmonic coefficient.

The means of claim 1 1 , wherein the means for calculating at least one harmonic coefficient further comprise means for performing a discrete Fourier transformation.

The means of any one of claims 10 to 12, wherein the minimum correlation finding means further comprises means for performing an optimal control.

An electrical inverter comprising

means for outputting an alternating current (AC) power to an electric grid or busbar, wherein the grid or busbar represents an unknown impedance (Z) to the inverter output, the unknown impedance (Z) comprising a resistance (R) and a reactance (X) with Z = R + X, and wherein the unknown grid impedance angle is determined by arctan R/X; and

means for adapting an output power angle of the inverter to the unknown grid impedance angle of the electric grid or busbar at a point of connection with the inverter output in accordance with one of the preceding claims. The electrical inverter of claim 14, wherein inverter is of the voltage source or the current source type.

A method for adapting the power rotation angle (3) of an electrical inverter, the method comprising the steps of

outputting an alternating current (AC) power to an electric grid or busbar, the electric grid or busbar representing at the inverter output an unknown impedance (Z) comprising a resistance (R) and a reactance (X), wherein an unknown grid impedance angle Θ is determined by arctan R/X;

detecting a momentary active power (P(t)) of the inverter- means for detecting a momentary reactive power (Q(t)) of the inverter; generating a rotated momentary active power (Prot( and a rotated momentary reactive power (QTOt(t)) by carrying out a rotation operation of the momentary active power (P(t)) and the momentary reactive power (Q(t)) by a power angle (6>); and

providing a reference frequency ( (t));

wherein the step of providing a reference frequency (fK ) further comprises applying a momentary jitter signal (jitter(t)) to a nominal frequency f0 for generating the reference frequency, so that the reference frequency (t) applied to the inverter output frequency controlling means varies in time; and the method further comprises the step of adapting the power rotation angle (#) in response to the rotated momentary active power (Prot( and the rotated momentary reactive power (£?rot( ))-

The method of claim 12, further comprising means for detecting fluctuations in time in one or both of the rotated momentary active power and the rotated momentary reactive power.

The means of claim 16 or 17, further comprising means for adapting the power rotation angle (§) until it approaches, or optionally matches, the grid impedance angle (Θ).

The means of claim 18, further comprising means for judging the quality of the match between the power rotation angle ($) and the grid impedance angle (Θ).

20. The means of claim 19, wherein the means for judging the quality of the match further comprise means for observing the rotated momentary active and reactive power which are altered through a jittering inverter frequency.

The method of any one of claims 16 to 20, further comprising controlling a momentary inverter output frequency ( vsi( ) m response to the rotated momentary active power (PTOt(t) and the reference frequency (f(t));

The method of claim 17, wherein the step of controlling the momentary inverter output frequency ( vsi( further comprises applying a frequency droop control.

The method of any one of claims 16 to 22, further comprising controlling a momentary inverter output voltage (wvsi( in response to the rotated momentary reactive power (Qmt(t)',

The method of claim 23, wherein the step of controlling a momentary inverter output voltage («vsi( further comprises applying a voltage droop control.

The method of any one of claims 16 to 24, wherein the step of adapting the power rotation angle (#) comprises finding a minimum correlation in response to only the rotated momentary reactive power (Qrot(t))-

The method of any one of claims 16 to 25, wherein the step of finding a minimum correlation further comprises calculating at least one harmonic coefficient.

The method of any one of claims 16 to 26, wherein the step of calculating at least one harmonic coefficient further comprises performing a discrete Fourier transformation.

The method of any one of claims 16 to 27, wherein the step of finding a minimum correlation further comprises performing an optimal control.

29. The method of any one of claims 16 to 28, wherein the inverter is of the voltage source or the current source type.

Description:
ELECTRICAL INVERTER AND METHOD OF OPERATION

Supported by the Fonds National de la Recherche, Luxembourg (PHD-09-183). Technical Field

The present invention relates to means for adapting an power rotation angle of an electrical inverter to match the grid impedance angle at the point of interconnection. The invention further relates to an electrical inverter comprising such means. The present invention also relates to a method of adapting power rotation angle of an electrical inverter to match the grid impedance angle at the point of interconnection.

Background The increasing penetration with distributed renewable energy sources currently characterizes the transformation of the electrical power grid into a sustainable, environmentally friendly and efficient power supply of the future. Renewable energy sources are mostly injecting electrical energy with voltage source inverters (in the following referred to as VSI). When VSIs are operated in droop control mode, mimicking the behaviour of classical rotating synchronous generators, they are able to sustain and equally share all loads of a - potentially islanded - power grid (Alfred Engler, Regelung von Batteriestromrichtern in modularen und erweiterbaren Inselnetzen, PhD thesis, Universitat Gesamthochschule Kassel, May 2001, incorporated herein by reference in its entirety).

The classic droop control concept relates active power load to the grid frequency and the reactive power load to the grid voltage magnitude by a static head line as indicated in Fig 1 which illustrates the known static droop by relating Q with voltage and P with frequency.. The VSI droop control works with a characteristic head line. The VSI adapts its injection frequency to the momentary rate of injected active power with

Wherein ^vsi is the inverter output frequency, fo is the base frequency, P is the momentary active power, P ra ted is the inverter's rated power, and droop is the maximum frequency deviation at Prated expressed in per cent.

Similarly, it adapts the voltage amplitude according to the actually injected reactive power. When injecting full rated active power, it changes the injection frequency to maximum droop (of e.g. 1% in Fig 1). Usually the base frequency fo is a fixed reference value for the droop controller. The working principle of the droop control of voltage source inverters assures that each inverter adapts its frequency independently, according to its actual power injection. When the injection frequency differs from the current grid frequency, the power angle Θ of the VSI voltage drifts with respect to the voltage angle at the end of the connected branch and induces a change in the power flow. This changed power flow provokes frequency changes in the other VSIs in the grid - as these adapt their respective injection frequencies according to the changed power flow - and the whole system finds a new equilibrium point at a new frequency.

This concept allows for decoupled control of frequency and voltage in power grids with inductive line characteristics, i.e. generally the medium and high voltage grids. Distributed renewable energy sources mostly inject their power in the low voltage grid. Due to the resistive nature of the low voltage lines, independent control of frequency and voltage is no longer possible with classic droop control. As a consequence, a change in the grid's active power load not only has an effect on frequency, but also on the grid voltage. Vice versa, a change in reactive power load not only influences the grid voltage, but also the grid frequency. This circumstance is illustrated in Fig 2 showing the known effect of rotated active and reactive power on frequency and voltage.

In an inductive grid, the load S 1 = P + ]Q 1 would have an effect of f \ and u \ on the power grid. But in a low voltage grid, the resistive part R of the line impedance Z = R + jX may be even larger than the reactance X. Thus the impedance relation R/X≠ 0 cannot be neglected and the grid impedance angle Θ = arctan /7? is less than 90°. The active/reactive power coordinates are rotated by the angle φ = — Θ with respect to the frequency/voltage coordinates, as shown in Fig 2. The very same load situation S \ has a rotated effect on the low voltage grid and causes a change of f and w rot on frequency and voltage. One solution to this problem is to not apply the droop control directly on the measured values of S = P + ]Q, but to first rotate the measured power backwards by the angle

φ = -— Θ instead and exert the control on the transformed power coordinates ro t, Qmt-

-frot cos ψ — sin φ ' P ' sin # — COS 0 P

Qiot sin ψ cos y Q cos # sin (9 Q

As the grid impedance angle Θ is usually not precisely known at the point of connection PoC, the rotation may be performed with an estimated rotation angle $ with Prot = sin $ · P - cos $ · Q and

This corrective control concept has been described e.g. by:

J. Van den Keybus, A. Woyte, J. Driesen, K. De Brabandere, B. Bolsens: "A voltage and frequency droop control method for parallel inverters", in 35th Annual IEEE Power Electronics Specialists Conference, pages 2501-2507, 2004, incorporated herein by reference in its entirety;

Frede Blaabjerg, Uffe Borup, Lucian Asiminoaei, Remus Teodorescu: "A digital controlled pv-inverter with grid impedance estimation for ens detection", in IEEE Transactions on Power Electronics, 20(6): 1480-1490, November 2005, incorporated herein by reference in its entirety; and

A. Luna, P. Rodiguez, R. Teodorescu, J Vasquez, J. Guerrero: "Adaptive droop control applied to distributed generation inverters connected to the grid", in IEEE Transactions on Industrial Electronics, 56(10):4088^4096, October 2009, incorporated herein by reference in its entirety.

Fig 3 illustrates how the droop control is applied to the rotated values of active and reactive power. In the block diagram of Fig 3 the currently injected active power P and reactive power Q are measured and/or calculated. The rotation block 31 rotates P and Q by an angle Θ into P mt and grot, on which the droop control block 32 for frequency / and the droop control block 33 for voltage u are exerted. Frequency β si and voltage «vsi of the AC source are applied to the power grid 34, resulting in active and reactive power. Since the grid impedance angle Θ at the point of connection is not known, it has to be detected in order to find the optimal power rotation angle $ for the VSI, which ideally should match the grid impedance angle Θ.

Prior Art

There are several known possibilities to obtain the power angle at the point of connection, mostly by determining the exact grid impedance (magnitude and angle) over a wide frequency range at the PoC [E. Gunther R. Bergeron A. Robert, T. Deflandre: "Guide for assessing the network harmonic impedance", in Proc. Inst. Elect. Eng., 2:301-310, 1997, incorporated herein by reference in its entirety]. Possible methods are:

By estimation The impedance angle could be estimated by general knowledge of the impedance relation R/X of low voltage grids. This assumed angle will most certainly be wrong, but the non-optimal control and the remaining coupling effect between frequency and voltage are accepted and might be considered "good enough".

By calculation

The grid utility company can provide the impedance angle at the point of connection. This value might be the result of calculations based on the grid model of the utility and thus deviate from the real value at the point of connection. The utility might not have a possibility to measure or calculate values for every point of connection. Depending on the number of VSIs installed by different suppliers, inquiring the impedance angle from the utility might take long and be uneconomic.

By specific hardware

The impedance can be measured at the PoC with a particular impedance measuring device. This approach is likely to be inefficient due to the necessary costly equipment and time and effort consumption.

By the VSI

The CPU or DSP of the VSI can be used to implement a measuring method. While numerous publications have proposed different schemes, not all methods can be implemented with the given limited capacities and capabilities of the VSI. Is has been proposed by D. Thomas, B. Palethorpe, M. Summer: "System impedance measurement for use with active filter control" in Power Electronics and Variable Speed Drives, volume 475, pages 24-28. IEE, September 2000, incorporated herein by reference in its entirety; and Frede Blaabjerg et al., loc. cit. that a measurement could be done e.g. by injection of harmonic currents. A. Luna et al., loc. cit. have suggested high precision voltage and current frequency measurement. Most commonly the impedance is calculated based on a measurement of voltage and current (transient or steady state) at a specific (inter)harmonic frequency. Some of these methods interpolate the impedance at fundamental frequency, calculating a value deviating from the real situation. Also, harmonic injection leads to waveform distortion during the measurement.

Disadvantages of the prior art

At least some drawbacks of most methods for determining the impedance angle at the point of connection are

• increased computation efforts; • for some methods specific and expensive hardware is required;

• often demand change of / interference with the specific current control algorithm;

• both impedance magnitude and angle are determined, although only the angle is needed for decoupled control of frequency and voltage.

Without knowledge of the grid power angle at the point of connection, a decoupled control of frequency and voltage cannot be exerted, leading to cross-influence of active power to voltage and reactive power to frequency and thus to undesired effects up to grid instability. Technical problem

It is therefore an object of the present invention to overcome or at least alleviate the disadvantages of the prior art. In particular it is an object of the present invention to propose an electrical inverter and a method which allow determining the exact impedance angle of the power grid at the inverter point of connection (PoC).

Summary of the invention

A first aspect of the present invention relates to means for adapting a power rotation angle ϋ of an electrical inverter to an unknown grid impedance angle Θ of the electric grid or busbar at a point of connection of the electrical grid or busbar with the inverter output. The power rotation angle adapting means comprise means for detecting a momentary active power of the inverter; means for detecting a momentary reactive power of the inverter; means for generating a rotated momentary active power and a rotated momentary reactive power by carrying out a rotation operation of the momentary active power and the momentary reactive power by a rotation angle; means for providing a reference frequency, means for controlling a momentary inverter output frequency in response to the rotated momentary active power and the reference frequency; and means for controlling a momentary inverter output voltage. The reference frequency providing means further comprise means for providing a nominal frequency, wherein the nominal frequency corresponds to the nominal grid or busbar frequency, and means for applying a momentary jitter signal to the nominal frequency, so that the reference frequency varies in time. The inverter further comprises means for adapting the power rotation angle in response to the rotated momentary active power and the rotated momentary reactive power.

The means for adapting the power rotation angle in response to the detected rotated momentary active power and the detected rotated momentary reactive power may further comprise means for detecting fluctuations in time in one or both of the rotated momentary active power and the rotated momentary reactive power.

The means for adapting the power rotation angle may further comprise means for adapting the power rotation angle & until it approaches, or optionally matches, the grid impedance angle Θ.

The means for adapting the power rotation angle may further comprise means for judging the quality of the match.

The means for judging the quality of the match may further comprise means for observing the rotated momentary active and reactive power which are altered through a jittering inverter frequency. The means may further comprise meaiis for controlling a momentary inverter output frequency in response to the rotated momentary active power and the reference frequency.

The means for controlling the momentary inverter output frequency may comprise frequency droop control means.

The means of any one of the preceding claims may further comprise means for controlling a momentary inverter output voltage in response to the rotated momentary reactive power.

The means for controlling a momentary inverter output voltage may comprise voltage droop control means.

The means for manipulating the power angle may comprise a minimum correlation finding means in response to only the rotated momentary reactive power.

The minimum correlation finding means further may comprise means for calculating at least one harmonic coefficient. The means for calculating at least one harmonic coefficient may further comprise means for performing a discrete Fourier transformation. The minimum correlation finding means may further comprise means for performing an optimal control, optionally e.g. a linear quadratic controller or a H-infinity controller, such as known to the skilled person (http://en.wikipedia.org/wiki/Optimal_control). A second aspect of the present invention relates to an electrical inverter comprising means for outputting an alternating current (AC) power to an electric grid or busbar, wherein the grid or busbar represents an unknown impedance to the inverter output, the unknown impedance comprising a resistance and a reactance with Z = R + X, and wherein the unknown grid impedance angle Θ is determined by arctan R/X. the inverter also comprises means for adapting an power rotation angle ϋ of the inverter to the unknown grid impedance angle Θ of the electric grid or busbar at a point of connection with the inverter output in accordance with one of the configurations of the first aspect.

The inverter may be of the voltage source or the current source type.

A third aspect of the present invention relates to a method for adapting the power angle of an electrical inverter, the method comprising the steps of: outputting an alternating current (AC) power to an electric grid or busbar, the electric grid or busbar representing at the inverter output an unknown impedance comprising a resistance and a reactance, wherein an unknown grid impedance angle Θ is determined by arctan R/X; detecting a momentary active power of the inverter; means for detecting a momentary reactive power of the inverter; generating a rotated momentary active power and a rotated momentary reactive power by carrying out a rotation operation of the momentary active power and the momentary reactive power by a rotation angle ϋ; providing a reference frequency, wherein the reference frequency corresponds to the nominal inverter output frequency. The step of providing a reference frequency further comprises applying a momentary jitter signal to a nominal frequency, so that the reference frequency applied to the inverter output frequency controlling means varies in time. The method further comprises the step of adapting the power rotation angle ϋ in response to the detected rotated momentary active power and the detected rotated momentary reactive power.

The method may further comprise means for detecting fluctuations in time in one or both of the rotated momentary active power and the rotated momentary reactive power. The method may further comprise adapting the power rotation angle & until it approaches, or optionally matches, the grid impedance angle Θ.

The method may further comprise judging the quality of the match. The step of judging the quality of the match may further comprise observing the rotated momentary active and reactive power which are altered through a jittering inverter frequency.

The method may further comprise controlling a momentary inverter output frequency in response to the rotated momentary active power and the reference frequency.

The step of controlling the momentary inverter output frequency may further comprise applying a frequency droop control.

The method may further comprise controlling a momentary inverter output voltage in response to the rotated momentary reactive power. The step of controlling a momentary inverter output voltage may further comprise applying a voltage droop control.

The step of manipulating the power rotation angle may further comprise finding a minimum correlation in response to only the rotated momentary reactive power.

The step of finding a minimum correlation further comprises calculating at least one harmonic coefficient.

The the step of calculating at least one harmonic coefficient further comprises performing a discrete Fourier transformation.

The step of finding a minimum correlation further comprises performing an optimal control. The may be applied to an inverter of the voltage source or the current source type.

Further advantageous embodiments are described in the dependent claims and the following description. Brief description of the drawings

These aspects and further embodiments of the present invention will be described in the following with reference to the enclosed non-limiting drawings. Fig 1 is a diagram illustrating the known static droop by relating Q with voltage and P with frequency. Fig 2 is a diagram showing the known effect of rotated active and reactive power on frequency and voltage.

Fig 3 is a diagram illustrating the known droop control applied to the rotated values of active and reactive power.

Fig 4 is a diagram showing the droop head line in accordance with the present invention.

Fig 5 is a block diagram of a control structure of an electrical inverter in accordance with the present invention.

Fig 5a is a block diagram of a control structure of an electrical inverter and angle quality evaluation and angle manipulation in accordance with the present invention. Fig 6 is a flowchart illustrating a method in accordance with the present invention.

Fig 7 is a block diagram of an embodiment with an optimal controller to minimise jitter Q. Fig 8 is a power over time diagram showing the reaction of P, P mt , Q and Q I0t . Detailed technical description of the invention

For decoupled control of frequency and voltage, only an adaptation to the grid impedance angle is required, the exact grid impedance magnitude is not measured. This is particularly done by exploiting the following VSI characteristic:

Different from classic electric generators with rotating masses, a VSI can adjust its output frequency arbitrarily, which the present invention proposes to exploit for impedance detection.

Based on the working principle explained on page 1, if a VSI slightly alters its injection frequency, it will detect a change in its injected power rates P TO t/P and Q Q. The present invention suggests using this effect in order to detect the grid impedance angle Θ for optimal decoupled droop control: When the output frequency ^vsi is jittered, the rotation angle θ of power rotation is altered until only P jitters and Q remains untouched from changes in frequency This means that the rotation angle & then matches the unknown grid impedance angle Θ at the inverter's point of connection and the rotation angle & is considered as optimal.

Further, the present invention proposes to combine an adaptive change of rotation angle until optimal adaptation to the impedance angle is reached, and a slowly jittering base frequency for the droop controller as an indicator of the rotation effect.

The exact signal form chosen for the jitter signal may be chosen from any suitable waveform, e.g. sinusoidal, sawtooth, triangle, binary, rectangular, or any other oscillating signal forms.

The present invention proposes to detect the optimal rotation angle 9, which is dependent on the grid impedance angle Θ, by the VSI applying a slowly oscillating base frequency fo + jitter(i) to the subordinate power or current controller. For the further discussion it is assumed that the example of a subordinate droop control has been chosen. As shown in Fig 4, the head line swings vertically with an amplitude A around the base value in the grey band. The jittering of the inverter's injection frequency is independent from the subordinate controller. It can be applied to droop controllers with different types of droops or other types of controllers, like e.g. maximum power point controllers.

If the base frequency is superposed with a jitter function jitterif), the jitter will also be noticeable at the power measurement. The measured active and reactive power are rotated by the rotation angle $ into rotated coordinates of active and reactive power Ρ τοι and Q TOt With a non-optimal rotation angle the swing will be visible in both rotated active power rot and rotated reactive power Q rot . This effect is exploited to adjust the rotation angle in accordance with the present invention: For detection, the power rotation angle ϋ is changed until the jitter in P rot reaches its maximum and/or the jitter in Qrot reaches its minimum.

Fig 5 a is a block diagram of a control structure of an electrical inverter and angle quality evaluation and angle manipulation in accordance with the present invention. The evaluation of the maximum jitter in rot or the minimum jitter in Q TOt may be separated from the adaptation of the rotation angle 9·. In block 57a, the evaluation of the rotation effect quality can be done automatically, e.g. in block 55a with a detection of the jitter signal in the rotated reactive power Q TO t(t) or e.g. by the calculation of harmonic Fourier coefficients, or it can be done manually by a human user. In block 58a the adaptation of the power rotation angle $ can be done e.g. in a linear sweep from 0 to π/2 (i.e. 90°), or e.g. with an optimal controller freely incrementing or decrementing θ between 0 and π/2 as the quality evaluation indicates, thus approaching the optimal rotation angle $ 0 pt in circular trajectories or e.g. with non-linear binary search algorithms or by a human user with a manual control.

Preferably the frequency jitters slowly around a base value fo with a swing period significantly larger than the grid period, i.e. T osc » 1/ f 0 , e.g. T osc ~ 0.5 seconds, and with an amplitude considerably smaller than the grid frequency A « f 0 , e.g. A ~ 0.00025 -fo .

Description of an embodiment

The frequency set point is overlaid with a sinusoidal jitter signal, e.g.

jitter(t) = A · sin (^i)

The optimization can e.g. be done by detecting the exciter signal jitterif) in the measured reactive power signal Q(t) . Since the oscillation period T osc is known, only one complex Fourier coefficient needs to be calculated with very low computing resources from the measured (discrete) reactive power measurements Q[k]. With the sample time r sam pie the number of samples per oscillation is

and for a sampling buffer of s oscillations the Fourier coefficient of the 5 th fundamental is calculated as

X S QW (cos {2nff k) - j sin (2π£*) ) . Usually we can set s = 1. In order to find the minimum correlation, the absolute value of the correlation signal

can be fed to an optimal controller.. The controller output is the rotation angle & which is fed into the rotation matrix. While the droop control is continuously enacted on the active and reactive power, the influence of the slow frequency oscillations on the reactive power will be influenced by the actual rotation angle

Fig 5 is a block diagram depicting a control structure of an electrical inverter in accordance with the present invention. The adaptation of the rotation angle 3 through a minimum correlation finder 56, is based on the measurement of the reactive power Q. The currently injected active power P and reactive power Q are measured and/or calculated. The rotation block 51 rotates P and Q by an angle Θ into P I0t and Q mt . In block 55, a jitter signal jitterif) is added to the reference frequency fo.

The droop control block 52 therefore applies a modulated headline corresponding to fo + jitterif) for droop control of frequency Jvsi- By contrast, the droop control block 53 applies a non-modulated headline for droop control of voltage «vsi-

Frequency J si and voltage MVSI of the AC source are applied to the power grid 54, resulting in active and reactive power. The jitter signal jitterif) which has been applied to the reference frequency f becomes detectable in the measured values of P(t) and Q(t), and in the rotated values P ro t(t) and Q TO t(t).

Via a minimum correlation finder 56 the rotation angle Θ is changed, until correlation coefficient reaches its minimum in Q mt . Thus the frequency variation reaches a minimum influence on Q, and simultaneously reaches a maximum influence on active power P when the rotation angle Θ matches the grid impedance angle at the point of connection.

At the point of minimal correlation between / and Q, the optimum power rotation angle & 0 p t is determined as the grid impedance angle Θ. With this optimum power rotation angle 9 opt the droop control can be applied to the rotated values of P and Q for independent control of frequency and voltage.

The general method is illustrated as a flowchart in Fig 6. When, e.g. in normal operation of the inverter, an adaptation of the power rotation angle 9 is requested, see step 61 , step 62 foresees that a jitter signal jitterif) is applied to the reference frequency fo. In step 63 the power rotation angle & is changed, and reactive power is monitored for a decrease in Q m and/or active power is monitored for an increase in rot . If this is the case, step 64 foresees a further change of the power rotation angle 9·, and step 63 is executed again. If no further change in reactive power Q mX and/or active power P mt is observed, step 65 provides that the output power rotation angle 9 is stored. Finally, in step 66 the jitter signal is stopped and normal operation is resumed.

The initial output power rotation angle 9 depends on the application environment and the search algorithm. In a medium voltage grid the initial output power rotation angle 9 may e.g. be chosen at 0, and then increase; in a low voltage grid it may be preferable to start at 45° and then let the minimum finder control operate. It may be advantageous if the minimum search finder can also apply negative increments, i.e. decrease the output power rotation angle In the embodiment described here, the rotation angle search began at an angle of $ = 0.

The search of the optimal rotation angle can be done arbitrarily in the range of 9 from 0 to π/2. The search can done e.g. linearly in a single sweep or e.g. in a circular search or e.g. in a non-linear binary search.

The increment of the power rotation angle θ depends on the sample rate, i.e. on the clock rate of execution of the program, respectively how fast the power can be measured and calculated, as well as on setting of the minimum finder control. The person skilled in the art is aware of the classical control problem to tune a controller so that it quickly reaches its set point, perhaps including an overshoot, but without becoming instable. The increment thus depends strongly on the implementation chosen. In the implemented embodiment the step size was variable due to the type of controllers used: the bigger the DFT coefficient, the bigger the step size. Fig 7 shows a block diagram of an implementation with an optimal controller to minimise jitter in Q TOt . Fig. 7 shows the adaptation of the rotation angle through a minimum correlation finder, based on the measurement of the reactive power Q mt . The currently injected active power P and reactive power Q are measured and/or calculated. The rotation block 71 rotates P and Q by an angle $ into P mt and Q mX . In block 75, a jitter signal jitterif) is added to the reference frequency fo.

The droop control block 72 therefore applies a modulated headline corresponding to fo + jitter{f) for droop control of frequency f s \ - By contrast, the droop control block 73 applies a non-modulated headline for droop control of voltage wvsi-

Frequency ^vsi and voltage wvsi of the AC source are applied to the power grid 74, resulting in active and reactive power. The jitter signal jitterif) which has been applied to the reference frequency f becomes detectable in the measured values of Pit) and Qif), and in the rotated values P ro t( and Q TO iif)-

Although the jitter signal jitterif) ma Y be chosen sinusoidal, sawtooth, triangle, binary or any other oscillating signal forms, in the present embodiment it is preferred that the jitter signal jitterif) be a sine wave, hence the minimum correlation finder 76 comprises a Digital Fourier Transformation (DFT) 77 and an optimal controller 78. The DFT 77 calculates preferably the coefficient of first sine fundamental of the jitter signal jitter{t) from the value of £? ro t( - Of course, one or more alternative fundamentals might also be used. Via a controller 78 the power rotation angle & is changed, until the DFT coefficient reaches its minimum in Q l0t . Thus the frequency variation reaches a minimum influence on Q, and simultaneously reaches a maximum influence on active power P when the power rotation angle & matches the grid impedance angle Θ at the point of connection. At the point of minimal correlation between and Q, the impedance angle is determined as $ opt . With this optimal angle $ opt the droop control can be applied to the rotated values of P and Q for independent control of frequency and voltage. As can be seen in Fig 8, in the chosen implementation the swing in Q minimises only two seconds after the algorithm was started.

Other implementations are, of course, possible. Instead of using an optimum controller on the Fourier coefficient in Q(f), a simple angle sweep from 0 to π/2 could be performed manually. In fact, it is obvious to the person skilled in the art there are several ways to implement the present invention in order to find the optimal power angle. The present invention is based on the concept to exploit a frequency oscillation and read back power (or current) to find the power rotation angle $ matching grid impedance angle Θ.

A subordinate droop control is not mandatory for the present invention to work. The power angle detection can be performed in inverters with different control concepts.

Advantages of the invention

There are several favourable properties of this solution compared to the state of the art: · An adaptation is performed only on the impedance angle. No need to calculate the exact impedance value.

• No particular measurement equipment is necessary.

• No interference with (and thus no dependency on) the underlying control mechanism (be it droop controller or other).

· No low-level interference with current controller (be it PI control on d/q coordinates, THD optimal controllers, etc.), as no additional current is injected.

• No voltage or current waveform distortion.

• Very low computational effort for VSI DSP or CPU as only one DFT coefficient is calculated (for chosen implementation).

• Rather low sample rate possible for measurement of power.

• Low power disturbance in grid, as oscillation amplitude can be chosen according to actual power injection.

· Works for single phase systems as well as for 3 phase systems.

• The present invention works also for other droop modes like Q(P) or φ(Δ/) on a stiff grid, such as those mentioned in J. Van den Keybus, loc. cit.

• Can be run superposed on normal operation and repeated as deemed necessary.

• Simultaneous adaptation of several VSI possible by choice of arbitrary oscillation periods T osc (longer duration for adaptation).

• VSI equipped with the present invention are able to control voltage and frequency in an independent manner at any point of connection.