Modeling for Underground Cable Water Tree Growth Dynamics

Underground cable faults, whether transient or permanent, are traceable to insulation failure problems, most of which are water tree initiated. Insulation breakdown, which usually leads to costly power outages, may be prevented by taking pre-emptive actions. The most decisive pre-emptive action is one in which real-time tracking of water tree advancement within the cable insulation system is possible. Such pre-emptive actions, however, depend on accurate modeling of the phenomenon. Earlier water tree models are static in that they focused on the cable insulation property change at a time segment. Thus, they lack the properties needed for tracking water tree progress and for determining the onset of transient and permanent faults. This paper presents a new approach to water tree modeling, focused on insulation degradation geometry in the form of parabolic expansion of water tree. We developed a dynamic model centered on the computation of the capacitance of a vented water tree as a function of time. The dynamic model accounts for the time-dependence of the radial growth of the water tree to track insulation degradation. The model was tested in predicting cross-linked polyethylene (XLPE) cable’s insulation lifespan. The result was found to be within the range of the recorded lifespan of field aged cables in the literature. Also, performance comparison with an earlier analytical model validated with COMSOL Hyperphysics software shows a significant correlation between them.

critical nature of the problem to power system reliability, this subject has been vigorously researched for well over three decades. The researches were majorly focused on understanding the dynamics of water tree formation, the sequence of events leading to its morphology into full-blown insulation degradation, and the circumstances of insulation breakdown emerging from it. Mostly, the presence of micron-size particulate impurities within the body of insulation or at any of the interfaces, microvoids due to a manufacturing defect, irregularities at the internal interfaces leading to protrusions and local electric field enhancement, and cracks within the insulation system arising from aging are some of the significant causes identified in the literature [1] [2]. However, another common source of water treeing, which usually results in insulation breakdown over time, stems from external injuries to a cable's insulation during transport, installation, or maintenance. This occurrence, aided by the powerful electric field created when the cable is energized, encourages the seepage of moisture into the insulation system, thus commencing the water treeing phenomenon. Over time, as water tree growth advances towards the conductor, the thickness of sound insulation between it and the outer screen reduces. At this point, a transient condition, such as a switching surge generating a momentarily considerable electric stress can result in a breakdown, leading to an earth fault. The intricacies of the identified causes of water tree formation and growth may be summarized as the result of interplaying electromechanical and electrochemical forces within the polymeric insulation system. These forces are set in motion when the twin factors of moisture and a strong, time-varying electric field are present within the insulation systems.
As the demand for a very high degree of reliability and quality electric power continues to take central focus in the management of electric power systems, it becomes imperative to look more closely at the causes of failure in critical power infrastructures such as underground cable systems and how such failure can be prevented. Detecting early warning signs of failure, and being able to track its advancement are means of preventing such failures and the power outage that would have resulted from it. In this paper, we provide a focus on a new water tree modeling methodology that is supportive of cable insulation failure tracking and power outage prevention.

XLPE Cable Construction and the Water Treeing Phenomenon
The Cross-linked Polyethylene (XLPE) cable is popular in medium-high voltage power transmission and distribution applications. Figure 1 shows the cross-sectional feature of a single core XLPE cable. A detailed description of the components and their functions are provided by [1]. The strength of the cable's polymeric insulation lies in the cross-linking process, which modifies the polymer's basic structure and gives it much more improved properties. The end-product is stronger, more tolerant of excessive temperature resulting from overloading in properties that are supportive of applications in a wide range of operational voltage. These improvements, therefore, place it at an advantage, especially in applications requiring operational voltages that are up to 500 kV. In spite of its many advantages, XLPE cables are still susceptible to the water treeing problem which is a major cause of insulation failure. There are three main types of water trees, as illustrated in Figure 2(a). The first is the bow-tie water trees, which are so named because of its resemblance to a bow-tie. They are usually formed within the insulation and constitutes no threat of insulation breakdown. The second and most troublesome class of all is the vented water tree, which may be formed at any of the interfaces. It is the main culprit in most insulation breakdown occurrences, mainly when formed at the outer (insulation-screen) interface [3]. Figure 2(b) illustrates how a vented water tree can result in insulation breakdown. Bush trees, which is an aggregation of many vented water trees, is the third class. It may be formed at any of the interfaces and constitute no risk of insulation breakdown.

