CALCULATION APPARATUS AND CALCULATION METHOD

Information

  • Patent Application
  • 20220393465
  • Publication Number
    20220393465
  • Date Filed
    November 06, 2020
    4 years ago
  • Date Published
    December 08, 2022
    2 years ago
Abstract
A calculation device 20 is linked with an electric power system 1 and optimizes energy utilization efficiency of a grid S1 including a power storage device 15. The calculation device 20 determines the range of a final charge state of a prediction target section of the power storage device 15 on the basis of power demand prediction for the next and following sections, and assesses, through the optimization calculation for the prediction target section, whether or not the final charge state of the section falls within the range.
Description
BACKGROUND
Technical Field

The present invention relates to a technique for optimizing energy use efficiency of a grid.


Description of Related Art

Introduction of a distributed power supply typified by a photovoltaic (PV) power generation system has been promoted from the viewpoint of reduction of dependence on a fossil fuel and an environmental problem. A PV system converts power generated by a solar power generation panel from a direct current (DC) to an alternating current (AC) using an inverter circuit of a power control device, and outputs the AC.


Japanese Patent No. 4765162 discloses a technique for controlling second power conversion means such that received power detected by received power detection means does not fall below predetermined power during output of both the generated power of the distributed power supply and the power from power storage means in the power storage type solar power generation system. Thus, the power from the power storage means is prevented from flowing back to a power system.


BRIEF SUMMARY

When an energy storage apparatus has sufficient reserve power, the grid including the energy storage apparatus can absorb excess or deficiency of the power by charge-discharge to compensate for the excess or deficiency. The reserve power is, when a surplus of power is charged and stored, a remaining capacity that can be charged. When a shortage of power is compensated for by the discharge, the reserve power is a capacity that can be discharged.


In the case of performing a calculation for optimizing the energy use efficiency of the grid, when the state of charge of the energy storage apparatus is determined only by the supply and demand prediction of the power in the section, sometimes the energy storage apparatus cannot secure the reserve power, which should absorb the excess or deficiency of the power caused in the next section or later, at the start of the next section.


In a calculation apparatus for optimizing energy use efficiency of a grid that is interconnected to a power system and includes an energy storage apparatus, a range of a final charge state in a prediction target section of the energy storage apparatus is determined based on a supply and demand prediction of power in a next section or later, and whether the final charge state of the section is matched with the range is evaluated by the optimization calculation of the prediction target section.


The above aspect can be applied to a calculation method for optimizing the energy use efficiency of the grid interconnected to the power system and including the energy storage apparatus. The present invention can be applied to a program performing the optimization calculation.


At the start of the next section, it is possible to secure the reserve power, which should absorb the excess or deficiency of the power caused in the next section or later, for the energy storage apparatus.





BRIEF DESCRIPTION OF THE FIGURES


FIG. 1 is a block diagram of a microgrid.



FIG. 2 is a block diagram of the microgrid.



FIG. 3 is a graph illustrating a correlation characteristic between a solar radiation amount and generated power.



FIG. 4 is a graph illustrating transitions of a prediction value of generated power and a prediction value of power consumption.



FIG. 5 is a simplified block diagram of the microgrid.



FIG. 6 is a graph illustrating a transition of a state of charge of a storage battery.



FIG. 7 is a graph illustrating the transition of the state of charge of the storage battery.



FIG. 8A is a view illustrating the number of candidates of received power.



FIG. 8B is a view illustrating the number of candidates of the received power.



FIG. 9 is a view illustrating an example of division of a prediction target period.



FIG. 10 is a graph illustrating a supply and demand amount of power in each section.



FIG. 11 is a flowchart of calculation processing.



FIG. 12 is an enlarged view of a part of FIG. 10.



FIG. 13 is a block diagram of the microgrid.





DETAILED DESCRIPTION OF VARIOUS EMBODIMENTS

In a calculation apparatus for optimizing energy use efficiency of a grid that is interconnected to a power system and includes an energy storage apparatus, a range of a final charge state in a prediction target section of the energy storage apparatus is determined based on a supply and demand prediction of power in a next section or later, and whether the final charge state of the section is matched with the range is evaluated by the optimization calculation of the prediction target section.


In this configuration, the range of the final charge state in the section is determined based on the supply and demand prediction of the power after the next section. Because the final charge state of the section is the starting charge state of the next section, the final charge state of the section is matched with the determined range, so that the reserve power, which should absorb the excess or deficiency of the power caused in the next section or later, can be secured for the energy storage apparatus at the start of the next section.


A surplus power amount in a first period, in which supply of the power exceeds demand, in the next section or later may be calculated based on the supply and demand prediction of the power in the next section or later, and an upper limit value of the final charge state in the section may be determined to be a value obtained by subtracting a numerical value corresponding to the surplus power amount from a use range upper limit of the energy storage apparatus. In this configuration, at the starting time point of the next section, the reserve power (free capacity) for charging the surplus power of the grid can be secured in the energy storage apparatus. Thus, energy loss can be prevented.


A minimum surplus power amount in the first period may be calculated based on the supply and demand prediction of the power in the next section or later, and the upper limit value of the final charge state in the section may be determined to be a value obtained by subtracting a numerical value corresponding to the minimum surplus power amount from the use range upper limit of the energy storage apparatus. In this configuration, at the starting time point of the next section, the reserve power (free capacity) for charging at least the minimum surplus power amount can be secured in the energy storage apparatus.


A deficient power amount, by which power supply falls below the demand to cause power shortage, may be calculated for a second period from a start of the next section to the first period, and a lower limit value of the final charge state in the section may be determined to be a value obtained by adding a numerical value corresponding to the deficient power amount to a use range lower limit of the energy storage apparatus. In this configuration, at the starting time point of the next section, the reserve power for discharging the shortage power of the grid can be secured in the energy storage apparatus. Therefore, the energy can be efficiently used.


A maximum deficient power amount in the second period may be calculated based on the supply and demand prediction of the power in the next section or later, and the lower limit value of the final charge state in the section may be determined to be a value obtained by adding a numerical value corresponding to the maximum deficient power amount to the use range lower limit of the energy storage apparatus. In this configuration, at the starting time point of the next section, the reserve power for discharging the maximum deficient power amount can be secured in the energy storage apparatus.


First Embodiment

1. Description of Microgrid S1


