Configuration Optimization of Cables in Ductbank Based on Their Ampacity

Electrical power companies are using more underground cables rather than overhead lines to distribute power to their customers. In practice, cables are generally installed in some compact ductbanks. Since the cost of underground cables is very expensive, using the entire space of a ductbank is extremely important. But such usage is limited due to the overheating of cables. Overheating is generally caused by overload, which means the carrying current exceeds the ampacity of a cable. The ampacity of a cable depends on not only the material and design of a cable but also the distance between different cables. Thus the configuration of cables determines the total ampacity value and the potential use of a ductbank. In this paper, the best configuration based on ampacity is achieved for a three-row, five-column ductbank that is buried at a depth of one meter below the earth’s surface. Both balanced and unbalanced scenarios are considered, and all cables have two available types to be selected.


Introduction
Underground cables have more advantages than overhead lines since cables offer better protection and are not as unsightly in appearance in urban areas.In practice, cables are generally installed in some compact ductbanks in order to provide easier installation of multiple cables in a concrete space [1], as shown in Figure 1.However, installation and maintenance of underground cables are a lot more expensive compared with overhead lines [2].According to [3], the cost of laying the cables is $125 -$200 per meter.Thus it is extremely critical to use the full potential of the ductbank.However such use is limited by the overheating of Figure 1.Cables in a ductbank for installation [9].cables.Overheating is the most significant factor in decreasing cable service life [4].Since cables are surrounded by soil instead of air, the speed of temperature diffusion is much slower than in air [4].The high currents carried by cable conductors are usually the cause of high temperatures.Overheating of cables generally results by overloading them [5].Each cable has a current limitation, called ampacity [4], that allows the cable to operate without problems.When the carrying current exceeds its ampacity, cable damage results, followed by failure that may be difficult to fix.Ampacity depends on the strength of the heat source, the material of cables, and the surrounding environment, including the ductbank and soil [5].When the thermal resistance of cable layers and soil is low, the heat can spread faster, and the ampacity of the cable is higher.Conversely, higher thermal resistance can mean lower ampacity value.
However, a cable's ampacity value is decided not only by its own characteristics but also by neighboring cables.The heat generated by one cable can influence the maximum value of the current of one nearby.This influence is called the mutual heating effect [1].In a ductbank, there are lots of available ducts that can be selected.So various cable configurations are possible.Different configurations cause different total ampacity value.The distance between two cables significantly influences ampacity value due to the mutual heating effect.So proper design of cable layout, i.e. using the entire space of a ductbank, can lead to maximum total current carrying capacity.Similarly, one cable configuration can offer only minimum total ampacity.This worst-case scenario is useful when a power system is being analyzed without knowing the exact layout of cables.
Thus, the configuration optimization of cables in a ductbank is hugely crucial.
Although some researchers have studied cable configuration optimization [1] [5] [6] [7] [8], they covered only the three-phase balanced condition and only one type of cable.However, most distribution systems are dealing with unbalanced loading, and the selection of cables must consider various types and designs.The objective of this paper is to figure out the best configuration for deli-Journal of Power and Energy Engineering vering more current, if needed, in one ductbank and avoiding overheating of the cables under both balanced and unbalanced conditions based on types of cables.