Water Tree Modeling
The study of water trees for preemptive actions to forestall insulation breakdown requires adequate modeling of the phenomenon. Three basic approaches have been used in the literature. They are analytical approaches, use of modeling and simulation software such as COMSOL Hyperphysics, and laboratory-moderated accelerated insulation aging processes [4] [5]. The analytical approach to water tree modeling requires the identification of a representative physical geometry of   [7]. Using such geometries, researchers have sought to model the problem as a changing permittivity or conductivity of insulation material with the growth of water trees or both.
The nature of modern power systems and the demand for a very high degree of reliability requires proactive approaches to dealing with faults to ensure system protection, safety, and service continuity. Consequently, of late, research activities aimed at modeling water trees for practical applications such as fault diagnosis in underground power systems were carried out by [8] [9] [10], among many others. However, the earlier water tree models are unsuitable for tracking water tree growth, identifying potential signs of impending insulation breakdown, and preventing it because of their staticity. This is because those models can only provide information about the state or behavior of failing insulation at a time segment. Because of this shortcoming, we introduce in this paper, the concept of water tree dynamic model. The dynamic model allows the computation of water tree capacitance as a function of time. Hence, it provides an insight into the state of an underground cable's insulation at every time segment during its lifetime. The steps taken to develop this model are presented in this paper.

Nature and View of Water Trees
Water trees are dendritic growths of micro-defects within the insulation of polymeric cables. [11] proposed and experimentally verified a mechanical damage theory of polymeric insulation. His result shows that the dendrites are made up of ellipsoidal pockets of water of microscopic sizes connected by tubular channels, which, under a microscope, looks like strings of pearls ( Figure 3). Each microvoid is separated from the next by about 3 μm of tubular channels with an approximated radius of 40 nm. Both the channels and the ellipsoidal microvoids are usually filled with water. Each of the pearly string of water ellipsoids is surrounded by XLPE insulation whose varying degrees of degradation determines the overall electrical permittivity of the area. A degraded area of insulation that contains these dendritic strings of water micro-pockets is referred to as a water tree. The overall approximate geometrical structure depends on the magnitude of the electric field, the system operating frequency, and the water content [12]. While the geometrical feature of water trees can be quite random, it has been discovered that vented water trees that grew in field aged cables which were operated at power frequency (50 -60 Hz) and high voltages usually maintain approximately parabolic geometry as illustrated in Figure 4(a) [12] [13]. This shape illustrates how an inward-growing vented water tree can be a potential source of insulation breakdown. The approximate geometry of the vented water tree which is reproduced in Figure 4(b) has three major parts, namely the root, the body, and the front. The modeling activities presented in the following sections of this paper are motivated in part by this geometrical feature.

Conceptualization of a Water Tree Geometry for Dynamic Modeling
As earlier proposed, the cross-sectional view of a vented water tree, as seen in Figure 3 and Figure 4, is approximately parabolic. According to [14], vented water tree growth progresses in cyclic layers. Following this report, each new layer of growth would be approximately parabolic in shape resulting in an axial growth towards the conductor with each layer. In building the vented water tree dynamic model ( Figure 5), which illustrates a simplified and enlarged sector of the  presented. To simplify the computations, the water treed section is approximated to the section ABCD while the section DEFC is undegraded XLPE between the water treed section and the conductor O. The radial distance NO and PO of the root and the front of the water tree from the center of the conductor is denoted by w r and pe r respectively, while the angle subtended at O by the water treed section is o θ . The water treed region averagely has a higher dielectric permittivity w ε than that of the surrounding insulation material, which has a constant dielectric permittivity pe ε [15]. Under the influence of a time-varying electric field, the water treed region exhibits a dynamic permittivity behavior, which is dependent on the density of the network of dendritic water inclusions in the region. This difference in permittivity between the water treed region and pure XLPE confers on each of them a unique capacitive property.
Hence, the two regions may be considered as two different capacitors w C and pe C connected in series. The water tree dielectric permittivity w ε , radial growth w l , and the angular expansion θ all changes with time as the region grows. The changes in these three quantities with respect to time accounts for the dynamic shape change of the water tree. It is on this shape change principle that the derivation and calculation of water tree dynamic capacitance are based.

Capacitance of Vented Water Tree
For insulation failure studies, it is customary to model insulation degradation as some measurable quantity whose value during the degradation process, may be S. Famakin, C. Kim Journal of Power and Energy Engineering computed. Since the insulation failure process results in a progressive change in the permittivity of its material, the capacitance of a failing insulation has been used as a measurable index of the extent of the failure process. The overall capacitance of a cable during the failure process is determined by the capacitance of the water treed region [7]. This capacitance and its change with respect to time is a function of the number and advancement of water trees within the insulation system. Consequently, it's successful modeling as a function of time is a means of computing and tracking the insulation degradation process.
Regarding the geometry of Figure 5, when a water treed cable is energized, the conductor's surface, water tree/insulation interface, and water tree root assume potentials c V , pe V , and w V , respectively. The potential difference d V between pe V and w V is given by: The electric field E at any point r within the water treed space having a dielectric permittivity w ε and bounded by the potentials pe V , and w V is given by: where λ is the charge per unit length of the region. Substituting (2) into (1), and integrating with respect to r in the limit w pe r r r ≤ ≤ , the potential difference d V is given by where NO Hence, the capacitance w C of the water treed region is given by As shown in Figure 5, , , c w l r l , and 0 ε are insulation thickness, conductor radius, water tree length, and permittivity of free space = 8.85 × 10 −12 F/m, respectively. Also, it may be seen that As may be observed from (6) we proceed to derive explicit expressions for each in the following sections. Figure 6 shows the simplified geometry of a vented water tree growing within a section of XLPE insulation and the variation of the permittivity with the growth.