A microgrid S is a small-scale power system interconnected to a power system 1, and includes at least a distributed power supply, an energy storage apparatus, and a load. The power system 1 may be a power system of a power provider or an independent power system constructed by a self-sustained operation output of a large power conditioner.



FIG. 1 is a block diagram of the microgrid S1. The microgrid S1 includes a solar power generation panel 10 that is the distributed power supply, a storage battery 15 that is the energy storage apparatus, a power conditioner 20 that is the power control device, and a load L.


The power conditioner 20 includes a first converter circuit 21, a second converter circuit 23, a DC link unit 25, a bidirectional inverter circuit 31, a relay 37, a control device 50, a DC voltage detection unit 27, an output current detection unit 33, and an output voltage detection unit 35.


The solar power generation panel 10 is connected to the first converter circuit 21. The first converter circuit 21 is a DC/DC converter, and boosts and outputs an output voltage (DC) of the solar power generation panel 10. The first converter circuit 21 may be a chopper.


The storage battery 15 is connected to the second converter circuit 23. For example, the storage battery 15 is a secondary battery. The second converter circuit 23 is a bidirectional DC/DC converter that discharges and charges the storage battery 15. The second converter circuit 23 may be a bidirectional chopper.


The solar power generation panel 10 and the storage battery 15 are connected in parallel to the DC link unit 25 through the first converter circuit 21 and the second converter circuit 23.


The DC link unit 25 is located between a connection point 24 of the converter circuits and the bidirectional inverter circuit 31. An electrolytic capacitor C1 is provided in the DC link unit 25. The electrolytic capacitor C1 is provided so as to stabilize voltage Vdc of the DC link unit 25.


The DC voltage detection unit 27 detects the voltage Vdc of the DC link unit 25. The voltage Vdc of the DC link unit 25 detected by the DC voltage detection unit 27 is input to the control device 50.


The bidirectional inverter circuit 31 is a bidirectional conversion circuit that selectively performs an inverse conversion (inverter) converting the DC into the AC and a forward conversion (converter) converting the AC into the DC. The bidirectional inverter circuit 31 is connected to the DC link unit 25, and converts DC power input from the DC link unit 25 into AC power and outputs the AC power during reverse conversion operation. Particularly, power corresponding to voltage increased from a reference value in the DC link unit 25 by power generation of the solar power generation panel 10 is input to the bidirectional inverter circuit 31. Consequently, the power corresponding to the voltage increased from the reference value is converted from the DC to the AC and output from the bidirectional inverter circuit 31.


The storage battery 15 can store surplus power of the solar power generation panel 10 through the second converter circuit 23. When the generated power of the solar power generation panel 10 is insufficient, the storage battery 15 can compensate for the shortage of the generated power by the discharge through the second converter circuit 23.


The bidirectional inverter circuit 31 is connected to the power system 1 using the system power supply 2 as an AC power supply through a relay 37.


The relay 37 is installed for interconnection with the power system 1. The microgrid S1 can be interconnected to the power system 1 by closing the relay 37.


The output current detection unit 33 detects output current Iinv of the bidirectional inverter circuit 31. The output voltage detection unit 35 is located on an output side of the bidirectional inverter circuit 31, and detects output voltage Vinv of the bidirectional inverter circuit 31.


The output current Iinv of the bidirectional inverter circuit 31 detected by the output current detection unit 33 and the output voltage Vinv of the bidirectional inverter circuit 31 detected by the output voltage detection unit 35 are input to the control device 50. The control device 50 calculates output power (active power) Pinv of the bidirectional inverter circuit 31 based on the output current Iinv and the output voltage Vinv of the bidirectional inverter circuit 31. The output power Pinv is “positive” during the reverse conversion and “negative” during the forward conversion.


The load L that is a demand facility is connected to a power line (trunk line) 5 connecting the bidirectional inverter circuit 31 and the power system 1 through a branch line 4. The power can be supplied to the load L from both the power conditioner 20 and the power system 1.


A power receiving point 3 is a power supply point by the power system 1, and is a boundary portion between the power system 1 and premises (microgrid S1) as illustrated in FIG. 1.


An external measuring instrument 40 such as an external transducer is provided in the power system 1 as an instrument detecting the power at the power receiving point 3.


The external measuring instrument 40 includes a received current detection unit 41 and a system voltage detection unit 43. The external measuring instrument 40 is installed corresponding to the power receiving point 3, and the received current detection unit 41 detects a received current at the power receiving point 3. The system voltage detection unit 43 detects the system voltage of the power system 1.


The external measuring instrument 40 calculates the received power (effective power) PRCV based on the received current and the system voltage. The received power PRCV detected by the external measuring instrument 40 is input to the control device 50. The received power PRCV can be used to determine the state of a power flow (hereinafter, simply referred to as a flow). The external measuring instrument 40 is a measuring instrument that measures the received power PRCV of the power receiving point 3.


In the received power PRCV, a forward flow (FIG. 1: the flow of electricity from the power system 1 to the microgrid S1) is set to “positive”, and a reverse flow (FIG. 2: the flow of electricity from the microgrid S1 to the power system 1) is set to “negative”.


The control device 50 includes a CPU 51 and a memory 53. The memory 53 stores a program predicting supply and demand of the power of the microgrid S1 and a program performing the optimization calculation optimizing the energy use efficiency of the microgrid S1. The control device 50 is an example of the “calculation apparatus” of the present invention.


The control device 50 can give a command to the bidirectional inverter circuit 31 to control switching between the forward conversion operation and the reverse conversion operation. The output of the bidirectional inverter circuit 31, namely, the output power Pinv can be controlled. The output power Pinv can be controlled by adjusting the output current Iinv.


The control device 50 turns on and off the first converter circuit 21 to control connection and non-connection of the solar power generation panel 10 to and from the DC link unit 25. The connection and disconnection of the storage battery 15 with respect to the DC link unit 25 can be controlled by turning on and off the second converter circuit 23. The control device 50 can control switching between charge and discharge of the storage battery 15 through the second converter circuit 23. The charge current and the discharge current of the storage battery 15 can be controlled through the second converter circuit 23.


2. Power Generation Prediction and Load Prediction of Microgrid S1


The control device 50 predicts the generated power [kW] of the solar power generation panel 10 and the power consumption [kW] of the load L.


