
==== Front
Heliyon
Heliyon
Heliyon
2405-8440
Elsevier

S2405-8440(24)12085-3
10.1016/j.heliyon.2024.e36054
e36054
Research Article
Expansion planning of hybrid electrical and thermal systems using reconfiguration and adaptive bat algorithm
Abbasi Ali Reza a
Zadehbagheri Mahmoud Ma.zadehbagheri@iau.ac.ir
b⁎
a Department of Electrical, Faculty of Engineering, Fasa University, Fasa, Iran
b Department of Electrical Engineering, Yasuj Branch, Islamic Azad University, Yasuj, Iran
⁎ Corresponding author. Ma.zadehbagheri@iau.ac.ir
11 8 2024
30 8 2024
11 8 2024
10 16 e3605421 6 2024
30 7 2024
8 8 2024
© 2024 The Author(s)
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).
_ This study introduces a comprehensive model for the concurrent expansion planning of various energy systems and their associated equipment. The need to reduce network costs, emissions, losses, and feeder loading, as well as to enhance network reliability and voltage profile, mandates the utilization of proper multi-objective planning models that respect all network constraints. The introduced framework includes units for generating both electrical and thermal energies. The model leverages conventional expansion alternatives such as the installation of new lines, network reconfiguration, rewiring, and the addition of new thermal and electrical generating units to the network. Expansion planning involves determining the optimal time, location, and type of new installations to meet future energy demands while minimizing costs and emissions. Reconfiguration refers to altering the network topology to improve reliability and reduce losses. The proposed expansion planning is formulated as a discrete, nonlinear, and non-convex optimization problem, which is solved using the Self Adaptive Learning Bat Algorithm (SALBA). This algorithm improves convergence speed and increases the diversity of the search population, enhancing the likelihood of finding the global optimum. Numerical simulations of the proposed methodology on two modified standard IEEE test systems corroborate the efficacy and feasibility of the suggested approach. Key innovations include the comprehensive modeling for concurrent expansion planning, the use of an advanced optimization algorithm, and a focus on reducing costs, emissions, losses, and feeder loading while enhancing network reliability and voltage profile.

Index Terms

Expansion planning
Emission
Reliability
Self-adaptive learning bat algorithm (SALBA)
==== Body
pmc Nomenclature

Symbol and description	
CInv	Cost of Investment	
CRel.	Cost of reliability	
CLoss	Cost of energy loss	
CSub	Cost of Generation	
NP	Number of power units	
NC	Number of CHP units	
Nh	Number of heat units	
NL	Number of lines	
Pjh	Electric power equivalent thermal power at the jth bus	
PiP	Net injected active power of the power unit at the ith bus	
PjC	Net injected active power of the CHP unit at the jth bus	
NCus	Total number of customers served	
λi	Average failure rate of the ith component	
ρi	Current compensation coefficient	
|Ii|	Absolute current amplitude of the ith feeder	
|Inewi|	Current magnitude of ith branch after reconfiguration	
|Ioldi|	Current magnitude of ith branch before reconfiguration	
Ui	Annual outage time of the ith component	
δi	Angle of the voltage at the ith bus	
Yij	Amplitude of admittance between the ith and jth buses	
θij	Angle of branch between the ith and jth buses	
CiElec.	Cost of a power unit that can be added at the ith bus	
CiCHP	Cost of a CHP unit that can be added at the ith bus	
CiTher.	Cost of a heat unit that can be added at the ith bus	
Cline	Cost of a circuit that can be added or upgraded between ith and jth bus	
T	Planning horizon time	
Pij,minLine	Minimum active power between ith and jth bus	
PijLine	Expected value of active power between ith and jth bus	
Pij,maxLine	Maximum active power between ith and jth bus	
Pd	Electric power demand of system	
HCj	Injected heat of the CHP unit at the ith bus	
Hkh	Injected heat of the heat unit at the kth bus	
Hd	Thermal demand of the system	
Hkhmin	Minimum thermal outputs of the kth unit	
Hkhmax	Maximum thermal outputs of the kth unit	
τi,kj,δk	CO2 emissions coefficients	
αi,βi,γi,ζi,λi	Emissions coefficients for the power-only units	
θj,ηj,πk,ρk	Emissions coefficients for the CHP units and heat-only units	
PjCmin(HjC)	Minimum power limit of CHP unit j which are functions of generated heat (HjC)	
PjCmax(HjC)	Maximum power limit of CHP unit j which are functions of generated heat (HjC)	
Vmax	Maximum voltage magnitude	
Vmin	Minimum voltage magnitude	
fi(X)	The ith objective function	
gi(X)	ith inequality constraint of MOP	
hi(X)	ith equality constraint of MOP	
Neq,Nueq	Number of equality and inequality constraints	
fimin/fimax	Lowest/Highest limit of ith objective function	
Pipmin	Minimum injected active power of the power unit at the ith bus	
Pipmax	Maximum injected active power of the power unit at the ith bus	
Swi	The ith sectionalizing switch	
List of abbreviations	
EEP E	nergy Expansion Planning	
CHP C	ombined Heat and Power	
DFR D	istribution Feeder Reconfiguration	
MOP M	ulti-objective Optimization Problem	
SALBA S	elf-Adaptive Learning Bat Algorithm	

1 Introduction

1.1 AProblem description and motivation

Increased efficiency due to the integration of combined thermal and electrical networks has motivated the implementation of these energy systems [1]. As a result, new operative tools need to be devised to address the Energy Expansion Planning (EEP) problem in systems containing combined thermal and electrical energy generating units. EEP involves answering four key questions about generating units and other network elements: the time, location, type, and method of integrating an element into the system over the entire planning horizon [2]. This integration must respect all network limits and satisfy the total load demand. Moreover, recent concerns related to fossil fuel resources, such as greenhouse gas emissions, suggest that planning goals should include the reduction of SO2, CO2, and NOx emissions. The introduction of Renewable Energy Resources (RERs) has also raised optimism about decreasing total environmental emissions compared to conventional power plants [3]. Consequently, the secure integration of Distributed Generating (DG) units into the network, considering economic and environmental concerns, has gained significant attention. Among various types of DGs, this paper focuses on three: power-only units, heat-only units, and Combined Heat and Power (CHP) units.

Through the expansion of energy units in future networks, enhancing system reliability should be considered, as reliability plays a vital role in keeping costs within a reasonable budget. As a result, this paper employs the Distribution Feeder Reconfiguration (DFR) strategy as a means of enhancing reliability. DFR aims to improve conditions by increasing reliability and reducing losses in the network through changes in network topology [4,5]. Within this context, the motivation of this research is to introduce a tool for the concurrent development of thermal and electrical energy systems. This tool facilitates decision-making for the installation of new units that generate electrical and thermal power, the addition of new lines, and the reconfiguration of the network over a specified timeframe. The goal is to ensure compliance with both equality and inequality constraints while reducing overall costs and emissions, and simultaneously enhancing the system's dependability.

