MICROGRID SYSTEM COMPRISING ENERGY MANAGEMENT SYSTEM OF ENERGY STORAGE SYSTEM (ESS)-CONNECTED PHOTOVOLTAIC POWER SYSTEM

Abstract

An improved EMS (Energy Management System) of ESS (Energy Storage System)—connected photovoltaic power system is provided, where the economic efficiency of the microgrid power transaction is maximized by minimizing the amount paid to the power system as a result of optimal operation as to the energy supply and demand in the process of transacting power surplus/shortage with the power system, the responsiveness to passive resource energy forecasting of supply and demand is improved by resolving the uncertainty of solar power generation forecasting and load forecasting, the deterioration of the available storage capacity of ESS is minimized, and contribution to solving the nation's power supply shortage is made by the operation based on the detailed identification of energy storage status of ESS.

Claims

1. A microgrid system comprising an energy generation (EG) system; one or more electrical load coupled to the EG system; an ESS (energy storage system) coupled to the EG system and the electrical load; an EMS (energy management system) for managing energy of the microgrid including the EG, the one or more electrical load, the ESS, and power transaction between the microgrid and a system power (grid); wherein the EMS comprises: a first module which forecasts power supply and demand, does operation scheduling, and controls ESS that stores EG electricity; a second module for checking and managing state of charge of the ESS; a third module for calculating and managing the discharge rate of the ESS; a fourth module for system power flattening control to reduce a difference between maximum value and minimum value of system power by checking the load, the EG system generation amount, charge/discharge rate of the ESS, and the system power (grid) power usage by the microgrid; a fifth module for system power smoothing control reducing a variation in the usage of system power by time periods; a sixth module for controlling power smoothing of the ESS; a seventh module for controlling the peak of power usage; an eighth module for controlling net zero energy operation so that power demand and supply are balanced; a ninth module for controlling a power demand response based on the amount of ESS discharge; and a control unit for controlling each module, wherein controlling operations of the modules are conducted so as to minimize the amount paid to the grid.

2. The microgrid system of claim 1 wherein the controlling operations of the modules are further conducted to reflect a difference between a maximum system power (grid) usage and a minimum system power (grid).

3. The microgrid system of claim 1 wherein the controlling operations of the modules are further conducted to reflect the peak power usage.

4. The microgrid system of claim 1 wherein the controlling operation of the modules are further conducted to reflect power smoothing of the ESS.

5. The microgrid system of claim 1 wherein the controlling operation of the modules are further conducted to reflect the net zero energy operation.

6. The microgrid system of claim 1 wherein the controlling operation of the modules are further conducted to reflect the power demand response based on the amount of ESS discharge.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0014] FIG. 1 shows a system diagram of a microgrid system according to certain embodiment.

[0015] FIG. 2 is an example that shows how system power usage and ESS charging/discharging power change depending on the external conditions (peak control, net zero operation, flattening and demand response, etc.) applied.

[0016] FIGS. 3-5 are exemplary views showing examples wherein information to be forecasted and monitored by an EMS of ESS-connected photovoltaic power system.

DESCRIPTION OF EXEMPLARY EMBODIMENTS

[0017] Hereinafter, preferred embodiments according to the present invention will be described in more detail with reference to the accompanying drawings. The present invention is related to a microgrid (MG) 400 which comprises an energy storage system (ESS) 410, and an energy management system (EMS) 100. First, a microgrid (MG) is a system composed of distributed power sources or MG's, various power loads, and measuring instruments, and is operated independently or in connection with the power grid of a power company, and is classified as shown in Table 1.

TABLE-US-00001 TABLE 1 Classification Application Characteristic Linked type MG for Customer-owned Electricity transaction + self-generation .fwdarw. electricity MG customers facilities rate↓, reliability↑ Always connected operation with distribution system Uninterrupted supply by switching to single operation in case of emergency Composed of one or more MGs in campus, hospital, building, military base, etc. Example) Smart E-Campus MG for Power company Power line operation efficiency + self-generation .fwdarw. distribution (By D/L section) purchase cost↓, reliability↑ Renewable E+ battery connection, operating all or part of the distribution system as MG .fwdarw. Power line operation efficiency, increasing new and renewable energy e.g.) Linked MG (Shinan), Canadian distribution-class MG project Standalone Energy Small and Self-generation + MG EMS .fwdarw. fuel cost↓, MG Independent Medium Sized quality/reliability↑ MG Islands Inverter controls voltage/frequency: improves power (Electricity quality Projects) High new and renewable energy introduction rate (50% or more) Ex) Kasa Island energy independence island Hybrid Large Sized Self-generation + MG EMS .fwdarw. fuel cost↓, MG Islands, Remote quality/reliability↑ Areas (Electricity Diesel generator controls voltage/frequency Projects) Diesel generator + renewable energy (less than 50%) Ex) Geochado, Ulleungdo, Deokjeokdo, etc.

[0018] FIG. 1 shows a system diagram of a microgrid system according to certain embodiment. Here, solid line arrows indicate the supply of power, and dotted line arrows the flow of information. Although both energy management system (EMS) 100 and data storage 200 belong to the Microgrid 400, the Microgrid Sub-system 430, in this case, is one that consists of new and renewable energy sources or EG sources, loads and energy storage system (ESS) 410. There are a number of new and renewable energy sources such as photovoltaic (PV) and wind turbine (WT) but only PV 420 was indicated for convenience. The power produced by PV will be stored in the ESS or consumed by the load. Also, PV/Load data are used to predict the future PV/load value after being stored in data storage. The Power Grid (Grid 300) supply electricity to the ESS 410 or load 430, whereas the unit cost of power is provided at the power exchange 500. The event server 600 assumes the role of notifying the situation wherein DR (Demand Response) or net zero operation is required. There will be no information about the external operating conditions from the event server in a scenario that does not include any special conditions; otherwise, the constraints and objective function will be changed after receiving external operating condition information. And the conditions may include peak control, power use flattening, power demand response and operation of net zero energy, etc.

[0019] FIG. 2 is a diagram that shows, as an example, how system power usage and ESS 410 charging/discharging power change depending on the external conditions applied. The picture on the upper left is a basic setting and the rest of the pictures in order of bottom left, bottom center, upper left and bottom right show the change when peak control, net zero operation, flattening and demand response have been applied, respectively.

[0020] An ESS 410 is a device for improving efficiency in terms of energy use or cost reduction by storing electricity when it is produced by renewable energy generation (solar power generation) or when the power cost in the system is low and supplying it when necessary.

[0021] An EMS 100 refers to a system that can monitor, analyze and control remotely and in real time by linking sensors, measurement equipment, and analysis S/W. An EMS like this has four main functions: energy monitoring, energy management, analysis and statistics, and energy control. Energy monitoring is to monitor energy usage and real-time loads, and energy management does various management function including such as energy operation policy, energy operation efficiency targeting, overload reserve, and operation improvement. In the analysis and statistics function, reports submitted by the government and management of reduction projects are analyzed, and finally, in energy control, peak management (load remote control) and EMS operation data are controlled. In other words, EMS collects and analyzes energy-related data (data provided by a demand manager or public information/utility). Then, based on the analyzed data, an energy operation schedule in consideration of a saving method, etc. is derived and transmitted to the control side, and control management is performed based on the derived schedule.

[0022] The present invention achieves utilizes these functions and achieves objectives such as maximizing the economic efficiency of the microgrid power transaction by minimizing the amount paid to the power system as a result of optimal operation as to the energy supply and demand in the process of transacting power surplus/shortage with the power system, improving the responsiveness to passive resource energy forecasting of supply and demand by resolving the uncertainty of solar power generation forecasting and load forecasting, minimizing the deterioration of the available storage capacity of ESS 410, and contribution to solving the nation's power supply shortage by the operation based on the detailed identification of energy storage status of ESS, thereby providing improved EMS 100 of ESS-connected photovoltaic power system. To this end, the present invention provides an energy management system equipped with an algorithm capable of optimizing the microgrid energy demand and supply, and also with visualization thereof.