Generated power PPV of the solar power generation panel 10 has a correlation with a solar radiation amount X as illustrated in FIG. 3. The prediction value of the generated power PPV can be obtained from a prediction data provider 70 through a network NW. The prediction data provider 70 may be a provider by a provider of the power conditioner 20 or a provider of a power generation company. The generated power PPV of the solar power generation panel 10 is converted into the AC by the bidirectional inverter circuit 31. Conversion efficiency η is conversion efficiency of the inverter circuit 31 when the DC is converted into the AC.


The power consumption PLOAD of the load L can be predicted from past data. For example, the power consumption of the next day can be predicted by statistically processing the data of the power consumption PLOAD for several days.


The power consumption PLOAD of the load L can be obtained from the received power PRCV of the power receiving point 3 and the output power Pinv of the bidirectional inverter circuit 31. In the case of the forward current (PRCV>0), the power consumption PLOAD of the load L is a sum of the output power Pinv and the received power PRCV. In the case of the reverse power flow (PRCV<0), the power consumption PLOAD of the load L is a difference between the output power Pinv and the received power PRCV.






P
LOAD=Pinv+PRCV   (A)



FIG. 4 is a graph illustrating results of power generation prediction and load prediction of the microgrid S1. A broken line is the prediction value of the generated power of the microgrid S1, and a solid line is the prediction value of the power consumption of the microgrid S1. In this example, the power generation prediction and the load prediction are performed every hour while setting a prediction target period T as one day, and the prediction value is a stepwise waveform in which the value changes every hour. The power generation prediction and the load prediction of the microgrid S1 may be performed by the control device 50, or the data predicted by another device may be obtained.


3. Optimization Calculation Optimizing Microgrid Energy Efficiency



FIG. 5 is a simplified block diagram of the microgrid S1.


The control device 50 performs the optimization calculation that optimizes the energy use efficiency of the microgrid S1. The optimization calculation is to obtain an optimum value of a target value (hereinafter, a received power target value PRCVref) of the received power PRCV that minimizes an objective function F of Mathematical formula 1 based on the power generation prediction and the load prediction (FIG. 4).


An equation 1 is an objective function evaluating the energy use efficiency of the microgrid S1.









F
=





t

T



{



k
1

·


T
MAX

(


P
RCVreft

,

P
PVt

,

P
LOADt


)


+


k
2

·


T
MIN

(


P
RCVreft

,


P

PVt
,




P
LOADt



)


+


k
3

·


P
RCV

(


P
RCVreft

,

P
PVt

,

P
LOADt


)



}


+


k
4

·

P
PEAK







[

Mathematical


formula


1

]







A first term and a second term of the objective function F are terms evaluating the use-restricted period of the storage battery 15. A third term and a fourth term are terms evaluating the electricity rate of the microgrid S1.


The use-restricted period includes two periods including a charge-restricted period TMAX and a discharge-restricted period TMIN. The charge-restricted period TMAX is a period during which the charge is restricted, and for example, is the case where a SOC of the storage battery 15 is located at the upper limit of the use range (in the case of full charge). The discharge-restricted period TMIN is a period during which the discharge is restricted, and for example, is the case where the SOC of the storage battery 15 is located at a lower limit of the use range.


In this example, k1>k2>k3>k4 is satisfied, and the evaluation of the use-restricted period of the storage battery 15 is prioritized in the evaluation of the energy use efficiency of the microgrid S1. k1=k2 and k3=k4 may be satisfied.


When the arithmetic calculation optimizing a received power target value PRCVref is performed, there are two conditions of an upper limit and a lower limit (an equation 2) of the received power target value PRCVref and an upper limit and a lower limit (an equation 3) of the output power [kW] of the storage battery 15 as constraint conditions.


(Upper and Lower Limit Constraint on Received Power Target Value)





PPCSCAP≤PRCVref t≤PPCSCAP (t∈T)  [Mathematical formula 2]


(Upper and Lower Limit Constraint on Output Power of Storage Battery)





PBATCAP≤(PLOAD t−PRCVref t−PPV t·η)PBATCAP(t∈T)   [Mathematical formula 3]










TABLE 1







T
Set of total time section t (prediction target period)


TMAX(PRCVref t, PPV t, PLOAD t)
Charge-restricted period in time section t


TMIN(PRCVref t, PPV t, PLOAD t)
Discharge-restricted period in time section t


PRCV(PRCVref t, PPV t, PLOAD t)
Received power estimated value in time section t [kW]


PPEAK
Received power peak value of prediction target period T [kW]


k1~k4
Weighting factor


PRCVref t
Received power target value in time section t [kW]


PPV t
Generation power prediction value of PV panel in time section t [kW]


PLOAD t
Power consumption prediction value of load in time section t [kW]


PPCSCAP
Power conditioner rated capacity [kW]


η
Conversion efficiency [%]


PBATCAP
Storage battery rated charge-discharge power [kW]









The SOC (state of charge) is a ratio of a charge amount to rated capacity of the storage battery 15. The SOC is an example of the “charge state” of the storage battery 15. In order to calculate the objective function F, the SOC of the storage battery 15 is required to be estimated. A method of simulating the SOC of the storage battery 15 and the estimated value of received power PRCVt will be described below.


A method for calculating the charge-restricted period TMAX and discharge-restricted period TMIN in a time section t will be described.


<Step 1>


A provisional output power prediction value PBATtmpt of the storage battery 15 at the time section t is calculated from an equation 4 and an equation 5.













0


(


P
LOADt

-

P
RCVreft

-


P
PVt

·
η


)




P
BATCAP












P
BATtmpt

=



(


P
LOADt

-

P
RCVreft


)

η

-

P
PVt








(

during


discharge

)





[

Mathematical


formula


4

]

















-

P
BATCAP




(


P
LOADt

-

P
RCVreft

-


P
PVt

·
η


)

<

0











P
BATtmpt

=



(


P
LOADt

-

P
RCVreft


)

·
η

-

P
Pvt








(

during


charge

)





[

Mathematical


formula


5

]







<Step 2>


From an equation 6 to an equation 9, a provisional state-of-charge prediction value SOCtmpt of the storage battery 15 in the time section t is calculated.


(a) In the case of t=0





WhBAT 0=SOCINI×BATCAP   [Mathematical formula 6]





WhBATtmp t=WhBAT 0−PBATtmp t×hour   [Mathematical formula 7]


(b) In the case of t≠0





WhBATtmp t=WhBAT t-1−PBATtmp t×hour   [Mathematical formula 8]










