|تعداد مشاهده مقاله||25,816,991|
|تعداد دریافت فایل اصل مقاله||10,659,935|
یک راه حل مصالحه ای بر اساس تصمیم گیری فازی به منظور برنامه ریزی ساعتی چند هدفه درریز شبکههای خوشه بندی شده با در نظر گرفتن عدم قطیعت منابع انرژی تجدیدپذیر
|هوش محاسباتی در مهندسی برق|
|مقاله 6، دوره 12، شماره 4، دی 1400، صفحه 57-72 اصل مقاله (807.83 K)|
|نوع مقاله: مقاله پژوهشی انگلیسی|
|شناسه دیجیتال (DOI): 10.22108/isee.2020.122174.1349|
|رضا ساکی1؛ اسماعیل رک رک* 1؛ میثم دوستی زاده1؛ محمد عابدینی2|
|1دانشکده فنی مهندسی -دانشگاه لرستان- خرم آباد- ایران|
|2استادیار، دانشکده مهندسی، دانشگاه آیت ا... بروجردی، بروجرد، ایران|
|Video Clustered microgrids؛ Fuzzy decision making؛ Information Gap Decision Theory؛ Renewable energy resources|
The growing development of energy consumption has led to a global energy crisis in recent years. Therefore, energy systems have an important role in modern societies. The several profits of smart microgrids, make them an actual remedy for improving the operation of distribution systems. Microgrids (MGs) can contribute to energy transactions and deliver the anticipation of the distribution system hub. The traditional methods of electrical energy transmission have undertaken some revisions. The traditional structure of the distribution network leads to some historical problems . Wall Street Journal reported, a cyber-attack on only nine electrical substations in the United States will cause main blackouts in the country . By applying the energy resources in distribution systems, smart MGs are formed as an approach to deal with those challenges. However, applying MGs requires advanced technology for controlling and managing several experiments such as technical restrictions and coordination to distribution system parameters . Clustered MGs (CMGs) are collections of MGs and a division of a distribution system. In present distribution systems, the distributed energy resources (DERs) seem vital to supply system loads under several conditions. The energy storage system (ESS) has an important role in DER to improve network performance. It can be used both as a load and as energy storage in the network . Limitation of renewable energy resources (RESs) and ESSs impact on network efficiency . They have a number of restrictions, such as uncertainty of RESs, feeder capacity limits, increasing feeder losses and power flow constraints. For overcoming these constraints, CMGs have been presented in this paper.
Feeder capacity constraint is considered for MGs which are connected via several independent links and lead to a reduction in energy losses via distribution feeders. Hourly robust energy management approach (REMA) is presented for distribution systems under the penetration of RESs. In addition to wind turbines (WTs) and photovoltaics (PVs) as RESs, the MGs are equipped with an energy storage system (ESS) and micro-turbines (MTs).
1.2. Literature Survey
By reviewing recent works on energy management in distribution systems, respected results can be seen; nevertheless, they have some drawbacks in comparison to this study. In  multi-vector energy management approach is proposed for islanded microgrid clusters, taking into account the role of resources and loads in optimally managing the energy flow through the CMGs. In , a bi-level optimization problem is presented in for energy management of a distribution system consisting of a community of MGs, with the economical, operational, and environmental objectives. In , the communication between MGs and distribution network is considered as the main aspect for improving the effective performances of the distribution systems. Reference  motivated on improving the faults in the energy planning of CMGs by a real-time distribution compensation method, which is used to the local MGs. Predictive control model based on energy management framework, is presented in . In this reference authors have tried to minimize the interaction between the distribution system and community of MGs. The suggested algorithm was responsible for optimal organization of different energy resources with the network of MGs, taking into consideration diverse operational, economical, and environmental limitations. A chance-constrained control strategy is proposed in  for improving the local operation of networked microgrids in the presence of renewable-based distributed generations and storage devices. In , a transactive control strategy is introduced for managing the charging/discharging energy pattern of energy storage systems to decline the peak load condition caused by electric vehicles (EVs). Reference  has concentrated on voltage control issues of CMGs considering the medium and low voltage levels individually. The work existing in  has taken into consideration the frequency regulation of networked MGs by regulating a frequency reference point for MGs locally; then it applied the tertiary control approach to manage the flow of power through the point of common of each MG. A three-level multi-layer control strategy is introduced in  which coordinated the operation strategies of energy storage devices and sub-MGs. In , a three-step analysis is performed to investigate the role of MGs as resilience resources in boosting the resilience of the distribution system, where a new MG formation is introduced to explore the resilience of multi-energy MGs. Finally, the distribution network reconfiguration is adopted in  to connect islanded MGs to decline load curtailment as much as possible. In terms of stability, different voltage and current stabilizer controllers are used in  to regulate voltage and current of CMGs, where the proposed control algorithm is defined as a plug-and-play controller. In  the chance of communicating different points of MGs, that are equipped with a drop controller, to form a larger grid is examined. Reference  applied the small-signal analysis for the management and control of large grids. In , the MGs are connected via distribution system arrangements to form CMGs that are strong in face of uncertainty of renewable-based generation units, improvement the entire system resilience, and reservation it from instability which can result in shutdowns. Table 1 presents the taxonomy of existing literature on scheduling and improving the operation of the distribution system.
Table (1): Taxonomy of existing literature on planning and improving the operation of the distribution system
This paper presents CMGs with RESs technologies, the demand response program (DRP) and ESSs. Optimal power transactions among CMGs and the grid are determined. Therefore, this article has the following innovations:
1.4. Paper organization
The structure of this paper is organized as follows. Section 2 discusses the modeling of the proposed method. Section 3 gives the mathematical formulation of the proposed method and the solution method. Section 4 presents the test case and its clustering structure. Section 5 shows and discusses the numerical results and Section 6 draws the conclusions.
2. Modeling and formulation
2.1. Objective Functions
The objective functions in this problem are defined in the form of three objective functions including the cost function of the generation units in each MG, the voltage deviation index (VDI), and the cost of active power losses as described in Eqs (1) -(3).
The active power that flows between MGs should be limited in distribution systems. Power balance must be met in each MG, while MGs can import or export energy from /to other MGs or distribution systems.
AC load flow is mentioned in Eqs (4)-(9).
MGs limitations containing load flow constraints and DRP, and restraining loads in certain hours have been deliberated in the REMA. So, the amount of active and reactive power is limited by the systems agent. In addition, limitations are imposed on the ESSs, where the most important limitations on their charging and discharging status are the corruption of these resources. RESs limitations are other constraints that are straightly related to the energy consumption of microgrids. The modeling and equations for each constraint of the MGs are discussed.
The limitations on active and reactive power balances in each MGs must be considered in the hourly programming. These constraints are shown in the equations which include the balancing of the active and reactive generation power in the MGs. The active and reactive power imported/exported among CMGs and of other clusters must be equal to the sum of the power generated by their DERs. In Eqs (10) and (11) these subjects are seen.
The DRP constraints are clear as restraining the amount of active and reactive power loads to a definite value. According to (14) and (15), the variations forms in the demand will be equal to the sum of base loads after applying the DRP. The DRP program is demonstrated in Eqs (12) - (16).
ESSs constraints contain charging and discharging can be done by applying the charging and discharging factor at the charging and discharging power (17), (18). The charge state parameter limitation is also shown in Eqs (20), (21).
Limitation of MTs reactive power is presented by its generation between zero and maximum under normal operating conditions. Eqs. (23), and (24) illustrate these restrictions.
The predicted power of the WTs is determined by the prediction coefficient. Hence, the active power generated is restricted by this coefficient and it is restricted by the lag or lead angle of reactive power. Limitation of WTs is demonstrated in Eqs (25), (26).
The predicted power of the PVs is determined by the prediction coefficient which is shown in Eq (27).
2.2. Branches Constraints
Branches constraints contain the state of the switches between MGs and clusters. The state of switches is presented by a decision variable that it’s designed in the optimization problem.
There are several methods to multi-objective optimization, the most important of which are comparative as follows.
2.3.1 Weighted sum
In this method, each objective is allocated a scalar weight that indicates its relative importance to other objectives. The optimization problem is then turned into optimizing the weighted sum of different objective functions. The multiple objective optimization problem is transformed into a single objective problem by using a weighted sum of the original multiple objectives. Although the weighted sum method is the simplest and most straightforward way of obtaining multiple points on the Pareto-optimal front, it does not work in finding the non-convex optimal cost. It is hard to select the weighting to ensure that the points are spread evenly on the Pareto front and it is difficult as well in finding all the optimal cost by changing the weights.
Constraint method has been commonly used as an interactive decision-making tool owing to its easy implementation. This method is also known as a trade-off method, which means that the decision-maker specifies a trade-off among the multiple objectives. In this method, one of the objectives is optimized while the others are treated as constraints. This method works in finding the non-convex optimal cost .
In this paper optimization problem is a non-convex problem, so the constraint method is the most efficient method compared to other methods. This is noteworthy that heuristics methods have a problem with the uniqueness of the answer, due to the long convergence, and they are not efficient for hourly planning.
2.4. Min-max Fuzzy decision-maker (FDM)
In this section, we have studied the optimization of objective functions through the constraint method in order to find the best compromise solution through the Pareto Front and Fuzzy Satisfying methods. In the proposed REMA problem, is optimized while and are considered as the constraints as in Eqs (30)-(32).
It is observed from Fig. 2 and Eqs (30)-(32) that and are constrained by parameter . This parameter varies from the minimum value to the maximum value of and gradually, and for any value of and , the modified single objective optimization problem, i.e., (1-29), is solved, and the optimal solutions like in Fig. 1 are obtained. It is noteworthy to mention that in (30) the constraints of the original multi-objective optimization problem, i.e., (1) - (29), which are described in Section 2, are also included. The set of all obtained solutions for the entire variations of is Pareto front of the multi-objective optimization problem.
By solving the REMA problem, a Pareto front is derived and it is required to select the best solution from this Pareto optimal set. The fuzzy decision-maker is used in this paper for this purpose. In this method, a fuzzy membership function is apportioned to each solution in the Pareto front. The fuzzy sets are defined by membership functions. These functions show the degree of membership in fuzzy sets using values from 0 to 1. The membership value ‘0’ indicates inconsistency with the sets, while ‘1’ means compatibility. Fuzzy membership functions that assign a degree of approval to each objective are defined based on which the best solution can be found out of the available Pareto-optimal solutions.
The best compromise solution which can be selected is designated by using the min-max method. In the min-max method the minimum value of , and for each solution is determined, and the solution with the maximum value of is selected as the best compromise solution. The membership functions are calculated as follow:
where refers to the th solution of the th objective function. in this study has 10 members. is the membership objective function for the fuzzy decision, with the maximum value of calculated for the minimum value of (i.e., ). This means that the maximum value of is obtained when is optimally minimized. is also assigned to the smallest value of the membership objective function . The best solution can then be selected using the fuzzy min-max proposition as follow:
Fig. 1. Description of Pareto front in -constraint method.
The best compromise solution which can be selected is described in  by using the min-max method. In the min-max method the minimum value of , and for each solution is determined, and the solution with the maximum value of min ( , , ) is selected as the best compromise solution. In this paper, , and are calculated as (33)-(35).
By solving the MO-HEM problem, a Pareto front is derived and it is required to select the best solution from this Pareto optimal set. The fuzzy decision-maker is used in this paper for this purpose. In this method, a fuzzy membership function is assigned to each solution in the Pareto front. The fuzzy membership is in the interval (0,1).
Table 4 shows the values for the objective function before CMGs and REMA and after CMGs and REMA in 24 hours. It is clear that by applying REMA in the presence of CMG, objective function of losses (OF loss), VDI and Generation cost (G cost) are decreased. This point is important that the numerical value of every objective function has been achieved based on every 24 hours.
Table 4 Values for objective functions
3. IGDT Approach
Information gap decision theory (IGDT) is a famous decision-making tool which was introduced by Ben-Haim in 2006 . IGDT has superiority over other techniques such as scenario-based, and Monte Carlo approaches that have been demonstrated in . Also, the accuracy of the results obtained from IGDT over other methodologies has been examined in . Lastly, the general merits of the IGDT technique are highlighted in  and compared with available uncertainty modeling techniques. Commonly, the IGDT technique is formulated as follows:
where is the uncertain input parameter, while x is the set of decision variables; the behavior of uncertain parameter is described by . The uncertainty set can be mathematically described as:
In (5), is the predicted value of the uncertain parameter, while α denotes the radius of uncertainty. Generally, the IGDT technique is a bi-level optimization problem that is solved in two levels. The first level, called the base case (BC) is solved without consideration for uncertainty, and there is no deviation between the uncertain parameter and its forecasted value. To describe this level, the following equations are used.
The amount of objective function in (6) is obtained in a deterministic optimization problem and considered as the BC for the next level. The second level can be solved from two viewpoints namely, (a) RAS: this strategy is chosen by risk-averse decision-maker which tries to increase the robustness of the system in face of uncertainty, and (b) opportunity seeker strategy: this strategy is followed by optimistic decision-maker which accepts the level of chance to decrease the value of the objective function from its predicted value in the BC. Since this study aims at maximizing the robustness of the whole distribution network including CMGs, it uses the RAS. Accordingly, the objective function, in this case, will be increased by a tolerable value which is defined by decision-maker, while it tries to maximize the uncertainty radius. Consequently, the result obtained from this strategy has a degree of robustness, even if the uncertain parameter would fluctuate.
4. Test Study and Simulation Parameters
The test study and simulation parameters are presented in this section. The effects of CMGs on the distribution system are also studied. The IEEE 94-bus distribution test system , as shown in Fig. 2, is used for simulation purposes. It consists of 37 branches, 32 sectionalizing switches, and 5 tie switches. The proposed complete framework of this study can be easily adapted to any other large-scale network. Construction is implemented on the standard 94 bus distribution test system in Portuguese. It is supposed that 10 MGs are connected to this network, and the MGs are clustered in three groups according to table 2. The location of each MG and the structure of networked MGs, as well as the interconnection paths, are shown on the one-line diagram of the test system in Fig. 2. The data of this system is available in . The proposed problem is simulated as a mixed-integer non-linear programming (MINLP) in GAMS  environment using SBB solver. The distribution system load demand pattern, electricity price, and predicted output power of PVs and WTs are given in Table 3. The rated power of WTs, PVs, and MTs that have been installed in each MG, and the load of MGs are given in .
Table (2): MGs clustered
The proposed model is solved for a 24 hours’ performance perspective.
Table (3): Forecasted hourly demand, output power of WTs & PVs, energy price.
5. Numerical Results
In this section, the results obtained for the proposed REMA are discussed. Then the effects of clustering on the distribution system are also examined. The hourly output power from PV, WT, and MT in each MG is demonstrated in Figs. (3), (4) and (5). These figures show the hourly planning power output for each DG in CMGs; this planning is done by the distributing system agent based on REMA. Output power from PV before 7 and after 17 is zero. It is noteworthy to mention that according to unequal constraints for power generation of DGs (Eqs 23 - 27) if power generation by DGs is more than power generation that the agent has decided for that DG in a special time, that extra power is curtailed to the upstream network, or otherwise it is saved by ESSs or charging stations of EVs. Curtailment power is the difference between the maximum of DG generation power and the maximum power that the agent defines for DGs based on the optimal condition. According to Figs. (3), (4) and (5) for example, it is illustrated that applying REMS, OF loss in 24 hours is 2.85 MW. Thereby in 24 hours distribution system and MGs losses, 280 KW is reduced.
Fig. 2. Single line diagram of 94 bus actual Portuguese distribution test system.
The SOC of ESS in each cluster is shown in Figs. 7, 8, and 9. Cluster number 3 has the maximum SOC and this phenomenon is illustrated by applying HEMS. ESSs are active sources in MGs. For example, in clusters at the time between 2 to 7, ESSs are charged, and at the time between 17 to 22, they are discharged. Fig. 6 shows the power transaction among clusters and the distribution system. For example, cluster#1, in half of the 24 hours, delivers energy to the distribution system in optimal conditions.
Fig.3.Hourly output power from PV in MGs
Fig.4.Hourly output power from MT in MGs
Fig. 5.Hourly output power from WT in MGs
Fig. 6. Hourly power transaction among clusters and the distribution system
Fig. 7. Hourly SOC of ESS in cluster 1
Fig. 8. Hourly SOC of ESS in cluster 2
Hourly voltage constraint condition in buses that have heavy loads compared with and without applying REMA is shown in Figs. 10, 11. It is clear that after applying this framework, the voltage profile is improved. In other words, in this paper, technical constraints are considered and they are improved after this plane. The active power losses of distribution network in this two condition is illustrated in Fig. 11. According to this figure, power losses of the system with REMA is 1,038 KW lower than that of without these technologies, which demonstrates the importance of CMGs in declining the power losses. The status of switches per hour among clusters and between MGs to the network is illustrated in table 5. It is clear that many loads that are settled in MGs, in many hours, are provided by the clusters. Therefore, switches that connect clusters are often closed.
Fig. 9. Hourly SOC of ESS in cluster 3
Fig.10. Hourly voltage profile before optimization
Fig .11. Active power loses of the distribution system in two scenarios
The variation of the DRP index in MGs of different clusters is shown in Fig. 7. It is evident from this figure that the behavior of responsive loads is different in each cluster. Besides, it can be seen that the customers increased their consumption rate in off-peak hours, while they consumed lower load in on-peak hours; which means the DSO can benefit from the DRP of CMGs.
Fig .12. Variation of DR index in different clusters
Table (5): The status of switches per hour among clusters and between MGs to network
In this paper, the REMA strategy is proposed for MGs in clusters. To improve the techno-economic characteristics of the distribution system, each MG is equipped with ESSs, WTs, MTs, PVs, and responsive loads, while the MGs and clusters connected via switch links. The proposed model aims at minimizing the operation cost, losses and, VDI of the distribution system, which can receive energy from upstream networks and MG clusters. The simulation result shows that the distribution system can benefit from CMGs and MGs. In summary, the main conclusions are as follows:
1. The IGDT method leads to the optimal design of switching statues in the REMA framework.
2. The distribution system should inject more active power to the CMGs to increase the strength of the entire system, including the network of CMGs, against uncertainty.
3. Case studies demonstrated that up to 40% of MGs total loads is provided by ESSs in some hours. Hence, ESS charging/discharging pattern plays a key role in REMA.
4. Both DRP and CMGs improve distribution system characteristics and prevent load shedding during daily operation of the distribution system.
5. Having applied REMA decreases total systems power losses and indicate the importance of CMGs in declining the power losses.
 Submission date: 04, 04, 2020
