Optimization of Power Balance Transaction Based on Renewable Energy Sources Using Artificial Salmon Tracking Algorithm for Modeling the Interconnected Grid Development

Article history: Received: 21 August, 2019 Accepted: 09 October, 2019 Online: 22 November, 2019 Since environmental requirements penetrate engineering processes to keep global warming and to reduce pollutant emissions, the system operation should be designed based on environmentally approach friendly. Operationally, the system processes are supported by energy suppliers to meet continuously energy transaction while clean and green energies are also targeted to keep the environmental conditions. In other word, the renewable energy source becomes an opportunity integration of a potential natural source inclusion into conventional energy producers. Technically, the energy balance is also important to optimize to get an optimal power portion during the operating period. These works are prioritized to search the balanced combination of the integrated energy mix composition. The works also present the novel computational intelligence to find out the energy portion using Artificial Salmon Tracking Algorithm. By considering technical requirements and tested on the integrated renewable energy plants, this algorithm is applied to optimize the economic power production. The works show dynamically the total power for feeding the energy consumption through local loads. The power production is also balanced in various combination portions of energy sources in accordance with the power demand as similar as discharged pollutants. Moreover, the computation has been obtained clearly for optimal solutions within 24 hours using a proposed algorithm.


Introduction
Electricity is one of the energy types which is used to be an important factor to support many technical processes. This energy has been affecting various activities caused by an easy conversion to other types for supplying many appliances [1], [2]. It uses not only industries but also in daily activities. Moreover, stable energy producers are needed to cover continuously an energy dispatching quota (EDQ) where energy users are located far away from the potential supplier [1], [3], [4]. In addition, this condition should be also supported by the reliable infrastructure to transmit and distribute the produced power to wide spreading loads at different locations. To cover the EDQ, the energy producers should consist of all activated power plants [5], [6]. Technically, the EDQ is also faced with the load demand changes in the period of time operation related to scheduled loads. Operationally, the EDQ is also used to divide each contribution in the integrated energy mix composition (IEMC) through a scheduling power output of energy plants associated with technical limitations, which is frequently operated based on an economic strategy [5], [7], [8]. In general, this strategy is supposed to the financial consideration based on the reducing fuel procurement fee and environmental compensation as the impact of the emission. Practically, this way concern in an integrated structure for the IEMC deal with generating sites for the generation systems, transmission systems, and distribution systems [4], [9], [10]. These conditions are associated with an energy stock to meet the total load and scheduled power outputs for the EDQ. Recently, a dynamic problem also becomes one of the main problems caused by demand changes within 24 hours considered reasonable composition energy sources for determining economically the operating cost and measuring energy producer contributions [7], [11]. Since global warming becomes an important aspect of the system operation, the application of air and environmental quality standards has increased, and pollution restrictions have also increased. Therefore, the operation of the system is increasingly fixed in implementing restrictions on pollution emissions due to the burning process of fossil fuels at energy production sites [12], [13]. Furthermore, to maintain the allowable emission limit, dynamically, the operation of the system is also increasingly considering releasing emissions as part of efforts to reduce its emission effects. Thus, the operation of the power plant must modify to find a more economical operating strategy by reducing pollutants, and reducing the cost of operating the electricity process. [5], [14]. The fossil fuel firing has contributed to the air quality condition through the disposal of contaminants. In other words, these combustions at thermal plants have also contributed to producing emissions, for example, CO, CO2, SOx, and NOx [5]. Moreover, the pollutant emission quota (PEQ) is also presented as a dynamic problem to meet hourly demand changes whereas the IEMC is subjected to the optimal composition for energy producers. To avoids complexity problems of both dispatching types, the EDQ and PEQ are determined together at the same time throughout single quota objective function (SQOF) to cover EDQ and PEQ problems [7], [15], [16]. Moreover, the SQOF becomes the main objective of the optimization problem under various technical constraints.
At present, to overcome these problems many ways are applied to obtain optimal conditions through optimization problems and by means of approaches that include mathematical programming and optimization techniques [10], [17]. In detail, these methods belong to traditional and evolutionary approaches. Nowadays, the optimization problem under various technical constraints becomes complex and huge models with excluding non-affecting small parameters for the system. To cover this condition, classical approaches have suffered to find out the solution where smart ways have been widely applied and are increasingly in demand as a way of calculating to solve various problems [18]. These approaches have an opportunity to wide implementations on numerical targets. Moreover, computational intelligence has many types associated with its procedures and inspirations to get the optimal solution. Furthermore, in its development, evolutionary computation is increasingly being developed and algorithms are arranged by imitating the behavior and mechanisms of flocks in nature. This was presented using optimization principles in accordance with natural mechanisms and structures to improve the performance of the classical approach [16], [19]. This paper presents a novel computational intelligence, artificial salmon tracking algorithm (ASTA), for determining the optimal solution of the SQOF based on the EDQ and PEQ problems. Technical limitations and environmental requirements are also applied to the conventional and hybrid energy systems to locate suitable solutions on hourly desired portions for the IEMC. In these works, the IEMC also covers an integrated renewable energy source (IRES) presented in a wind energy source (WES) and solar energy source (SES) for the 24 hours operating period. Both potential sources are installed at selected buses of the system to present the infrastructure development of the conventional system connection.

