A system for dynamic optimisation and stabilisation of an electrical grid by using a multilevel method of Additive-Increase/Multiplicative-Decrease-AIMD method, wherein said system is adjusted to regulate power in the electrical grid, which comprises at least two levels, as well as smart meters, characterized in that the system comprising at least two AIMD clusters on adjacent grid levels, wherein each AIMD cluster comprises an AIMD server and at least one AIMD consumer, wherein the grid level of the AIMD server defines the level of the entire AIMD cluster, and the AIMD consumer in this AIMD cluster is on the next lower level, wherein each AIMD server executes a first part of the AIMD algorithm and each AIMD consumer executes a second part of the AIMD algorithm,
- wherein the first part of the AIMD algorithm is executed as follows and wherein N represents the level of an individual AIMD server:
- if a net power PN on level Nis smaller than a dynamic maximum power PDNmax, that is PN<PDNmax, then the AIMD server on level N assigns a value 0 to a control parameter cN (cN=0);
- if PN≥PDmax, then the AIMD server on level N assigns a value 1 to the control parameter cN (cN=1),
- wherein the value of the control parameter cN in the relevant time window tis sent to all AIMD consumers within this AIMD cluster,
- wherein the second part of the AIMD algorithm is executed as follows:
- a) if cN=0 and the dynamic maximum power of the AIMD consumers on level N−1, which are connected to the AIMD server on level N, in the relevant time window PD(N-1)max(t) is lower than a fixed maximum power of these AIMD consumers PS(N-1)max, that is PD(N-1)max(t)<PS(N-1)max, then the dynamic maximum power of these AIMD consumers in a subsequent time window PD(N-1)max(t+Δt) is increased by an additive constant α(N-1), that is PD(N-1)max(t+Δt)=PD(N-1)max(t)+α(N-1), otherwise the power PD(N-1)max does not change;
- b) if cN=1, the dynamic maximum power in the subsequent time window PD(N-1)max(t+Δt) is decreased by using a randomly generated parameter γ(N-1), of which the random value is between 0 and 1, where the generated parameter γ(N-1) is compared to a predefined probability space to determine which predefined probability space is relevant in this time window and wherein each probability space has a corresponding predefined multiplicative decrease factor β(N-1)n, where n stands for the corresponding probability space, which is used to decrease the dynamic maximum power in the subsequent time window PD(N-1)max(t+Δt), and the probability of decrease for a specific multiplicative decrease factor (β(N-1)n) is calculated using a following equation:
- P
- D
-
- (
- N
- -
- 1
- )
-
- max
- (
- t
- +
- Δ
-
- t
- )
- =
- P
- D
-
- (
- N
- -
- 1
- )
-
- max
- (
- t
- )
- ·
- β
- (
- N
- -
- 1
- )
-
- n
- ,
- where β(N-1)n is defined in relation to the probability space of the randomly generated parameter γ(N-1),
- wherein the value of the calculated dynamic maximum power PD(N-1)max(t+Δt), which is calculated by the AIMD consumer, which belongs to the AIMD cluster on level N, represents input data for the AIMD server, which belongs to the AIMD cluster on a lower level N−1 to execute the first part of the AIMD algorithm when PD(N-1)max is compared to the net power P(N-1),
- where the net power PN, PN-1 on a specific level is measured using the smart meter on this level.