1.2 BLiterature review

Earlier studies have primarily focused on enhancing the expansion of systems that utilize a single type of energy. For instance, Ref. [6] presents a development planning approach to analyze the distribution of heat-only energy demand and supply across Europe. By using a comprehensive thermal atlas, it identifies high-demand areas and aligns them with potential supply zones. Ref. [7] examines the potential risks associated with the expansion of liquefied natural gas infrastructure in Germany, highlighting that such expansion could increase dependency on gas and slow down the transition to renewable energies. Ref. [8] employs an adaptive robust optimization framework for planning the expansion of natural gas infrastructure in Europe, aiming to identify projects that can withstand fluctuations in gas demand and supply, thereby maintaining system stability. Ref. [9] introduces a multi-objective approach to planning the expansion of power generation and transmission infrastructure, incorporating capacitor bank allocation and demand response programs. Additionally, references [10,11] explore expansion planning in electricity distribution and transmission networks. The wide-ranging benefits of planning for energy expansion have garnered the focus of many investigators. Consequently, a variety of objectives have been specified. Numerous research efforts have been undertaken with the goal of reducing costs [12,13]. A novel optimization method based on squirrel and cat algorithms is introduced for optimal expansion planning of distribution systems. The approach aims to reduce costs and improve the efficiency of the power distribution system [12]. Several strategies have been formulated with the aim of improving reliability [14,15]. Ref. [14] evaluates the reliability and economic aspects of generation expansion planning by incorporating the variability of wind energy sources. It aims to optimize the balance between cost and reliability in power systems that integrate significant amounts of wind energy. Ref. [16] introduces an index to enhance the voltage stability via planning for network expansion. In addition to the aforementioned planning concerns, various studies have focused on the reduction of losses in power grids as an objective function, as indicated in Refs. [17,18]. A stochastic mixed-integer convex programming model is presented for long-term distribution system expansion planning, focusing on mitigating greenhouse gas emissions in Ref. [19]. The model aims to optimize the reinforcement of the distribution network while minimizing environmental impact. Moreover, the capability of DFR in making the network situation better through EEP is referenced in Ref. [20]. It should be noted that some studies solve the problem by considering only one of the mentioned criteria, while others use multiple objective functions. These studies have been developed for systems that are specific to a single type of energy, such as district heating, natural gas, and electricity systems.

The interdependence of electricity, gas, and heat networks has recently led to increased attention on multi-carrier energy systems. In Ref. [21], an optimal operation model is presented for coordinated multi-carrier energy hubs within integrated electricity and gas networks. The model aims to enhance efficiency and reliability by optimizing the interactions between these interconnected energy systems. Ref. [22] presents a mixed-integer linear optimization model for the integrated planning of municipal multi-energy carrier systems, focusing on achieving CO2 reduction targets. The model aims to design cost-effective energy systems that incorporate grid infrastructures and energy conversion units while minimizing greenhouse gas emissions. Ref. [23] presents a robust model for the expansion planning of integrated energy systems that include energy hubs, with a focus on ensuring reliability and efficiency in the planning process. Refs. [24,25] explores the combined expansion planning of gas and electricity networks, aiming to optimize the integration and coordination of both systems. Reference [26] reviews advanced models for the expansion planning of integrated power, natural gas, and hydrogen systems, with a focus on cost minimization and optimal facility placement.

Reference [27] discusses the coordinated expansion planning of natural gas and electric power systems, highlighting the importance of detailed representation for making optimal investment decisions. Refs. [26,27] both aim to improve system reliability and efficiency while managing uncertainties and reducing losses. Refs. [[28], [29], [30]] discuss optimal expansion planning of electricity and natural gas distribution systems with a focus on energy hubs and reconfiguration. Ref. [28] addresses the co-planning of reconfigurable electricity and natural gas distribution systems incorporating energy hubs to enhance operational performance. Ref. [29] focuses on resilience-constrained expansion planning of integrated power, gas, and heat distribution networks with reconfiguration to handle rare events. Ref. [30] presents a hierarchical robust expansion co-planning model for multi-energy distribution systems considering reconfiguration. Table 1 illustrates the classification of the mentioned literature review.Table 1 Review of literature.

Table 1#
Ref	Main concentrations in the objective function	Reconfiguration	Adding new	Rewiring	Objective Function	Optimization Method	
	Cost	Loss	Reliability	Feeder loading	Emission	Units	Load points	Line	Single	Multi	
[21]	✓	✓			✓		✓	✓	✓		✓		MILP		
[22]	✓		✓		✓		✓		✓			✓	MILP		
[23]	✓		✓				✓		✓		✓		MILP		
[24]	✓				✓		✓		✓		✓		MILP		
[25]	✓		✓		✓		✓	✓	✓			✓	MILP		
[26]	✓						✓		✓			✓	MILP		
[27]	✓		✓				✓		✓		✓		MILP		
[28]	✓		✓		✓	✓	✓		✓			✓	MIQP		
[29]	✓				✓	✓	✓		✓			✓	MILP		
[30]	✓		✓			✓	✓		✓			✓	MILP		
PM*	✓	✓	✓	✓	✓	✓	✓	✓	✓	✓		✓	SALBA		

1.3 CResearch gap

The main research gaps and shortcomings are listed below.• Separate Consideration of Parameters: Previous studies have considered the parameters and objectives discussed in this paper, but they have done so separately rather than in a unified manner.

• Limited Diversity in Optimization Methods: Most studies have utilized MILP methods, showing little diversity in optimization approaches. Exploring other methods, such as evolutionary algorithms, could improve results.

• Focus on Cost and Reliability: The majority of studies concentrate on cost reduction and reliability enhancement. There is a need for more research on other aspects, such as system sustainability and flexibility.

• Inadequate Multi-Carrier Energy Studies: Studies on multi-carrier energy systems have not adequately addressed critical aspects such as reconfiguration, the creation of new load points for network expansion, rewiring, feeder loading, reliability, and environmental issues.

• Integrated Model Requirement: There is a need for an integrated model that addresses the parallel expansion planning of energy networks, including both thermal and electrical energy systems.

• Insufficient Attention to Environmental Impacts: Some studies do not adequately address environmental impacts. Further investigation in this area could contribute to the development of more sustainable systems.

1.4 DProcedure

As explained previously, the suggested EEP is a sophisticated problem with multiple objectives, requiring strong and robust algorithms to determine the best overall solution. For this reason, an improved Bat Algorithm (BA) is used for solving the suggested EEP. The BA emulates the bats behavior in chasing prey or searching for food. This algorithm is conceptually straightforward, making its implementation easy. However, the initial version of the algorithm is limited by the homogeneity of its population. To address this, various self-adjusting enhancement measures were incorporated into the original algorithm to improve its convergence speed and expand the diversity of its search population. Distributing the population throughout the entire search area increases the likelihood of discovering the global optimum. The strength and effectiveness of the proposed method were demonstrated by testing it on two distinct IEEE standard test systems.

1.5 Contributions