Energy Mixed Approach
Dynamically, the energy production due to all possible combinations of the hybrid system based on conventional energy sources (CES) and the IRES [20]- [23]. The hybrid system is commonly constructed using many types of generating units and it also uses the various voltage level systems to cover installed power plants depending on the CES and the IRES. Technically, this condition steers up to maintain a power output combination from joined energy producers based on the CES and the IRES which is decided to cover the total hourly demand. Moreover, this operation faces decreasing an emission and the operating fee. To integrate two dynamic problems with different targets for reducing pollutant discharge and decreasing running charge of the CES, a penalty factor is one of the important variables [5], [14]. Many previous studies reported that the emission has corresponded to fossil fuelbased power plants to meet the given load [7]. In these studies, the integration of the PEQ into the EDQ problem is not combined using an ascending order method of the penalty factor but it is approached using a new technique as given in a dynamic penalty factor (DPF) approach associated with the allowed emission discharge (AED), the total produced emission (TPE), and the over rate emission coefficient (OREC) [6], [7]. This method is an alternative approach for the penalty factor based on a dynamic computation in the processes. Mathematically, the DPF and OREC is discussed clearly in [24] as presented in Equations (1) to Equation (3).
where ORECz is the over rate emission coefficient of the z th iteration, TPEzs is the total produced emission of the s th generating unit of the z th iteration (kg/h), TAEzs is the total allowed emission of the s th generating unit of the z th iteration (kg/h), nGz is the number of generating units of the z th iteration exceeded the allowed emission, hz is a penalty factor set of the z th iteration ($/kg), hGzs is the individual penalty factor of the s th generating unit exceeded the allowed emission of the z th iteration step ($/kg), dhz is the dominant penalty factor of the z th iteration ($/kg), and rhz is the selected hGzs of the z th iteration for the highest TPEzs.
In particular, the penalty factor covers environmental effects from the conventional system operation using fossil fuel-based. Recently, the energy system operation (ESO) is developed using interconnected structures for connecting the CES and the SES which is located at different sites of the IRES [20], [25]. The ESO is also used to deliver power outputs at generation sites to the energy usages at faraway locations using interconnected topologies of the transmission and distribution systems. Operationally, the ESO is operated within 24 hours to provide a highly reliable electric energy based on the optimal energy dispatching [3], [26]. In general, the system operation is carried out within 24 hours, and this period considers the power output scheduled for 24 hours to meet energy consumption where the consumption changes load demand from this hour to the next according to technical and environmental requirements. In detail, this relation covers all possible combinations of suitable energy producers associated with reasonable power production [27], [28], which is presented as given in Equation (4) to Equation (6).
where t is period intervals of time (t=1, 2, 3, …, e), e is the total period, ng is the iterating number, PP is a power production, PCES is a power production of the CES, PWES is the power production of the WES, PSES is the power production of the SES, TD is the total demand, TLoss is the total loss, EP is the energy production, and T is a duration of present and past hours.
Moreover, to obtain the optimal solution, the SQOF is operated by entering several constraints as technical criteria during operation, including power balance, power load flow with shrinkage losses for lines, power capacity limits, fluctuating voltages, power transfer capability limits, ramp limits, and emission standards. The ESO covers the CES, WES, and SES which is optimized together to reach the optimal portion of the power producer combination

Optimizing Procedure
In this section, ASTA is compiled using its procedures based on the exploring and surviving steps. Algorithm procedures and parameters are very important to present computational abilities while searching the optimal solution [29], [30]. This algorithm is inspired by migration of Salmon fish. In general, the salmon migrating history is illustrated in Figure 1 in terms of spawning fish in freshwater, migrating to the ocean, and returning home from the ocean. During these phases, the Salmon will face many predators and obstacles. By considering Salmon's behavior, ASTA is constructed as given in Figure 2.
Computationally, these steps are illustrated in Figure 3 covered a transformation of the SQOF into programming sequences of the hybrid system. Moreover, the system is also evaluated using the Newton Raphson method for determining balanced energy performances as referred to in [31]. Furthermore, Figure 3 is also used to guide the procedures and hierarchies for optimizing the SQOF. This figure consists of two parts as given in the left and right sides. The left side is used to transform the mathematical model of the problem into a programming procedure as the sequencing processed flowchart. One other is used to optimize the problem based on the sequencing algorithm to search for the best solution. As illustrated in Figure 3, ASTA is programmed using pseudo-codes for searching the best solution using parameters.   In particular, Figure 4 shows the energy system model for applying ASTA and for optimizing the SQOF. The developed and standard models are very helpful to meet the problem caused by completed data of the real system and difficulty data collection on the existing operation [22], [32]. In these studies, the system consists of the CES, WES, and SES centers installed at different locations. In detail, the system covers five loads which are connected to Bus 1, Bus 2, Bus 3, Bus 4, and Bus 5. Moreover, conventional energy sources are installed at Bus 1, Bus 2, and Bus 3 whereas the WES is centered at Bus 1 and the SES is integrated to Bus 6. By considering this model and conditions, the system is optimized based on demand changes for 24 hours as provided in Table 1 with the period time operation is 24 hours. These demands also cover for the day and night peak loads.