SOC
tmpt

=



W


h
BATtmpt


BATCAP

×
100





[

Mathematical


formula


9

]







<Step 3>


The presence or absence of deviation from the upper and lower limits the presence or absence of the matching of the upper and lower limits in the time section t is determined, and the state-of-charge prediction value SOCt of the storage battery 15, the output power PBATt of the storage battery 15, the charge-restricted period TMAX, and the discharge-restricted period TMIN are determined.


(a) In the case of SOCtmpt≥SOCMAX






T
MAX(PRCVref t,PPV t,PLOAD t)=1,TMIN(PRCVref t,PPV t,PLOAD t)=0   [Mathematical formula 10]






P
BAT t=(SOCtmp t-1−SOCMAX)×BATCAP   [Mathematical formula 11]


(b) In the case of SOCtmpt≤SOCMIN






T
MAX(PRCVref t,PPV t,PLOAD t)=0,TMIN(PRCVref t,PPV t,PLOAD t)=1   [Mathematical formula 12]






P
BAT t=(SOCtmnp t-1−SOCMIN)×BATCAP  [Mathematical formula 13]


(c) In the case of SOCMIN<SOCtmpt<SOCMAX






T
MAX(PRCVref t,PPV t,PLOAD t)=0,TMIN(PRCVref t,PLOAD t)=0   [Mathematical formula 14]





PBAT t=PBATtmp t   [Mathematical formula 15]


Subsequently, an initial value (in the case of t=0) of the charge amount prediction value WhBATt of the storage battery 15 can be obtained from equations 16, 17.





WhBAT 0=SOCINI×BATCAP   [Mathematical formula 16]





WhBAT t=WhBAT 0−PBAT t×hour   [Mathematical formula 17]


Subsequently, the charge amount prediction value WhBATt of the storage battery 15 can be obtained for each time section (when t≠0) by an equation 18. The prediction value SOCt of the SOC of the storage battery 15 can be obtained from an equation 19.





WhBAT t=WhBAT t-1−PBAT t×hour   [Mathematical formula 18]






[

Mathematical


formula


19

]







SOC
t

=



W


h
BATt


BATCAP

×
100





The received power prediction value PRCV at the time section t can be calculated from equations 20, 21.





(PPV t+PBAT t)≥0 . . . (during power conditioner inverse conversion)






P
RCV(PRCVref t,PPV t,PLOAD t)=PLOAD t−(PPV t+PBAT t)·η  [Mathematical formula 20]











(


P
PVt

+

P
BATt


)

<

0






(

during


power


conditioner


conversion

)








P
RCV

(



P

RCVreft
,




P
PVt


,

P
LOADt


)

=


P
LOADt

-


(


P
PVt

+

P
BATt


)

η







[

Mathematical


formula


21

]







A receiving power peak prediction value PPEAK in the prediction target period T can be calculated from an equation 22.






P
PEAK
=P
RCV(PRCVref t,PPV t, PLOAD t)





PPEAK_old≤PRCV(PRCVref t,PPV t, PLOAD t)  [Mathematical formula 22]

    • (here, at t=0, PPEAK_old=PPEAK INI)










TABLE 2







PBATtmp t
Provisional output power prediction value of storage battery in time section t [kW]


WhBAT t−1
Storage battery charge amount in previous time section t1 [kWh]


WhBATtmp t
Provisional charge amount prediction value of storage battery in time section t [kWh]


SOCINI
Storage battery charge state initial value in prediction target period T [%]


BATCAP
Storage battery rated capacity [kWh]


hour
Pitch width [h]


SOCtmp t
Provisional charge state prediction value of storage battery in time section t [%]


SOCMAX
Storage battery charge state upper limit [%]


SOCMIN
Storage battery charge state lower limit [%]


PBAT t
Storage battery output power prediction value in time section t [kw]


WhBAT t
Storage battery charge amount prediction value in time section t [kwh]


SOCt
Storage battery charge state prediction value in time section t [%]


PPEAKold
Peak power until previous time section t − 1 in prediction target period T [kw]


PPEAK INI
Peak power initial value in prediction target period T [kw]









From equations 4 to 19, when the received power target value PRCVreft is determined with respect to a generated power prediction value PPVt×η and a power consumption prediction value PLOADt of the load L, a SOCt of the storage battery 15 can be estimated.


In each time section t, the control device 50 estimates the SOCt of the storage battery 15 using the received power target value PRCVreft as a variable, and calculates four terms of the objective function F from the result.


Such an operation is performed (twenty-four cycles are executed when the prediction target period T is one day and when the arithmetic cycle is one hour) in the prediction target period T. Then, when the values of the objective function F are compared to each other with respect to the pattern of the combination of the received power target values PRCVreft, the combination of the received power target values minimizing the objective function F, namely, the received power target value PRCVreft of each time section t can be determined for the prediction target period T.



FIGS. 6 and 7 are graphs in which a horizontal axis represents time [h], a left vertical axis represents power [kW], and a right vertical axis represents the SOC [%]. Y1 (bold line) indicates a transition of the received power [kW], and Y2 indicates a transition of the generated power [kW]. Y3 represents a transition of the power consumption [kW] of the load, Y4 represents a transition of the output power [kW] of the storage battery, and Y5 (broken line) represents a transition of the SOC [%] of the storage battery 15.



FIG. 6 illustrates the case where the received power target value PRCVref is determined by an empirical rule, and FIG. 7 illustrates the case where the received power target value PRCVref is determined based on the objective function F.


When the received power target value PRCVref is determined by the empirical rule (FIG. 6), the SOC of the storage battery 15 is maintained around 100 [%] in the time zone of 16:00 to 18:00. For this reason, the storage battery 15 is required to restrict acceptance of the charge in the time zone of 16:00 to 18:00.


When the received power target value PRCVref is determined based on the objective function F (FIG. 7), the SOC of the storage battery 15 changes with a margin with respect to an SOC upper limit line Lim 1 (SOCMAX=100 [%]) and an SOC lower limit line Lim 2 (SOCMIN=10 [%]), and both the charge and the discharge are always possible. That is, when the received power target value PRCVref is compared with the case of the empirical rule, the SOC is suppressed to be low in the time zone (part A) of 16:00 to 18:00, and the point that the charge in this time zone is not limited is improved.


4. Power Control of Microgrid