The innovations and main highlight of the paper can be summarized as follows.• Comprehensive Modeling for Concurrent Expansion Planning of Various Energy Systems: This paper introduces an all-encompassing model for the concurrent expansion planning of various energy systems and their associated equipment. The model includes units for generating both electrical and thermal energies and utilizes conventional expansion alternatives like the installation of new lines, network reconfiguration, rewiring, and adding new thermal and electrical generating units and new load points to the network.

• Use of Modified Advanced Optimization Algorithm: The expansion planning problem is formulated as a discrete, nonlinear, and non-convex optimization problem, solved using the Self Adaptive Learning Bat Algorithm (SALBA). This algorithm improves convergence speed and increases the diversity of the search population, enhancing the likelihood of finding the global optimum.

• Focus on Reducing Costs, Emissions, Losses, and Feeder Loading: The proposed model aims to reduce network costs, emissions, losses, and feeder loading while enhancing network reliability and voltage profile.

• Utilization of Distribution Feeder Reconfiguration (DFR): The paper employs DFR as a strategy to improve system reliability by determining the best network topology to increase reliability indices. DFR provides better conditions in terms of increasing reliability and reducing losses by changing the network topology.

• Validation of the Proposed Model: Numerical simulations of the proposed methodology on two modified standard IEEE test systems corroborate the efficacy and feasibility of the suggested approach.

These innovations contribute to a comprehensive and efficient approach to energy system expansion planning, addressing both economic and environmental aspects.

1.6 Paper organization

The subsequent sections of the paper are structured as follows: Section 2 formulates the objective functions and constraints relevant to the issue at hand. Section 3 discusses the capability of DFR in improving reliability indices. Sections 4, 5 describe the original BA with proposed modifications and the application procedure, respectively. Section 6 delineates the simulated test systems, accompanied by the numerical outcomes and their corresponding discussions. The final section summarizes the key points and contributions of the paper.

2 Problem formulation

2.1 AMathematical representation of EEP

The proposed EEP is a non-linear, large-scale optimization problem characterized by non-convex and non-smooth solutions. The problem exhibits multiple modes and several local optima. As discussed earlier, the proposed EEP incorporates a range of objective functions, including minimizing total cost and minimizing emissions. These objectives are formulated as follows.

2.2 Minimizing total cost

Several sources impose costs on the network and must be considered when calculating the total network cost. These sources include system investment costs, the cost of the expected value of undistributed energy, energy loss costs, and the cost of energy imported from the upstream network. The total cost is the sum of these costs as follows:(1) FCost=CInv+CRel+CLoss+CSub

As obvious the total cost consists of four different terms that are explained in the following.• Investment Cost: As mentioned previously, three types of generating units are considered: heat-only, power-only, and CHP units. The investment cost (CInv) includes the costs associated with the installation of these units, as well as the costs related to the installation and/or reinforcement of lines, as follows:

NC, NP, and Nh represent the number of CHP, power-only, and heat-only units, with indices i, j, and k referring to them, respectively. Additionally, the amount of generated heat and power is denoted by H and P, and NL represents the number of lines. The terms of the above formulation can be calculated using the following equations:(4) Ci(PiP)=CiElec.*T*PiP(3)Cj(PjC,HjC)=CiCHP*T*(PjC+Pjh)Ck(Hkh)=CiTher.*T*Pkh(5)CL(Lij)=Cline*nij*Lij(6)

• Cost of Expected Value of Non-Distributed Energy: ECOST is considered the system reliability index to specify the satisfactory level of dependability for consumers. It determines the economic evaluation of redundancy allocation and network reinforcement, as well as the identification of weak points in a system. Once the weak points in the system are detected, the operator should devise proper operation policies and maintenance schedules. Calculating ECOST requires evaluating load-point reliability indices. These indices consist of average outage time, annual system outage time, and average failure rate, which are computed as follows:

(7) rS=UsλS;US=∑iλiri;λS=∑iλi

For the ith bus, the calculation of ECOST is as follows [20]:(8) ECOSTi=LiCiλi

where Ci denotes the interruption cost, which is determined by a nonlinear function based on the duration of the interruption. The Ci value, related to the time of interruption, is determined using the CCDF. To assess the costs associated with consumer interruptions, it is necessary to compute the system's load-point reliability indices for each switching configuration separately. In this case, a standard CCDF is employed to convert the energy losses experienced by customers into their monetary equivalents. Fig. 1 shows a typical CCDF [31].Fig. 1 Typical CCDF.

Fig. 1

For nbus load points in the distribution feeder, the calculation of total ECOST is as follows:(9) CRel=∑i=1nbusECOSTi=∑i=1nbusLiCiλi

It is noted that CCDF specifies the value of bus interruption cost [31].• Cost of Energy loss: To emphasize the importance of reducing energy loss in the system, a cost corresponding to its quantity is assigned according to the equation below:

(10) CLoss=T×CPloss×∑i=1NbrRi|Ii|2

In which resistance and current of ith branch are denoted by Ri and |Ii| respectively and Nbr shows number of branches.• Cost of Generation: Csub represents the yearly cost of imported energy. The network cost can be ascertained in the following way:

(11) Csub=price×Psub×T

(12) Psub=∑j=1n|Vsub||Vj||Ysub,j|cos(θsub,j−δsub−δj)

2.2.1 Minimizing emission

This paper considers three sources of emissions: NOx, SO2, and CO2. The first two gas emissions are functions of the output power of the units. For power-only units, the function is exponential, while for heat-only and CHP units, it is linear, as follows:(13) EmTypeofDGNox,So2=(∑i=1Np{10−2(αt+βtPt+γiPi2)+ζie(λi,Pi)}+∑j=1Nc{(θj+ηj)Oj}+∑k=1Nh{(πK+ρK)THK})*T

In the above formulation, the output power of CHP and power-only units are denoted by Pt and Oj, respectively. The emission factors for units that only generate power are represented by αi, βi, γi, ζi, and λi. The emission coefficients for CHP units and heat-only units are denoted by θj, ηj and πk, and ρk, respectively. Similarly, CO2 emissions are represented by a linear function of the output power of the units as follows:(14) EmTypeofDGCO2=(∑i=1NpτiPi+∑i=1NckjOj+∑k=1NhσkTHk)*T

CO2 emissions coefficients are shown by τi, kj, σk. The average emissions due to the power plants of the main grid (Emgrid) are calculated as follows:(15) Emgrid=(ENoXgrid+ESo2grid+ECO2grid)*Psub

where ENoXgrid, ESo2grid and ECO2grid represent the average rates of NOx, SO2, and CO2 emissions of all committed power plants. The total emission of the network is then calculated by summing the above functions as follows:(15a) FEmission=EmTypeofDGNox,So2+EmTypeofDGCO2+Emgrid

• Constraints

There are both equality and inequality constrains in the proposed formulation that need to be respected during EEP.

2.2.2 Equality constraints

• Preserving radial structure of Distribution system