Ampacity Calculation
It is known that a cable's ampacity is based on the highest allowable temperature that cable can hold without overheating, and it is influenced by the mutual heating effect of nearby cables.To properly design a cable system and optimize cable configuration, calculating the ampacity value of various cables with different cross-sections and sizes is extremely important.Typically, an underground cable consists of four layers, including a conductor layer, insulation layer, shield layer, and jacket layer [10].This paper used COMSOL [11], which is a powerful multi-physics simulation software, to model a shielded cable, as shown in Figure 2.
Several publications proposed different methods to calculate cables' parameters and their ampacities for both single and multiple cable configurations [12] [13] [14] [15].Among these publications, two of them are widely used: the Neher and McGrath method [16] and IEC Standards 287-3-2 [17].These two methods are similar.They summarize all existing principles and equations to calculate cable ampacity in different conditions, including single cable, multiple cables without ductbank, and multiple cables with ductbank.These two methods are then classified and summarized by Dr. George J. Anders [4], as shown below.
In order to calculate the ampacity of cable i, a thermal circuit includes heating sources, and the thermal resistance of different layers is built based on the highest allowable cable temperature.
( ) where i I is the ampacity of cable i, j I is the carrying current of cable j, i R is the AC resistance of cable j, 1 λ is cable j's shield loss factor, 1 λ is cable jacket loss factor, n is the conductors number, j µ is the loss factor, dj W is the dielectric loss of cable j, T 1 , T 2 , T 3 , T 4 are the thermal resistance of different layers including ductbank and soil, max θ is the highest temperature that allows the cable to operate without problems, amb θ is the ambient temperature, and the reduction factor of conductor temperature.All these parameters depend only on the material and design of the cable and of the surrounding conditions [18].
The influence of mutual heating from nearby cables is corrected by

Optimization of Cables in a Ductbank
From the Equations ( 1)- (6) shown in Section 2.1, it can be noticed that to find the ampacity of cable i, carrying currents of all other cables should be pre-known, given the mutual heating effect.So if these equations are applied to all cables, a set of mutually interconnected equations is obtained.However, a set of interrelated equations is challenging to solve, and frequently, the iteration method can be used.But it is time-consuming, and not convergent in some conditions.So a more efficient method is the optimization method, which is recommended by Dr. Moutassem [5].Finding the ampacity value of each cable for a specific configuration could be described as an optimization problem.The objective function is the sum of all carrying currents.The constraints are that the temperatures of all cables are smaller than the highest allowable temperatures.In this paper, the same equations are used to decide the best configuration based on ampacity for cables in a ductbank assumed a constant frequency system.The detailed transformation steps of this optimization problem are summarized in Appendix B.
The optimization problem for multiple cables installed in a ductbank for a specific configuration can be summarized as follows: Journal of Power and Energy Engineering The objective function is The constraint for cable 1 is: Similarly, the constraints for the other cables are: Using MATLAB, the constraints can be acquired in a matrix form.
( ) where all elements in matrices c and d are calculated based on Equations ( 13)-( 14) in Appendix B and matrix c has one on its diagonal terms.
The procedure for finding the total ampacity value for a specific configuration of cables in a ductbank is completed.The next step is to find the configuration that leads to the maximum total ampacity value and minimum total ampacity value.The method applied in this paper includes three steps.Firstly, assume all ducts have their own cables with some initial guess as to current values.Secondly, randomly choose some of these cables to have current equal to zero, which means these ducts don't have cables installed in them.Thirdly, find the best and worst configuration that produces the maximum total ampacity value and minimum total ampacity value.But since, in this program, the types of cables should be selected automatically, one more step is added that introduces additional ducts for different cable types selection.The steps of configuration optimization of cables in a ductbank are shown in Figure 3.
In this paper, a three-row, five-column ductbank is selected.It is buried at a depth of one meter below the earth's surface.The distance between two ducts in the same row is 0.3 meter, and the distance between each row is 0.5 meter, which is shown in Figure 4.Both balanced and unbalanced conditions are considered.
In a balanced scenario, all cables are equally loaded.For an unbalanced scenario,

Configuration Optimization for a Balanced Condition
For two cables per phase, the second type of cable is selected, and the maximum ampacity of each cable is 655 A, as shown in Figure 5(a).The minimum ampacity of each cable is 559 A, as shown in Figure 5(b).This configuration is reasonable since all cables are located near each other.So the heat generated by one  carefully, an opposite conclusion can be reached.For example, the distance between the two cables of phase "a" in Figure 5(a) is longer than the distance between the two cables of phase "a" in Figure 5(c).Similar conclusions can be obtained in these two figures by comparing the distances between different cables.
The conclusion is that the configuration of Figure 5(a) produces larger total ampacity than the configuration in Figure 5(c).The detailed results for maximum ampacity value are shown in Table 1.
The temperature limitation of type two cable is 90˚C and the resulting temperatures of all cables are below 90˚C.Moreover, the temperatures of the cables are symmetric in Table 1 since the arrangement of cables is symmetric in the ductbank and the cables are in balanced condition.
For three cables per phase, the second type of cable is selected, and the maximum ampacity of every cable is 566 A, which is shown in Figure 6

