Microgrid power flow monitoring and control
A microgrid power flow monitoring and control system is described herein. The control system may determine active and reactive power sharing shortage on the electric power delivery system. The control system may utilize the control strategies of generation units, such as ISO control, droop control and constant power control to estimate power flow within a microgrid or other isolated system. A control strategy of one or more generators may be modified based on the determined power flow.
1. A method for evaluating a power flow on a bus, i, in at least a portion of an electric power delivery system to perform a control operation, comprising:
obtaining power characteristic values of power on the bus, i, of the electric power delivery system, including at least a frequency value and a voltage value of the power on the bus, i;
calculating a power characteristic-dependent load value based on at least one of the power characteristic values and a percentage of a load on the bus, i, that is power characteristic-dependent;
obtaining data identifying at least one of an operating mode and a control mode of the electric power delivery system;
calculating a control-operating value based on a product of a control-mode coefficient, β i , and an operating value, X, of at least one generator connected to the bus, i;
obtaining data identifying an angle of the power on the bus, i;
calculating a power generation reference value, P i Ref , based on:
a sum of the calculated power characteristic-dependent load value, the calculated control-operating value, and the obtained voltage value and
the angle of the power on the bus, i; and,
a shared vale S v calculated to satisfy S v =|V i |Σ n=1 N |Y in ∥V n |cos (δ i −δ n −θ in ),
where:
|V i | is the voltage magnitude on the bus, i,
N represents a total number of busses in the electric power delivery system, such that n represents bus n, of the total N busses in the electric power delivery system during mathematical summation,
|Y in | is the magnitude of a Y bus matrix element related to the bus, i, and bus n,
θ in is an angle of a Y bus matrix element related to the bus, i, and bus n,
δ i is an angle of the bus, i, and
δ n is an angle of bus n; and
performing a control operation based on the calculated power generation reference value, P i Ref .
2. The method of claim 1 , wherein the data identifying at least one of the operating mode and the control mode of the electric power delivery system identifies an operating mode of the electric power delivery system as either grid-connected or isolated.
3. The method of claim 1 , wherein the data identifying at least one of the operating mode and the control mode of the electric power delivery system identifies an operating mode of the electric power delivery system as isolated and identifies a control mode as at least one of isochronous mode, droop mode, and constant power mode.
4. The method of claim 1 , wherein the data identifying at least one of the operating mode and the control mode of the electric power delivery system identifies the electric power delivery system as at least one grid-connected, isolated with isochronous mode generation, isolated with droop mode generation, and isolated with constant power mode generation.
5. The method of claim 1 , wherein the identified voltage and angle of the power on the bus are used to calculate a Y vector.
6. The method of claim 1 , further comprising generating a modified Newton Raphson expression.
7. The method of claim 1 , wherein calculating the power generation reference value, P i Ref , is performed by summing the calculated power characteristic-dependent load value, the calculated control-operating value, and the shared value, S v .
8. The method of claim 1 , wherein the power characteristic-dependent load value is calculated by:
determining a percentage of the load on the bus, i, that varies quadratically with respect to the obtained power characteristic;
determining a percentage of the load on the bus, i, that varies linearly with respect to the obtained power characteristic;
determining a percentage of the load on the bus, i, that is power characteristic independent; and
multiplying the active power load on the bus, i, at a nominal power characteristic multiplied by a function of (i) the percentage of the load on the bus, i, that varies quadratically with respect to the obtained power characteristic, (ii) the percentage of the load on the bus, i, that varies linearly with respect to the obtained power characteristic, and (iii) the percentage of the load on the bus, i, that is power characteristic independent.
9. The method of claim 8 , wherein the power characteristic-dependent load value comprises a voltage-dependent load value.
10. The method of claim 8 , wherein the power characteristic-dependent load value comprises a frequency-dependent load value.
11. The method of claim 10 , wherein the frequency-dependent load value is calculated to satisfy the expression:
( c 1i ƒ pu 2 +c 2i ƒ pu +c 3i ) P Load_i N ,
where
c 1i represents the percentage of the load on the bus, i, that varies quadratically with respect to the obtained frequency value,
c 2i represents the percentage of the load on the bus, i, that varies linearly with respect to the obtained frequency value,
c 3i represents the percentage of the load on the bus, i, that is frequency independent,
ƒ pu represents the frequency value in per unit of the power on the bus, i, and
P Load_i N represents the active power value of the load on the bus, i, at nominal frequency.
12. The method of claim 10 , wherein
the frequency-based load value is calculated to satisfy the expression:
( c 1i ƒ pu 2 +c 2i ƒ pu +c 3i ) P Load_i N ,
where
c 1i represents a percentage of a load on the bus, i, that varies quadratically with respect to the measured frequency,
c 2i represents a percentage of the load on the bus, i, that varies linearly with respect to the measured frequency,
c 3i represents a percentage of the load on the bus, i, that is frequency independent,
ƒ pu represents the frequency value in per unit of the power on the bus, i, and
P Load_i N represents the active power value of the load on the bus, i, at nominal frequency; and
such that active the power generation reference value, P i Ref , is calculated to satisfy the expression:
P i Ref =( c 1i ƒ pu 2 +c 2i ƒ pu +c 3i ) P Load_i N +β i X+|V i |Σ n=1 N |Y in ∥V n |cos(δ i −δ n −θ in ).
13. A method for evaluating a power flow on a bus, i, in at least a portion of an electric power delivery system to perform a control operation, comprising:
obtaining power characteristic values of power on the bus, i, of the electric power delivery system, including at least a frequency value and a voltage value of the power on the bus, i;
calculating a power characteristic-dependent load value by:
determining a percentage of the load on the bus, i, that varies quadratically with respect to the obtained power characteristic;
determining a percentage of the load on the bus, i, that varies linearly with respect to the obtained power characteristic;
determining a percentage of the load on the bus, i, that is power characteristic independent; and
multiplying the active power load on the bus, i, at a nominal power characteristic multiplied by a function of (i) the percentage of the load on the bus, i, that varies quadratically with respect to the obtained power characteristic, (ii) the percentage of the load on the bus, i, that varies linearly with respect to the obtained power characteristic, and (iii) the percentage of the load on the bus, i, that is power characteristic independent;
obtaining data identifying at least one of an operating mode and a control mode of the electric power delivery system;
calculating a control-operating value based on a product of a control-mode coefficient, β i , and an operating value, X, of at least one generator connected to the bus, i;
obtaining data identifying an angle of the power on the bus, i;
calculating a power generation reference value, P i Ref , based on:
a sum of the calculated power characteristic-dependent load value, the calculated control-operating value, and the obtained voltage value and
the angle of the power on the bus, I; and
performing a control operation based on the calculated power generation reference value, P i Ref .
14. The method of claim 13 , wherein the power characteristic-dependent load value comprises a voltage-dependent load value.
15. The method of claim 13 , wherein the power characteristic-dependent load value comprises a frequency-dependent load value.
16. The method of claim 15 , wherein the frequency-dependent load value is calculated to satisfy the expression:
( c 1i ƒ pu 2 +c 2i ƒ pu +c 3i ) P Load_i N
where:
c 1i represents the percentage of the load on the bus, i, that varies quadratically with respect to the obtained frequency value,
c 2i represents the percentage of the load on the bus, i, that varies linearly with respect to the obtained frequency value,
c 3i represents the percentage of the load on the bus, i, that is frequency independent,
ƒ pu represents the frequency value in per unit of the power on the bus, i, and
P Load_i N represents the active power value of the load on the bus, i, at nominal frequency.
17. The method of claim 15 , wherein:
the frequency-based load value is calculated to satisfy the expression:
( c 1i ƒ pu 2 +c 2i ƒ pu +c 3i ) P Load_i N
where:
c 1i represents a percentage of a load on the bus, i, that varies quadratically with respect to the measured frequency,
c 2i represents a percentage of the load on the bus, i, that varies linearly with respect to the measured frequency,
c 3i represents a percentage of the load on the bus, i, that is frequency independent, and,
ƒ pu represents the frequency value in per unit of the power on the bus, i, and
P Load_i N represents the active power value of the load on the bus, i, at nominal frequency.
18. The method of claim 17 , wherein the identified voltage and angle of the power on the bus are used to calculate a shared value, wherein the shared value S v is calculated to satisfy the expression:
S v =|V i |Σ n=1 N |Y in ∥V n |cos(δ i −δ n −θ in )
where:
|V i | is the voltage magnitude on the bus, i,
N represents a total number of busses in the electric power delivery system, such that n represents bus n, of the total N busses in the electric power delivery system during mathematical summation,
|Y in | is the magnitude of a Y bus matrix element related to the bus, i, and bus n,
θ in is an angle of a Y bus matrix element related to the bus, i, and bus n,
δ i is an angle of the bus, i, and
δ n is an angle of bus n.
19. The method of claim 18 , wherein the shared value is calculated such that the power generation reference value, P i Ref , is calculated to satisfy the expression:
P i Ref =( c 1i ƒ pu 2 +c 2i ƒ pu +c 3i ) P Load i N +β i X+|V i |Σ n=1 N |Y in ∥V n |cos(δ i −δ n −θ in ).