Generally, distribution networks are configured in a radial pattern, which offers advantages such as robust protection, ease of design, and affordability. Therefore, it is crucial to maintain the radial structure of these distribution systems. During optimization, if a loop emerges, a sectionalizing switch within the loop is activated to preserve the system's radiality. In this research, the 'depth-first search algorithm' is employed to identify any loops present in the system [32]. After identifying a loop, a branch within that loop is selected and opened at random. Subsequently, the 'depth-first search algorithm' is applied repeatedly until the network achieves a fully radial state.• Power flow constraints:

(16) PK=VK∑j∈NVj(gjkcosθjk+bjksinθjk)

(17) QK=VK∑j∈NVj(gjksinθjk−bjkcosθjk)

In these equations, the conductance and susceptance between the jth and kth buses are denoted by gjk and bjk respectively and the impedance angle of the respective branch is shown by θjk.

2.2.3 Inequality constraints

• Limits on heat production and consumption:

(18) ∑j=1NCHCj+∑k=1NhHkh≥Hd

Hd signifies the system's demand for heat.• Limits on capacity of conventional units

(19) Pipmin≤Pip≤Pipmaxi=1,…,NP

Here Pipmin and Pipmax show the min and max output powers of the ith unit.• Limits on capacity of CHP units:

(20) PjCmin(HjC)≤PjC≤PjCmax(HjC)j=1,…,NC

(21) HjCmin(PjC)≤HjC≤HjCmax(PjC)j=1,…,NC

Limits on generation of heat units:(22) Hkhmin≤Hkh≤Hkhmaxk=1,…,Nh

Hkhminand Hkhmax show the lower and upper bounds of thermal outputs for the kth unit respectively.• Limits on bus voltage:

(23) Vmin≤Vk≤Vmaxk=1,2,…,Nbus

Vk is the voltage magnitude at the kth bus. Vmin and Vmax represent the lower and upper bounds on permissible bus voltages respectively.• Limits on Line loading:

(24) PFk≤PFkmax

here PFk and PFKmax are the amount of power flow in the kth branch and its maximum level.

3 Effect of Dfr on reliability

An analysis of failure statistics reveals that a significant number of interruptions in distribution systems are due to overhead and underground cable failures [33]. The resistive ohmic loss in underground cables generates heat, leading to insulation failure, which is primarily made of cross-linked ethylene-propylene or polyethylene rubber [34]. This failure occurs when the generated heat exceeds the tolerable levels defined for the cable. The temperature increases as the ohmic loss rises, which is a direct function of the square of the current magnitude through the cable. As the temperature in the cable rises, moisture can infiltrate the cable, causing damage to the dielectric insulation. The sensitivity of these insulations to moisture is the primary cause of breakdowns under overload conditions [33]. When the current flow through the cable increases and the temperature rises in overhead lines, the ground clearance is reduced. This reduction in ground clearance leads to increased electrical breakdowns and a higher failure rate [35]. All these conditions result in decreased system reliability. DFR has been shown to be capable of changing the flow and amount of current through the feeders of a system by opening or closing some remote switches located in the network. DFR alters the flow of power across feeders at almost no cost. It also provides the optimal structure for the system to minimize active loss. Reductions in active loss lead to decreases in the system's maximum loadability. Consequently, DFR can effectively enhance the system's reliability by reducing the current in overloaded lines and decreasing losses.

Altogether, DFR can be implemented as a strategy to reduce the failure rate and subsequently enhance reliability. To represent the mathematical formulation of the problem, the initial failure rate of a system line prior to reconfiguration can be indicated as λiinit, and the optimal failure rate achieved by complete active/reactive current compensation can be represented as λibest. The current failure rate is computed as follows [36]:(25) λinew=ρi(λiinit−λibest)+λibest;ρi=|Inewi||Ioldi|

It should be particularly noted that the aforementioned method for reducing failure rates is generally most beneficial for feeders whose current levels are near their rated capacities. Therefore, feeders with low currents were ignored while performing this strategy. The criteria for selecting lines for this strategy rely on their electrical current; thus, only lines that carry a current greater than 80 % of their nominal capacity are chosen.

4 Optimization procedure

The basics of multi-objective optimization are discussed here, and an interactive fuzzy satisfying approach is introduced to solve the suggested EEP. In this method, the objective functions are converted into a min-max problem and solved as a single objective problem. Additionally, the proposed evolutionary algorithms are explained at the end.

4.1 AMulti-objective optimization framework

In a Multi-objective Optimization Problem (MOP), multiple objectives compete with one another. The problem should be solved to find a solution that satisfies all the objectives to some extent. A typical MOP problem can be represented as follows:(26) minF=[f1(X),f2(X),…,fn(X)]Ts.t.gi(X)≺0i=1,2,…,Nueqhi(X)=0i=1,2,…,Neq

Here, fi(X) denotes the ith objective function, while equality and inequality constraints are represented by gi(X) and hi(X), respectively. The number of equality constraints is Neq and the number of inequality constraints is Nueq. Fig. 2 illustrates the flow chart of the multi-objective optimization model.Fig. 2 Multi-objective optimization planning model.

Fig. 2

4.2 BFuzzy Modelling for normalization of objective functions