Result
In this section, the IEMC is presented dynamically within 24 hours and it is optimized using ASTA. The 24 hours operation common approaches for the existing operation is based on all integrated energy producers [7], [27]. By considering ASTA's parameters detailed in Section 3, the energy mixed producer is determined optimally as given in Table 2 for the optimal power production. Technically, this table shows the individual power commitment of the energy sources of the CES while the WES and SES are considered free for the natural energy sources as given in Table 3. From this table, it is known that this table informs the scheduled power production for 24 hours determined totally in 9,990.21 MW with the power fluctuation and contribution are illustrated in Figure 5 and Figure 6. These aspects are inlined with many previous works that the system has been delivered in variable portions associated with demand changes at the energy consumers [33].
According to Table 2 and Table 3, it is also known that all conventional generations have produced power outputs in different capacities. These different capacities show the generated participation in the system for supporting the provided energy for the load demand as reported in [14], [15]. The energy contributors cover all centers as detailed in both tables. In total, Center 1 takes a role in the highest power production of around 5,273.98 MW. Furthermore, the lowest contributor is supported by the CES Center 3 in 1,523.93 MW. By considering the IRES, Table 3 presents all penetrations within 24 hours for the SES and the WES. This penetration is very important to measure a renewable energy inclusion into the existing system with a certain portion of the participants to control the total energy production [10]. These penetrations cover hourly operations, totally, the system is penetrated in 1,450.00 MW of the SES and 2,522.93 MW of the WES. In detail, the system provides around 13,870.21 MW with the discharged pollutant of the CES is listed in Table 4 for all operating times. This emission is categorized in the produced emission, permitted emission, and over emission. As given in Table 4, the total emission is produced in 19,662.34 kg where the pollution is allowed around 8,491.65 kg based on standard emission. In addition, the system has 11,170.69 kg of the over emission during existing energy sources. It means that the system should be maintained to keep the emission level. In particular, an energy balance is one of the important parameters in the system operation while the energy should be used in the final users at various types of appliances based on the power capacities [11]. In this case, the energy consumption has corresponded to individual power productions. This aspect is also used to present individual contributors to the power unit commitment. In these works, individual performances of the energy producer unit are presented in Figure 7 and Figure 8 covered hourly power production changes and hourly power fluctuations. In total, the IEMC produces 9.990.21 MW with the power produced performances as given in Figure 7 and Figure 8. Figure 7 shows the fluctuation of power productions for the CES from the present to the next hours associated with the own scheduled power production. The detailed contributor to the power plants is presented in Figure 7. This figure informs that the highest contributor comes from G1 of the CES Center 1. Moreover, the IEMC also emits the total over emission of around 19,662.34 kg whereas the IRES penetrates are given in Figure 8. This figure illustrates the hourly penetration to the system covered in 24 hours for the WES and the SES. In addition, the CES is performed in Figure 7 for the hourly power production during an existing system to supply the load center.   By considering the 24 hours operation, the power balance is detailed in Figure 9. This figure shows all participation of the energy producers balanced to the total hourly demand. It is also known that the main power producers are covered in conventional energy producers. In addition, the IRES is operated to support this optimal condition which is linked to the hourly composition. In these works, another problem is given in the power loss as depicted in Figure 10 designed for the 24 hours operation. The system has fluctuated power losses for the 24 hours operation associated with day and night loads.

Conclusion
As stated earlier, the composition of the energy mix is integrated by considering integrated renewable energy sources and various technical requirements, as well as environmental constraints. The formula is outlined in a computational problem that is searched using the Artificial Salmon Tracking Algorithm to get the optimal composition for 24 hours. The calculation results show that dynamically the part of the energy produced is affected by IRES as given in solar and wind energy sources for 24-hour operations. In addition, ongoing conventional energy producers contribute with different capacities to support individual commitments to the production of power units. Furthermore, from these works, the implementation of real algorithms for large systems is suggested for future work in line with the placement of distributed renewable energy.