Optimal Bidding Strategy for Renewable Microgrid with Active Network Management

Active Network Management (ANM) enables a microgrid to optimally dispatch the active/reactive power of its Renewable Distributed Generation (RDG) and Battery Energy Storage System (BESS) units in real time. Thus, a microgrid with high penetration of RDGs can handle their uncertainties and variabilities to achieve the stable operation using ANM. However, the actual power flow in the line connecting the main grid and microgrid may deviate significantly from the day-ahead bids if the bids are determined without consideration of the real-time adjustment through ANM, which will lead to a substantial imbalance cost. Therefore, this study proposes a formulation for obtaining an optimal bidding which reflects the change of power flow in the connecting line by real-time adjustment using ANM. The proposed formulation maximizes the expected profit of the microgrid considering various network and physical constraints. The effectiveness of the proposed bidding strategy is verified through the simulations with a 33-bus test microgrid. The simulation results show that the proposed bidding strategy improves the expected operating profit by reducing the imbalance cost to a greater degree compared to the basic bidding strategy without consideration of ANM.


Introduction
Microgrids (MGs) are medium or low voltage distribution networks that comprise distributed generations, electrical loads, and energy storage systems [1,2].A Microgrid Operator (MGO) aims to achieve stable and efficient operation of the MG through the proper use of the available resources [3,4].In particular, in the grid-connected mode, a MGO can participate in the electricity market to buy or sell energy to maintain the balance between supply and demand of the MG in a most economical way [5,6].Recently, much research has been conducted on establishing an optimal bidding strategy for the MGO in charge of operating MGs with high penetration of Renewable Distributed Generation (RDG) systems, the so-called renewable MGs [7][8][9][10][11].Stochastic approaches have been commonly used in the relevant research [7][8][9][10][11] to obtain an optimal bidding strategy for renewable MGs, because RDGs produce significantly uncertain and variable output compared to conventional generators.Optimal bidding and operation strategies are derived in [7][8][9][10][11] to minimize expected operating costs through a stochastic approach considering the uncertainties in the outputs of Wind Turbines (WTs) or Photovoltaics (PVs), and demand.In particular, the previous approaches in [7][8][9] focus on the scheduling of dispatchable Energies 2016, 9, 48 2 of 15 resources in the MG rather than strategic bidding of the MGO in the electricity market.Further in [10,11], the MGO considers the imbalance cost caused by the difference between the submitted bid and the actual power flow in the line connecting the main grid and MG.In [11], a tradeoff between maximizing the expected operating profit of the MG and reducing the risk of low expected operating profit in unfavorable scenarios is considered to develop a bidding strategy for the MGO.Furthermore, in [11] it is shown that the MGO can increase its profit by utilizing demand response as an additional controllable resource.
The bidding strategies proposed in [7][8][9][10][11] do not however consider network conditions and physical constraints, such as lower and upper voltage limits for each bus and line capacities.Furthermore, although the reactive power balance in the MG is as important as the active power balance for stable operation of a MG network, none of previous approaches described in [7][8][9][10][11] consider the reactive power demands in their formulation of an optimal bidding strategy.For these reasons, there is a high risk of violating the physical constraints in real-time operation.As a solution for resolving this problem, the MGO can curtail the active power outputs of RDGs, adjust their power factors, or schedule active/reactive powers of the Battery Energy Storage Systems (BESSs) in real-time operation.
The above real-time adjustments of the MGO can be realized through Active Network Management (ANM) [12][13][14][15].Owing to the real-time adjustments using ANM, the power flow in the line connecting the main grid and MG can deviate significantly from the submitted bid, which is obtained according to the bidding strategies proposed in [10,11].For this reason, the MGO is likely to be charged for a significant imbalance cost.Therefore, this study proposes an optimal bidding strategy that enables the MGO to consider the change of the power flow on the connected line by real-time adjustments using ANM, through a scenario-based stochastic optimization method.By using the proposed bidding strategy, the imbalance cost can be reduced and the expected operating profit can be improved.
The remainder of the paper is organized as follows: in Section 2, the assumptions for the market participation of the MGO are summarized.Section 3 describes the basic bidding strategy of the MGO without consideration of ANM.Section 4 introduces the features of ANM in MG operation, and proposes the optimization approach to achieve an optimal bidding strategy considering ANM for day-ahead bidding.Section 5 presents the limitations of the basic bidding strategy described in Section 3 through simulation, and verifies the improvement in the daily expected operating profit when the proposed bidding strategy is used.Finally, Section 6 contains the conclusions of this study.