[0023] The EMS 100 according to the present invention comprises: a first module which forecasts power supply and demand, does operation scheduling, and controls ESS that stores photovoltaic (PV) electricity; a second module for checking and managing state of charge of the ESS; a third module for calculating and managing the discharge rate of the ESS; a fourth module for system power leveling control to reduce a difference between maximum value and minimum value of system power by checking a load, photovoltaic power generation amount, charge/discharge rate of ESS, and system power usage by the microgrid; a fifth module for system power smoothing control reducing a variation in the usage of system power by time slots; a sixth module for controlling power smoothing of the ESS; a seventh module for controlling the peak of power usage; an eighth module for controlling net zero energy operation so that power demand and supply are balanced; a ninth module for controlling a power demand response based on the amount of ESS discharge; and a control unit for controlling each module. Here, each module is an element having an operation processing function, and can be calculated and processed in the following manner.

[First Module]

[0024] 1) The following equation was used to express the maximum value of p.sup.d.sup.b(i) which varies according to the situation. p.sup.d.sup.b(i) is the discharge rate of ESS in the i-th time period, which can offset the amount purchased from the grid. p.sub.d,min(i)=p.sub.d(i)−Δp.sub.d(i) is the minimum (expected) value of total demand, Δp.sub.d(i) is the expected error of the total demand, p.sub.ESS.sup.dis,max(i) is the maximum discharge rate of the ESS in the i-th time period.

P.sub.max.sup.d.sup.b(i)=min(max(0,p.sub.d,min(i)),p.sub.ESS.sup.dis,max(i))

[0025] 2) The following equation was used to express the maximum value of p.sup.c.sup.s(i) which varies according to the situation. p.sup.c.sup.s(i) is the charging rate of ESS in the i-th time period, which can offset the amount sold to the grid. p.sub.d,max(i)=p.sub.d(i)+Δp.sub.d(i) is the maximum (expected) value of total demand, p.sub.ESS.sup.chg,max(i) is the maximum charging rate of the ESS in the i-th time period.

P.sub.max.sup.c.sup.s(i)=min(max(0,−p.sub.d,max(i)),p.sub.ESS.sup.chg,max(i))

[0026] 3) The following equation was used to express the maximum value of p.sup.d.sup.u(i) which varies according to the situation. p.sup.d.sup.u(i) is the discharge rate of ESS in the i-th time period, which is the amount to cope with uncertainty.

P.sub.max.sup.d.sup.u(i)=min(max(0,p.sub.d,max(i))−max(0,p.sub.d,min(i)),p.sub.ESS.sup.dis,max(i))

[0027] 4) The following equation was used to express the maximum value of p.sup.c.sup.u(i) which varies according to the situation. p.sup.c.sup.u(i) is the charging rate of ESS in the i-th time period, which is the amount to cope with uncertainty.

P.sub.max.sup.c.sup.u(i)=min(min(0,p.sub.d,max(i))−min(0,p.sub.d,min(i)),p.sub.ESS.sup.chg,max(i))

[0028] 5) The following expression is used to express the range of p.sup.d.sup.b(i) that varies according to the situation.

P.sub.max.sup.d.sup.b(i)(1−δ.sub.bs2(i))≤p.sup.d.sup.b(i)≤p.sub.max.sup.d.sup.b(i) [0029] i) 0≤p.sup.d.sup.b(i)≤P.sub.max.sup.d.sup.b(i) if δ.sub.bs2(i)=1 [0030] ii) p.sup.d.sup.b(i)=P.sub.max.sup.d.sup.b(i) if δ.sub.bs2(i)=0

[0031] 0≤δ.sub.bs1(i)≤1, 0≤δ.sub.bs2(i)≤1

This condition is that δ.sub.bs1(i), δ.sub.bs2(i) each has a value between 0 and 1, and it is a condition that each has a value of 0 or 1 because each of the variable is integer. δ.sub.bs1(i) is a binary variable indicating whether there is purchase of power form the grid, and δ.sub.bs2(i) a binary variable indicating whether there is a sales of power to the grid. (1: yes, 0: no).

[0032] 6) The following expression is used to express the range of p.sup.d.sup.s(i) that varies according to the situation. p.sup.d.sup.s(i) is the discharge rate of ESS in the i-th time period, which is the amount sold to the grid.

p.sup.d.sup.s(i)≤p.sub.ESS.sup.dis,max(i)(1−δ.sub.bs2(i)) [0033] i) p.sup.d.sup.s(i)≤0 if δ.sub.bs2(i)=1 [0034] ii) p.sup.d.sup.s(i)≤p.sub.ESS.sup.dis,max(i) if δ.sub.bs2(i)=0

[0035] 7) The following expression is used to express the range of p.sup.c.sup.s(i) which varies according to the situation.

P.sub.max.sup.c.sup.s(i)δ.sub.bs1(i)≤p.sup.c.sup.s(i)≤P.sub.max.sup.c.sup.s(i) [0036] i) p.sup.c.sup.s(i)=P.sub.max.sup.c.sup.s(i) if δ.sub.bs1(i)=1 [0037] ii) 0≤p.sup.c.sup.s(i)≤P.sub.max.sup.c.sup.s(i) if δ.sub.bs1(i)=0

[0038] 8) The following expression is used to express the range of p.sup.c.sup.b(i) that varies according to the situation. p.sup.c.sup.b(i) is the charging rate of ESS in the i-th time period, which is the amount purchased from the grid.

p.sup.c.sup.b(i)≤p.sub.ESS.sup.chg,max(i)δ.sub.bs1(i) [0039] i) p.sup.c.sup.b(i)≤p.sub.ESS.sup.chg,max(i) if δ.sub.bs1(i)=1 [0040] ii) p.sup.c.sup.b(i)≤0 if δ.sub.bs1(i)=0

[0041] 9) The following expression is used to express the range of p.sup.d.sup.b(i) that varies according to the situation.

p.sub.max.sup.d.sup.b(i)(1−δ.sub.bs1(i))≤p.sup.d.sup.b(i) [0042] i) 0≤p.sup.d.sup.b(i) if (δ.sub.bs1(i),δ.sub.bs2(i))=(1,1), (δ.sub.bs1(i),δ.sub.bs2(i))=(1,0) [0043] ii) P.sub.max.sup.d.sup.b(i)≤p.sup.d.sup.b(i) if (δ.sub.bs1(i),δ.sub.bs2(i))=(0,1), (δ.sub.bs1(i),δ.sub.bs2(i))=(0,0)

[0044] 10) The following expression is used to express the range of p.sup.d.sup.b(i) that varies according to the situation.

P.sub.max.sup.d.sup.u(i)(1−δ.sub.bs2(i))≤p.sup.d.sup.u(i)≤P.sub.max.sup.d.sup.u(i)(1−δ.sub.bs1(i)) [0045] i) p.sup.d.sup.u(i)=0 if (δ.sub.bs1(i),δ.sub.bs2(i)) [0046] ii) P.sub.max.sup.d.sup.u(i)≤p.sup.d.sup.u(i)≤0 if (δ.sub.bs1(i),δ.sub.bs2(i))=(1,0) [0047] iii) 0≤p.sup.d.sup.u(i)≤P.sub.max.sup.d.sup.u(i) if (δ.sub.bs1(i),δ.sub.bs2(i))=(0,1) [0048] iv) p.sup.d.sup.u(i)=P.sub.max.sup.d.sup.u(i) if (δ.sub.bs1(i),δ.sub.bs2(i))=(0,0)

[0049] 11) The following expression is used to express the range of p.sup.c.sup.s(i) that varies according to the situation.

P.sub.max.sup.c.sup.s(i)δ.sub.bs2(i)≤p.sup.c.sup.s(i) [0050] i) P.sub.max.sup.c.sup.s(i)≤p.sup.c.sup.s(i) if (δ.sub.bs1(i),δ.sub.bs2(i))=(1,1) [0051] ii) 0≤p.sup.c.sup.s(i) if (δ.sub.bs1(i),δ.sub.bs2(i))=(1,0) [0052] iii) P.sub.max.sup.c.sup.s(i)≤p.sup.c.sup.s(i) if (δ.sub.bs1(i),δ.sub.bs2(i))=(0,1) [0053] iv) 0≤p.sup.c.sup.s(i) if (δ.sub.bs1(i),δ.sub.bs2(i))=(0,0)