Configuration Optimization for a Special Unbalanced Example
For unbalanced condition, a particular example:   For three unbalanced cables per phase, where c is the highest loaded phase, the best configuration based on ampacity is shown in Figure 8, and the detailed results are shown in Table 4.

Configuration Optimization for General Unbalanced Condition
If the highest loaded phase is changed from phase c to phase b and then to phase a, a general pattern for the unbalanced condition is observed, where H means highest loaded phase; L means lowest loaded phase; M means medium loaded phase.It is noticed that the best configuration for balanced condition and unbalanced condition based on ampacity is different according to Figure 5(a), Figure 6(a), Figure 9 and Figure 10.But the worst configuration is always arranging all cables near each other.

Conclusion
Installing cables in ductbanks occurs more and more frequently nowadays since their installation is easy.Use of the full potential of a ductbank is extremely important for reasons of economy.This paper proposes using the optimization method to find the best and worst configuration for cables in a ductbank based on ampacity.The best and worst configurations are decided for both balanced and unbalanced scenarios, which are different from common sense without optimization.For an unbalanced condition, a particular example is studied and   then extended to a general pattern for unbalanced cables in a ductbank.In the future, the impacts of optimization during abnormal condition will be discussed.
The study will include different faulted phases and loading conditions.Based on the results and conclusion, best case and cable configuration under abnormal condition will be presented.
Figure 10.The best configuration for general three unbalanced cables per phase.
( ) For all other cables, similar result equations can be obtained as well.
Let  (14) so that the ampacity calculation can be solved as an optimization problem.

Appendix C. Data of Two Types of Cables
Table C1.Data of first type cable [5].Table C2.Data of second type cable [5].

Figure 2 .
Figure 2. Common layers arrangement of cables simulated in COMSOL.
detail as a particular example of unbalanced cases.Then the general patterns for unbalanced conditions are also obtained.In this paper, cables have two available types that can be selected.The detailed data of these two types of cables are listed in Appendix C.

Figure 3 .
Figure 3. Procedure for configuration optimization of cables in a ductbank.
(a).The minimum ampacity of each cable is 495 A, which is shown in Figure 6(b).Similar to Figure 5(b), the configuration in which all cables are arranged together leads to the smallest total ampacity value.According to common sense without optimization, usually the configuration in Figure 6(c) is applied.For the configuration in Figure 6(c), the ampacity value of each cable is 541 A, which is smaller than the optimization result.The detailed results for maximum ampacity value are shown in

Figure 6 .
Figure 6.The optimization result compared with common sense for three balanced cables per phase.(a) Best configuration by optimization; (b) Worst configuration by optimization; (c) Common sense without optimization.
cables per phase, the best configuration is shown in Figure 7 when phase c is the highest loaded phase and phase b is the medium loaded phase.The second type of cable is selected, and the detailed results are shown in Table 3. Journal of Power and Energy Engineering

Figure 7 .
Figure 7.The best configuration for two unbalanced cables per phasein a particular example.

Figure 8 .
Figure 8.The best configuration for three unbalanced cables per phasein a particular example.

Figure 9 .
Figure 9.The best configuration for general two unbalanced cables per phase.

Table 2 .Table 1 .
Detailed results for two balanced cables per phase.

Table 2 .
Detailed results for three balanced cables per phase.

Table 3 .
Detailed results for two unbalanced cables per phase.

Table 4 .
Detailed results for three unbalanced cables per phase. 1