Water Tree Permittivity
According to [7] [16], and [17], water inclusions and the interlinking channels making up the water tree are densest at the root. That density decreases towards the front of the water tree. Similarly, the relative permittivity has been found to follow the same trend. According to [16], it decreases exponentially from a value Figure 6. Water tree permittivity variation with length.
S. Famakin, C. Kim Journal of Power and Energy Engineering of about 6.9 at the root to 2.3 at the front. This value has been adopted in this paper. The argument for a value of 6.9 appears to be more plausible in view of the fact that water treed insulation only contains a microscopic quantity of moisture and remains predominantly insulation.
In modeling the relationship between water tree capacitance and water tree length, [7] also used COMSOL Hyperphysics software to develop a graphical illustration that shows that water tree capacitance indeed varies with the length of water trees as a decreasing exponential function. We mathematically present this as: Then, thus, we write By applying (14) to (13), we obtain 6.9 A = and 1.10 k = . Consequently, where w l′ is the normalized water tree length given by

Water Tree Radial Growth
As earlier proposed, a vented water tree is approximately parabolic. The general expression for such a parabola takes the form The constants a and c controls the width along the x-axis and depth along the y-axis of the parabola, respectively. Applying this fact to water tree growth, it implies that a controls the expansion of the root of the water tree while c controls its radial growth. Since each of these growths takes place with respect to time, (17) may be recast as Based on the description and images of real-time vented water tree growth seen under the microscope and described by [14] and [18], radial growth was found to be faster than the angular expansion. Guided by this description, we approximated functions of best fit for the two coefficients as: and ( ) ( ) where R is the radius of the circle described by the cable's circular cross-section.
Hence, (18) Equation (21) is the vented water tree shape-changing equation describing water tree growth with time. Using (21), water tree length at a given time may be computed from y R x = − made by the cable's cross-section is equated to (18), which describes the water tree at any given time. The intersection points ( )

Water Tree Circumferential Growth
We obtain solutions 1 x and 2 x to (23) as expressed in (24) where 1 x is positive, and 2 x is negative.
Considering only 1 x , then 1 y of the intersection, ( )   breaches the insulation entirely and reaches the conductor (with its length the same as the insulation thickness), which is the condition for insulation breakdown. Disregarding the stochasticity that is peculiar to the water treeing problem, our shape-changing equation was able to predict an inception-insulation breakdown time, which falls within the range reported in the literature by [19].
A further performance test of the dynamic models was done by comparing it graphically with two other models (Chen's and Stucki's) [7] [20], that were found in the literature. Figure 10 illustrates this.
While Chen's model was built analytically on the properties of XLPE cables just like our dynamic model (DM), Stucki's model was developed from measurements of water tree capacitance and its length. It was carried out on the field aged samples of Polyethylene (PE) insulation from a cable retired from service.
Our model follows the same trend as Chen's even though there is a difference of an order of magnitude between them. This may be attributed to the fact that the authors used a relative permittivity of 88 for the computation of water tree capacitance, while a value of 6.9 was used in this paper for reasons earlier stated.  Apart from this difference, we can say that there is a significant correlation between the two models. However, our model is uncorrelated with Stucki's. This may be because PE is known to be prone to abrasion and is less stable in the presence of chemical agents than XLPE. This fact may affect water treeing dynamics in the material and the value of the measured capacitances. However, the operating environment of the cable before its retirement from service, the circumstances of its retirement, and other factors related to its years of service is unknown. Hence, no firm conclusion can be drawn to justify its behavior when compared with our dynamic model and Chen's model.

Conclusion
This paper has presented a detailed methodology for developing a dynamic growth model of water tree for insulation degradation tracking and breakdown anticipation. The model formulates a water tree shaping function of time, which led to two time-varying predictive variables, namely, water tree length and water tree capacitance. A significant strength of the model is in the computation of a cable's insulation end of life. In the future, the randomness in the growth and shape of water trees may be factored into its modeling and considered as an area of further studies to improve it. The development of this model is a precursor to additional studies and further modeling activities that are expected to result in the development and testing of methods and systems for anticipation and detection of self-clearing fault conditions in underground cable systems. This condition, which usually precedes final insulation failure is highly evasive. Identifying it and tracking its onset is useful in preventing unannounced underground feeder outages that usually arise from insulation's final failure.

Conflicts of Interest
The authors declare no conflicts of interest regarding the publication of this paper.