The control device 50 controls the power of the microgrid S1 such that the received power PRCV of the microgrid S1 follows the received power target value PRCVref calculated using the objective function F.


For example, when a measurement value of the received power PRCV is lower than the received power target value PRCVref during the forward flow in the power receiving point 3 and the reverse conversion operation in the inverter circuit 31, the output power of the storage battery 15 is narrowed to reduce the difference between the measurement value of the received power PRCV and the received power target value PRCVref. When the measurement value of the received power PRCV is higher than the received power target value PRCVref, the difference between the measurement value of the received power PRCV and the received power target value PRCVref is reduced by increasing the output power of the storage battery 15.


As described above, the received power PRCV of the microgrid S1 can follow the received power target value PRCVref calculated using the objective function F by adjusting the output of the storage battery 15 according to the difference between the measurement value of the received power PRCV and the target value.


5. Calculation Load Reduction Method


In order to determine the optimum value of the received power target value PRCVref, the SOC and the like of the storage battery 15 are required to be estimated by the number of candidates of the received power target value PRCVref for each time section t of the prediction target period T, and each term of the objective function F is required to be calculated. For example, in the case where the prediction target period T is one day and the number of time sections is twenty-four, when the received power target value PRCVref is changed in the range of 0 to 10 [kW] in increments of 1 [kW], the total number of candidates of the received power target value PRCVref is 1124=9.85×1024, and the calculation load increases in proportion to the total number of candidates.


The calculation load of the optimization calculation can be reduced using the following two methods.


(A) Constraint on change width of received power target value


(B) Division of prediction target period T


The method (A) will be described.


Because there is a possibility that the change in the received power PRCV affects the power quality of the power system 1, desirably the change width is as small as possible. For this reason, as illustrated in an equation 23, the constraint is provided on the change width of the received power target value PRCVref per time section t.





ΔPRCVref≤(PRCVref t−PRCVref t-1)≤ΔPRCVref(t∈T)   [Mathematical formula 23]

    • <Q PRCVref:Change width maximum value of received power target value [kW]



FIGS. 8A and 8B are views illustrating a combination of candidates of the received power target value PRCVreft. FIG. 8A illustrates the case where the constraint is not provided on the change width of the received power target value PRCVref, and FIG. 8B illustrates the case where the maximum value of the change width of the received power target value PRCVref is set to ±2 [kW].


The number of candidates of the received power target value PRCVref per time section t is “11” when the constraint is not provided on the change width of the received power target value PRCVref (in the case of FIG. 8A). On the other hand, the number of candidates of the received power target value PRCVref per time section t is “5” when the maximum value of the change width of the received power target value PRCVref is set to ±2 [kW] (in the case of FIG. 8B).


In this way, the constraint is provided on the change width of the received power target value PRCVref, so that the number of candidates of the received power target value PRCVref per time section t can be reduced to greatly reduce the calculation load. The influence on the power system 1 can be reduced, and the quality of electricity can be maintained.


The method (B) will be described.


B1. Division of Sections


The prediction target period T is divided into a plurality of sections ΔT. The control device 50 performs the optimization calculation optimizing the energy efficiency of the microgrid S1 using each of the divided sections ΔT as one prediction target period. That is, the optimum value of the received power target value PRCVref of each time cross section t is determined using the objective function F.


For example, when the prediction target period T is one day, the prediction target period T is divided into four sections ΔT1 to ΔT4 as illustrated in FIG. 9. Then, the control device 50 obtains the optimization value of the received power target value PRCVref using the objective function F with each section ΔT as one prediction target period.


In this case, the received power target value PRCVref in the entire prediction target period T is obtained by connecting the received power target values PRCVref of the respective time sections t of the respective sections ΔT determined as the optimum values in time series.


When the prediction target section T is divided, the received power target value PRCVref of the entire prediction target period T is calculated in a plurality of times. Thus, the number of combinations of candidates of the received power target value PRCVref per one time is reduced, so that the calculation load can be greatly reduced.


The number of divisions N of the prediction target period T may be other than four, such as three divisions. The number of time sections of each section ΔT may or may not be the same. Each section ΔT is preferably longer than the cycles of the power generation prediction and the load prediction. In the example of FIG. 9, ΔT=6 hours and the cycle of the power generation prediction and the load prediction=1 hour.


B2. Reserve Power of Storage Battery at Start of Next Section



FIG. 10 is a graph illustrating the supply and demand prediction of the power for the prediction target section T. The supply and demand of the power can be predicted from the prediction value of the generated power of the solar power generation panel 10, the prediction value of the power consumption of the load L, and the received power target value from the power system 1. In FIG. 10, Y1 represents a received power target value (minimum), Y2 represents a generated power prediction value, and Y3 represents a prediction value of the power consumption.


In the case where the prediction target section T is divided to perform the optimization calculation for each section ΔT, when the calculation is performed in consideration only of the supply and demand prediction in the section ΔTn, sometimes the reserve power that should absorb the excess or deficiency of the power in the next section ΔTn+1 cannot be secured in the storage battery 15 at the start of the next section ΔTn+1. For this reason, there is a problem in how to secure the reserve power of the storage battery 15 in anticipation of the supply and demand prediction of the power in the next section ΔTn+1 at the final time point of the section ΔTn. The reserve power is an empty capacity that can be charged when a surplus of power is charged and stored. When a shortage of power is compensated for by the discharge, the capacity is a capacity that can be discharged.


In order to solve the above problem, the control device 50 of the power conditioner 20 calculates the range of the final SOC of the storage battery 15 in the section ΔTn based on the supply and demand prediction of the power in next section ΔTn+1.


Because the final SOC of the section ΔTn=the start SOC of the next section ΔTn+1, the reserve power of the storage battery 15 can be secured at the starting time point of the section ΔTn+1 by matching the final SOC of the storage battery 15 of the section ΔTn with the calculated range. That is, the storage battery 15 can be operated with a margin with respect to the use range upper limit SOCMAX and the use range lower limit SOCMIN, and the excess or deficiency of the power can be compensated for by being absorbed by the charge-discharge. FIG. 11 is a flowchart of the calculation processing for calculating the SOC range.


(S10) Calculation of Upper Limit Value SOCUPPER of final SOC


The control device 50 calculates the upper limit value SOCUPPER of the final SOC of the storage battery 15 in the section ΔTn from the surplus power amount WC in the next section ΔTn+1.