To solve the multi-objective problem, the objective functions are normalized using a fuzzy technique. This method converts the impreciseness of the objectives into exact values through fuzzy sets and a membership function (μfi), as shown in Eq. (28). The membership function, which takes values between 0 and 1, depicts the level of adaptability of each objective. Due to the varying characteristics and behaviors of the examined objective functions, a fuzzy membership function is used to express all the objective functions uniformly [42].(27) μfi(X)={0fi≻fimaxfimax−fifimax−fiminfimin≤fi(X)≤fimax1fi≺fimin

The minimum acceptable level (fimin) and the maximum level (fimax) of the ith objective function are determined through the individual optimization of that function.

4.3 CInteractive fuzzy satisfying technique

The preferred solution is extracted from the collection of optimal solutions using an interactive fuzzy satisfaction technique [37]. First, separate membership functions are assigned to all objectives, and desired levels of preference (reference membership values) are specified for each function. The reference value (μri) determines the achievement level of each objective and takes a value between 0 and 1. If an objective is more important, higher values of μri are assigned to it.

Once the membership functions and reference values are set, the following min-max problem is solved to obtain the desired solution:(28) minX∈φ(maxi|μri−μfi(X)|)

Ω is the set of non-inferior solutions. It should be noted that if there is no biased preference among the objectives, the following problem will be solved instead:(29) maxP∈Ω(minn=1,…,n0|μfn|)

4.4 Self-adaptive Learning Bat algorithm

4.4.1 Original BA

Bats locate nourishment or prey using echolocation. Echolocation involves determining the location of objects by interpreting their responses to specific subsonic signals. The bat dispatches a signal in the surrounding environment and remains on the lookout for the echoed signal to return from its prey or food. This concept has been modeled to create a metaheuristic optimization technique known as the Bat Algorithm (BA) [38].

To emulate the algorithm, it is important to understand that every bat uses echolocation to measure distances. A bat located at position Xm or moving at a velocity of Velm emits a signal. This signal has a frequency, denoted as fmin, a wavelength, represented by λ, and an amplitude, which indicates its loudness, denoted as A0. The echoes returned from various objects differ, allowing bats to distinguish between prey and other objects. Adjustments are made to the wavelength and its emission rate, denoted as r. Simultaneously, the amplitude of the signal for the ith bat, indicated by Ai, is gradually reduced from a high initial value, A0, to the minimum value, Amin. The process begins with the creation of an initial set of bats, randomly distributed within the search domain. The evolutionary process involves modifying the positions of the bats based on subsequent adjustments.(30) Velmk+1=Velmk+(Gbestk‐Xmk)fmm=1,…,NBatXmk+1=Xmk+Velmk+1m=1,…,NBat

(31) fm=fmmin+rand(.)(fmmax−fmmin)m=1,…,NBat

To imitate the random patterns of bat flight, a random number (rand) is generated and compared to rmk. If rand exceeds this threshold, a new solution is created near the bat's current location in the following way:(32) Xnew=Xold+εAk

After the bats' locations are improved through the changes mentioned earlier, a new random entity, Xnew, is created if its signal rate, ri, exceeds a randomly determined number. This new solution will be added to the existing population if it meets the following condition.(33) [rand≺Ai]&[f(Xi)≺f(Xgbest)]

As previously stated, the amplitude values of signals produced by bats diminish progressively, as expressed by the following formulation:(35) Aik+1=αAik(34)rik+1=ri0[1−exp(−γk)]

4.4.2 Self-adaptive learning (SAL) mechanism

The initial version of the BA has certain limitations, such as the risk of getting trapped in local optima and a slow convergence rate towards the best solution. Adaptability is a key feature of any heuristic optimization technique. Therefore, a modification phase is added to the original BA to enhance this feature. The SAL method increases diversity within the population, helping the algorithm escape local optima. During this phase, each entity encounters four different modification strategies, each with varying probabilities, and must adaptively choose the most suitable modification based on its probability. The likelihood of selecting a strategy depends on how well it optimizes the fitness function. The strategies are as follows:

Mutation strategy 1: This approach is intended to explore various directions near the optimal solution, thereby enhancing the algorithm's overall search efficiency. Using the optimal solution, denoted as Gbestk, a new entity is created as follows:(36) Xm,1k=Xmk+round(1+rand(.))(Gbestk−Meank)m=1,…,N1k

Mutation strategy 2: This approach, inspired by natural processes, aims to enhance the algorithm's local search capabilities, leading to a faster. The core idea of this method is to seek the optimal solution while avoiding the least favorable one. This is achieved by creating a unique individual as follows:(37) Xm,2k=Xr1k+rand(.)(Gbestk−Worstk)m=1,…,N2k

Mutation strategy 3: The third method enhances the variety among the bat population through the creation of an individual based on mutation in the following way:(38) Xm,3k=Xr1k+rand(.)1(Xr2k−Xr3k)+rand(.)2(Xr4k−Xr5k)m=1,…,N3k

Mutation strategy 4: Analogous considerations that applied to the former approach are pertinent to this one as well. In this method, the newly devised individual, which relies on mutation, is formulated as follows:(39) Xm,4k=(Xr1k+Xr2k+Xr3k3)+(R1−R2)(Xr2k−Xr1k)+(R3−R2)(Xr2k−Xr3k)+(R1−R3)(Xr3k−Xr1k)m=1,…,N4k

where R1=F(Xr1K)/R′,R2==F(Xr2K)/R′,R3==F(Xr3K)/R′,R′=|F(Xr1K)|+|F(Xr2K)|+|F(Xr3K)|,N1it,N2it,N3it, and N4it represent the quantities of bats that migrate in accordance with mutations 1, 2, 3, and 4 during the k-th iteration. Vectors Xr1k,Xr2k,Xr3k,Xr4k,andXr5k are selected at random from the population. The strategy that has a greater likelihood will be chosen for execution from the given options If the new test individuals created by this strategy improve the fitness function, they will replace the ith bat. Initially, each strategy has an equal probability of 0.25. This probability is updated after LP iterations based on the success rate of the trial solutions generated by the ath mutation strategy (SRak), as described below:(40) Probak=SRak∑a=14SRak;a=1,…,4

The rate of success is determined by the nsait count of trial solutions yielded by the a-th mutation strategy and the nfait count of trials that did not succeed, as calculated in the subsequent manner:(41) SRak=∑it=k−LPk−1nsait∑it=k−LPk−1(nsait+nfait)+δ;a=1,…,4;k≻LP

The fixed parameter δ is established at 0.01 to prevent a zero value for the success rate. As per equation (40), the total probabilities for choosing strategies sum up to one. Every member of the bat group might execute the strategy deemed most likely using a Roulette Wheel Mechanism (RWM) during each cycle to attain the optimal solution. The procedure for selecting the method is outlined in the following description.

5 Application of optimization technique to Eep

To implement the suggested approach for the multi-objective EEP, follow these steps.Step 1. Initialize the Algorithm:

• Define the input data and generate the initial bat population randomly. It is important to note that the bats represent feasible solutions for the problem.

• The required data includes:

Network Data: Sectionalizing switches, tie switches, nominal voltage, branch data, bus data, etc.

Reliability Data: Repair rate and failure rate of network components, etc.

Algorithm Parameters: Population size, termination criterion, adjusting parameters, etc.

Step 2. Evaluate All Objective Functions:

• Execute the power flow analysis by considering the decision variables and the constraints associated with the problem for each solution.

• This process determines the values of the objective functions and their corresponding membership functions.

Step 3. Apply the Interactive Fuzzy Satisfaction Technique:

• Implement the interactive fuzzy satisfaction technique across the entire population.

• Store the results in the repository.

Step 4. Identify the Best and Worst Solutions:

• Select the best and worst solutions from the population based on the evaluation criteria.

Step 5. Update the Bat Population:

• Move the bat population to new positions as described in section IV.D.1.

Step 6. Execute the Self-Adaptive Modification Phase:

• Perform the proposed self-adaptive modification phase as outlined in section IV.D.2.

Step 7. Update the Repository with New Non-Dominated Solutions:

• Use the new non-dominated solutions to update the repository.

• Check the repository size as described in Section IV.B–C.

Step 8. Repeat the Process if Necessary:

• If the termination criterion is not satisfied, return to step 4 and continue the process.

6 Numerical simulations

6.1 Case description

The proposed methodology is applied to two different modified IEEE test systems to demonstrate its feasibility and robustness. The reliability parameters are defined consistently for both cases and are explained in this section. The line with the highest impedance is estimated to have a failure rate of 0.4 failures per year, while the line with the lowest impedance has a failure rate of 0.1 failures per year. The failure rates for the other lines are calculated relative to these values, based on the length of each line [39]. Additionally, the required repair time for the lines is set to 6 h. A main breaker is installed in the main feeder, and a sectionalizer is planned to be placed at the beginning of each branch or line. The emission coefficients for the units are provided in Table 2.Table 2 Emission coefficients.

Table 2Coefficient	αi	βi	γi	ζi	ηj	θj	
Value	0.02543	−0.06047	0.05638	0.0005	0.0015	0.00015	
Coefficient	πk	ρk	τi	kj	σk	-	
Value	0.008	0.001	0.28	0.2	0.4	–	

Regarding the optimization techniques, the population size for SALBA is set to 30, and the maximum number of iterations is set to 50. Additionally, both αBat and γBat are set to 0.9. The variable (rg) is assigned a random value, and Ag is set to 1. Furthermore, to confirm the validity of the results, they are benchmarked against three widely-used evolutionary algorithms: Honey Bee Mating Optimization (HBMO), Particle Swarm Optimization (PSO), and Genetic Algorithm (GA) [[39], [40], [41],43]. The settings for these algorithms are as follows.• PSO: The cognitive and social parameters are both set to 2. The initial and final inertia weight factors are 0.9 and 0.4, respectively. The population size (Npop) is 30, and the maximum number of iterations (Itmax) is 50.

• GA: The maximum number of iterations (Itmax) is 50, the population size (Npop) is 30, the crossover rate is 0.3, and the mutation probability is 0.09.

• HBMO: The parameters Smax, αHBMO, NWorker, NDreone, NSperm, and NBrood are set to 1, 0.93, 10, 16, 15, and 15, respectively.

These parameters are consistently applied to all test systems. Additionally, the error criterion for the iterative method is set at 0.001. The simulations were carried out in a Matlab 7.4 environment on a Pentium P4, Core 2 Duo 2.2 GHz PC with 2 GB of RAM.

6.2 Case I

This 13.8 kV distribution network, as shown in Fig. 3, includes 16 load points, two distribution substations, 24 branches, and five looping branches [35]. Each branch is isolated from the network by an isolator switch. The network has five branches that can be upgraded, eight new load points, and 16 paths for installing additional load points. The original system does not include any thermal loads; therefore, arbitrary thermal loads are considered in this paper. The capacities of these thermal loads are 1 MW, 800 kW, and 650 kW, located at buses 2, 11, and 14, respectively.Fig. 3 Schematic representation of the network with 18 nodes.

Fig. 3

Other parameters and specifications of the system are provided in Table 3. Table 4 shows the new connection points and load points. It is assumed that the cost of electrical and thermal units is not affected by their technology or type, so a standard cost is assigned to each unit. Units with a capacity of 1 MVA, which only supply power, can be installed at any node that handles loads. In contrast, the 1 MW CHP units and the 150 kW heating units can only be connected to nodes with a thermal demand, and their connection order is determined by a predefined priority. The optimal values for these parameters are obtained by running the HBMO, PSO, GA, and SALBA algorithms 50 times each.Table 3 Technical/Cost parameters.

Table 3

Table 4 New points of connection and load and their code.

Table 4Load Point	Connection Points	Load Point	Connection Points	
7	S/S-2, 6, 8	11	S/S-2,10, 15	
8	4, 8, 12	13	S/S-1, 9, 14	
9	S/S-1, 10, 13	14	13, 15	
10	5, 9, 11	15	11, 14, 16	

To demonstrate the effectiveness of the SALBA algorithm, a comparative display of the convergence plots for the objective function is presented. Fig. 4 shows that the algorithms require different numbers of iterations to achieve convergence.Fig. 4 Evaluating the convergence trends for the objective function through ITLO, HBMO, PSO, and GA algorithms for 18-bus test system.

Fig. 4

The figure confirms the advantage of the proposed algorithm over the PSO, GA, and HBMO methods. Since the final results of these algorithms are equivalent, only the outcomes from the proposed method are detailed. The suggested EEP issue has been resolved, and the corresponding results are presented in Table 5, Table 6.Table 5 Findings for every individual objective function.

Table 5Objective Function	Values	
Cost	28271.52	
Emission	7362.31	

Table 6 Outcomes associated with the newly incorporated generation units.

Table 6Bus	No. of the units	
Power units	Heat units	CHP units	
2	0	0	2	
10	2	0	0	
11	0	2	1	
14	1	1	1	

Table 5 presents the values of the objective functions, while Table 6 shows the locations, types, and numbers of newly installed generating units. Fig. 5 illustrates the single-line schematic of the distribution system, representing the optimal compromise solution.Fig. 5 The best compromised solution for the reconfiguration network.

Fig. 5

In this figure, the branch labeled S/S-1-5 indicates that the branch between substation 1 and bus 5 needs to be rewired. Other branches that need rewiring are S/S-2-12 and 12–16. Additionally, the newly installed lines are S/S-1-9, 4–8, 5–10, 7–8, 8–12, 9–13, 11–15, 13–14, and 15–16. Switches are labeled according to the buses they connect to; for example, 'switch3-4′ refers to the switch located on the segment connecting buses 3 and 4, which needs to be opened to maintain the radial configuration of the network.

6.3 Case II

Fig. 6 [36] illustrates the single-line schematic of the secondary test setup. This setup consists of an 11 kV radial distribution network that includes two substations and four feeders, with a total of 70 nodes and 78 branches, including tie branches.Fig. 6 Schematic representation of the network with 70 nodes.

Fig. 6

Similar to the previous case, arbitrary thermal loads are considered for this system as specified in Table 7. Additionally, there are 8 existing branch sections that can be upgraded, 3 new load points, and 9 paths for the installation of new branch sections, as detailed in Table 8.Table 7 Information the thermal loads and the system.

Table 7Bus No.	9	27	44	57	62	
Thermal load (kW)	200	250	150	300	250	

Table 8 Information of new candidate branches.

Table 8No. of new branch	From	To	R(Ω)	X(Ω)	
1	71	9	0.2	0.3	
2	71	15	0.1	0.15	
3	71	30	0.15	0.25	
4	72	45	0.1	0.15	
5	72	59	0.1	0.15	
6	72	60	0.15	0.25	
7	73	21	0.2	0.3	
8	73	27	0.15	0.25	
9	73	28	0.1	0.15	

The selection of candidate locations for the installation of various units is based on the same criteria as the previous case. However, the maximum size of power units is 200 kW, the size of CHP units is 300 kW, and the capacity of heat units is 50 kW. Fig. 7 displays the convergence graphs for the objective function using the HBMO, PSO, GA, and SALBA algorithms. It shows that the proposed algorithm reaches the global optimum solution in just 22 iterations, while the other algorithms require over 30 iterations to converge to global solutions. In other words, SALBA achieves convergence approximately 1.4 times faster than the others. The figure clearly indicates that the SALBA algorithm is more effective for the EEP task than its counterparts, achieving convergence to the global optimum in fewer iterations, which is crucial in the computational operations of the power systems field. The capability of DFR is examined in this case by performing the proposed EEP. Analyzing the results given in Table 9, Table 10, Table 11 reveals the efficacy and superiority of SALBA. It is noted that the SALBA parameters were obtained by running the algorithm 100 times.Fig. 7 Evaluating the convergence trends for the objective function through ITLO, HBMO, PSO, and GA algorithms for 70-bus test system.

Fig. 7

Table 9 Results of network reconfiguration and objective function.

Table 9Objective Functions Values	Open sectionalizing switches	Closed tie switches	
Cost	Emission	
13442.97	2628.12	S48–49, S14-15,S27-28,S58–59, S42–48, S65-66	tie1, tie4, tie5, tie9, tie10 and tie11	

Table 10 Result for new added units.

Table 10Bus	No. of units	
Power Units	Heat Units	CHP units	
9	1	0	1	
27	0	1	1	
40	2	0	0	
44	0	0	1	
54	1	0	0	
57	1	2	1	
62	0	1	1	

Table 11 Result for new upgraded or added lines.

Table 11New added lines	Upgraded lines	
(9, 71), (45, 72), (27, 73)	(22,66), (66,67)	

SALBA has been successful in promptly reaching the best solutions for the objective functions. Fig. 8 illustrates the single-line representation of the test system after optimization. The objective function value associated with the positions of sectionalizing and tie switches is reported in Table 9.Fig. 8 The best compromised solution of the distribution network.

Fig. 8

According to this table, the following switches are opened: switch 14–15, switch 27–28, switch 42–48, switch 48–49, switch 58–59, and switch 65–66. The following switches are closed: tie1, tie4, tie5, tie9, tie10, and tie11. Altering the open/close status of these switches ensures the radial topology of the system. It is also noted that achieving an optimal value for the objective functions required the installation of new generating units, as detailed in Table 10.

Additionally, new lines 9–71, 45–72, and 27–73 are installed. To meet the problem requirements, lines 22–66 and 66–67 are rewired, as detailed in Table 11.

7 Conclusion

In this study, a comprehensive static Energy Expansion Planning (EEP) model was proposed to optimally plan the future expansion of multi-type energy generating systems, including both thermal and electrical systems, along with new planning alternatives. The considered generating units include power-only units, heat-only units, and Combined Heat and Power (CHP) units. The installation and advancement of network elements were fulfilled through the proposed EEP to guarantee a secure, cost-effective, and reliable configuration of the network capable of satisfying the predicted load demand.

The implementation of the proposed EEP model resulted in significant improvements across several key metrics. Firstly, there was a substantial reduction in overall network costs. The optimized expansion planning minimized the need for costly infrastructure investments while ensuring efficient energy distribution. Additionally, the integration of RERs and the optimized placement of new generating units led to a notable decrease in greenhouse gas emissions, successfully reducing SO2, CO2, and NOx emissions compared to conventional power plants. Furthermore, the reconfiguration of the network topology through DFR significantly reduced energy losses. The optimized network configuration ensured minimal transmission and distribution losses, enhancing overall system efficiency. The proposed methodology also improved network reliability by optimizing the placement and operation of generating units and reconfiguring the network. The reliability indices showed marked improvement, indicating a more robust and dependable energy system. Moreover, the optimized expansion planning and network reconfiguration resulted in a more stable voltage profile across the network, significantly enhancing the voltage stability index and ensuring consistent and reliable power delivery. The problem was solved using the SALBA algorithm, which takes advantage of adaptive modifications to diversify the searching population and achieve faster convergence to the global optimum solution. Application of the proposed modeling and solving methodology to two different test systems corroborated the veracity and robustness of the procedure. The methodology demonstrated its applicability to realistic large-scale systems and successfully minimized cost, minimized loss, enhanced reliability, and reduced emissions securely.

In conclusion, the results indicate that the proposed concurrent expansion planning model, combined with the advanced optimization capabilities of SALBA, provides a comprehensive and effective solution for the expansion of hybrid electrical and thermal systems. For future work, it is recommended that other aspects of EEP, such as dynamic modeling and stochastic planning of the system along with other resources such as batteries or demand respond, be analyzed.

CRediT authorship contribution statement

Ali Reza Abbasi: Writing – review & editing, Writing – original draft, Visualization, Validation, Supervision, Software, Resources, Project administration, Methodology, Investigation, Funding acquisition, Formal analysis, Data curation, Conceptualization. Mahmoud Zadehbagheri: Writing – original draft, Visualization, Validation, Supervision, Software, Resources, Project administration, Methodology, Investigation, Funding acquisition, Formal analysis, Data curation, Conceptualization.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
==== Refs
References

1 Liu Z. Nazari-Heris M. Robust bidding strategy of interconnected multi-carrier systems in the electricity markets under the uncertainty in electricity load Sustain. Energy Technol. Assessments 57 2023 103245 10.1016/j.seta.2023.103245
2 Mahdavi M. Integrated generation-transmission expansion planning considering power system reliability and optimal maintenance activities Int. J. Electr. Power Energy Syst. 145 2023 108688 10.1016/j.ijepes.2022.108688
3 Sohrabi Tabar V. Risk-based day-ahead planning of a renewable multi-carrier system integrated with multi-level electric vehicle charging station, cryptocurrency mining farm and flexible loads J. Clean. Prod. 380 Part 1 2022 134898 10.1016/j.jclepro.2022.134898
4 Mahdavi M. An accurate evaluation of consumption pattern in reconfiguration of electrical energy distribution systems 2021 IEEE Industry Applications Society Annual Meeting (IAS) 2021 1 7 10.1109/IAS48185.2021.9677155 Vancouver, BC, Canada
5 Mahdavi Robust distribution network reconfiguration in the presence of distributed generation under uncertainty in demand and load variations IEEE Trans. Power Deliv. 38 5 2023 3480 3495 10.1109/TPWRD.2023.3277816
6 Möller B. Heat Roadmap Europe: identifying local heat demand and supply areas with a European thermal atlas Energy 158 2018 281 292 10.1016/j.energy.2018.06.025
7 Brauers H. Liquefied natural gas expansion plans in Germany: the risk of gas lock-in under energy transitions Energy Res. Social Sci. 76 2021 10.1016/j.erss.2021.102059
8 Riepin I. Adaptive robust optimization for European strategic gas infrastructure planning Appl. Energy 324 2022 119686 10.1016/j.apenergy.2022.119686
9 Davoodi A. Multi-objective dynamic generation and transmission expansion planning considering capacitor bank allocation and demand response program constrained to flexible-securable clean energy Sustain. Energy Technol. Assessments 47 2021 101469 10.1016/j.seta.2021.101469
10 Pawar A. A novel squirrel-cat optimization based optimal expansion planning for distribution system Sustainable Computing: Informatics and Systems 43 2024 10.1016/j.suscom.2024.101017
11 Goodarzi S. Tight convex relaxation for TEP problem: a multiparametric disaggregation approach‏ IET Generation, Transmission & Distribution 14 14 2022 2810 2817‏ 10.1049/iet-gtd.2019.1270
12 Pawar A. A novel squirrel-cat optimization based optimal expansion planning for distribution system Sustainable Computing: Informatics and Systems 43 2024 10.1016/j.suscom.2024.101017
13 Rezasi A. Comparison parametric and non-parametric methods in probabilistic load flow studies for power distribution networks Electr. Eng. 104 2022 3943 3954 10.1007/s00202-022-01590-9
14 Kumar R. Reliability and economics evaluation for generation expansion planning incorporating variability in wind energy sources Elec. Power Syst. Res. 224 2023 109720 10.1016/j.epsr.2023.109720
15 Mohammadi M. Probabilistic load flow in distribution networks: an updated and comprehensive review with a new classification proposal Elec. Power Syst. Res. 222 2023 109497 10.1016/j.epsr.2023.109497
16 El-Meligy M. Voltage stability-constrained robust transmission expansion planning using binary discretization technique Elec. Power Syst. Res. 231 2024 110303 10.1016/j.epsr.2024.110303
17 Lekan Gbadamosi S. Nwulu N. A comparative analysis of generation and transmission expansion planning models for power loss minimization Sustainable Energy, Grids and Networks 26 2021 100456 10.1016/j.segan.2021.100456
18 Baleanu D. Recent developments of energy management strategies in microgrids: an updated and comprehensive review and classification Energy Convers. Manag. 297 2023 117723 10.1016/j.enconman.2023.117723
19 Juan M. A stochastic mixed-integer convex programming model for long-term distribution system expansion planning considering greenhouse gas emission mitigation Int. J. Electr. Power Energy Syst. 108 2019 86 95 10.1016/j.ijepes.2018.12.042
20 Tejada D. Transmission network expansion planning considering repowering and reconfiguration Int. J. Electr. Power Energy Syst. 69 2015 213 221 10.1016/j.ijepes.2015.01.008
21 Ghasemi H. MILP model for integrated expansion planning of multi-carrier active energy systems IET Generation, Transmission & Distribution ume13 7 2019 1177 1189 10.1049/iet-gtd.2018.6328
22 Bakhshi H. A multi-objective mixed integer linear programming model for integrated electricity-gas network expansion planning considering the impact of photovoltaic generation Energy 222 2021 119933 10.1016/j.energy.2021.119933
23 Zhang X. Reliability-based optimal planning of electricity and natural gas interconnections for multiple energy hubs IEEE Trans. Smart Grid 8 4 2017 10.1109/TSG.2015.2498166 1658–167
24 Chaudry M. Combined gas and electricity network expansion planning Appl. Energy 113 2014 1171 1187 10.1016/j.apenergy.2013.08.071
25 Khaligh V. Anvari-Moghaddam A. Stochastic expansion planning of gas and electricity networks: a decentralized-based approach Energy 186 2019 115889 10.1016/j.energy.2019.115889
26 Klatzer T. State-of-the-art expansion planning of integrated power, natural gas, and hydrogen systems Int. J. Hydrogen Energy 47 47 2022 20585 20603 10.1016/j.ijhydene.2022.04.293
27 Zadehbagheri M. Sutikno T. Kiani M. A new method of virtual direct torque control of doubly fed induction generator for grid connection Int. J. Electr. Comput. Eng. 13 1 2023 1201 1214 10.11591/ijece.v13i1.pp1201-1214
28 Zhou Xianzheng Guo Chuangxin Wang Yifei Li Wanqi Optimal expansion Co-planning of reconfigurable electricity and natural gas distribution systems incorporating energy hubs Energies 10 1 2017 124 10.3390/en10010124
29 Sabzpoosh Saravi V. Resilience-constrained expansion planning of integrated power–gas–heat distribution networks Appl. Energy 323 2022 119315 10.1016/j.apenergy.2022.119315
30 Zadehbagheri M. Pishavaei M. Ildarabadi R. Sutikno T. The coordinated control of FACTS and HVDC using H-infinity robust method to stabilize the inter-regional oscillations in power systems Int. J. Power Electron. Drive Syst. 8 3 2017 1274 1284
31 Goel L. Billinton R. Evaluation of interrupted energy assessment rates in distribution systems IEEE Trans. Power Del. 6 4 1991 1876 1882 http://www.algs4.cs.princeton.edu/41undirected
32 Kavousi-Fard A. A novel adaptive modified harmony search algorithm to solve multi-objective environmental/economic dispatch J. Intell. Fuzzy Syst. 26 6 2014 2817 2823 10.3233/IFS-130949
33 Lewin P.L. Theed J.E. Davies A.E. Larsen S.T. Method for rating power cables buried in surface troughs Proc. Inst. Elect. Eng., Gen. Transm. Distrib. 146 4 1999 360 364 10.1049/ip-gtd:19990360
34 Khoramnia R. An smart stochastic approach to model plug-in hybrid electric vehicles charging effect in the optimal operation of micro-grids J. Intell. Fuzzy Syst. 28 2 2015 835 842 10.3233/IFS-141365
35 Kavousi-Fard Abdollah Reliability enhancement using optimal distribution feeder reconfiguration Neurocomputing 106 2013 1 11 10.1016/j.neucom.2012.08.033
36 Sakawa M. Yano H. An interactive fuzzy satisficing method for multiobjective nonlinear programming problems with fuzzy parameters Fuzzy Set Syst. 30 3 1989 221 238 10.1016/0165-0114(89)90017-1
37 Akbari K. Optimal placement of distributed generation in radial networks considering reliability and cost indices J. Intell. Fuzzy Syst. 30 2 2016 1077 1086 10.3233/IFS-151883
38 Veena Vani B. An efficient battery swapping and charging mechanism for electric vehicles using bat algorithm Comput. Electr. Eng. 118 Part A 2024 109357 10.1016/j.compeleceng.2024.109357
39 Niknam T. “An efficient multi-objective HBMO algorithm for distribution feeder reconfiguration Expert Syst. Appl. 38 3 2011 2878 2887
40 Ansari J. Load frequency control in power systems by a robust backstepping sliding mode controller design Energy Rep. 10 2023 1287 1298 10.1016/j.egyr.2023.08.008
41 Gandhi C.P. A novel hyperbolic fuzzy entropy measure for discrimination and taxonomy of transformer winding faults IEEE Trans. Instrum. Meas. 71 2022 1 8 10.1109/TIM.2022.3212522 Art no. 3526408
42 Movahedpour M. Kiani M.J. Zadehbagheri M. Mohammadi S. Microgrids Operation by Considering Demand Response and Supply Programs in the Presence of IGDT-Based Reverse Risk 10 2022 IEEE Access 48681 48700 10.1109/ACCESS.2022.3172422
43 Zadehbagheri M. Ma'arif A. Ildarabadi R. Ansarifard M. Suwarno I. Design of multivariate PID controller for power networks using GEA and PSO Journal of Robotics and Control (JRC) 4 1 2023 108 117 10.18196/jrc.v4i1.15682