[0054] 12) The following expression is used to express the range of p.sup.c.sup.u(i) that varies according to the situation.

P.sub.max.sup.c.sup.u(i)δ.sub.bs1(i)≤p.sup.c.sup.u(i)≤P.sub.max.sup.c.sup.u(i)δ.sub.bs2(i) [0055] i) p.sup.c.sup.u(i)=P.sub.max.sup.c.sup.u(i) if (δ.sub.bs1(i),δ.sub.bs2(i))=(1,1) [0056] ii) P.sub.max.sup.c.sup.u(i)≤p.sup.c.sup.u(i)≤0 if (δ.sub.bs1(i),δ.sub.bs2(i))=(1,0) [0057] iii) 0≤p.sup.c.sup.u(i)≤P.sub.max.sup.c.sup.u(i) if (δ.sub.bs1(i),δ.sub.bs2(i))=(0,1) [0058] iv) p.sup.c.sup.u(i)=0 if (δ.sub.bs1(i),δ.sub.bs2(i))=(0,0)

[0059] 13) The following expression is used to reflect the range of p.sup.d.sup.u(i) in the objective function. w.sub.d(i) is a variable introduced to cope with the uncertainty in the i-th time period, and z.sub.d is a variable introduced to cope with to the uncertainty bearing no physical meaning. c.sub.buy(i) is the price of purchasing a unit amount of electricity from the grid in the i-th time period, c.sub.sell(i) is the price at which the unit amount of electricity is sold to the grid in the i-th time period, and dt is the unit time.

[00001]$0.5.Math.c_{buy}\left(i\right)-c_{sell}\left(i\right).Math.p^{d_u}\left(i\right)dt\le z_d+w_d\left(i\right)i)0\le z_d+w_d\left(i\right)\mathrm{if}.Math.c_{buy}\left(i\right)-c_{sell}\left(i\right).Math.=0\mathrm{ii})p^{d_u}\left(i\right)\le \frac{z_d+w_d\left(i\right)}{0.5.Math.c_{buy}\left(i\right)-c_{sell}\left(i\right).Math.\mathrm{dt}}\mathrm{if}.Math.c_{buy}\left(i\right)-c_{sell}\left(i\right).Math.\ne 0$

[0060] 14) The following expression is used to reflect the range of p.sup.c.sup.u(i) in the objective function. w.sub.c(i) is a variable introduced to cope with the uncertainty in the i-th time period, and z.sub.c is a variable introduced to cope with to the uncertainty bearing no physical meaning.