Acceptance date: 08 , 09, 2020
Corresponding author: Esmaeel Rokrok, Department of Electrical Engineering - Lorestan University - Khorramabad - Iran
 Saleh. M, Althaibani. A, Esa. Y, A. Mohamed, "Impact of clustering MGs on their stability and resilience during blackouts," International Conference on Smart Grid and Clean Energy Technologies (ICSGCE) pp.195 – 200, 2015
 Tweed. K, "Attack on Nine Substations Could Take Down U.S. Grid," IEEE Spectrum, 2014.
 Shuaia. Z, YingyunSun. B, JohnShen. Z, et al, "MG stability:Classification and a review," Renewable and Sustainable Energy Reviews, Vol.58, pp.167–179, 2016
 Anurag. C, Saini. R, "A review on Integrated Renewable Energy System based power generation for stand-alone applications Configurations,storage options,sizing methodologies and control," Renewable and Sustainable Energy Reviews, Vol.38, pp.99–120, 2014.
 Kianmehr. E, Nikkhah. S, Vahidinasab.V, et al, "A Resilience-based Architecture for Joint Distributed Energy Resources Allocation and Hourly Network Reconfiguration," IEEE Transactions on Industrial Informatics, Vol. 15, pp.1-12, 2019.
 Yang. M, W. Yanwei, Yang. F, et al," Hierarchical control for DC microgrid clusters with high penetration of distributed energy resources," Electric Power Systems Research , Vol.148 , pp.14-20, 2017.
 Xiong. X, wang. J, Jing. T, et al,"Power optimization control of microgrid cluster," Electric Power Automation Equipment," Vol. 37, pp. 10-17, 2009
 Mehta. R , Srinivasan. D , Trivedi. A , et al,"Hybrid Planning Method Based on Cost-Benefit Analysis for Smart Charging of Plug-In Electric Vehicles in Distribution Systems," IEEE Transactions on Smart Grid, Vol. 10 , pp. 523 – 534, 2019.
 Lu. T, Zhaoyu. W , Qian. A , Lee. W," Interactive Model for Energy Management of Clustered Microgrids," IEEE Transactions on Industry Applications , Vol. 53 , pp. 1739 – 1750, 2017.
 Palizban. O, Kauhaniemi. K,"Energy storage systems in modern grids—Matrix of technologies and applications," Journal of Energy Storage, vol. 6, pp. 248-259, 2016.
 Osama. R. A, Zobaa. A. F and Abdelaziz. A, "A Planning Framework for Optimal Partitioning of Distribution Networks Into Microgrids," in IEEE Systems Journal, Vol. 14, No. 1, pp. 916-926, 2020.
 Khavari. F, Badri. A, Zangeneh. A and Shafiekhani. M ,"Energy management in multi-microgrids via an aggregator to override point of common coupling congestion," IET Generation, Transmission & Distribution, Vol. 13 , pp. 634 - 642 , 2019.
 Asimakopoulou. G. E, Dimeas. A. L and Hatziargyriou. N. D, "Leader-Follower Strategies for Energy Management of Multi-Microgrids," IEEE Transactions on Smart Grid, Vol. 4, pp. 1909-1916, 2013.
 Nunna. H and Doolla. S, "Demand Response in Smart Distribution System With Multiple Microgrids," IEEE Transactions on Smart Grid, Vol. 3, pp. 1641-1649, 2012.
 Nick. M, Cherkaoui. R and Paolone. M, "Optimal Planning of Distributed Energy Storage Systems in Active Distribution Networks Embedding Grid Reconfiguration," IEEE Transactions on Power Systems, Vol. 33, pp. 1577-1590, 2018.
 Arasteh. H, Sepasian. M. S, and Vahidinasab. V, "An aggregated model for coordinated planning and reconfiguration of electric distribution networks," Energy, Vol. 94, pp. 786–798, 2016.
 Wang. Y, Mao. S and Nelms. R, "On Hierarchical Power Scheduling for the Macrogrid and Cooperative Microgrids," IEEE Transactions on Industrial Informatics, Vol. 11, pp. 1574-1584, 2015.
 Liu. X, Gao. B, Zhu. Z, et al, "Non-cooperative and cooperative optimization of battery energy storage system for energy management in multi-microgrid, "IET Gener. Trans. Distrib, Vol 25, pp. 2369–2377, 2018.
 Wang. X, Gong. Y and Jiang. C, "Regional Carbon Emission Management Based on Probabilistic Power Flow With Correlated Stochastic Variables," IEEE Transactions on Power Systems, Vol. 30, pp. 1094-1103, 2015.
 Jafari. A, Shahbazian. A, Fereidunian. A , Nikoofard. A, "Improving self-healing of smart distribution network by allocating switches and distributed generation resources using soft computing," Computational Intelligence in Electrical Engineering, Vol. 11, , 2020.
 Chiandussi. C, Codegone. M, Ferrero. M , "Comparison of multi-objective optimization methodologies for engineering applications," Computers & Mathematics with applications, Vol. 63, pp.912-942, 2012.
 Nasr. M, Nikkhah. S, Gevork. B, et al, " A multi-objective voltage stability constrained energy management system for isolated microgrids," Electrical Power and Energy Systems, Vol. 117, 2020.
 Soroudi. A, Rabiee. A, Keane. S, " Information gap decision theory approach to deal with wind power uncertainty in unit commitment," Electr Power System Research, Vol. 145, pp.137-148, 2017.
 Rabiee. A, Soroudi. A, Mohammadi-ivatloo. B and et al, "Corrective voltage control scheme considering demand response and stochastic wind power," IEEE Transaction on Power System, Vol.29, pp.2965-2973, 2014.
 Wang. X, Gong. Y, Jiang. C," Regional carbon emission management based on probabilistic power flow with correlated stochastic variables," IEEE Transaction on Power System, Vol.30, pp. 1094-1103, 2014.
 Soroudi. A, Siano. P, Keane. A, "Optimal DR and ESS scheduling for distribution losses payments minimization under electricity price uncertainty," IEEE Transaction on Smart Grid, Vol.7, pp. 261-272, 2016.
 Parvania. M,Fotuhi-Firuzabad. M, Shahidehpour. M, "Optimal demand response aggregation in wholesale electricity markets, "IEEE Transaction on Smart Grid, Vol.4, pp.1957-1965, 2013.
تعداد مشاهده مقاله: 379
تعداد دریافت فایل اصل مقاله: 249