The surplus power amount WC is a value obtained by integrating the surplus power in the first period (A1 in FIG. 10) in which the difference between the generated power prediction value PPVt of the solar power generation panel 10 that is the distributed power supply and the power consumption prediction value PLOADt of the load L becomes continuously positive in the next section ΔTn+1.


The surplus power is a value obtained by adding the difference between the power consumption PLOAD and the received power PRCV to the generated power PPV. The surplus power depends on the magnitude of the received power PRCV supplied from the power system 1, and is smaller as the received power PRCV is smaller.


When the storage battery 15 does not have the reserve power capable of charging at least the minimum surplus power amount WC, the energy loss is caused. For this reason, in this example, the minimum surplus power amount WC in the next section ΔTn+1 is calculated based on the generated power prediction value PPVt, the power consumption prediction value PLOADt, and the minimum received power target value PRCVrefmin t.


A minimum received power target value PRCVrefmin t in the next section ΔTn+1 can be obtained by subtracting the maximum value of the change width set in (A) for each time section from a final received power prediction value PRCV tend in the section ΔTn.


The surplus power amount WC may be calculated over a plurality of sections such as ΔTn+1 and ΔTn+2 (as an example, there may be 6 hours to 8 hours).


The upper limit value SOCUPPER of the final SOC of the storage battery 15 in the section ΔTn is obtained by subtracting the SOC corresponding to the surplus power amount WC from the use range upper limit SOCMAX of the storage battery 15 such that the free capacity that can charge the surplus power amount WC in the first period A1 can be secured. The surplus power amount WC may be a minimum value or other values.


In FIG. 12, “Lim 1” indicates the use range upper limit SOCMAX, and “Lim 3” indicates the upper limit SOCUPPER of the final SOC. “Y5” is the SOC of the storage battery 15. The use range upper limit SOCMAX is a limit value (a value that cannot be charged any more) on the charge side of the SOC.














In


the


case


of



P
PVt


-

P
LAODt



0





SOC
UPPER

=

SOCMAX
-





t


Δ


T

n
+
1






{


P
PVt

-

(


P
PVt

-


(


P
LAODt

-

P
RCVrefmint


)

·
η




}


BATCAP







Here
,


SOC
UPPER

:

P
RCVvefmint

:
upper


limit


value


of


final


soc


of


storage


battery


in


divided


section


Δ



T
n

[
%
]



minimum


received


power


target


value


of


time


section


t


in


in


divided


section


Δ



T

n
+
1


[
kW
]








[

Mathematical


formula


24

]







The microgrid S1 is not only in the case where the power is excessive but also in the case where the power is insufficient due to a supply and demand balance of the power. When the final SOC of the storage battery 15 is in the vicinity of the use range lower limit SOCMIN, sometimes there is the case where it is difficult to cope with the discharge from the section start to the time section tcross at which the power supply becomes excessive in the next section. Therefore, the control device 50 calculates a lower limit value SOCLOWER of the final SOC so as to leave the remaining power that can be discharged.


(S20) Calculation of Lower Limit Value SOCLOWER of Final SOC


The lower limit value SOCLOWER of the final SOC of the storage battery 15 in the section ΔTn is obtained by calculating the deficient power amount WD for the second period A2 from the start of the next section ΔTn+1 to the start of the first period A1 (to the time section tcross in which the supply is excessive). The deficient power amount WD is a value obtained by integrating the deficient power in the second period A2.


The deficient power is a value obtained by subtracting the generated power PPV and the received power PRCV from the power consumption PLOAD. The deficient power depends on the magnitude of the received power PRCV supplied from the power system, and becomes larger as the received power PRCV becomes smaller.


When the storage battery 15 has the reserve power with which the maximum deficient power amount WD can be discharged, the minimum received power can be maintained, and the energy use efficiency is high. Consequently, in this example, the maximum deficient power amount WD is calculated for the second period A2 of the next section ΔTn+1 based on the generated power prediction value PPVt, the power consumption prediction value PLOADt of the load L, and the minimum received power target value PRCVrefmint.


The reverse power flow is generated when the discharge amount of the storage battery 15 exceeds the deficient power amount WD, so that it is not allowable to exceed the deficient power amount WD because the reverse power flow from the storage battery 15 to the system is generated.


A value obtained by adding the SOC corresponding to the deficient power amount WD to the use range lower limit SOCMIN of the storage battery 15 so as to secure the capacity capable of discharging the deficient power amount WD is the lower limit value SOCLOWER of the final SOC of the section ΔTn. The deficient power amount WD may be a maximum value or other values. In FIG. 12, “Lim 2” indicates the use range lower limit SOCMIN, and “Lim 4” indicates the lower limit value SOCLOWER of the final SOC. Y5 is the SOC of the storage battery 15. The use range lower limit SOCMIN is a limit value (a value that cannot be discharged any more) on the discharge side of the SOC.














In


the


case


of



P
LAODt


-

P
PVt



0





SOC
LOWER

=

SOCMIN
+









t
=
0


t
=

t
cross




(


P
LAODt

-











P
PVt

·
η

-

P
RCVrefmint


)





BATCAP
·
η




(

t


Δ


T

n
+
1




)








Here
,


t
cross

:
excessive


supply


reaching


time


section


of


dividing


section







SOC
LOWER

:
lower


limit


value


of


final


soc


of


storage


battery


of


divided


section


Δ



T
n

[
%
]







[

Mathematical


formula


25

]







(S30) Penalty During Deviation


As described above, a range Z (see FIG. 12) having the lower limit of SOCLOWER and the upper limit of SOCUPPER can be obtained for the final SOC of the storage battery 15 in the section ΔTn.


When the final SOC of the storage battery 15 deviates from the range Z, namely, when the final SOC is not matched with the range Z, the final SOC is calculated as a penalty. When the final SOC does not deviates from the range Z, namely, when the final SOC is matched with the range Z, there is no penalty.





[Mathematical formula 26]





A) SOCtend−SOCUPPER≥0





PENALTYOVERn=SOCtend−SOCUPPER   (26A)





B) SOCLOWER−SOCtend≥0





PENALTYSHORTn=SOCLOWER−SOCtend  (26B)





C) SOCLOWER<SOCtend<SOCUPPER





PENALTYOVERn=PENALTYSHORTn=0   (26C)

    • Here, SOCtend: final soc of storage battery of divided section ΔTn [%]
    • PENALTYOVERn: upper limit value deviation penalty OD final SOC of storage battery of divided section ΔTn
    • PENALTYSHORTn: lower limit value deviation penalty OD final SOC of storage battery of divided section ΔTn