Microgrid Operators as Market Participant
In this study, it is assumed that a MGO is allowed to participate in the day-ahead energy market and submit bids as either a supplier or customer, according to the expected direction of the power flow in the line connecting the main grid and MG.In order to forecast the quantity and direction of power flow in the connecting line, it is necessary for the MGO to know forecasted values for electricity demands of buses and RDG outputs in the MG network.Various artificial intelligence-based forecasting techniques can be effectively used to obtain more accurate forecasting data in a MG environment [16,17].
However, the MGO in charge of operating a renewable MG is unwillingly to violate the submitted bid in the day-ahead bidding because of the large uncertainties in the forecasted outputs of RDGs.The difference between the submitted bid in the day-ahead bidding and the actual power flow in the line connecting in the real-time operation is defined as an imbalance [18].In order to keep the main grid reliable and stable in the presence of this imbalance, the System Operator (SO) of the main grid has to pay extra costs to secure reserves, or prepare flexible resources.Therefore, the SO penalizes the MGO who causes the imbalance.Figure 1 shows the principle of the charging imbalance cost for market participants.If power flow in the connecting line is greater than the upper limit of the tolerance band, the violation is charged at the penalty price for under-scheduling imbalance.Under-scheduling imbalance causes a power shortage in the main grid.Moreover, if the power flow on the connecting Energies 2016, 9, 48 3 of 15 line is less than the lower limit of the tolerance band, the violation is charged at the penalty price for over-scheduling imbalance.Over-scheduling imbalance causes a surplus of power in the main grid.However, no imbalance cost is charged if the power flow in the connecting line is within the tolerance band of the submitted bid.Table 1 summarizes the market rules for the calculation of imbalance cost used in the previous researches [10,11] as well as in the electricity market [19].As shown in Table 1, [10,11] use a simplified rule to calculate the imbalance cost.They do not consider the tolerance band is which the penalty price is not charged.On the contrary, the tolerance band is considered in the practical market rule for imbalance cost calculation in [19].Meanwhile, the amount of energy exchanged with the MGO is very small compared to the whole amount traded in the market.Thus, in this study it is assumed that the MGO is a price-taker that does not affect market price.
Energies 2016, 9, 48 3 of 14 in the main grid.However, no imbalance cost is charged if the power flow in the connecting line is within the tolerance band of the submitted bid.Table 1 summarizes the market rules for the calculation of imbalance cost used in the previous researches [10,11] as well as in the electricity market [19].As shown in Table 1, [10,11] use a simplified rule to calculate the imbalance cost.They do not consider the tolerance band is which the penalty price is not charged.On the contrary, the tolerance band is considered in the practical market rule for imbalance cost calculation in [19].
Meanwhile, the amount of energy exchanged with the MGO is very small compared to the whole amount traded in the market.Thus, in this study it is assumed that the MGO is a price-taker that does not affect market price.

Basic Bidding Strategy Without Active Network Management
This section describes the basic bidding strategies based on the previous studies without considering ANM [10,11].In the basic approach, the MGO determines optimal bids for the next 24 h in the day-ahead market through scenario-based stochastic optimization in order to maximize the expected operating profit, represented as Equation (1), subject to the constraint given by Equations ( 2)-( 6): subject to:

Basic Bidding Strategy Without Active Network Management
This section describes the basic bidding strategies based on the previous studies without considering ANM [10,11].In the basic approach, the MGO determines optimal bids for the next 24 h in the day-ahead market through scenario-based stochastic optimization in order to maximize the expected operating profit, represented as Equation (1), subject to the constraint given by Equations ( 2)-( 6 subject to: Energies 2016, 9, 48 4 of 15 The functions of revenue and other costs in Equation ( 1) are explained in Section 4 below.The power flow on the connecting line P s CON,t should not exceed the capacity of the corresponding line as in Equation ( 2).Further, the charging/discharging power P s BESS,i,t is limited by the capacity of the power conditioning system as formulated in Equation ( 3).The state-of-charge SoC s i,t of the i-th BESS should remain within the range of Equation ( 4) to avoid damage to the battery due to deep discharging or overcharging, and SoC s i,t can be calculated using the equality constraint in Equation ( 5).For simplicity, the charging/discharging efficiencies of BESSs are assumed to be 100% in this study.Finally, the constraint on the active power balance of the supply and demand in the MG is represented in Equation ( 6).

Features of Active Network Management in Microgrid Operation
A renewable MG suffers from problems because of the intermittent outputs of RDGs.For example, an under-voltage problem may occur at certain buses in the MG if there is no wind or sunshine at certain times.On the other hand, an over-voltage problem may occur at certain buses in the MG, if RDGs generate excessive output according to the weather conditions at certain times.ANM can be a solution to these operational problems.ANM is a smart-grid technology that enables a MGO to measure the state of each bus, such as the voltage and angle, in real-time; calculate the optimal dispatch of resources considering the network conditions and physical constraints; and then deliver the dispatched signals to resources through a communication infrastructure [12][13][14][15].Therefore, a MGO with ANM is able to stably and efficiently operate a renewable MG by optimally curtailing the outputs of RDGs and loads and scheduling the active/reactive powers of BESSs without violating physical constraints.

Formulation of Optimal Bidding Strategy with Active Network Management
Similar to the basic bidding strategy described in Section 3, the proposed bidding strategy aims to maximize the daily expected operating profit.However, the proposed bidding strategy is significantly different from the basic bidding strategy in that additional control variables, constraints and extra cost factors related to ANM are included in its formulation.

Objective function
The objective function of the optimization problem for an optimal bidding strategy with ANM is as follows: A MGO with ANM is able to adjust the forecasted demand Ps D,i,t to the dispatched load P s L,i,t during the real-time operation.Therefore, unlike the basic bidding strategy that uses Ps D,i,t , the revenue in the proposed bidding strategy can be calculated with P s L,i,t as: R s " The MGO needs to pay an incentive or a compensation for load curtailment to consumers who have contracted curtailable loads.This study assumes that all of the loads in the MG are curtailable.The cost of load curtailment can be modeled as follows: where the incentive price for load curtailment can be set by multiplying λ t by γ LCURT .The form of the compensation price γ LCURT λ t varies according to the policy, tariff, or contracts with consumers.The generation cost paid for RDGs should be calculated with the dispatched active power output P s G,i,t , rather than the forecasted output Ps G,i,t because the MGO with ANM is able to curtail the outputs of RDGs.In this study, it is assumed that the curtailed quantity Ps G,i,t ´Ps G,i,t is compensated at the compensation price γ GCURT λ t , where γ GCURT can be determined according to the various policies of the electricity company, or regulator.The costs incurred from these two factors are represented as: Because the lines of the distribution network in the MG have larger R/X ratios than a high voltage transmission line, it is not desirable to neglect the cost of the line loss due to the line resistance.In particular, the MGO with ANM is able to calculate this cost by monitoring state information such as the voltage and angle, at each bus.The cost of the line loss can be calculated at λ t as follows: The MGO participates in the day-ahead energy market to buy or sell energy for balancing the supply and demand of the MG.The exchange cost C s EX charged for line flow P s CON,t which is delivered through the line connecting the main grid and MG, can be represented as follows: Note that P s CON,t is positive when the MG draws power from the main grid, and negative when the MG provides power to the main grid.
In this study, the imbalance cost is modeled according to the principle described in Figure 1 as follows: Ipxq " Further, it is assumed that the MGO has BESS(s) for economic and stable MG operation.Accordingly, the MGO has to consider the cost of battery wear in order to avoid excessive operation of the BESS.The cost of battery wear can be adopted from [20] as:

Load constraints
By considering the load curtailment capability of ANM, the ranges of the dispatched active/reactive loads can be represented as: Energies 2016, 9, 48 6 of 15 However, the load power factor of each bus should be maintained even in a curtailed schedule because the power factor reflects the characteristics of the load [21].Hence, it is necessary to consider the constraints on the power factors of the loads as follows:

Generation constraints
In general, the maximum output of a conventional generator is specified at its nominal capacity.On the contrary, the maximum outputs of RDGs, such as WTs and PV systems, depend on the weather conditions at the time and are less than the given nominal capacity.Thus, the MGO can adjust the output of an RDG within a certain range given by: The active/reactive powers can be freely adjusted within the given capability curve, because an RDG is typically an inverter-based machine.For example, Figure 2a shows the capability curve for both WT and PV with the power factor that ranges from the 0.9 lagging to the 0.9 leading.This capability curve in Figure 2a can be represented as inequality constraints as follows: Energies 2016, 9, 48 6 of 14

Generation constraints
In general, the maximum output of a conventional generator is specified at its nominal capacity.On the contrary, the maximum outputs of RDGs, such as WTs and PV systems, depend on the weather conditions at the time and are less than the given nominal capacity.Thus, the MGO can adjust the output of an RDG within a certain range given by: , , , , , The active/reactive powers can be freely adjusted within the given capability curve, because an RDG is typically an inverter-based machine.For example, Figure 2a shows the capability curve for both WT and PV with the power factor that ranges from the 0.9 lagging to the 0.9 leading.This capability curve in Figure 2a can be represented as inequality constraints as follows: )) tan(cos ( )) (a) (b)