[00002]$\mathrm{0.5`}.Math.c_{buy}\left(i\right)-c_{sell}\left(i\right).Math.p^{c_u}\left(i\right)dt\le z_c+w_c\left(i\right)i)0\le z_c+w_c\left(i\right)\mathrm{if}.Math.c_{buy}\left(i\right)-c_{sell}\left(i\right).Math.=0\mathrm{ii})p^{c_u}\left(i\right)\le \frac{z_c+w_c\left(i\right)}{0.5.Math.c_{buy}\left(i\right)-c_{sell}\left(i\right).Math.\mathrm{dt}}\mathrm{if}.Math.c_{buy}\left(i\right)-c_{sell}\left(i\right).Math.\ne 0$

[0061] 15) Set the objective function as J, and set it so that J tends to decrease according as the situation fits the purpose.

J≤t

It is a condition that J has a value oft or less. Under this condition, The objective is to find a controlling status which minimizes J. J is minimized by minimizing t, where t is a variable needed to derive an operation schedule to cope with uncertainty.

[00003]$J=.Math.\left(\begin{array}{c}[c_{\mathrm{buy}}\left(i\right)\left\{p^{c_b}\left(i\right)-p^{d_b}\left(i\right)\right\}-c_{\mathrm{sell}}\left(i\right)\left\{p^{d_s}\left(i\right)-p^{c_s}\left(i\right)\right\}+\\ 0.5\left\{c_{\mathrm{sell}}\left(i\right)+c_{\mathrm{buy}}\left(i\right)\right\}\left\{p^{c_u}\left(i\right)-p^{d_u}\left(i\right)\right\}]\mathrm{dt}+w_c\left(i\right)+w_d\left(i\right))\end{array}\right)+\Gamma z_c+\Gamma z_d$

When the customer charges the storage device, the value of J increases, and when the storage device is discharged, the value of J decreases. The underlined part is to cope with the uncertainty, and the rest that is not underlined is: (electricity cost)−(saving/additional gain).

0≤δ.sub.bs1(i)≤1, 0≤δ.sub.bs2(i)≤1

This condition is that δ.sub.bs1(i), δ.sub.bs2(i) each has a value between 0 and 1, and it is a condition that each has a value of 0 or 1 because each of the variable is integer. δ.sub.bs1(i) is a binary variable indicating whether there is purchase of power form the grid, and δ.sub.bs2(i) a binary variable indicating whether there is a sales of power to the grid. (1: yes, 0: no).

z.sub.d≥0, z.sub.c≥0, w.sub.d(i)≥0, w.sub.c(i)≥0

This condition means that each of the variables z.sub.d, z.sub.c, w.sub.d(i), w.sub.c(i) has a value greater than or equal to zero.

[Second Module]

[0062] 1) When δ.sub.dis(i) and δ.sub.chg(i) are added, the condition that it has a value of 0 or more is expressed by the following expression. δ.sub.dis(i) is a binary variable indicating whether the energy storage device is discharged in the i-th time period, and δ.sub.chg(i) is a binary variable indicating whether the energy storage device is charged in the i-th time period.

−δ.sub.dis(i)−δ.sub.chg(i)≤0 (δ.sub.dis(i)+δ.sub.chg(i)≥0)

[0063] 2) When δ.sub.dis(i) and δ.sub.chg(i) are added, the condition that it has a value of 1 or less is expressed by the following expression.

δ.sub.dis(i)+δ.sub.chg(i)≤1

[0064] 3) The minimum value of the ESS discharge rate is expressed using the following expression. (p.sup.d.sup.b(i)+p.sup.d.sup.s(i)+p.sup.d.sup.u(i)) is the discharge rate of the ESS.

P.sub.ESS.sup.dis,min(i).Math.δ.sub.dis(i)≤p.sup.d.sup.b(i)+p.sup.d.sup.s(i)+p.sup.d.sup.u(i)

[0065] 4) The maximum value of the ESS discharge rate is expressed using the following expression. P.sub.ESS.sup.dis,max(i) is the maximum value of the ESS discharge rate in the i-th time period.

p.sup.d.sup.b(i)+p.sup.d.sup.s(i)+p.sup.d.sup.u(i)≤P.sub.ESS.sup.dis,max(i).Math.δ.sub.dis(i)

[0066] 5) The minimum value of the Ess charge rate is expressed using the following expression. (p.sup.c.sup.s(i)+p.sup.c.sup.b(i)+p.sup.c.sup.u(i)) is the charging rate of the ESS.

p.sub.ESS.sup.chg,min(i).Math.δ.sub.chg(i)≤p.sup.c.sup.s(i)+p.sup.c.sup.b(i)+p.sup.c.sup.u(i)

[0067] 6) The maximum value of the Ess charge rate is expressed using the following expression. p.sub.ESS.sup.chg,max(i) is the maximum value of the Ess charge rate in the i-th time period.

p.sup.c.sup.s(i)+p.sup.c.sup.b(i)+p.sup.c.sup.u(i)≤p.sub.ESS.sup.chg,max(i).Math.δ.sub.chg(i)

[0068] 7) 0≤δ.sub.dis(i)≤1, 0≤δ.sub.chg(i)≤1

[0069] It is the condition that each of δ.sub.dis(i) and δ.sub.chg(i) has values between 0 and 1. Since each of the variables is integer variables, as a result, each will have a value of 0 or 1.

[0070] 8) p.sup.d.sup.b(i)≥0, p.sup.d.sup.s(i)≥0, p.sup.d.sup.u(i)≥0, p.sup.c.sup.s(i)≥0, p.sup.c.sup.b(i)≥0, p.sup.c.sup.u(i)≥0

[0071] It is a condition that each of the variables p.sup.d.sup.b(i), p.sup.d.sup.s(i), p.sup.d.sup.u(i), p.sup.c.sup.s(i), p.sup.c.sup.b(i), p.sup.c.sup.u(i) has a value greater than or equal to zero.

[Third Module]

[0072] 1) The sum of the SOC increased by charging from the first time period to the i-th time period can be expressed by the following expression. η.sub.chg is the charging efficiency of the ESS, and E.sub.cap is the capacity of the ESS.

[00004]$\underset{k=1}{\overset{i}{.Math.}}\frac{{\eta }_{chg}dt}{E_{cap}}\left(p^{c_s}\left(k\right)+p^{c_b}\left(k\right)+p^{c_u}\left(k\right)\right)$

The sum of SOC reduced by discharge from the first time period to the i-th time period can be expressed by the following expression. η.sub.dis is the discharge efficiency of the ESS.

[00005]$-\stackrel{i}{\underset{k=1}{.Math.}}\frac{dt}{{\eta }_{\mathrm{dis}}{E}_{cap}}\left(p^{d_b}\left(k\right)+p^{d_s}\left(k\right)+p^{d_u}\left(k\right)\right)$

The sum of the SOC reduced by natural discharge from the first time period to the i-th time period can be expressed by the following equation. P.sub.ESS.sup.self(i) is the natural discharge rate of the ESS.

[00006]$-\stackrel{i}{\underset{k=1}{.Math.}}\frac{dt}{E_{cap}}{p}_{ESS}^{self}\left(k\right)$

Using these, the condition that the state of charge determined by charging and discharging up to the i-th time period is greater than or equal to the minimum value of the state of charge in the i+1 th time period can be expressed by the following expression. SOC(i) is the state of charge in the i-th time period, and SOC.sub.min(i) is the minimum value of the state of charge in the i-th time period.

[00007]${\mathrm{SOC}}_{\min }\left(i+1\right)\le \mathrm{SOC}\left(1\right)+\underset{k=1}{\overset{i}{.Math.}}\frac{{\eta }_{chg}dt}{E_{\mathrm{cap}}}\left(p^{c_s}\left(k\right)+p^{c_b}\left(k\right)+p^{c_u}\left(k\right)\right)-\stackrel{i}{\underset{k=1}{.Math.}}\frac{dt}{{\eta }_{\mathrm{dis}}{E}_{\mathrm{cap}}}\left(p^{d_b}\left(k\right)+p^{d_s}\left(k\right)+p^{d_u}\left(k\right)\right)-\stackrel{i}{\underset{k=1}{.Math.}}\frac{dt}{E_{\mathrm{cap}}}{p}_{ESS}^{\mathrm{self}}\left(k\right)$

[0073] 2) The sum of the SOC increased by charging from the first time period to the i-th time period can be expressed by the following expression.

[00008]$\underset{k=1}{\overset{i}{.Math.}}\frac{{\eta }_{chg}dt}{E_{cap}}\left(p^{c_s}\left(k\right)+p^{c_b}\left(k\right)+p^{c_u}\left(k\right)\right)$

[0074] The sum of SOC reduced by discharge from the first time period to the i-th time period can be expressed by the following expression.

[00009]$-\stackrel{i}{\underset{k=1}{.Math.}}\frac{dt}{{\eta }_{\mathrm{dis}}{E}_{cap}}\left(p^{d_b}\left(k\right)+p^{d_s}\left(k\right)+p^{d_u}\left(k\right)\right)$

[0075] The sum of SOC reduced by natural discharge from the first time period to the i-th time period can be expressed by the following expression.

[00010]$-\stackrel{i}{\underset{k=1}{.Math.}}\frac{dt}{E_{cap}}{p}_{ESS}^{\mathrm{self}}\left(k\right)$

[0076] Therefore, the condition that the state of charge determined by charging and discharging up to the i-th time period is less than or equal to the maximum value of the state of charge in the i+l-th time period can be expressed by the following expression. SOC.sub.max(i) is the maximum value of the state of charge in the i-th time period.

[00011]$\mathrm{SOC}\left(1\right)+\stackrel{i}{\underset{k=1}{.Math.}}\frac{{\eta }_{chg}dt}{E_{\mathrm{cap}}}\left(p^{c_s}\left(k\right)+p^{c_b}\left(k\right)+p^{c_u}\left(k\right)\right)-\stackrel{i}{\underset{k=1}{.Math.}}\frac{dt}{{\eta }_{\mathrm{dis}}{E}_{\mathrm{cap}}}\left(p^{d_b}\left(k\right)+p^{d_s}\left(k\right)+p^{d_u}\left(k\right)\right)-\stackrel{i}{\underset{k=1}{.Math.}}\frac{dt}{E_{\mathrm{cap}}}{p}_{ESS}^{\mathrm{self}}\left(k\right)\le {\mathrm{SOC}}_{\max }\left(i+1\right)$

[Fourth Module]

[0077] 1) The condition that the power demand and power supply are balanced can be expressed by the following equation:

Load−PV−ESS−Grid=NetDemand−ESS−Grid=0

[0078] where Load=load [0079] PV=electricity generated by the photovoltaic system [0080] ESS=ESS charge/discharge rate [0081] Grid=consumption of grid power [0082] NetDemand=Total Demand
Therefore, NetDemand−ESS=Grid. Here, for NetDemandMin which is the minimum value of the total demand, the following expression is used which is equivalent to p.sub.g.sup.min≤NetDemandMin−ESS.

p.sub.g.sup.min≤p.sub.d,min(i)+p.sup.c.sup.b(i)+p.sup.c.sup.u(i)+p.sup.c.sup.s(i)−p.sup.d.sup.b(i)−p.sup.d.sup.u(i)−p.sup.d.sup.s(i)

[0083] 2) The condition that the power demand and power supply are balanced can be expressed by the following equation:

Load−PV−ESS−Grid=NetDemand−ESS−Grid=0

[0084] where Load=load [0085] PV=electricity generated by the photovoltaic system [0086] ESS=ESS charge/discharge rate [0087] Grid=consumption of grid power [0088] NetDemand=Total Demand
Therefore, NetDemand−ESS=Grid. Here, for NetDemandMax which is the maximum value of the total demand, the following expression is used which is equivalent to NetDemandMax−ESS≤p.sub.g.sup.max.

p.sub.d,max(i)+p.sup.c.sup.b(i)+p.sup.c.sup.u(i)+p.sup.c.sup.s(i)−p.sup.d.sup.b(i)−p.sup.d.sup.u(i)−p.sup.d.sup.s(i)≤p.sub.g.sup.max [0089] 3) c.sub.g.sup.flat(p.sub.g.sup.max−p.sub.g.sup.min)T
By adding this equation to J, the difference between the maximum grid power value and the minimum grid power value (p.sub.g.sup.max−p.sub.g.sup.min) is reflected in the objective function. To reduce the value of J, the value of is (p.sub.g.sup.max−p.sub.g.sup.min) reduced, so the use of system power is flattened. If the value of the penalty constant c.sub.g.sup.flat for flattening of grid power is larger, the effect of flattening is reflected more. T is the control period.

[Fifth Module]

[0090] 1)

p.sub.g.sup.dev(i)≤[p.sub.d,min(i)−p.sup.d.sup.b(i)−p.sup.d.sup.u(i)−p.sup.d.sup.s(i)+p.sup.c.sup.b(i)+p.sup.c.sup.u(i)+p.sup.c.sup.s(i)]−[p.sub.d,max(i−1)−p.sup.d.sup.b(i−1)−p.sup.d.sup.u(i−1)−p.sup.c.sup.b(i−1)+p.sup.c.sup.b(i−1)+p.sup.c.sup.u(i−1)+p.sup.c.sup.s(i−1)],

[0091] If i=1,

−p.sub.g.sup.dev(i)≤[p.sub.d,min(i)−p.sup.d.sup.b(i)−p.sup.d.sup.u(i)−p.sup.d.sup.s(i)+p.sup.c.sup.b(i)+p.sup.c.sup.u(i)+p.sup.c.sup.s(i)]−[p.sub.g.sup.buy,prev−p.sub.g.sup.sell,prev]

The condition that the power demand and power supply are balanced can be expressed by the following equation:

Load−PV−ESS−Grid=NetDemand−ESS−Grid=0

[0092] where Load=load [0093] PV=electricity generated by the photovoltaic system [0094] ESS=ESS charge/discharge rate [0095] Grid=consumption of grid power [0096] NetDemand=Total Demand

[0097] Therefore, NetDemand−ESS=Grid.

[0098] Utilizing this, [p.sub.d,min(i)−p.sup.d.sup.b(i)−p.sup.d.sup.u(i)−p.sup.d.sup.s(i)+p.sup.c.sup.b(i)+p.sup.c.sup.u(i)+p.sup.c.sup.s(i)] is the minimum value of the grid power consumption in the i-th time period. Also, [p.sub.d,max(i−1)−p.sup.d.sup.b(i−1)−p.sup.d.sup.u(i−1)−p.sup.d.sup.s(i−1)+p.sup.c.sup.b(i−1)+p.sup.c.sup.u(i−1)+p.sup.c.sup.s(i−1)] is the maximum value of the grid power usage in the i−1th time period. Therefore, the given expression is a condition that (minimum value of grid power consumption in the i-th time period)−(maximum value of grid power consumption in the i−1 th time period) has a value greater than or equal to −p.sub.g.sup.dev(i). When i=1, use [p.sub.g.sup.buy,prev−p.sub.g.sup.sell,prev] instead of (maximum value of grid power consumption in the i−1 th time period).

[0099] 2)

[p.sub.d,max(i)−p.sup.d.sup.b(i)−p.sup.d.sup.u(i)−p.sup.d.sup.s(i)+p.sup.c.sup.b(i)+p.sup.c.sup.u(i)+p.sup.c.sup.s(i)]−[p.sub.d,min(i−1)−p.sup.d.sup.b(i−1)−p.sup.d.sup.u(i−1)−p.sup.d.sup.s(i−1)+p.sup.c.sup.b(i−1)+p.sup.c.sup.u(i−1)+p.sup.c.sup.s(i−1)]≤p.sub.g.sup.dev(i),

[0100] if i=1,

[p.sub.d,max(i)−p.sup.d.sup.b(i)−p.sup.d.sup.u(i)−p.sup.d.sup.s(i)+p.sup.c.sup.b(i)+p.sup.c.sup.u(i)+p.sup.c.sup.s(i)]−[p.sub.g.sup.buy,prev−p.sub.g.sup.sell,prev]≤p.sub.g.sup.dev(i)

The condition that the power demand and power supply are balanced can be expressed by the following equation:

Load−PV−ESS−Grid=NetDemand−ESS−Grid=0

[0101] where Load=load [0102] PV=electricity generated by the photovoltaic system [0103] ESS=ESS charge/discharge rate [0104] Grid=consumption of grid power [0105] NetDemand=Total Demand

[0106] Therefore, NetDemand−ESS=Grid.

Utilizing this, [p.sub.d,max(i)−p.sup.d.sup.b(i)−p.sup.d.sup.u(i)−p.sup.d.sup.s(i)+p.sup.c.sup.b(i)+p.sup.c.sup.u(i)+p.sup.c.sup.s(i)] is the maximum value of the grid power consumption in the i-th time period. Also, [p.sub.d,min(i−1)−p.sup.d.sup.b(i−1)−p.sup.d.sup.u(i−1)−p.sup.d.sup.s(i−1)+p.sup.c.sup.b(i−1)+p.sup.c.sup.u(i−1)+p.sup.c.sup.s(i−1)] is the minimum value of the grid power usage in the i−1th time period. Therefore, the given expression is a condition that (maximum value of grid power consumption in the i-th time period)−(minimum value of grid power consumption in the i−1 th time period) has a value less than or equal to p.sub.g.sup.dev(i). When i=1, use [p.sub.g.sup.buy,prev−p.sub.g.sup.sell,prev] instead of (minimum value of grid power consumption in the i−1 th time period).

[0107] 3) 0≤p.sub.g.sup.dev(i)

[0108] It is a condition that p.sub.g.sup.dev(i) is greater than or equal to zero.

[00012]$\begin{array}{cc}4)\stackrel{N}{\underset{i=1}{.Math.}}{c}_{sm,g}.Math.\mathrm{dt}.Math.p_g^{\mathrm{dev}}\left(i\right)& \end{array}$

If this expression is added to J, p.sub.g.sup.dev(i) which is a variable for smoothing the power consumption is reflected in the objective function. To reduce the value of J, the value of p.sub.g.sup.dev(i) is reduced, so that the difference of power consumption from the previous time period is reduced for each time period, thereby grid power smoothing is performed. If the value of c.sub.sm,g is larger, the smoothing effect is reflected more.

[Sixth Module]

[0109] 1)

−P.sub.ESS.sup.dev(i)≤[p.sup.d.sup.b(i)+p.sup.d.sup.b(i)+p.sup.d.sup.s(i)−p.sup.c.sup.b(i)−p.sup.c.sup.b(i)−p.sup.c.sup.s(i)]−[p.sup.d.sup.b(i−1)+p.sup.d.sup.u(i−1)+p.sup.d.sup.s(i−1)−p.sup.c.sup.b(i−1)−p.sup.c.sup.u(i−1)−p.sup.cis s(i−1)]

[0110] if i=1,

−p.sub.ESS.sup.dev(i)≤[p.sup.d.sup.b(i)+p.sup.d.sup.u(i)+p.sup.d.sup.s(i)−p.sup.c.sup.b(i)−p.sup.c.sup.u(i)−p.sup.c.sup.s(i)]−[p.sub.ESS.sup.dis,prev−p.sub.ESS.sup.chg,prev]

Here, [p.sup.d.sup.b(i)+p.sup.dis u(i)+p.sup.dis s(i)−p.sup.c.sup.b(i)−p.sup.c.sup.u(i)−p.sup.c.sup.s(i)] is ESS charge/discharge rate in the i-th time period. Also, [p.sup.d.sup.b(i−1)+p.sup.d.sup.u(i−1)+p.sup.d.sup.s(i−1)−p.sup.c.sup.b(i−1)−p.sup.c.sup.u(i−1)−p.sup.c.sup.s(i−1)] is the ESS charge/discharge rate in the i−1th time period. Therefore, the given expression is a condition that (ESS charge/discharge rate in the i-th time period)−(ESS charge/discharge rate in the i−1th time period) has a value greater than or equal to −p.sub.ESS.sup.dev(i). When i=1, use [p.sub.ESS.sup.dis,prev−p.sub.ESS.sup.chg,prev] instead of (ESS charge/discharge rate in the i−1 time period).

[0111] 2)

[p.sup.d.sup.b(i)+p.sup.d.sup.u(i)+p.sup.d.sup.s(i)−p.sup.c.sup.b(i)−p.sup.c.sup.u(i)−p.sup.cis s(i)]−[p.sup.d.sup.b(i−1)+p.sup.d.sup.u(i−1)+p.sup.d.sup.s(i−1)−p.sup.c.sup.b(i−1)−p.sup.c.sup.u(i−1)−p.sup.c.sup.s(i−1)]≤p.sub.ESS.sup.dev(i)

[0112] If i=1,

[p.sup.d.sup.b(i)+p.sup.d.sup.u(i)+p.sup.d.sup.s(i)−p.sup.c.sup.b(i)−p.sup.c.sup.u(i)−p.sup.c.sup.s(i)]−[p.sub.ESS.sup.dis,prev−p.sub.ESS.sup.chg,prev]≤p.sub.ESS.sup.dev(i)

Here [p.sup.d.sup.b(i)+p.sup.d.sup.u(i)+p.sup.d.sup.s(i)−p.sup.c.sup.b(i)−p.sup.c.sup.u(i)−p.sup.c.sup.s(i)] is the ESS charge/discharge rate in the i-th time period. Also, [p.sup.d.sup.b(i−1)+p.sup.d.sup.u(i−1)+p.sup.d.sup.s(i−1)−p.sup.c.sup.b(i−1)−p.sup.c.sup.u(i−1)−p.sup.c.sup.s(i−1)] is the ESS charge/discharge rate in the i−1th time period. Therefore the given expression is a condition that (ESS charge/discharge rate in the i-th time period)−(ESS charge/discharge rate in the i−1th time period) has a value less than or equal to p.sub.ESS.sup.dev(i). When i=1, [p.sub.ESS.sup.dis,prev−p.sub.ESS.sup.chg,prev] is used instead of (ESS charge/discharge rate in the i−1th time period).

[0113] 3) 0≤p.sub.ESS.sup.dev(i)

It is a condition that p.sub.ESS.sup.dev(i) has a value greater than or equal to zero.

[00013]$\begin{array}{cc}4)\stackrel{N}{\underset{i=1}{.Math.}}{c}_{sm,ESS}.Math.\mathrm{dt}.Math.p_{ESS}^{dev}\left(i\right)& \end{array}$

If this term is added to J, p.sub.ESS.sup.dev(i) which is a variable for smoothing the ESS charge/discharge rate is reflected in the objective function. To reduce the value of J, the value of p.sub.ESS.sup.dev(i) is reduced, so that the difference of ESS charge/discharge rate from the previous time period is reduced for each time period, thereby ESS charge/discharge rate smoothing is performed. If the value of C.sub.sm,ESS is larger, the smoothing effect is reflected more.

[Seventh Module]

[0114] 1)) p.sup.c.sup.b(i)+p.sup.c.sup.u(i)+p.sup.c.sup.s(i) is the ESS charge rate, and p.sup.d.sup.b(i)+p.sup.d.sup.u(i)+p.sup.d.sup.s(i) is the ESS discharge rate, therefore p.sup.c.sup.b(i)+p.sup.c.sup.u(i)+p.sup.c.sup.s(i)−p.sup.d.sup.b(i)−p.sup.d.sup.u(i)−p.sup.d.sup.s(i) is the ESS charge/discharge rate. The condition that the power demand and power supply are balanced can be expressed by the following equation:

Load−PV−ESS−Grid=NetDemand−ESS−Grid=0

[0115] where Load=load [0116] PV=electricity generated by the photovoltaic system [0117] ESS=ESS charge/discharge rate [0118] Grid=consumption of grid power [0119] NetDemand=Total Demand

[0120] Therefore, NetDemand−ESS=Grid.

And the equation NetDemandMax−ESS=GridMax may be formulated when NetDemandMax is the maximum value of total demand, and GridMax is the maximum value of the grid power consumption, the following inequality can be formulated to express a condition that the maximum value of grid power consumption is less than the peak size.

[00014]$-p^{d_b}\left(i\right)-p^{d_u}\left(i\right)-p^{d_s}\left(i\right)+p^{c_b}\left(i\right)+p^{c_u}\left(i\right)+p^{c_s}\left(i\right)+{\stackrel{\_}{p}}_{d,\max }\left(i\right)\le p_{PC}^{g_{buy}}\left(i\right)+p_{\mathrm{PC}}^{g_b^{\mathrm{peak}}}\left(i\right)$

If the maximum value of the grid power consumption is large and the inequality does not hold with the set peak size, the condition is relaxed by increasing the relaxation variable p.sub.PC.sup.g.sup.buy(i).

[00015]$2)p_{\mathrm{PC}}^{g_{\mathrm{buy}}}\left(i\right)\le \left(1-{\delta }_{\mathrm{PC}}^{g_{\mathrm{buy}}}\left(i\right)\right)\left(p_{\mathrm{ESS}}^{\mathrm{chg},\max }\left(i\right)+{\stackrel{\_}{p}}_{d,\max }\left(i\right)-p_{\mathrm{PC}}^{g_{\mathrm{buy}}^{\mathrm{peak}}}\left(i\right)\right)i)p_{\mathrm{PC}}^{g_{\mathrm{buy}}}\left(i\right)\le \left(1-1\right)\left(p_{\mathrm{ESS}}^{\mathrm{chg},\max }\left(i\right)+{\stackrel{\_}{p}}_{d,\max }\left(i\right)-p_{\mathrm{PC}}^{g_{\mathrm{buy}}^{\mathrm{peak}}}\left(i\right)\right)$$\mathrm{if}{\delta }_{\mathrm{PC}}^{g_{\mathrm{buy}}}\left(i\right)=1$

(if the peak control is successful), therefore p.sub.PC.sup.g.sup.buy(i)≤0

[00016]$\begin{array}{cc}\mathrm{ii})p_{PC}^{g_{\mathrm{buy}}}\left(i\right)+p_{\mathrm{PC}}^{g_{\mathrm{buy}}^{\mathrm{peak}}}\left(i\right)\le p_{ESS}^{\mathrm{chg},\max }\left(i\right)+{\stackrel{\_}{p}}_{d,\max }\left(i\right)& \end{array}$

if δ.sub.PC.sup.g.sup.buy(i)=0 (if the peak control fails), therefore

[00017]$\begin{array}{cc}{p}_{\mathrm{PC}}^{g_{\mathrm{buy}}}\left(i\right)\le p_{\mathrm{ESS}}^{\mathrm{chg},\max }\left(i\right)+{\stackrel{\_}{p}}_{d,\max }\left(i\right)-p_{\mathrm{PC}}^{g_{\mathrm{buy}}^{\mathrm{peak}}}\left(i\right),& \end{array}$

which is equivalent to

[00018]$\begin{array}{cc}{p}_{PC}^{g_{\mathrm{buy}}}\left(i\right)+p_{\mathrm{PC}}^{g_{\mathrm{buy}}^{\mathrm{peak}}}\left(i\right)\le p_{ESS}^{\mathrm{chg},\max }\left(i\right)+{\stackrel{\_}{p}}_{d,\max }\left(i\right)& \end{array}$

It is a condition that even if the peak size, which is the upper limit of the grid power consumption, and the value of the relaxation variable are added, it does not have to be greater than the sum of the ESS maximum charging rate and total demand.

[0121] 3) 0≤p.sub.PC.sup.g.sup.buy(i)≤1

δ.sub.PC.sup.g.sup.buy(i) has a value between zero and 1, therefore 0 or 1 because it is an integer variable.

[0122] 4) 0≤p.sub.PC.sup.g.sup.buy(i)

p.sub.PC.sup.g.sup.buy(i) has a nonnegative value.

[00019]$\begin{array}{cc}5)\underset{i\in {\mathrm{PC}}_{g,\mathrm{buy}}}{.Math.}\left[-c_{\mathrm{PC},1}^{g_{\mathrm{buy}}}\left\{1+{\left(c_{\mathrm{PC},2}^{g_{\mathrm{buy}}}\right)}^i\right\}\mathrm{dt}.Math.{\delta }_{\mathrm{PC}}^{g_{\mathrm{buy}}}\left(i\right)+c_{\mathrm{PC},3}^{g_{\mathrm{buy}}}\mathrm{dt}.Math.p_{\mathrm{PC}}^{g_{\mathrm{buy}}}\left(i\right)\right]& \end{array}$

By adding this term to J, whether the peak control is successful will be reflected in the objective function. i∈P.sub.g,buy means that the i-th time period is included in the time periods for performing peak control. If the peak control is successful, δ.sub.PC.sup.g.sup.buy(i)=1, therefore the value of j is decreased by the term −c.sub.PC,1.sup.g.sup.buy{1+(c.sub.PC,2.sup.g.sup.buy).sup.i}dt.Math.1. If the value of c.sub.PC,1.sup.g.sup.buy is larger, the success of peak control is reflected more. If the peak control fails, δ.sub.PC.sup.g.sup.buy(i)=0, therefore −c.sub.PC,1.sup.g.sup.buy{1+(c.sub.PC,2.sup.g.sup.buy).sup.i}dt.Math.0=0, and the value of J will not decrease. Since 0≤p.sub.PC.sup.g.sup.buy(i), the value of J increases by c.sub.PC,3.sup.g.sup.buydt.Math.p.sub.PC.sup.g.sup.buy(i). If c.sub.PC,3.sup.g.sup.buy is larger, then whether the peak control fails is reflected more.

[Eighth Module]

[0123] 1)

{p.sub.d,min(i)+p.sub.ESS.sup.chg,max(i)}δ.sub.IO.sup.g(i)−p.sup.d.sup.b(i)−p.sup.d.sup.u(i)−p.sup.d.sup.s(i)+p.sup.c.sup.b(i)+p.sup.c.sup.u(i)+p.sup.c.sup.s(i)≤p.sub.ESS.sup.chg,max(i) [0124] i) The condition that the power demand and power supply are balanced can be expressed by the following equation:

Load−PV−ESS−Grid=NetDemand−ESS−Grid=0

[0125] where Load=load [0126] PV=electricity generated by the photovoltaic system [0127] ESS=ESS charge/discharge rate [0128] Grid=consumption of grid power [0129] NetDemand=Total Demand

[0130] Therefore, NetDemand−ESS=Grid.

When net zero operation is conducted, Grid=0. Therefore NetDemand−ESS=0. For NetDemand, NetDemandMin≤NetDemand holds, and applying NetDemand=ESS, the condition is expressed by the following inequality.

p.sub.d,min(i)≤p.sup.d.sup.b(i)+p.sup.d.sup.u(i)+p.sup.d.sup.s(i)−p.sup.c.sup.b(i)−p.sup.c.sup.u(i)−p.sup.c.sup.s(i)

This is an inequality that can be obtained by applying δ.sub.IO.sup.g(i)=1 in the following inequality and is a condition for successful net zero energy operation.

{p.sub.d,min(i)+p.sub.ESS.sup.chg,max(i)}δ.sub.IO.sup.g(i)−p.sup.d.sup.b(i)−p.sup.d.sup.u(i)−p.sup.d.sup.s(i)+p.sup.c.sup.b(i)+p.sup.c.sup.u(i)+p.sup.c.sup.s(i)≤p.sub.ESS.sup.chg,max(i) [0131] ii) If δ.sub.IO.sup.g(i)=0 (when net zero energy operation fails),

{p.sub.d,min(i)+p.sub.ESS.sup.chg,max(i)}.Math.0−p.sup.d.sup.b(i)−p.sup.d.sup.u(i)−p.sup.d.sup.s(i)+p.sup.c.sup.b(i)+p.sup.c.sup.u(i)+p.sup.c.sup.s(i)≤p.sub.ESS.sup.chg,max(i)

[0132] Therefore,

−p.sup.d.sup.b(i)−p.sup.d.sup.u(i)−p.sup.d.sup.s(i)+p.sup.c.sup.b(i)+p.sup.c.sup.u(i)+p.sup.c.sup.s(i)≤p.sub.ESS.sup.chg,max(i)

[0133] p.sup.c.sup.b(i)+p.sup.c.sup.u(i)+p.sup.c.sup.s(i) is the Ess charge rate, and p.sup.d.sup.b(i)−p.sup.dis u(i)−p.sup.d.sup.s(i) the ESS discharge rate, therefore p.sup.c.sup.b(i)+p.sup.c.sup.u(i)+p.sup.c.sup.s(i)−p.sup.d.sup.b(i)−p.sup.d.sup.u(i)−p.sup.d.sup.s(i) is the ESS charge/discharge rate, and the above inequality is a condition for the maximum Ess charge rate.

[0134] 2)

−{p.sub.d,max(i)−p.sub.ESS.sup.dis,max(i)}δ.sub.IO.sup.g(i)+p.sup.d.sup.b(i)+p.sup.d.sup.u(i)+p.sup.d.sup.s(i)−p.sup.c.sup.b(i)−p.sup.c.sup.u(i)−p.sup.c.sup.s(i)≤p.sub.ESS.sup.dis,max(i) [0135] i) The condition that the power demand and power supply are balanced can be expressed by the following equation:

Load−PV−ESS−Grid=NetDemand−ESS−Grid=0

[0136] where Load=load [0137] PV=electricity generated by the photovoltaic system [0138] ESS=ESS charge/discharge rate [0139] Grid=consumption of grid power [0140] NetDemand=Total Demand

[0141] Therefore, NetDemand−ESS=Grid.

When net zero operation is conducted, Grid=0. Therefore NetDemand−ESS=0. For NetDemand, NetDemand≤NetDemandMax holds, and applying NetDemand=ESS, the condition is expressed by the following inequality.

p.sup.d.sup.b(i)+p.sup.d.sup.u(i)+p.sup.d.sup.s(i)−p.sup.c.sup.b(i)−p.sup.c.sup.u(i)−p.sup.c.sup.s(i)≤p.sub.d,max(i)

This is an inequality that can be obtained by applying in the following inequality and is a condition for successful net zero energy operation.

−{p.sub.d,max(i)−p.sub.ESS.sup.dis,max(i)}δ.sub.IO.sup.g(i)−p.sup.d.sup.b(i)+p.sup.d.sup.u(i)+p.sup.d.sup.s(i)−p.sup.c.sup.b(i)−p.sup.c.sup.u(i)−p.sup.c.sup.s(i)≤p.sub.ESS.sup.dis,max(i) [0142] ii) δ.sub.IO.sup.g(i)=0 (when net zero energy operation fails),

−{p.sub.d,max(i)−p.sub.ESS.sup.dis,max(i)}.Math.0+p.sup.d.sup.b(i)+p.sup.d.sup.u(i)+p.sup.d.sup.s(i)−p.sup.c.sup.b(i)−p.sup.c.sup.u(i)−p.sup.c.sup.s(i)≤p.sub.ESS.sup.dis,max(i)

[0143] Therefore,

p.sup.d.sup.b(i)+p.sup.d.sup.u(i)+p.sup.d.sup.s(i)−p.sup.c.sup.b(i)−p.sup.c.sup.u(i)−p.sup.c.sup.s(i)≤p.sub.ESS.sup.dis,max(i)

[0144] p.sup.c.sup.b(i)+p.sup.c.sup.u(i)+p.sup.c.sup.s(i) is the Ess charge rate, and p.sup.d.sup.b(i)+p.sup.d.sup.u(i)+p.sup.d.sup.s(i) the ESS discharge rate, therefore p.sup.d.sup.b(i)+p.sup.d.sup.u(i)+p.sup.d.sup.s(i)−p.sup.c.sup.b(i)−p.sup.c.sup.u(i)−p.sup.c.sup.s(i) is the ESS charge/discharge rate, and the above inequality is a condition for the maximum ESS discharge rate.

[0145] 3) The condition that the ESS charge/discharge rate is less than or equal to the sum of the “minimum value of total demand” and the “relaxation variable when condition is not satisfied” is expressed by the following inequality.

p.sup.d.sup.b(i)+p.sup.d.sup.u(i)+p.sup.d.sup.s(i)−p.sup.c.sup.b(i)−p.sup.c.sup.u(i)−p.sup.c.sup.s(i)≤p.sub.d,min(i)+p.sub.IO.sup.g(i)

[0146] 4) The condition that the ESS charge/discharge rate is greater than or equal to the sum of the “maximum value of total demand” and the “relaxation variable when condition is not satisfied” is expressed by the following inequality.

p.sub.d,max(i)−p.sub.IO.sup.g(i)−p.sup.d.sup.b(i)+p.sup.d.sup.u(i)+p.sup.d.sup.s(i)−p.sup.c.sup.b(i)−p.sup.c.sup.u(i)−p.sup.c.sup.s(i)

[0147] 5) 0≤δ.sub.IO.sup.g(i)≤1

δ.sub.IO.sup.g(i) has a value between zero and 1, and therefore 0 or 1 because it is an integer variable.

[0148] 6) 0≤p.sub.IO.sup.g(i)

p.sub.IO.sup.g(i) has a nonnegative value.

[00020]$\begin{array}{cc}7)\underset{i\in {\mathrm{IO}}_g}{.Math.}\left[-c_{\mathrm{IO},1}^g\left\{1+{\left(c_{\mathrm{IO},2}^g\right)}^i\right\}\mathrm{dt}.Math.{\delta }_{\mathrm{IO}}^g\left(i\right)+c_{\mathrm{IO},3}^g\mathrm{dt}.Math.p_{\mathrm{IO}}^g\left(i\right)\right]& \end{array}$

By adding this formula to J, whether the net zero energy operation is successful will be reflected in the objective function. i∈IO.sub.g means that the i-th time period is included in the time periods for performing net zero energy operation. If the net zero energy operation is successful, δ.sub.IO.sup.g(i)=1, therefore the value of J is decreased by the term −c.sub.IO,1.sup.g{1+(c.sub.IO,2.sup.g).sup.i}dt.Math.1. If the value of c.sub.IO,1.sup.g is larger, whether the operation is successful is reflected more. If the net zero energy operation fails, δ.sub.IO.sup.g(i)=0, therefore −c.sub.IO,1.sup.g{1+(c.sub.IO,2.sup.g).sup.i}dt.Math.0=0, and the value of J will not decrease. Since 0≤p.sub.IO.sup.g(i), the value of J increases by c.sub.IO,3.sup.gdt.Math.p.sub.IO.sup.g(i). If c.sub.IO,3.sup.g is larger, then whether the peak control fails is reflected more.

[Ninth Module]

[0149] 1) The condition that the ESS discharge amount must have a value greater than or equal to the set amount of demand response power is expressed by the following inequality.

P.sub.save−P.sub.DR≤Σ{(p.sup.dis b(i)+p.sup.dis u(i)+p.sup.d.sup.s(i)−p.sup.c.sup.b(i)−p.sup.c.sup.u(i)−p.sup.c.sup.s(i))dt}

[0150] If it is not satisfied, the condition is relaxed by increasing the value of the supplementary variable subtracted from the demand response power quantity.

[0151] 2) P.sub.DR≤(1−δ.sub.DR)[P.sub.save+Σ{p.sub.ESS.sup.chg,max(i)dt}] [0152] i) If δ.sub.DR−1 (when demand response is successful)

P.sub.DR≤(1−1)[P.sub.save+Σ{p.sub.ESS.sup.chg,max(i)dt}]

[0153] Therefore, P.sub.DR≤0 [0154] ii) If δ.sub.DR=0 (when demand response fails)

P.sub.DR≤(1−0)[P.sub.save+Σ{p.sub.ESS.sup.chg,max(i))dt}]

Therefore, Σ{−p.sub.ESS.sup.chg,max(i)dt}≤P.sub.save−P.sub.DR
Utilizing the inequality above,

Σ{−p.sub.ESS.sup.chg,max(i)dt}≤P.sub.save−P.sub.DR≤Σ{(p.sup.d.sup.b(i)+p.sup.d.sup.u(i)+p.sup.d.sup.s(i)−p.sup.c.sup.b(i)−p.sup.c.sup.u(i)−p.sup.c.sup.s(i))dt}

Therefore,

[0155]
Σ{−p.sub.ESS.sup.chg,max(i)dt}≤Σ{(p.sup.d.sup.b(i)+p.sup.d.sup.u(i)+p.sup.d.sup.s(i)−p.sup.c.sup.b(i)−p.sup.c.sup.u(i)−p.sup.c.sup.s(i))dt}

Exchanging the terms, the condition for the ESS maximum charging rate is obtained as the following inequality.

{(−p.sup.d.sup.b(i)−p.sup.d.sup.u(i)−p.sup.d.sup.s(i)+p.sup.c.sup.b(i)+p.sup.c.sup.u(i)+p.sup.c.sup.s(i))dt}Σ{p.sub.ESS.sup.chg,max(i)dt}

[0156] 3) 0≤δ.sub.DR≤1

δ.sub.DR has a value between zero and 1, and therefore 0 or 1 because it is an integer variable.

[0157] 4) 0≤P.sub.DR

P.sub.DR has a nonnegative value.

[0158] 5) For each DR,

−c.sub.DR,1{1+(c.sub.DR,2).sup.i}.Math.δ.sub.DR+c.sub.DR,3.Math.P.sub.DR

By adding this formula to J, whether the demand response is successful will be reflected in the objective function. If the demand response is successful, δ.sub.DR=1, therefore the value of J is decreased by the term −c.sub.DR,1{1+(c.sub.DR,2).sup.i}dt.Math.1. If the value of c.sub.DR,1 is larger, whether it is successful is reflected more. If the demand response fails, δ.sub.DR=0, therefore −c.sub.DR,1{1+(c.sub.DR,2).sup.i}dt.Math.0=0, and the value of J will not decrease. Since 0≤P.sub.DR, the value of J increases by c.sub.DR,3.Math.P.sub.DR. If c.sub.DR,3 is larger, then whether it fails is reflected more.

Example

[0159] Using the modules described above, tests were conducted as shown in tables 2-4.

TABLE-US-00002 TABLE 2 Time Period 0~1 1~2 2~3 3~4 4~5 5~6 6~7 7~8 8~9 9~10 10~11 11~12 Demand 6.5 6.8 6.6 6.2 4.6 5 6.9 8.9 12.4 17.8 26.5 27.7 Time Period 12~13 13~14 14~15 15~16 16~17 17~18 18~19 19~20 20~21 21~22 22~23 23~0 Demand 25.7 24 24.8 24.5 24.6 26.6 23.7 17.6 14.3 12.9 9.8 9.5

TABLE-US-00003 TABLE 3 Time Period 0~1 1~2 2~3 3~4 4~5 5~6 6~7 7~8 8~9 9~10 10~11 11~12 Supply 0 0 0 0 0 0 0 3 6 9 15 15 Time Period 12~13 13~14 14~15 15~16 16~17 17~18 18~19 19~20 20~21 21~22 22~23 23~0 Supply 15 18 12 12 6 3 0 0 0 0 0 0

TABLE-US-00004 TABLE 4 Time Period 0~1 1~2 2~3 3~4 4~5 5~6 6~7 7~8 8~9 9~10 10~11 11~12 Buying 66.1 66.1 66.1 66.1 66.1 66.1 66.1 66.1 66.1 96.5 111.3 111.3 Price Per Unit Selling 66.1 66.1 66.1 66.1 66.1 66.1 66.1 66.1 66.1 96.5 111.3 111.3 Price per Unit Time Period 12~13 13~14 14~15 15~16 16~17 17~18 18~19 19~20 20~21 21~22 22~23 23~0 Buying 96.5 96.5 96.5 96.5 96.5 111.3 111.3 111.3 96.5 96.5 111.3 66.1 Price Per Unit Selling 96.5 96.5 96.5 96.5 96.5 111.3 111.3 111.3 96.5 96.5 111.3 66.1 Price per Unit

[0160] The settlement cost with the grid of the case when the EMS/ESS is not operated and that when the EMS/ESS is operated with the optimal operation schedule is shown in Table 5. Here no external conditions (peak control, net zero energy operation, demand response) were imposed

TABLE-US-00005 TABLE 5 Total Power Cost Cost Saving by the ESS With the Operation Operation of Capacity of EMS/ESS EMS/ESS 40 24370.8 215.5 50 24341.6 244.7 60 24313.0 273.3 70 24283.8 302.5 80 24255.3 331.0 90 24226.1 360.2 100 24197.6 388.7 110 24168.4 417.9 120 24139.8 446.5

[0161] Table 3 shows forecasted values (kW) of the power demand schedule, Table 4 shows forecasted values (kW) of the power supply schedule, and Table 4 shows the forecasted value (won/kWh) of the power unit price schedule. These forecasts may be monitored or displayed (visualized) externally as shown in FIGS. 3 to 5. Here, the power forecast, power supply and demand status, power transaction price, current weather, etc. can be displayed, and the power forecast can include peak time forecast, maximum demand forecast, operational reserve, and operational reserve ratio. And the power supply and demand status may include supply capacity, current demand, supply reserve power, supply reserve ratio, etc., and the power transaction price may include maximum, average, minimum, current price, and the current weather may include lowest/highest and the current temperature.

[0162] While the embodiments of the present disclosure and their advantages have been described in detail, it should be understood that various changes, substitutions and alterations may be made herein without departing from the scope of the present disclosure.