The equation 27 is an objective function E obtained by adding an evaluation term of the penalty to four items with respect to the objective function F. The optimization value of the received power target value PRCVref is calculated using the objective function E with the section ΔT as one prediction target period.


When the optimization calculation is performed for each section ΔTn, the final SOC of the storage battery 15 in the section ΔTn is set to the initial SOCINI of the storage battery 15 in the next division section ΔTn+1, and the steps of S10 to S30 are repeatedly performed by the number of sections. Then, at the completion of the entire section, it is possible to obtain an optimal solution of the received power target value PRCVref that optimizes the energy use efficiency of the microgrid S1 for the prediction target period T.


A solution having no penalty or a small penalty is obtained as the optimal solution by adding the penalty evaluation term to the objective function E, so that it is possible to prevent that the final SOC of the section ΔTn is not matched with the range Z.












[

Objective


function

]




E
=





n

N






t


Δ


T
n





{



k
1

·


T
MAX

(


P

RCVreft
n


,

P

PVt
n


,

P

LOADt
n



)


+


k
2

·


T
MIN

(


P

RCVreft
n


,

P

PVt
n


,

P

LOADt
n



)


+


k
3

·


P
RCV

(


P

RCVreft
n


,

P

PVt
n


,

P

LOADt
n



)


+

(



k
5

·

PENALTY

OVER
n



+


k
6

·

PENALTY

SHORT
n




)


}



+


k
4

·

P
PEAK








Here
,

k
5

,


k
6

:
weight


coefficient







[

Mathematical


formula


27

]







The magnitude relationship between the coefficients k5, k6 that determine a weight of the penalty and other coefficients may be k1>k2>k3>k4>k5>k6. k5>k6>k1>k2>k3>k4 may be satisfied. The magnitude relationship between the two coefficients may be k5>k6 or k5<k6.


6. Effects


The storage battery 15 can charge the surplus power of the microgrid S1 and can discharge the surplus power during the shortage of the power to compensate it. When the use of the storage battery 15 is restricted, the surplus power cannot be stored and the energy is lost, so that the energy use efficiency of the microgrid S1 decreases.


In this configuration, the use-restricted period of the storage battery 15 can be minimized by obtaining the optimum value of the received power target value PRCVreft based on the objective function E or the objective function F. Thus, the energy use efficiency of the microgrid S1 can be optimized from the viewpoint of the use efficiency of the storage battery 15.


In this configuration, the prediction target period T is divided into a plurality of sections ΔT to perform the optimization calculation optimizing the energy use efficiency of the microgrid S1. Consequently, the calculation load can be reduced, and the time for calculating the optimal solution can be shortened.


In this configuration, the range Z of the final SOC in the section ΔTn is determined based on the supply and demand prediction of the power in and after the next section ΔTn+1. Because the final SOC of the section ΔT is the start SOC of the next section ΔTn+1, the storage battery 15 can secure the reserve power that should absorb the excess or deficiency of the power in and after the next section ΔTn+1 at the start of the next section ΔTn+1 by matching the final SOC of the section ΔTn with the range Z.


In this configuration, the surplus power amount WC in the first period A1 in which the supply of power exceeds the demand is calculated based on the supply and demand prediction of the power in and after the next section ΔTn+1. The upper limit value SOCUPPER of the final SOC in the section ΔTn of the storage battery 15 is determined to be a value obtained by subtracting the SOC corresponding to the surplus power amount WC from the use range upper limit SOCMAX. In this configuration, at the starting time point of the next section ΔTn+1, the reserve power charging the surplus power of the microgrid S1 can be secured in the storage battery 15. Thus, energy loss can be prevented.


In this configuration, the deficient power amount WD is calculated for the second period A2 from the start of the next section ΔTn+1 to the first period A1. The lower limit value SOCLOWER of the final SOC in the section ΔTn of the storage battery 15 is determined to be a value obtained by adding the SOC corresponding to the deficient power amount WD to the use range lower limit SOCMIN. In this configuration, at the starting time point of the next section ΔTn+1, the reserve power discharging the deficient power of the microgrid S1 can be secured in the storage battery 15. Therefore, the energy can be efficiently used.


Other Embodiments

The present invention is not limited to the embodiments described with reference to the above description and the drawings. For example, the following embodiments are also included in the technical scope of the present invention.


(1) In the first embodiment, the grid having the linear power line (trunk line) 5 is illustrated as an example of the microgrid S1. Alternatively, a grid having a circular power line (trunk line) may be used. A microgrid S2 in FIG. 13 includes a power line 100 having a circular shape. A solar power generation panel 110 and a wind power generator 120 are connected to the power line 100 through power converters 115, 125. A load 130 and a storage battery 140 are connected to the power line 100. The power line 100 of the microgrid S2 is connected to the power system 1 through the interconnection line 105.


The microgrid S2 includes a control device 150. The control device 150 divides the prediction target period T into the plurality of sections to perform the optimization calculation optimizing the energy use efficiency of the microgrid S2, thereby calculating the optimal solution of the received power target value PRCVref. The control device 150 performs power control of the microgrid S2 such that the received power of the power receiving point 3 follows the calculated received power target value PRCVref. Specifically, the received power PRCV of the power receiving point 3 is monitored based on the output of the meter 160 provided at the power receiving point 3. When the difference exists with respect to the received power target value PRCVref, the control device 150 charges or discharges the storage battery 140 through the power converter 145 to reduce the difference. In this way, the received power PRCV can follow the calculated target value PRCVref, and the energy use efficiency of the microgrid S2 can be optimized. The control device 150 is an example of the “calculation apparatus” of the present invention.


The distributed power supply is a generic term for all small-scale power generation facilities dispersedly disposed adjacent to a demand site. For example, the distributed power supply may be a biomass power generator or the like in addition to the solar power generation panel 10 and the wind power generator 120. The distributed power supply may be a power source using renewable energy or a power source using fossil fuel.


(2) In the first embodiment, the calculation optimizing the energy use efficiency of the microgrid S1 is performed using the received power target value as the variable. The calculation optimizing the energy use efficiency of the microgrid may be performed using the SOC of the storage battery 15 as the variable. The present technique can be applied not only to the microgrid (small-scale power system) but also to the grid (power system) that is interconnected to another power system and includes the energy storage apparatus. In the first embodiment, the storage battery 15 is exemplified as the energy storage apparatus. Alternatively, the energy storage apparatus may be a capacitor or the like.


