iconOpen Access

ARTICLE

A Multi-Stage Expansion Planning Method for Rural Distribution Networks with Flexible Interconnection

Yueyang Ji1, Yaohui Peng1, Haoran Ji1,*, Xinran Na1, Yuxuan Chen1, Wei Li2, Shengbin Chen2

1 State Key Laboratory of Smart Power Distribution Equipment and System, Tianjin University, Tianjin, 300072, China
2 Electric Power Research Institute, China Southern Power Grid (CSG), Guangzhou, 510663, China

* Corresponding Author: Haoran Ji. Email: email

(This article belongs to the Special Issue: Digital and Intelligent Planning and Operation Technologies for Flexible Distribution Network)

Energy Engineering 2026, 123(8), 2 https://doi.org/10.32604/ee.2025.074599

Abstract

With the increasing penetration of distributed generations and continuous growth of loads, traditional rural distribution networks face severe challenges in both hosting capacity and reliability. Addressing these issues requires planning approaches that strike a balance between economic efficiency in infrastructure development and resilience in operation. Considering the dynamic growth of distributed generations and rural loads over the planning horizon, this paper presents a multi-stage expansion planning approach that coordinates flexible interconnection devices (FIDs) with substation and line construction to improve both economic performance and system reliability. The proposed method account for the time-varying growth of DGs and loads, as well as the declining investment cost of power electronic devices across multiple planning stages. The model holistically considers both economic efficiency and operational reliability, formulating the problem as a mixed-integer second-order cone programming (MISOCP) model to ensure computational efficiency. Case studies conducted on a practical 138-node rural distribution network in Guangxi, China, demonstrate the effectiveness of the proposed method. Compared to traditional single-stage or single-resource planning strategies, results indicate that the proposed multi-stage coordinated strategy achieves a significant reduction in total annualized cost while simultaneously enhancing system reliability, effectively mitigating voltage violations, and achieving a 100% PV accommodation rate without curtailment. This work provides a practical and adaptive planning framework for rural distribution networks, offering valuable insights for achieving cost-effective and resilient network development under rural energy transition.

Keywords

Rural distribution networks; coordinated planning; flexible interconnection device; multi-stage planning; reliability

1  Introduction

Driven by dual-carbon goals and rural revitalization strategy, rural areas in China are undergoing rapid energy transition and electrification [1]. The extensive integration of distributed renewable energy, particularly photovoltaics [2], is transforming rural distribution networks from traditional passive and unidirectional systems into active networks with bidirectional power flows and complex spatiotemporal characteristics [3]. While this transformation facilitates the development and consumption of renewable energy, it also poses significant challenges to planning and operation of rural distribution networks [4]. In remote regions, historical underinvestment and geographic constraints have led to persistently weak network structures and long power supply radii [5]. These issues result in severe voltage fluctuations, high power losses [6], and poor fault recovery capability [7]. As distributed generator (DG) penetration continues to rise, issues such as reverse power flow and voltage violations become increasingly prominent. These developments underscore the inadequacy of traditional planning approaches in meeting future demands for safe, efficient, and flexible rural distribution networks [8].

To enhance power supply capacity, traditional distribution network planning has primarily focused on substation siting optimization [9] and feeder expansion [10]. Although such solutions can improve network capacity, they suffer from long investment cycles, diminishing marginal returns, and limited adaptability to the dynamic hosting capacity requirements driven by rapid changes in generation and load patterns [11]. In recent years, advancements in flexible power electronics and energy storage technologies have offered new perspectives for distribution planning. A multi-criteria decision framework for optimal deployment of mobile solid-state hydrogen energy storage systems is proposed in [12] to enhance the reliability and sustainability of renewable-powered distribution networks. In particular, flexible interconnection devices (FIDs) can provide power regulation [13], voltage support [14] and load transformation [15] without altering the physical topology [16]. Research has demonstrated the multiple benefits of FIDs in active distribution networks. FIDs have been widely explored to enhance distribution network performance under high photovoltaic penetration. For instance, they have been used to improve hosting capacity [17], reduce operational losses [18], and support multi-objective robust optimization of network resources [19]. FIDs have also been leveraged to enhance fault recovery and self-healing capabilities, enabling rapid load restoration during power outages [20].

At the planning level, several studies have incorporated FIDs into expansion decision-making. A siting and sizing model of FIDs is established in [21], effectively enhancing the economy, flexibility, and reliability of distribution networks. Research in [22] demonstrates that FIDs can simultaneously mitigate harmonic distortion, voltage deviation, and power losses under renewable and load uncertainties through optimal allocation and coordinated active/reactive power control. Moreover, recent efforts have explored the coordinated planning of FIDs with conventional grid assets to exploit their complementary strengths in flexibility and capacity [23]. A mixed-integer linear programming model for collaborative expansion planning of multi-energy distribution networks is developed in [24], which adeptly manages the complexities involved in the construction and reinforcement of multiple resources.