BESS constraints
A BESS consists of a battery cell, which charges/discharges a direct current, and a power conditioning system, which converts a direct current into an alternating current: , , 0 , , , 1 The power conditioning system of the BESS can adjust the power factor of the output of the BESS within the four quadrants of the capability curve in Figure 2b [22].The capability curve can be modeled by Equation (23).Equations ( 24) and ( 25) are the same as those for the basic bidding strategy in Section 3.

Power flow constraints
The MGO with ANM is able to obtain state information, such as the voltage and angle, at each bus.Hence, the power flow related constraints can be included in the optimization:

BESS constraints
A BESS consists of a battery cell, which charges/discharges a direct current, and a power conditioning system, which converts a direct current into an alternating current: The power conditioning system of the BESS can adjust the power factor of the output of the BESS within the four quadrants of the capability curve in Figure 2b [22].The capability curve can be modeled by Equation (23).Equations ( 24) and ( 25) are the same as those for the basic bidding strategy in Section 3.

Power flow constraints
The MGO with ANM is able to obtain state information, such as the voltage and angle, at each bus.Hence, the power flow related constraints can be included in the optimization: where: At Bus 1, which is linked to the line connecting the main grid and MG, the active/reactive powers should be satisfied for the power balancing Equation (26).Similarly, the active/reactive power balancing Equation ( 27) must be satisfied at all buses except the first bus.The constraints for Bus 1 in Equation ( 26) additionally include P s CON,t and Q s CON,t .P t F,i,s and Q t F,i,s can be represented as in Equation (28).The power flows in the distribution lines and the connecting line should not exceed their respective line capacities, which can be represented as follows: Additionally, the power flow in the connecting line should be maintained within the power factor range predetermined by the contract with the SO: `Qs CON,t ˘2 ď ´Ps CON,t tanpcos ´1pPFR CON qq ¯2 Finally, the voltages of all buses in the network should be within the regulation range as follows:

Simulation Settings
The performance of the proposed bidding strategy is verified through simulations with a renewable MG based on the IEEE 33-bus distribution network in [23].The structure of the MG is shown in Figure 3. Detailed data for the branches are taken from [23].It is assumed that five WTs, five PVs, and five BESSs are installed in the MG.The capacities of the WTs and PVs are selected to render the MG self-sufficient, such that the sum of the nominal capacities is greater than the daily peak load of the MG.Further, the maximum charging/discharging power in MVA of a BESS is set to half of its capacity in MWh.The specific data related to the WTs, PVs, and BESSs are listed in Table 2.The base values of power and voltage are set to 1 MVA and 12.66 kV, respectively.The forecasted hourly demands in the MG as presented in Figure 4 are taken from the California Independent System Operator (CAISO) for 2 July 2015 [24].The active/reactive power demand of each bus is assumed to be distributed according to the data in [23].It is assumed that the ratio of the forecasted hourly output of each WT to its capacity is the same for all WTs, and the values for this ratio are taken from [25].The same assumption is applied to the forecasted hourly output of each PV, and the ratio from [7] is used for the all PVs.The forecasted hourly outputs during 24 h for WTs and PVs are presented in Figure 5.The forecasting errors for the demand, the WT output, and the PV output in each hour are all assumed to be normally distributed with zero mean.The standard deviations are set to 2% of the forecasted demand [26], 10% of the forecasted WT output [25], and 20% of the PV output [10], respectively.On the basis of these forecasting error models, 1000 scenarios are generated using Monte Carlo simulations, and 10 representative scenarios are selected from them using the scenario-reduction method in [10].
Energies 2016, 9, 48 8 of 14           The day-ahead market prices of CAISO for 2 July 2015 [20] are used for λt, which are presented in Figure 6.For the retail price ρt, the Time-of-Use (TOU) pricing scheme of the Pacific Gas and Electric Company (PG&E, Sacramento, CA, USA) with off-peak ($38.71/MWh),mid-peak ($77.42/MWh), and on-peak ($190.35/MWh)prices is used [27].The day-ahead market prices of CAISO for 2 July 2015 [20] are used for t , which are presented in Figure 6.For the retail price ρ t , the Time-of-Use (TOU) pricing scheme of the Pacific Gas and Electric Company (PG&E, Sacramento, CA, USA) with off-peak ($38.71/MWh),mid-peak ($77.42/MWh), and on-peak ($190.35/MWh)prices is used [27].Two simulation cases are designed to examine the effect of the penalty price by applying different market rule as listed in Table 3. Case I and Case II have the same tolerance band, but the penalty price factor for Case II is set to 1.5 times larger than that for Case I.The values of the other parameters in Equations ( 7)-(31) are listed in Table 4.The optimization problems for both the basic and proposed bidding strategies are solved using GAMS/BONMIN solver on a PC with a 2.30 GHz Intel ® Core i5 CPU, 8 GB RAM, and a 64-bit Windows ® 7 operating system.The relative gap of GAMS/BONMIN Solver is set to 0.1%.7)-(31).

Generation curtailment compensation factor
Fast charging/cycling cost parameters αB, βB 1/0.75 Minimum/maximum state-of-charge of the BESS 20%/95% ∆t 1 h Predetermined power factor range of RDGs 0.9 Lagging-0.9Leading Predetermined power factor range of the line connecting the main grid and MG 0.95 Lagging-0.95Leading Lower/upper voltage limit VLL, VUL 0.95 p.u./1.05 p.u. Two simulation cases are designed to examine the effect of the penalty price by applying different market rule as listed in 3. Case I and Case II have the same tolerance band, but the penalty price factor for Case II is set to 1.5 times larger than that for Case I.The values of the other parameters in Equations ( 7)-(31) are listed in Table 4.The optimization problems for both the basic and proposed bidding strategies are solved using GAMS/BONMIN solver on a PC with a 2.30 GHz Intel ® Core i5 CPU, 8 GB RAM, and a 64-bit Windows ® 7 operating system.The relative gap of GAMS/BONMIN Solver is set to 0.1%.7)-(31).

Simulation Results
If ANM is not used in real-time operation of renewable MG, the RDGs would simply generate the maximum outputs determined by the weather conditions with a fixed unity power factor.In addition, the BESSs would operate according to the active power schedules without considering the network conditions.Then, the voltages at some buses in the MG may not be maintained between the lower and upper voltage limits, and the power flows in some distribution lines may exceed their respective line capacities.Indeed, with the given simulation settings, the under-voltage problem is observed in Buses 10-18 and 28-33 at t = 17 for one representative scenario in Figure 7a.This is because the reactive power supplies from the RDGs or BESSs are not at all scheduled to meet the reactive power demand in the MG, when the basic bidding strategy without ANM being applied.On the other hand, the over-voltage problem is not observed for all representative scenarios, because the output patterns of RDGs are not sufficiently excessive to cause the voltage values of some buses to exceed the upper-voltage limit with the given simulation settings.Further, the power flow of each distribution line does not violate the line capacities with the given simulation settings.Because the line capacities of the test MG are already sufficiently designed to endure the power flow required to meet the demand, the power flow can be considerably reduced by the outputs of the RDGs.Although the under-voltage problem is only observed with the given simulation settings, the over-voltage problem and the violation of line capacity can be observed in other configurations, particularly with large outputs of the RDGs, less demand, or smaller line capacities.In order to resolve the under-voltage problem as shown in Figure 7a, the MGO needs real-time adjustments with ANM such as a curtailment of the outputs of the RDGs and loads and an optimal schedule the active/reactive powers of BESSs.Therefore, it is assumed that real-time adjustments with ANM is also applied, even if the basic bidding strategy is used.As result of real-time adjustments using ANM, the voltages of all buses in the MG can be maintained within the voltage limits as shown in Figure 7b.However, the power flow in the connecting line significantly deviate from the submitted In order to resolve the under-voltage problem as shown in Figure 7a, the MGO needs real-time adjustments with ANM such as a curtailment of the outputs of the RDGs and loads and an optimal schedule the active/reactive powers of BESSs.Therefore, it is assumed that real-time adjustments with ANM is also applied, even if the basic bidding strategy is used.As result of real-time adjustments using ANM, the voltages of all buses in the MG can be maintained within the voltage limits as shown in Figure 7b.However, the power flow in the connecting line significantly deviate from the submitted bid based on the basic bidding strategy.In other words, the basic bidding strategy can cause a large imbalance cost due to the modified power flow in the connecting line while addressing the physical constraints.On the contrary, the proposed bidding strategy considers such as the ANM capabilities when determining day-ahead bids.Therefore, the deviation from the submitted bid and the resulting imbalance cost can be small.
The simulated daily expected operating revenue, costs and profit for the basic and proposed bidding strategies are presented in Table 5.From Table 5, the daily expected operating profit of the proposed bidding strategy is improved compared to the basic bidding strategy by 1.5% ($41.01) for Case I and 2.34% ($63.48) for Case II.This is because the basic bidding strategy causes relatively larger imbalance costs than the proposed bidding strategy.These larger imbalance costs are due to the way in which the MGO determines the bids in the basic bidding strategy without consideration of the real-time adjustments using ANM. Figure 8 shows the hourly imbalance normalized to the bid for each hour in the basic and proposed strategies in both cases.It can be clearly seen that the proposed bidding strategy significantly reduces the imbalance.Consequently, as presented in Table 5, the MGO can increase the expected operating profit using the proposed bidding strategy considering ANM.By comparing Case I and Case II, the improvement in the profit in the proposed bidding strategy significantly increases as a higher penalty price factor is applied.Specifically, the improvement in the profit in the proposed bidding strategy increases by 54.8% more in Case II compared to Case I.This increase occurs because the imbalance cost of the basic and proposed bidding strategies change in different directions.Compared to the results of Case I, the imbalance cost of the basic bidding strategy in Case II increases from $54.04 to $78.04, but that of the proposed bidding strategy decreases from $13.10 to $9.55.Therefore, the difference in the imbalance cost between the two strategies in Case II has a larger value of $68.90 than the $40.94 in Case I.This is because the proposed bidding strategy can determine bids to avoid an excessive imbalance cost at the expense of the greater usage of the control capability of ANM, such as the output curtailment of RDG, the load curtailment, and the BESS.It can be deduced from the result in Table 5 that the differences in the generation curtailment cost, load curtailment cost, and battery wear cost are greater in Case II than in Case I.By comparing Case I and Case II, the improvement in the profit in the proposed bidding strategy significantly increases as a higher penalty price factor is applied.Specifically, the improvement in the profit in the proposed bidding strategy increases by 54.8% more in Case II compared to Case I.This increase occurs because the imbalance cost of the basic and proposed bidding strategies change in different directions.Compared to the results of Case I, the imbalance cost of the basic bidding strategy in Case II increases from $54.04 to $78.04, but that of the proposed bidding strategy decreases from $13.10 to $9.55.Therefore, the difference in the imbalance cost between the two strategies in Case II has a larger value of $68.90 than the $40.94 in Case I.This is because the proposed bidding strategy can determine bids to avoid an excessive imbalance cost at the expense of the greater usage of the control capability of ANM, such as the output curtailment of RDG, the load curtailment, and the BESS.It can be deduced from the result in Table 5 that the differences in the generation curtailment cost, load curtailment cost, and battery wear cost are greater in Case II than in Case I.

Conclusions
The performance of the proposed bidding strategy is verified through simulations with a renewable MG based on the IEEE 33-bus distribution network in [23].The structure of the MG is shown in Figure 3.This study focuses on an optimal bidding strategy for a renewable MG with ANM.It is found that the proposed bidding strategy can significantly reduce the expected imbalance cost and increase the expected operating profit compared to the basic bidding strategy.This is because the proposed bidding strategy can determine day-ahead bids considering the change in the power flow on the connecting line by the real-time adjustments using ANM.Additionally, it is found that the improvement in the profit increases as the penalty price increases.Nevertheless, large expected imbalances still exist at some time stages in Figure 8, even if the proposed bidding strategy is applied.This imbalance comes from the large uncertainties associated with RDG and can harm the stability and reliability of the main grid.Therefore, it is necessary to develop the method for MGO to effectively manage this imbalance in real-time operation in further research.

Figure 2 .
Figure 2. Capability curves of (a) a RDG and (b) a BESS.

Figure 2 .
Figure 2. Capability curves of (a) a RDG and (b) a BESS.

Figure 3 .
Figure 3. Renewable microgrid based on the IEEE 33-bus distribution network.

Figure 3 .
Figure 3. Renewable microgrid based on the IEEE 33-bus distribution network.

Figure 3 .
Figure 3. Renewable microgrid based on the IEEE 33-bus distribution network.

Figure 5 .
Figure 5. Forecasted hourly output of WT and PV normalized to their nominal capacity.

Figure 5 .
Figure 5. Forecasted hourly output of WT and PV normalized to their nominal capacity.

Figure 6 .
Figure 6.Day-ahead market price in CAISO on 2 July 2015.

Figure 6 .
Figure 6.Day-ahead market price in CAISO on 2 July 2015.

Figure 7 .
Figure 7. Simulated voltage profiles of the test MG at t = 17 for one representative scenario (a) without the real-time operation with ANM and (b) with the real-time operation with ANM for the basic bidding strategy.

Figure 7 .
Figure 7. Simulated voltage profiles of the test MG at t = 17 for one representative scenario (a) without the real-time operation with ANM and (b) with the real-time operation with ANM for the basic bidding strategy.

Figure 8 .
Figure 8. Magnitudes of the hourly expected imbalances normalized to the bids for each hour in the basic/proposed bidding strategies for (a) Case I and (b) Case II.

Figure 8 .
Figure 8. Magnitudes of the hourly expected imbalances normalized to the bids for each hour in the basic/proposed bidding strategies for (a) Case I and (b) Case II.

Table 1 .
Market rules for imbalance cost.

Table 1 .
Market rules for imbalance cost.

Table 2 .
Capacities and locations of WTs, PVs, and BESSs.

Table 2 .
Capacities and locations of WTs, PVs, and BESSs.

Table 2 .
Capacities and locations of WTs, PVs, and BESSs.

Table 3 .
Assumed market rules for imbalance cost in simulation.

Table 4 .
Values of other parameters in Equations (

Table 3 .
Assumed market rules for imbalance cost in simulation.

Table 4 .
Values of other parameters in Equations (

Table 5 .
Expected Operating Revenue, Costs and Profit of Two Bidding Strategies for the Simulation Cases.