(3) In the first embodiment, the objective function F includes four terms of the first term to the fourth term, and the four terms are multiplied by the weighting coefficients k1 to k4 and added. In the first embodiment, the magnitude relationship of the weighting factors is k1>k2>k3>k4. Alternatively, the magnitude relationship may be k3>k4>k1>k2. The electricity rate of the microgrid can be preferentially evaluated by reversing the magnitude relationship of the weighting factors, so that the energy use efficiency of the microgrid can be optimized from the viewpoint of the electricity rate. The objective function F may be only the first and second terms evaluating the use-restricted period of the storage battery 15, or may be only the third and fourth terms evaluating the electricity rate of the microgrid. The objective function F may be any other function as long as the function optimizes the energy utilization efficiency of the microgrid. The same applies to the objective function E.


(4) In the first embodiment, the optimization calculation optimizing the energy use efficiency of the microgrid is performed by the control device 50 of the power conditioner 20. A communication subject is not limited to the power conditioner 20, but may be an arithmetic device provided separately from the power conditioner 20.


(5) In the first embodiment, the calculation optimizing the energy use efficiency of the microgrid has been described using the objective functions E, F. Alternatively, the calculation may be performed using, for example, AI.


(6) In the first embodiment, the upper limit value SOCUPPER and the lower limit value SOCLOWER of the final SOC in the section ΔTn are obtained. For example, in the next section ΔTn+1, when the second period A2 in which the power is insufficient does not exist before the first period A1, the upper limit value SOCUPPER of the final SOC may be obtained based on the surplus power amount WC in the first period A1, and the lower limit value of the final SOC may be the use range lower limit SOCMIN. On the contrary, in the next section ΔTn+1, when only the period in which the power is insufficient exists, the lower limit value SOCLOWER of the final SOC may be obtained based on the deficient power amount WD, and the upper limit value of the final SOC may be the use range upper limit SOCMAX.


(7) In the first embodiment, the SOC [%] of the storage battery 15 was evaluated in the calculation optimizing the energy use efficiency of the microgrid. The residual capacity [Ah] of the storage battery 15 may be evaluated instead of the SOC [%]. That is, the range of the final residual capacity in the prediction target section ΔTn of the storage battery 15 may be determined based on the supply and demand prediction of the power after the next section, and whether the final residual capacity in the section ΔTn is matched with the determined range may be evaluated in the optimization calculation of the prediction target section ΔTn. The SOC and the residual capacity are examples of the “state of charge” of the storage battery.


(8) In the first embodiment, the prediction target period T is divided into a plurality of sections in order to reduce the load of optimization calculation. The present technique is not limited to dividing or not dividing the prediction target period, but can be applied to the case where the optimization calculation is performed for an arbitrary prediction target section.

Claims
  • 1. A calculation apparatus for optimizing energy use efficiency of a grid that is interconnected to a power system and includes an energy storage apparatus, wherein a range of a final charge state in a prediction target section of the energy storage apparatus is determined based on a supply and demand prediction of power in a next section or later, and whether the final charge state of the section is matched with the range is evaluated by the optimization calculation of the prediction target section.
  • 2. The calculation apparatus according to claim 1, wherein a surplus power amount in a first period, in which supply of the power exceeds demand in the next section or later, is calculated based on the supply and demand prediction of the power in the next section or later, andan upper limit value of the final charge state in the section is determined to be a value obtained by subtracting a numerical value corresponding to the surplus power amount from a use range upper limit of the energy storage apparatus.
  • 3. The calculation apparatus according to claim 2, wherein: a minimum surplus power amount in the first period is calculated based on the supply and demand prediction of the power in the next section or later, andthe upper limit value of the final charge state in the section is determined to be a value obtained by subtracting a numerical value corresponding to the minimum surplus power amount from the use range upper limit of the energy storage apparatus.
  • 4. The calculation apparatus according to claim 2, wherein: a deficient power amount, by which power supply falls below the demand to cause power shortage, is calculated for a second period from a start of the next section to the first period, anda lower limit value of the final charge state in the section is determined to be a value obtained by adding a numerical value corresponding to the deficient power amount to a use range lower limit of the energy storage apparatus.
  • 5. The calculation apparatus according to claim 4, wherein: a maximum deficient power amount in the second period is calculated based on the supply and demand prediction of the power in the next section or later, andthe lower limit value of the final charge state in the section is determined to be a value obtained by adding a numerical value corresponding to the maximum deficient power amount to the use range lower limit of the energy storage apparatus.
  • 6. A calculation method for optimizing energy use efficiency of a grid that is interconnected to a power system and includes an energy storage apparatus, the calculation method comprising: determining a range of a final charge state in a prediction target section of the energy storage apparatus based on a supply and demand prediction of power in a next section or later; andevaluating whether the final charge state of the section is matched with the range by optimization calculation of the prediction target section.
  • 7. The calculation apparatus according to claim 3, wherein: a deficient power amount, by which power supply falls below the demand to cause power shortage, is calculated for a second period from a start of the next section to the first period, anda lower limit value of the final charge state in the section is determined to be a value obtained by adding a numerical value corresponding to the deficient power amount to a use range lower limit of the energy storage apparatus.
  • 8. The calculation apparatus according to claim 7, wherein a maximum deficient power amount in the second period is calculated based on the supply and demand prediction of the power in the next section or later, andthe lower limit value of the final charge state in the section is determined to be a value obtained by adding a numerical value corresponding to the maximum deficient power amount to the use range lower limit of the energy storage apparatus.
Priority Claims (2)
Number Date Country Kind
2019-206995 Nov 2019 JP national
2020-177870 Oct 2020 JP national
CROSS REFERENCE TO RELATED APPLICATIONS

This application is a national stage application, filed under 35 U.S.C. § 371, of International Application No. PCT/JP2020/041446, filed Nov. 6, 2020, which international application claims priority to and the benefit of Japanese Application Nos._2019-206995, filed Nov. 15, 2019, and 2020-177870, filed Oct. 23, 2020; the contents of all of which as are hereby incorporated by reference in their entireties.

PCT Information
Filing Document Filing Date Country Kind
PCT/JP2020/041446 11/6/2020 WO