However, current research still exhibits notable limitations. First, most planning models overlook the functional complementarity and investment synergy potential between configurations of FIDs and traditional capacity expansion measures such as new substation construction and line reinforcement. Second, existing models are predominantly based on single-stage static assumptions, failing to adequately account for the dynamic growth of generations and loads over time in rural areas, as well as the time-varying characteristics of equipment investment costs and technical parameters across different planning stages.

To address the aforementioned issues, this paper proposes a multi-stage, multi-resource planning method tailored for rural distribution grids. It aims to achieve coordinated, dynamic, and cost-effective allocation of FIDs and traditional resources such as substations. The main contributions are summarized as follows: (1) A multi-stage expansion planning framework is developed to account for the growth of distributed generations and loads. This framework allows investment decisions to be adjusted dynamically across multiple planning stages, enabling the planning process to better accommodate the evolution of sources and load patterns, as well as the declining tendency of FID costs. (2) A multi-resource collaborative optimization model is established that considers the coordinated allocation of FIDs, substations, and line construction simultaneously. The proposed model incorporates a comprehensive objective function that integrates stage-dependent investment costs, operational losses, and reliability assessments considering FID-assisted fault restoration. Thereby, economic objectives and reliability objectives are unified within the multi-stage framework. (3) Case studies based on a real rural distribution network in Guangxi, China, are conducted to validate the effectiveness of the proposed method. The results demonstrate that the proposed method can significantly reduce the total cost over the planning horizon and substantially improve reliability compared to single-stage planning approaches.

The remainder of this paper is organized as follows. Section 2 details the mathematical model of the multi-stage coordinated planning problem. Section 3 presents the case parameters and discusses the calculation results. Section 4 summarizes the paper and outlines future research directions.

2  Multi-Stage Planning Model for Rural Distribution Networks

To accommodate the dynamic growth of generations and loads in rural distribution networks during the planning horizon, A multi-stage coordinated expansion planning model that integrates flexible interconnection devices with substation construction is established.

2.1 Planning Objectives Considering Both Economic Efficiency and Reliability

Considering the construction and operational costs of FIDs and substations, as well as the operational and outage losses of rural distribution networks, the comprehensive objective function can be formulated as follows:

minf=fFID+fSL+fE+fR(1)

2.1.1 Annual Construction and Operation Cost of FID

The FID annual construction and operating cost fFID is regarded as the sum of the construction cost fC and the operating cost fO. The multi-stage construction cost of FIDs is shown in (2b)(2d), and the multi-stage operational cost of FIDs is formulated as (2e) and (2f).

fFID=fC+fO(2a)

fC=r=1NrfrC(2b)

frC=d(1+d)y(1+d)y1ijΩFIDcr,ijFIDSr,ijFID,r=1(2c)

frC=d(1+d)y(1+d)y1ijΩFIDcr,ijFID(Sr,ijFIDSr1,ijFID),r2(2d)

fO=r=1NrfrO/Nr(2e)

frO=ηijΩFIDcr,ijFIDSr,ijFID(2f)

2.1.2 Annual Construction and Operation Costs of Substations and Lines

The annual investment, construction, and operational maintenance costs for new substations and transmission lines fSL is shown in (3a). It comprises the sum of the costs for new substations frSub, new constructed transmission lines frL1, and reconstructed transmission lines frL2.

fSL=r=1Nr(frSub+frL1+frL2)(3a)

frSub={(ξ+γ)iΩSub(crSub,0δr,iSub+crSubSr,iSub),r=1(ξ+γ)iΩSub(crSub,0(δr,iSubδr1,iSub)+crSub(Sr,iSubSr1,iSub)),r2(3b)

frL1={(ξ+γ)crL1ijΩL1δr,ijL1LijL1,r=1(ξ+γ)crL1ijΩL1(δr,ijL1δr1,ijL1)LijL1,r2(3c)

frL2={(ξ+γ)crL2ijΩL2δr,ijL2,r=1(ξ+γ)crL2ijΩL2(δr,ijL2δr1,ijL2),r2(3d)

2.1.3 Annual Operational Loss Cost of Rural Distribution Network

The average annual operating loss cost fE is shown in Eq. (4a). ErLOSS represents the expected daily operating losses of rural distribution network, encompassing both line losses and the operating losses of the FID installed on the distribution network.

fE=365cPr=1NrErLOSS/Nr(4a)

ErLOSS=h=1NHt=1NTρh(ijΩbrijIr,h,t,ij2+ijΩFIDPr,h,t,ijFID,L)Δt(4b)

2.1.4 Annual Power Outage Loss Cost of Rural Distribution Network

The annual power outage cost fR is shown in Eq. (5a). The power outages ErPSI is calculated based on the power outage durations, as shown in Eqs. (5b) and (5c). The cost of power outage losses ErLOSS,F is shown in Eq. (5d).

fR=r=1Nr(cRErPSI+cPErLOSS,F)/Nr(5a)

ErPSI=kΩfλkiΩkLTr,k,iPr,k,iLOAD(5b)

Tr,k,i=(1μr,k,i)trep+μk,i0tiso+(μr,k,iμk,i0)(tiso+tFID)(5c)

ErLOSS,F=kΩfλk(treptisotFID)×(ijΩbrijIr,k,ij2+ijΩFIDPr,k,ijFID,L)(5d)

2.2 Multi-Stage Configuration Constraints of FIDs and Substations

The multi-stage planning for rural distribution networks requires coordinated configuration of FIDs and substations, and construction of transmission lines.

2.2.1 Multi-Stage Configuration Constraints of Flexible Interconnection Devices

During the planning process, configuration constraints of FIDs must be considered. The minimum configuration capacity for FIDs should be established as shown in Eqs. (6a) and (6b). Eqs. (6c) and (6d) represent the installation quantity limits for FIDs. Additionally, the capacity of FID at the current stage is always greater than or equal to that of the previous stage, as illustrated in Eqs. (6e) and (6f).

mijδr,ijFIDMFID(6a)

Sr,ijFID=smmij(6b)

ijΩFIDδr,ijFIDNFID,g(6c)

ijΩFIDδr,ijFIDNFID,n,iΩn(6d)

δr1,ijFIDδr,ijFID,r2(6e)

Sr1,ijFIDSr,ijFID,r2(6f)

2.2.2 Multi-Stage Configuration Constraints of Substations and Lines

The capacity constraints of substations must be considered during the planning process, as shown in Eq. (7a). The configuration of new substations will result in constructions of new lines and new interconnections, as shown in Eqs. (7b) and (7c). Eqs. (7d) and (7e) represent multi-stage substation capacity constraints, indicating that the capacity will not be reduced in subsequent planning stages. Eqs. (7f) and (7g) represent multi-stage line construction constraints, allowing only line additions during the planning process.

0Sr,iSubM×δr,iSub(7a)

ijΩL1δr,ijL1=δr,iSub(7b)

ijΩL2δr,ijL2=δr,iSub(7c)

δr1,iSubδr,iSub,iΩSub,r2(7d)

Sr1,iSubSr,iSub,iΩSub,r2(7e)

δr1,ijL1δr,ijL1,ijΩL1,r2(7f)

δr1,ijL2δr,ijL2,ijΩL2,r2(7g)

2.3 Operational Constraints of Rural Distribution Networks

2.3.1 Operational Constraints under Normal Conditions

(1) Power Flow Constraints

The branch flow method is employed to model the fundamental power flow constraints of rural distribution networks. Eqs. (8a) and (8b) represent the active and reactive power balance constraints at nodes, respectively. Eqs. (8c)(8e) denote the branch flow constraints for general nodes, nodes of lines to be constructed, and nodes of lines to be reconstructed, respectively. The active and reactive power constraints for substation nodes are shown in Eqs. (8g)(8j).

jiΩb(Pr,h,t,jirjiIr,h,t,ji2)+Pr,h,t,i=ilΩbPr,h,t,il(8a)

jiΩb(Qr,h,t,jixjiIr,h,t,ji2)+Qr,h,t,i=ilΩbQr,h,t,il(8b)

Ur,h,t,i2Ur,h,t,j2+(rij2+xij2)Ir,h,t,ij22(rijPr,h,t,ij+xijQr,h,t,ij)=0,ijΩb(ΩL1ΩL2)(8c)

|Ur,h,t,i2Ur,h,t,j2+(rij2+xij2)Ir,h,t,ij22(rijPr,h,t,ij+xijQr,h,t,ij)|M×(1δr,ijL1),ijΩL1(8d)

|Ur,h,t,i2Ur,h,t,j2+(rij2+xij2)Ir,h,t,ij22(rijPr,h,t,ij+xijQr,h,t,ij)|M×δr,ijL2,ijΩL2(8e)

Pr,h,t,ij2+Qr,h,t,ij2=Ir,h,t,ij2Ur,h,t,i2(8f)

|Pr,h,t,iDG+Pr,h,t,iFIDPr,h,t,iPr,h,t,iLOAD|M×δr,iSub,iΩSub(8g)

Pr,h,t,i=Pr,h,t,iDG+Pr,h,t,iFIDPr,h,t,iLOAD,iΩnΩSub(8h)

|Qr,h,t,iDG+Qr,h,t,iFIDQr,h,t,iQr,h,t,iLOAD|M×δr,iSub,iΩSub(8i)

Qr,h,t,i=Qr,h,t,iDG+Qr,h,t,iFIDQr,h,t,iLOAD,iΩnΩSub(8j)

(2) Secure Operation Constraints

To ensure the safe operation of rural distribution networks, node voltage constraints and branch current constraints are specified as shown in Eqs. (9a) and (9b), respectively. Additionally, considering substation reverse power flow rate limitations, the reverse power flow capacity of 10 kV rural distribution networks must not exceed 80% of the capacity of the transformer of the superior 110 kV substations, as indicated in Eq. (9c).

(Umin)2Ur,h,t,i2(Umax)2(9a)

0Ir,h,t,ij2(Imax)2(9b)

0.8StranPr,h,t,ij,tranΩSup,ijΩtr(9c)

(3) Operational Constraints of FIDs

FID operation must satisfy the active power balance constraint, as shown in Eq. (10a), where equipment losses are given by Eq. (10b). Furthermore, the apparent power during FID operation must not exceed the equipment capacity, as indicated in Eq. (10c).

Pr,h,t,iFID+Pr,h,t,jFID+Pr,h,t,iFID,L+Pr,h,t,jFID,L=0(10a)

AiFID(Pr,h,t,iFID)2+(Qr,h,t,iFID)2=Pr,h,t,iFID,L(10b)

(Pr,h,t,iFID)2+(Qr,h,t,iFID)2Sr,ijFID(10c)

2.3.2 Operation Constraints under Fault Condition

(1) Radial Constraint of Rural Distribution Networks

Except for branches interconnected by FIDs, power restoration of rural distribution networks after faults should adhere to the principle of radial operation as follows.

αr,k,ij=βr,k,ij+βr,k,ji,ijΩbΩFID(11a)

αr,k,ij=δr,ijL1,ijΩL1Ωk(11b)

αr,k.ij=1δr,ijL2,ijΩL2(11c)

ijΩbΩFIDβr,k,ij=1θr,k,j,jΩnΩ0(11d)

ijΩbΩFIDβr,k,ij=0,jΩ0(11e)

(2) Power Restoration Constraints based on FIDs

Considering that the FID can switch to the Uacθ control mode during faults to provide voltage support for power outage areas, as shown in Eqs. (12a)(12c). Furthermore, it is assumed that nodes in the power outage areas can only be restored when supported by a stable voltage source, as depicted in Eqs. (12d)(12h).

(UVf,min)2MVf(1βr,k,ij)Ur,k,j2(UVf,max)2+MVf(1βr,k,ij),ijΩFID(12a)

αr,k,ijδr,ijFID,ijΩFID(12b)

βr,k,ij=0,ijΩFID,iΩkNS(12c)

sr,k,ijiΩbΩFID(βr,k,jisr,k,j),iΩn(Ω0ΩSub)(12d)

sr,k,i1θr,k,i(12e)

μr,k,iMSupsr,k,i,iΩkNS(12f)

sr,k,i=1,iΩ0(12g)

sr,k,i=δr,iSub,iΩSub(12h)

Additionally, in fault scenarios, the subscripts of variables subject to power flow constraints should be adjusted accordingly. For instance, Pr,k,i represents the active power at node i during the kth fault scenario in the rth stage, and Pr,k,iLOAD represents the load restoration amount at node i during the kth fault scenario in the rth stage. Due to page limitations, further elaboration is omitted here.

3  Case Study

The case study is based on a real 10 kV rural distribution network in Guangxi, China. The model is implemented in MATLAB R2022b using YALMIP and solved with the commercial solver CPLEX. The optimality tolerance is set to 1% to balance computational efficiency and solution accuracy. All simulations are run on a computer with an Intel Core Ultra 7265K processor (3.90 GHz), 32 GB RAM, and Windows 10. To capture the temporal characteristics of distributed generation and load, two representative days are selected, one from winter and the other from summer. Each typical day is divided into six time periods, with corresponding photovoltaic output and load power curves used as input data.

3.1 The 138-Node Rural Distribution Network

A practical rural distribution network in Guangxi, China, is used to validate the proposed model. The system comprises four 10 kV feeders. Feeders A and B exhibit high photovoltaic penetration rates. The distribution lines of Feeder C and D are long and heavily-loaded. The topology of the system is shown in Fig. 1. The case study comprises a total of 138 nodes and 134 branches, with photovoltaic (PV) connections at 24 nodes. The load and DG conditions for each area in the case study are summarized in Table 1.

images

Figure 1: Topology of the 138-node rural distribution network

images

The maximum installation capacity for FID is 5 MV·A, with a unit configuration capacity of 0.1 MV·A. The loss coefficient is 0.02, and the initial unit investment cost for FID is 800 CNY/(kV·A). The discount rate of FID is 0.08, with a service life of 20 years and an operation cost factor of 0.01. To ensure safe and stable distribution network operation, the current amplitude upper limit for 10 kV lines is set at 1000 A. The safe operating voltage range of rural distribution network is set to 0.95 to 1.05 p.u. Load and DG output curves are shown in Fig. 2. Two typical daily operating curves are considered as representative scenarios for planning, each comprising 6 time slots.

images

Figure 2: Typical scenarios for photovoltaic generation and load level

3.2 Scheme Settings

The proposed method is applied to plan the rural distribution network over a 15-year horizon. The line fault between node 27 and node 28 of feeder A is used as a fault scenario for reliability assessment. Table 2 lists the investment costs of different equipment or construction items in each planning stage. Table 3 shows the projected growth of load and distributed generation across the stages.

images

images

Based on the above parameters, the following four schemes are defined.

•   Scheme 1: The original operational state of the system without planning.

•   Scheme 2: Considering the economic efficiency, substations and FIDs are configured with a planning horizon of 15 years in order to evaluating the economic benefit of infrastructure investment.

•   Scheme 3: Considering both the economic efficiency and reliability, substations and FIDs are configured with a planning horizon of 15 years, evaluating the economic and reliability improvement.

•   Scheme 4: Considering both the economic efficiency and reliability, substations and FIDs are configured with a planning horizon of 15 years in three stages for the accommodation of the source and load variation.

3.3 Results Analysis

The detailed configurations of FIDs and substations in Schemes 2 and 3 are presented in Table 4.

images

The optimal configurations of FIDs and substations in each stage of Scheme 4 are shown in Fig. 3 and the detailed configurations are presented in Table 5.

images

Figure 3: The optimal configurations of FIDs and substations

images

As shown in Table 5, the locations that substations and FIDs are installed remain unchanged across the three stages, while the capacity of the installed FIDs is progressively expanded over the planning stages. This suggests that incrementally increasing FID capacity in a staged manner while meeting stage-specific supply capacity and reliability constraints can achieve better economic performance without compromising planning objectives, as the unit investment cost of FIDs decreases.

Table 6 lists the planning performance metrics under the four schemes.

images

As shown in Table 6, compared with Scheme 1, Scheme 2 significantly improves system performance by planning FIDs and substations. The node voltages are maintained within secure limits, and both operational losses and total cost are reduced, indicating that the coordinated deployment of FIDs and substations effectively enhances the network’s operational efficiency and economic performance. The PV accommodation rate also increases from 84.38% to 99.7% in the final year of the planning period.

However, Scheme 2 optimizes only for economic cost and ignores reliability under fault conditions, resulting in relatively high outage losses. In contrast, Scheme 3 jointly considers economic and reliability objectives, reducing the average annual outage loss to 1.3 thousand CNY and lowering the corresponding SAIDI index. Meanwhile, the PV accommodation rate reaches 100%.

Scheme 4 introduces a multi-stage planning framework that divides the 15-year horizon into three stages (each stage lasts for five years) and dynamically adjusts investment decisions by evolving source-load profiles and equipment costs. Results show that Scheme 4 also reduces outage losses by 99.2% compared with Scheme 2. The investment cost is reduced by 9.3% compared with Scheme 3. Scheme 4 achieves a total annualized cost of 2.094 million CNY. This represents a significant reduction of 29.7% compared to Scheme 1, and also outperforms the other planning scenarios with cost reductions of 6.2% against Scheme 2 and 2.2% against Scheme 3. Furthermore, voltage violations are fully eliminated, indicating improved voltage stability and operational security. This indicates that Scheme 4 outperforms all other schemes across key metrics.

The multi-stage planning results under different FID cost reduction rate are shown in Table 7.

images

As shown in Table 7, when the cost of FID devices in case 1 decreases rapidly, the installed FID capacity in the first stage is relatively low. In contrast, the capacity constructed in the final stage is higher than other cases. This indicates that the investment in FIDs is postponed due to the rapid decline in equipment cost. In case 3, where the FID unit cost remains constant, the FID installed capacity of the first stage is significantly higher than other cases. The results demonstrate that variations in FID cost trajectories can significantly influence the timing and scale of FID capacity expansion in the planning process. As the FID cost declines more rapidly, the proposed model tends to postpone FID installation to later stages with lower unit costs, thereby reducing the capacity planning in earlier, high-cost stages.

Table 8 presents the multi-stage planning results under different load growth rates. As shown in the table, a higher load growth rate leads to a larger total planned FID capacity, and the incremental FID capacity between consecutive stages also increases correspondingly. This indicates that as the load expands more rapidly, the model allocates additional FID installations to maintain network reliability and voltage stability across the planning horizon.

images

4  Conclusions

This paper presents a multi-stage expansion planning method for rural distribution networks with flexible interconnection. The proposed approach co-optimises economic and reliability costs accounting for the dynamic growth of load and distributed generation. By deploying FIDs and substations at critical locations through multiple planning stages, the method significantly improves system performance. Calculation results demonstrate effective reduction in both operational losses and outage costs. It effectively mitigates voltage violations and achieves a 100% PV accommodation rate without curtailment. The proposed framework provides a quantitative cost-benefit tool that enables planners to identify the most economical investment pathways. For grid operators, the proposed framework serves as a quantitative decision-support tool to determine the optimal timing and scale of FID, substation, and feeder investments under evolving demand and renewable penetration. For policymakers, the model enables cost–benefit assessment and investment prioritisation, supporting the formulation of regulatory and incentive mechanisms that promote flexible and resilient rural distribution networks. Beyond the specific case of Guangxi, the proposed approach can be readily adapted to different regional cases by adjusting local parameters such as investment costs, renewable generation characteristics, and reliability targets.

The proposed multi-stage coordinated planning framework enables adaptive investment pathways that align capacity expansion with evolving sources and loads and declining FID costs. However, the current model is based on deterministic load and PV generation curves. Future research will extend the proposed framework by integrating stochastic or scenario-based optimisation to capture uncertainties in renewable generation and demand evolution. In addition, multi-energy coordination such as energy storage system will be considered to further enhance system flexibility and economic performance. These extensions will strengthen the practical applicability of the proposed approach, and further enhance the hosting capacity and operational performance of rural distribution networks.

Nomenclature

ΩL1, ΩL2, Ωk Set of lines to be constructed, lines to be reconstructed and faulty lines in fault scenario k.
Ωb, ΩFID, Ωtr Set of branches, FID branches, branches connected to 110 kV substations.
Ωn, Ω0, ΩkL Set of nodes, source nodes and power-off nodes in fault scenario k.
Ωf, ΩkNS Set of fault scenarios and non-source nodes in fault scenario k.
ΩSub, ΩSup Set of newly constructed substation nodes and superior 110 kV substation nodes.
Pr,h,t,iFID, Qr,h,t,iFID Active and reactive power injected by FID at node i at time t in scenario h of stage r.
Pr,h,t,ij, Qr,h,t,ij Active and reactive power flow from node i to node j at time t in scenario h of stage r.
Pr,h,t,i, Qr,h,t,i Active and reactive power injected at node i at time t in scenario h of stage r.
Pr,h,t,iDG, Qr,h,t,iDG Active and reactive power injected by DG at node i at time t in scenario h of stage r.
Pr,h,t,iLOAD, Qr,h,t,iLOAD Active and reactive load at node i at time t in scenario h of stage r.
Pr,h,t,ijFID,L Power loss of FID installed at branch ij at time t in scenario h of stage r.
Ur,h,t,i2, Ir,h,t,ij2 Square of the voltage amplitude of node i and the current amplitude of branch ij at time t in scenario h of stage r.
δr,iSub, δr,ijFID Binary variables indicating whether a new substation is to be constructed at node i and whether an FID is to be constructed between node i and node j in stage r.
δr,ijL1, δr,ijL2 Binary variables indicating whether a new line ij is to be constructed and whether line ij is to be reconstructed in stage r.
αr,k,ij, βr,k,ij, θr,k,j Switch status of branch ij, relationship between node i and node j, and node state of node j in fault scenario k of stage r.
sr,k,i Binary variables indicating whether node i has power supply in fault scenario k of stage r.
Sr,iSub, Sr,ijFID Capacity of new substation i and FID installed at branch ij in stage r.
Tr,k,i Power outage duration of node i in failure scenario k of stage r.
μr,k,i, μk,i0 Restoration coefficient of node i with FID and without FID in fault scenario k of stage r.
λk Annual average occurrence frequency of failure scenario k.
sm, mij Unit configuration capacity of FID and number of FID modules installed on branch ij.
ρh Probability of scenario h.
Stran Capacity of the transformer of superior 110 kV substation.
crSub,0, crSub, cr,ijFID Fixed cost and unit capacity cost for new substation, and unit capacity cost of FID located at branch ij in stage r.
crL1, crL2 Unit cost for line construction and fixed cost for line reconstruction in stage r.
cP, cR Unit electricity price and unit cost of outage losses.
d, y, η Discount rate, lifetime, and coefficient of the annual operational costs of FID.
Umax, Umin Upper and lower limit of node voltage amplitude.
UVf,max, UVf,min Upper and lower limit of voltage amplitude on the fault side of FID.
Imax Upper limit of branch current.
NFID,g, NFID,n Maximum number of FIDs installed in rural distribution network and at single node.
Nr, NH, NT Number of planning stages, operational scenarios and time slots within a single scenario.
LijL1 Length of newly constructed line ij.
AiFID Loss coefficient of FID at node i.
MFID Installation limit of FID modules.
M, MSup, MVf Sufficiently large constants in planning constraints and power recovery constraints.
rij, xij Resistance and reactance of branch ij.
ξ, γ Capital recovery factor and equipment maintenance factor.
trep, tiso, tFID Required restoration time by repairing the faulty branch, by isolating the faulty branch using isolating switch and by FID.

Acknowledgement: Not applicable.

Funding Statement: This study was funded by the Smart Gird-National Science and Technology Major Project of China (2024ZD0800600).

Author Contributions: Conceptualization, Yueyang Ji and Haoran Ji; methodology, Yuxuan Chen; validation, Yaohui Peng and Xinran Na; Yaohui Peng and Xinran Na; writing—original draft preparation, Yueyang Ji, Xinran Na, Yuxuan Chen and Yaohui Peng; writing—review and editing, Wei Li and Shengbin Chen; supervision, Haoran Ji. All authors reviewed the results and approved the final version of the manuscript.

Availability of Data and Materials: Not applicable.

Ethics Approval: Not applicable.

Conflicts of Interest: The authors declare no conflicts of interest to report regarding the present study.

References

1. Zhao H, Yao M, Bai Z, Zhang Y, Ding Z, Hu Z. Study on comprehensive efficiency evaluation of rural power grid under rural revitalization strategy considering regional differences. Energy Eng. 2023;120(10):2211–31. doi:10.32604/ee.2023.029371. [Google Scholar] [CrossRef]

2. Li C, Gao K, Chen C, Li J, Wang X, Shao Y, et al. A distributed photovoltaics ordering grid-connected method for analyzing voltage impact in radial distribution networks. Energy Eng. 2024;121(10):2937–59. doi:10.32604/ee.2024.052167. [Google Scholar] [CrossRef]

3. Wang C, Song G, Li P, Ji H, Zhao J, Wu J. Optimal siting and sizing of soft open points in active electrical distribution networks. Appl Energy. 2017;189(3):301–9. doi:10.1016/j.apenergy.2016.12.075. [Google Scholar] [CrossRef]

4. Gupta R, Pena-Bello A, Streicher KN, Roduner C, Farhat Y, Thöni D, et al. Spatial analysis of distribution grid capacity and costs to enable massive deployment of PV, electric mobility and electric heating. Appl Energy. 2021;287:116504. doi:10.1016/j.apenergy.2021.116504. [Google Scholar] [CrossRef]

5. Chawda GS, Shaik AG. Power quality improvement in rural grid using adaptive control algorithm to enhance wind energy penetration levels. IEEE Trans Smart Grid. 2023;14(3):2075–84. doi:10.1109/tsg.2022.3215841. [Google Scholar] [CrossRef]

6. Rasool MI, Saeed MR, Mohammed MA, Suhry OA, Ihsan H. Voltage profile enhancement and power loss reduction with economic feasibility using small capacity distribution transformers. Energy Eng. 2022;119(6):2447–67. doi:10.32604/ee.2022.021871. [Google Scholar] [CrossRef]

7. Fu X, Zhou Y. Collaborative optimization of PV greenhouses and clean energy systems in rural areas. IEEE Trans Sustain Energy. 2023;14(1):642–56. doi:10.1109/tste.2022.3223684. [Google Scholar] [CrossRef]

8. Zhang J, Ji Y, Zhou J, Jia Y, Shi G, Wang H. Cooperative AC/DC voltage margin control for mitigating voltage violation of rural distribution networks with interconnected DC link. IEEE Trans Power Deliv. 2025;40(2):1014–29. doi:10.1109/tpwrd.2025.3535712. [Google Scholar] [CrossRef]

9. Sun DI, Farris DR, Cote PJ, Shoults RR, Chen MS. Optimal distribution substation and primary feeder planning via the fixed charge network formulation. IEEE Trans Power Appar Syst. 1982;101(3):602–9. doi:10.1109/tpas.1982.317273. [Google Scholar] [CrossRef]

10. Wang M, Yang M, Fang Z, Wang M, Wu Q. A practical feeder planning model for urban distribution system. IEEE Trans Power Syst. 2023;38(2):1297–308. doi:10.1109/tpwrs.2022.3170933. [Google Scholar] [CrossRef]

11. Shi M, Zhang J, Ge X, Fei J, Tan J. Power optimization cooperative control strategy for flexible fast interconnection device with energy storage. Energy Eng. 2023;120(8):1885–97. doi:10.32604/ee.2023.025788. [Google Scholar] [CrossRef]

12. Shayeghi H, Davoudkhani IF. Fuzzy-stochastic planning and optimization of mobile solid-state hydrogen storage in renewable-powered distribution networks. Int J Hydrogen Energy. 2025;168:151048. doi:10.1016/j.ijhydene.2025.151048. [Google Scholar] [CrossRef]

13. Liu J, Ma L, Wang Q. Energy management method of integrated energy system based on collaborative optimization of distributed flexible resources. Energy. 2023;264:125981. doi:10.1016/j.energy.2022.125981. [Google Scholar] [CrossRef]

14. Heidari S, Fotuhi-Firuzabad M. Integrated planning for distribution automation and network capacity expansion. IEEE Trans Smart Grid. 2019;10(4):4279–88. doi:10.1109/tsg.2018.2855218. [Google Scholar] [CrossRef]

15. Yang Y, Pei W, Huo Q, Sun J, Xu F. Coordinated planning method of multiple micro-grids and distribution network with flexible interconnection. Appl Energy. 2018;228:2361–74. doi:10.1016/j.apenergy.2018.07.047. [Google Scholar] [CrossRef]

16. Jiang X, Zhou Y, Ming W, Yang P, Wu J. An overview of soft open points in electricity distribution networks. IEEE Trans Smart Grid. 2022;13(3):1899–910. doi:10.1109/tsg.2022.3148599. [Google Scholar] [CrossRef]

17. Li J, Zhang L, Zhang B, Tang W. Coordinated planning for flexible interconnection and energy storage system in low-voltage distribution networks to improve the accommodation capacity of photovoltaic. Glob Energy Interconnect. 2023;6(6):700–13. doi:10.1016/j.gloei.2023.11.004. [Google Scholar] [CrossRef]

18. Yan X, Wang Q, Bu J. High penetration PV active distribution network power flow optimization and loss reduction based on flexible interconnection technology. Electr Power Syst Res. 2024;226:109839. doi:10.1016/j.epsr.2023.109839. [Google Scholar] [CrossRef]

19. Qi Q, Wu J, Long C. Multi-objective operation optimization of an electrical distribution network with soft open point. Appl Energy. 2017;208:734–44. doi:10.1016/j.apenergy.2017.09.075. [Google Scholar] [CrossRef]

20. Wang Y, Su X, Song M, Jiang W, Shahidehpour M, Xu Q. Sequential load restoration with soft open points and time-dependent cold load pickup for resilient distribution systems. IEEE Trans Smart Grid. 2023;14(5):3427–38. doi:10.1109/tsg.2023.3262199. [Google Scholar] [CrossRef]

21. Zhang J, Wang T, Liao Z, Tang Z, Wang H, Yue J, et al. Flexible interconnection strategy for distribution networks considering multiple soft open points siting and sizing. Electr Power Syst Res. 2025;241:111335. doi:10.1016/j.epsr.2024.111335. [Google Scholar] [CrossRef]

22. Ebrahimi H, Shahnia F, Nikdel N, Galvani S. Renewable energy and demand uncertainty-aware stochastic allocation and management of soft open points for simultaneous reduction of harmonic distortion, voltage deviations and losses. Comput Electr Eng. 2025;123:110208. doi:10.1016/j.compeleceng.2025.110208. [Google Scholar] [CrossRef]

23. Bai L, Jiang T, Li F, Chen H, Li X. Distributed energy storage planning in soft open point based active distribution networks incorporating network reconfiguration and DG reactive power capability. Appl Energy. 2018;210:1082–91. doi:10.1016/j.apenergy.2017.07.004. [Google Scholar] [CrossRef]

24. Zare P, Dejamkhooy A, Davoudkhani IF. Efficient expansion planning of modern multi-energy distribution networks with electric vehicle charging stations: a stochastic MILP model. Sustain Energy Grids Netw. 2024;38:101225. doi:10.1016/j.segan.2023.101225. [Google Scholar] [CrossRef]


Cite This Article

APA Style
Ji, Y., Peng, Y., Ji, H., Na, X., Chen, Y. et al. (2026). A Multi-Stage Expansion Planning Method for Rural Distribution Networks with Flexible Interconnection. Energy Engineering, 123(8), 2. https://doi.org/10.32604/ee.2025.074599
Vancouver Style
Ji Y, Peng Y, Ji H, Na X, Chen Y, Li W, et al. A Multi-Stage Expansion Planning Method for Rural Distribution Networks with Flexible Interconnection. Energ Eng. 2026;123(8):2. https://doi.org/10.32604/ee.2025.074599
IEEE Style
Y. Ji et al., “A Multi-Stage Expansion Planning Method for Rural Distribution Networks with Flexible Interconnection,” Energ. Eng., vol. 123, no. 8, pp. 2, 2026. https://doi.org/10.32604/ee.2025.074599


cc Copyright © 2026 The Author(s). Published by Tech Science Press.
This work is licensed under a Creative Commons Attribution 4.0 International License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • 1286

    View

  • 419

    Download

  • 0

    Like

Share Link