Next Article in Journal
Techno-Economic Assessment of a Standalone Hybrid System Using Various Solar Tracking Systems for Kalpeni Island, India
Previous Article in Journal
Review of the Hydrogen Evolution Reaction—A Basic Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Induced Voltage on Pipeline Located Close to Parallel Distribution System

1
Samsung Electronics Device Solutions, Pyeongtaek 17786, Korea
2
Department of Electrical and Computer Engineering, Sungkyunkwan University, Suwon 440-746, Korea
3
Smart Grid Research Center, Korea Electrotechnology Research Institute (KERI), Gwangju 61751, Korea
4
Department of Mechanical Engineering, University of Canterbury, Christchurch 8041, New Zealand
5
Department of Electronic and Electrical Engineering, Sungkyunkwan University, Suwon 440-746, Korea
*
Author to whom correspondence should be addressed.
Energies 2021, 14(24), 8536; https://doi.org/10.3390/en14248536
Submission received: 22 October 2021 / Revised: 25 November 2021 / Accepted: 3 December 2021 / Published: 17 December 2021
(This article belongs to the Topic Power Distribution Systems)

Abstract

:
Gas and oil pipelines are widely used to supply customers. They are often laid in parallel with high voltage power lines, sharing right of way. When the pipeline is located near overhead high-voltage power, corrosion caused by induced voltages from AC power lines can occur in utility pipelines. Therefore, the calculation of induced voltage is always required for both interference and maintenance workers as well as pipeline facility safety. For parallel distribution lines, the calculation method for the induced voltage is not suitable due to the excessive error caused by using the screening factors of the neutral current. For a more practical analysis of the induced voltage, a new analysis is needed using the actual neutral current and which also considers the overhead ground conductor and the neutral conductor. This paper analyzed the induced voltage from a parallel distribution system using Carson’s formula and vector analysis. Simulation analysis results are verified by separate Electromagnetic Transient Program (EMTP) simulation.

1. Introduction

There are important issues to consider by both electrical utilities and pipeline operators when an underground pipeline is located proximal to the ground level in overhead high-voltage power lines [1,2,3,4,5,6], namely safety and corrosion. Due to the corrosion of gas and oil pipelines caused by induced voltages from AC power lines [2,3,4,5], a pipeline AC mitigation system is needed to avoid costly repairs and maintenance. The primary focus of mitigation is to reduce the induced voltages at normal and fault conditions according to the NACE standard [6] for the safety of maintenance operators and damage in pipeline facilities. This corrosion can also lead to an equipment malfunction [7,8].
The interferences through which the AC power lines generate induced voltage and current on nearby pipes are inductive coupling and conductive coupling. The inductive coupling is related to the electromotive force (EMF) induced in the pipe–earth circuit, while the conductive coupling is caused by the injection of current into the soil by a transmission line during phase-to-ground fault conditions or subjected to currents [9].
The interference issue of induced voltage on pipelines parallel to power lines has been treated a lot in the literature [10,11,12]. Recent literature related to this issue was published by A. Popoli [13], who analyzed interference caused by the vicinity of a metallic pipeline buried in the soil to a transmission line using parameters calculated with the finite difference method (FEM). G. Lucca [9] analyzed the induced voltage and current from a transmission line in fault condition on a nearby buried pipeline with the EMF. However, these works focus on the interference of transmission lines and pipelines.
There are many reports on the electromagnetic interference (EM) of the induced voltage using the Electromagnetic Transient Program (EMTP) [14,15,16,17]. A practical calculation method using EMTP is presented in [14], which assumes multi-conductor power lines with a variety of tower types (complex approximation layouts), and multilayered soil models are parallel to a nearby pipeline. Heo et al. [15] analyzed the induced voltage generated by overhead transmission lines under normal cases using the line and cable constant (LCC) module. Isogai and Ametani [16] analyzed the effect of induced voltages to the pipeline on power lines. The results show the three-phase horizontal lines tend to induce higher voltages than vertical lines when the same system parameters for the transmission line and ground are used.
In order to reduce the electromagnetic interference between power lines and pipelines, the most practical solution consists in worsening the electromagnetic coupling. Then, “the use of screening conductors buried in the soil above the metallic pipelines appears as a feasible solution to attenuate the voltage generated on the pipeline” by power lines [18]. Especially for designing an optimized mitigation system, a proper calculation is required to accurately predict the induced voltage on metallic pipelines.
Previous research on induced voltage is mainly focused on overhead transmission lines [8,9,10,11,12,13,14,15,16]. The induced voltage on the pipeline is critical because the unbalanced current flow in the distribution system can increase the induced voltage on the pipeline. Recently, it has been relatively common to find double-circuit lines (DCL) where distribution lines are physically parallel due to their significant environmental and economic advantages over single-circuit lines (SCL). Thus, it is necessary to analyze induced voltage from DCL [19,20,21].
The calculation method in [22] is not suitable for parallel distribution lines due to the excessive error caused by using a constant screening factor of the neutral current. The constant value of the screening factors that are currently used in the calculation of the induced voltage in distribution lines is the main cause of error [23,24], as it cannot consider unbalanced loading conditions in the distribution line. Thus, a new practical method is needed for analyzing the induced voltage using the actual neutral current without using the screening factors from the overhead ground conductor and the neutral conductor.
Two types of overhead distribution lines, SCL and DCL, are considered in this paper. These two types of distribution line circuits are only available from the Korea Electric Power Corporation (KEPCO) in Korea. The modeling and calculation methods of induced voltages using the actual neutral current from SCL and DCL distribution line models are presented in Section 2. The simulation results, vector analysis based on the simulation results, and comparison of case studies are demonstrated in Section 3. The conclusions are presented in Section 4.

2. Voltage Calculation between Pipeline and Distribution System

2.1. Induced Voltage Calculation in SCL

Figure 1 depicts an SCL overhead distribution laid parallel to the pipeline line. The ground wire, G indicates the overhead ground conductor. Subscripts A, B, C, N, and P indicate the conductor phase A, B, C, neutral, and pipeline, respectively.
An induced voltage on the pipeline line can be calculated using Carson’s formula [9,23]. The relationship between voltage and current from Figure 1 is defined:
V G V A V B V C V N V P = Z GG Z GA Z GB Z GC Z GN Z GP Z AG Z AA Z AB Z AC Z AN Z AP Z BG Z BA Z BB Z BC Z BN Z BP Z CG Z CA Z CB Z CC Z CN Z CP Z NG Z NA Z NB Z NC Z NN Z NP Z PG Z PA Z PB Z PC Z PN Z PP I G I A I B I C I N I P
where the calculation of neutral current assumes that the neutral grounding point of the 3-phase lines is common [19].
The self and mutual impedance are obtained by Carson’s formula. The mutual impedance of the system can be calculated as a lump sum of impedances between each conductor and the pipeline [23,25].
Z i i = R + 1.588 e 3 f + 2.022 e 3 f ln 1 G M R + 7.6786 + 1 2 ln ρ f Z i j = 1.588 e 3 f + 2.022 e 3 f ln 1 D i j + 7.6786 + 1 2 ln ρ f
where
  • R = resistance of the conductor;
  • GMR = geometric mean radius of the conductor;
  • Ρ = earth resistivity;
  • f = system frequency;
  • Dij = distance from conductor i to conductor j;
  • Zii = self-impedance of the conductor I;
  • Zij = mutual impedance between conductors i and j.
The voltage of the overhead ground conductor (VG) is expressed:
V G = Z GG I G + Z GA I A + Z GB I B + Z GC I C + Z GN I N = Z GG I G + Z GABC I ABC + Z GN I N
where
Z GABC I ABC = Z GA I A + Z GB I B + Z GA I C I N = I A + I B + I C
Overhead ground conductors are on top of the distribution lines to avoid lightning strikes. VG can be assumed to be zero [15] because it is grounded. Accordingly, the current of the overhead ground conductor, IG can be obtained on the
I G = Z GABC I ABC + Z GN I N Z GG
Then, the induced voltage, VP on the pipeline, can be calculated:
V P = Z PG I G + Z PA I A + Z PB I B + Z PC I C + Z PN I N = Z PG I G + Z PABC I ABC + Z PN I N = V G + V ABC + V N
where
Z PABC I ABC = Z PA I A + Z PB I B + Z PC I c

2.2. Induced Voltage Calculation in DCL

Figure 1 depicts the SCL overhead distribution laid parallel to the pipeline line. The ground wire, G, indicates the overhead ground.
Figure 2 depicts the DCL overhead distribution system laid parallel to the pipeline line. A matrix extension [23] is applied to calculate the induced voltage in the DCL. The subscripts U and L indicate the upper and lower 3-phase lines, respectively.
Then, the pipeline laid parallel to the overhead distribution line with the DCL can be expressed:
Energies 14 08536 i001
Mutual impedance matrices (ZXY) in Equation (6) represent the parallel overhead distribution lines and the pipeline. These matrices are applied to calculate the induced voltage of the pipeline and take into account the coupling with adjacent parallel distribution lines. Thus, for example:
V G = Z GG I G + Z GUA I UA + Z GUB I UB + Z GUC I UC + Z GLA I LB + Z GLB I LB + Z GLC I LC + Z GN I N = Z GG I G + Z GU I U + Z GL I L + Z GN I N
where
Z GU I U = Z GUA I UA + Z GUB I UB + Z GUC I UC Z GL I L = Z GLA I LB + Z GLB I LB + Z GLC I LC
VG can be assumed to be zero in the DCL configuration [21], so the current of the overhead ground conductor can be determined:
I G = Z GU I U + Z GL I L + Z GN I N Z GG
For DCLs, we assume the upper and lower side use the common ground of the neutral line, which is equal to the total current of the upper and lower sides using the superposition method [19]. Then, the neutral current can be calculated:
I N = I UA + I UB + I UC + I LA + I LB + I LC
Then, the pipeline-induced voltage (VP) in the DCL can be calculated using this result:
V P = Z PG I G + Z PU I U + Z PL I L + Z PN I N = V G + V U + V L + V N
where
Z PU I U = Z PUA I UA + Z PUB I UB + Z PUC I UC Z PL I L = Z PLA I LA + Z PLB I LB + Z PLC I LC

2.3. System Modeling

Figure 3 illustrates the two types of overhead distribution line (ODL) configurations with the pipeline. In this paper, to analyze the induced voltage on a pipeline located close to overhead distribution lines, the ODL is implemented in both SCL and DCL configurations, as shown in Figure 3.
Both SCL and DCL types are classified by existing distribution line types of the Korea Electric Power Corporation (KEPCO) [19]. The data obtained from the pipelines and conductors are applied in this study. The key details are listed in Table 1 and shown in Figure 4 [26].
All case studies are analyzed based on a 1 km distance of the parallel distribution lines. The system frequency, soil resistivity, and power factors are 60 (Hz), 100 (Ω m), and 0.9 (pu), respectively.

3. Simulation and Results

The induced voltage on the pipeline is analyzed using EMTP based on various conditions: unbalanced loading conditions, different pole types, and the separation distance between the ODL and the pipeline. In the case studies, VP is a calculation value from (5), (10), VEMTP is the EMTP simulation result, and the difference will be calculated based on VEMTP. The unbalanced loading ratio between phases is kept to be less than 30% in the distribution.

3.1. Simulation Results in SCL

The EMTP simulation result of SCL cases is shown in Table 2. For the SCL case study shown in Table 2, Case 1 has a balanced load, Cases 2A–2C have unbalanced single-phase loads, and Case 3 has an unbalanced three-phase load.
Table 2 shows that the difference between the calculated value and the simulation result is less than 0.3% error. The maximum and minimum induced voltage to the pipeline are found in Case 2A and Case 2C, respectively. The induced voltage in Case 2C is nearly 50% of the induced voltage from Case 2A. The simulation results indicate that the induced voltages are dependent on the load condition, as shown in Figure 3a.
By using Equation (5), the induced voltage of VG, VABC, and VN are calculated and shown in Table 3. Due to the difference in direction between the VABC and VG + VN phase angles, VABC is defined as inducing side voltage caused by a three-phase current, and VG + VN is defined as a shielding side voltage caused by the current of an overhead ground conductor and a neutral conductor.
The SCL induced voltages in Table 3 are classified as shielding side, inducing side, and induced voltages, as shown in Table 4.
For single-phase unbalanced load cases, the induced voltage of Case 2C is nearly half of the induced voltage in Case 2A and 2B. This difference is due to the phase angle of the neutral current, which is closely arranged to the C-phase line and the pipeline. The shielding effect of the neutral current was verified by both calculation and vector analysis, as shown in Figure 5, where the shielding effect is multiplied by five for illustration purposes. The resultant vector VP, VP is much smaller than VABC and VG + VN in both Case 2A and Case 2C.

3.2. Simulation Results in DCL

For the case study of DCL, three-phase lines added from the existing SCL to the DCL were divided into the upper and the lower side, as shown in Figure 3b. The load conditions of the upper and lower sides were changed and are defined:
  • Case 4: Three-phase balanced load in the upper and lower sides.
  • Case 5: Single-phase (Uc, Lc-phase) unbalanced load in the upper and lower sides.
  • Case 6: Three-phase unbalanced load in the upper and lower sides.
  • Case 7A–7C: Single-phase unbalanced load on the upper side and balanced load on the lower side.
  • Case 8A–8C: Balanced load on the upper side and single-phase unbalanced load on the lower side.
The results for VP, which is the calculated induced voltage using Equation (10) and EMTP simulation, are compared in Table 5. The simulation results show the minimal differences between the calculated values and simulation results.
In Table 6, the induced voltage VG, VU, VL, and VN in DCL cases can be derived using Equation (10). However, the inducing side voltage is further categorized as an upper side voltage (VU) or a lower side voltage (VL). VU is caused by an upper three-phase current, VL is caused by a lower three-phase current, and VN is caused by neutral current (sum of the three-phase current on the upper and lower sides) [19].
The induced voltages in DCL configurations in Table 6 are classified as the shielding side, inducing side, and induced voltages, as with SCL cases. Due to the difference in direction between VU + VL and VG + VN phase angles, the induced voltages in the DCL case are classified as the shielding side and inducing side voltages, as shown in Table 7.
The vector analysis of the induced voltage VP is shown in Figure 6. As the VP vectors in Case 7B and Case 8B depict, they are very small compared to VU + VL and VG + VN.

3.3. Comparison of Case SCL and DCL Studies

Figure 7 shows the induced voltages in SCL and DCL configurations with the same load condition. The induced voltage in the SCL case is larger than the induced voltage in the DCL case, except for the C-phase unbalanced load case. The Case 2C and Case 5 results indicate that the induced voltage in the DCL case is larger than the inducing side voltage and VN (shielding side), as shown in Table 8.
Figure 8 shows that the induced voltage with the upper and lower sides of the DCL carry the same single-phase unbalanced load. The induced voltage in the DCL configuration is different than in the SCL configuration due to inductive coupling on the pipeline located at different single-phase unbalanced loads.
The induced voltage of the lower side in Case 8A–8C is larger than the induced voltage of the upper side in Case 7A–7C, due to the effect of lower side-induced voltages, which are closely arranged to the pipeline. However, the result indicates that the lower single-phase unbalanced load is smaller than the upper single-phase unbalanced load. Since the inducing and shielding side voltage magnitudes in Case 8 are smaller than the voltage magnitude in Case 7, the simulation results agree with the vector analysis shown in Table 7.
Table 9 and Table 10 show EMTP simulation results by separation distance in the SCL case. The separation distance is ±90 m between 1 km lengths of the parallel overhead distribution lines.
The induced voltage due to the inductive coupling on the pipeline located at a different unbalanced load of single-phase simulation cases, Case 2A–2C, is shown in Figure 9. The simulation results indicate that the induced voltage is negligible at the center of the ODL and the peak, where the pipeline is located at a separation distance near ±15 m, as shown in Figure 9 and Table 10.
Table 11 shows EMTP simulation results by separation distance in the DCL case. Induced voltage due to the inductive coupling on the pipeline located at a different unbalanced load of single-phase (Case 7A–7C) is shown in Figure 10.
The induced voltage has a peak value at the center of the ODL with DCL and gradually decreases as the transverse position of the pipeline increases. Therefore, when comparing SCL and DCL cases by the separation distance, SCL shows asymmetric characteristics about the y-axis, whereas DCL shows asymmetric characteristics.

4. Conclusions

This paper analyzed the induced voltage from a parallel distribution system using Carson’s formula and vector analysis. The calculation method for the induced voltage is not suitable due to the excessive errors caused by the screening factors of the neutral current.
For a more practical analysis of the induced voltage, a new analysis is needed to use the actual neutral current and to consider the overhead ground conductor and the neutral conductor in the distribution system. The shielding effect of the neutral current was verified by both calculation and vector analysis. Thus, the proposing method is demonstrated without using the screening factor.
The calculation results are verified by the EMTP simulation and vector composition method. Moreover, various case studies are analyzed according to the load conditions, separation distance, and pole types of the distribution system. In addition, the induced voltage of the DCLs has fewer EM effects than the induced voltage of the SCLs in the same load condition. From the results that did not generate very much of an error, it is expected that the proposed method can be useful to calculate the induced voltage on pipelines located close to the parallel distribution system.

Author Contributions

Conceptualization, methodology, C.-H.K. and H.-S.K.; software, validation, H.-S.K.; formal analysis, H.-S.K. and C.-H.K.; data curation, H.-S.K.; Writing—Original draft preparation, H.-S.K. and H.-Y.M. and J.G.C.; Writing—Review and editing, H.-S.K. and H.-Y.M. and J.G.C.; supervision, C.-H.K. and J.G.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Adedeji, K.; Ponnle, A.; Abe, B.; Jimoh, A.; Abu-Mahfouz, A.; Hamam, Y. A review of the effect of ac/dc interference on corrosion and cathodic protection potentials of pipelines. Int. Rev. Electr. Eng. 2018, 13, 495–508. [Google Scholar] [CrossRef] [Green Version]
  2. Christoforidis, G.C.; Labridis, D.P.; Dokopoulos, P.S. A hybrid method for calculating the inductive interference caused by faulted power lines to nearby buried pipelines. IEEE Trans. Power Deliv. 2005, 20, 1465–1473. [Google Scholar] [CrossRef]
  3. Gummow, R.A.; Segall, S.M.; Fieltsch, W. Pipeline AC mitigation misconceptions. In Proceedings of the Northern Area Western Conference, Calgary, AB, Canada, 15–18 February 2010. [Google Scholar]
  4. Abdel-Gawad, N.M.K.; El Dein, A.Z.; Magdy, M. Mitigation of induced voltages and AC corrosion effects on buried gas pipeline near to OHTL under normal and fault conditions. Electr. Power Syst. Res. 2015, 127, 297–306. [Google Scholar] [CrossRef]
  5. Kaboli, A.; Hr, S.; Al Hinai, A.; Al-Badi, A.; Charabi, Y.; Al Saifi, A. Prediction of Metallic Conductor Voltage Owing to Electromagnetic Coupling Via a Hybrid ANFIS and Backtracking Search Algorithm. Energies 2019, 12, 3651. [Google Scholar] [CrossRef] [Green Version]
  6. NACE SP0177-2014. Mitigation of Alternating Current and Lightning Effects on Metallic Structures and Corrosion Control Systems; Nace International: Houston, TX, USA, 2014. [Google Scholar]
  7. International Organization for Standardization. Corrosion of Metals and Alloys-Determination of AC Corrosion-Protection Criteria; Standard EN ISO 18086:2017; International Organization for Standardization: Geneva, Switzerland, 2017. [Google Scholar]
  8. Ouadah, M.; Touhami, O.; Ibtiouen, R.; Benlamnouar, M.F.; Zergoug, M. Corrosive effects of the electromagnetic induction caused by the high voltage power lines on buried X70 steel pipelines. Int. J. Electr. Power Energy Syst. 2017, 91, 34–41. [Google Scholar] [CrossRef]
  9. Lucca, G. AC interference from faulty power cables on buried pipelines: A two-step approach. IET Sci. Meas. Technol. 2021, 15, 25–34. [Google Scholar] [CrossRef]
  10. Martins-Britto, A.G.; Rondineau, S.R.M.J.; Lopes, F.V. Power Line Transient Interferences on a Nearby Pipeline Due to a Lightning Discharge. In Proceedings of the International Conference on Power Systems Transients (IPST 2019), Perpignan, France, 16–20 June 2019. [Google Scholar]
  11. Martins-Britto, A.; Moraes, C.M.; Lopes, F.V.; Rondineau, S. Low-frequency Electromagnetic Coupling Between a Traction Line and an Underground Pipeline in a Multilayered Soil. In Proceedings of the 5th Workshop on Communication Networks and Power Systems (WCNPS 2020), IEEE Xplore, Brasília, Brazil, 12–13 November 2020; pp. 1–6. [Google Scholar]
  12. Peppas, G.D.; Papagiannis, M.-P.; Koulouridis, S.; Pyrgioti, E.C. Induced Voltage on an Aboveground Natural Gas/Oil Pipeline Due to Lightning Strike on a Transmission Line. In Proceedings of the 2014 International Conference on Lightning Protection (ICLP), Shanghai, China, 11–18 October 2014; pp. 461–467. [Google Scholar]
  13. Popoli, A.; Cristofolini, A.; Sandrolini, L. A numerical model for the calculation of electromagnetic interference from power lines on nonparallel underground pipelines. Math. Comput. Simul. 2021, 183, 221–233. [Google Scholar] [CrossRef]
  14. Moraes, C.M.; Martins-Britto, A.G.; Felipe, V. Electromagnetic Interferences Between Power Lines and Pipelines Using EMTP Techniques. In Proceedings of the 2020 Workshop on Communication Networks and Power Systems (WCNPS) Communication Networks and Power Systems (WCNPS), 2020 Workshop, Brasilia, Brazil, 1–6 November 2020. [Google Scholar]
  15. Heo, J.Y.; Seo, H.C.; Lee, S.J.; Kim, Y.S.; Kim, C.H. A Simulator for Calculating Normal Induced Voltage on Communication Line. J. Electr. Eng. Technol. 2014, 9, 1394–1400. [Google Scholar] [CrossRef] [Green Version]
  16. Hiroshi, I.; Akihiro, A.; Yuji, H. An investigation of induced voltages to an underground gas pipeline from an overhead transmission line. Electr. Eng. 2008, 164, 43–50. [Google Scholar]
  17. Leuven EMTP Center. Alternative Transients Program ATP Rule Book; EEUG, Canadian/American EMTP User Group: Portland, OR, USA, 1992. [Google Scholar]
  18. Popoli, A.; Sandrolini, L.; Cristofolini, A. Reduction in the Electromagnetic Interference Generated by AC Overhead Power Lines on Buried Metallic Pipelines with Screening Conductors. Electricity 2021, 2, 316–329. [Google Scholar] [CrossRef]
  19. Park, K.W.; Seo, H.C.; Kim, C.H.; Kim, C.H.; Jung, C.S.; Yoo, Y.P.; Lim, Y.H. Analysis of the Neutral Current for Two-Step-Type Poles in Distribution Lines. IEEE Tran. Power Deliv. 2009, 24, 1483–1489. [Google Scholar] [CrossRef]
  20. Oka, K.; Koizumi, S.; Oishi, K.; Yokota, T.; Uemura, S. Analysis of a neutral grounding method for a three phase four-wire 11.4 kV distribution system. In Proceedings of the Transmission and Distribution Conference and Exhibition 2002: Asia Pacific. IEEE/PES, Yokohama, Japan, 6–10 October 2002; Volume 2, pp. 998–1003. [Google Scholar]
  21. Kim, H.S.; Rhee, S.B.; Yeo, S.M.; Kim, C.H.; Lyn, S.H.; Kim, S.A.; Weon, B.J. Calculation of an Induced Voltage on Telecommunication Lines in Parallel Distribution Lines. Trans. KIEE 2008, 57, 1688–1695. [Google Scholar]
  22. National Radio Research Agency. The Technical Regulation on Detailed Calculation Methods of Power Induction Voltages; RRA Notice No. 2014-11; National Radio Research Agency: Naju, Korea, 2014. [Google Scholar]
  23. Kersting, W.H. The Modeling and Analysis of Parallel Distribution Lines. IEEE Trans. Ind. Appl. 2006, 42, 1126–1132. [Google Scholar] [CrossRef]
  24. Electric Power Research Institute. EPRI AC Transmission Line Reference Book—200 kV and Above; Electric Power Research Institute: Washington, DC, USA, 2014. [Google Scholar]
  25. Papadopoulos, T.A.; Apostolidis, A.K.; Chrysochos, A.I.; Christoforidis, G.C. Frequency-Dependent Earth Impedance Formulas Between Overhead Conductors and Underground Pipelines. In Proceedings of the IEEE International Conference on Environment and Electrical Engineering and 2019 IEEE Industrial and Commercial Power Systems Europe (EEEIC/I&CPS Europe), Genova, Italy, 11–14 June 2019; pp. 11–14. [Google Scholar]
  26. Amtani, A.; Hosakawa, Y. EMTP Simulations and Theoretical Formulation of Induced Voltages to Pipelines from Power Lines; IPST: Lyon, France, 2007. [Google Scholar]
Figure 1. Induced voltage calculation of SCL.
Figure 1. Induced voltage calculation of SCL.
Energies 14 08536 g001
Figure 2. Induced voltage calculation of DCL.
Figure 2. Induced voltage calculation of DCL.
Energies 14 08536 g002
Figure 3. Configurations of distribution lines and pipelines. (a) SCL. (b) DCL.
Figure 3. Configurations of distribution lines and pipelines. (a) SCL. (b) DCL.
Energies 14 08536 g003
Figure 4. Configuration of the pipelines; rg1 = 191.3 (mm), rg2 = 203.2 (mm), rg3 = 206.4 (mm); ρg = 1.59 × 10−7 (Ω m), μg = 280, εg = 2.30.
Figure 4. Configuration of the pipelines; rg1 = 191.3 (mm), rg2 = 203.2 (mm), rg3 = 206.4 (mm); ρg = 1.59 × 10−7 (Ω m), μg = 280, εg = 2.30.
Energies 14 08536 g004
Figure 5. Vector diagram of SCL (Case 2A, Case 2C), where VP is multiplied by five for clarity, as shown in the legend.
Figure 5. Vector diagram of SCL (Case 2A, Case 2C), where VP is multiplied by five for clarity, as shown in the legend.
Energies 14 08536 g005
Figure 6. Vector diagram of DCL (Case 7B, Case 8B), where VP is multiplied by five for clarity and illustrative purposes.
Figure 6. Vector diagram of DCL (Case 7B, Case 8B), where VP is multiplied by five for clarity and illustrative purposes.
Energies 14 08536 g006
Figure 7. Comparison of SCL and DCL.
Figure 7. Comparison of SCL and DCL.
Energies 14 08536 g007
Figure 8. Comparison of single-phase unbalanced load in DCL.
Figure 8. Comparison of single-phase unbalanced load in DCL.
Energies 14 08536 g008
Figure 9. Induced voltage according to the separation distance in SCL.
Figure 9. Induced voltage according to the separation distance in SCL.
Energies 14 08536 g009
Figure 10. Induced voltage according to the separation distance in DCL.
Figure 10. Induced voltage according to the separation distance in DCL.
Energies 14 08536 g010
Table 1. Conductor type and data.
Table 1. Conductor type and data.
Line TypeCable Type
(mm2)
Radius
(cm)
Resistance
(Ω/km)
Overhead ground conductorACSR 320.390.898
Distribution linesACSR 1600.910.182
Neutral conductorACSR 950.6750.301
Table 2. Load condition results of SCL.
Table 2. Load condition results of SCL.
Case StudyLoad ConditionVP
(V)
VEMTP
(V)
VDiff.
(%)
A
(MVA)
B
(MVA)
C
(MVA)
Case 11110.5130.5150.2
Case 2A1.3110.7130.7160.3
Case 2B11.310.7000.6980.2
Case 2C111.30.3530.3560.3
Case 311.11.20.4650.4670.2
Table 3. SCL calculation of induced voltage using (5).
Table 3. SCL calculation of induced voltage using (5).
Case StudyVG (V)VABC (V)VN (V)
Case 10.228 ∠ −21.60°0.389 ∠ 48.60°8.6 × 10−15 ∠ 163.42°
Case 2A0.199 ∠ 14.32°6.236 ∠ −127.33°6.783 ∠ 52.70°
Case 2B0.430 ∠ −21.67°6.894 ∠ 109.67°6.783 ∠ −67.30°
Case 2C0.198 ∠ −57.54°6.976 ∠ −4.71°6.783 ∠ 172.7°
Case 30.262 ∠ −38.85°4.255 ∠ 24.80°3.919 ∠ −157.27°
Table 4. Inducing and shielding side analysis of SCL.
Table 4. Inducing and shielding side analysis of SCL.
Case StudyInducing SideShielding SideInduced Voltage
VABC (V)VG + VN (V)VP (V)
Case 10.389 ∠ 48.60°0.228 ∠ −21.60°0.513 ∠ 23.83°
Case 2A6.236 ∠ −127.33°6.940 ∠ 51.68°0.713 ∠ 42.99°
Case 2B6.894 ∠ 109.67°7.090 ∠ −64.82°0.700 ∠ 6.17°
Case 2C6.976 ∠ −4.7°6.658 ∠ 174.01°0.353 ∠ 20.23°
Case 34.255 ∠ 24.80°3.801 ∠ −153.79°0.465 ∠ 13.21°
Table 5. Load condition results of DCL.
Table 5. Load condition results of DCL.
Case StudyLoad ConditionVP
(V)
VEMTP
(V)
VDiff.
(%)
UA
(MVA)
UB
(MVA)
UC
(MVA)
LA
(MVA)
LB
(MVA)
LC
(MVA)
Case 41111110.3300.3290.1
Case 5111.3111.30.4010.3881.3
Case 611.11.211.11.20.4060.3961.0
Case 7A1.3111110.4300.4260.4
Case 7B11.311110.8190.8081.1
Case 7C111.31110.2840.2750.9
Case 8A1111.3110.4120.4100.2
Case 8B11111.310.5250.5160.9
Case 8C111111.30.2530.2520.1
Table 6. DCL calculation of induced voltage using (10).
Table 6. DCL calculation of induced voltage using (10).
Case StudyVG (V)VU (V)VL (V)VN (V)
Case 40.322 ∠ −21.9°0.354 ∠ −132.7°0.391 ∠ 48.4°8.6 × 10−15 ∠ 172.4°
Case 50.343 ∠ −88.0°6.421 ∠ −10.6°7.041 ∠ −5°13.697 ∠ 172.4°
Case 60.411 ∠ −52.6°3.500 ∠ 19.8°4.292 ∠ 24.5°7.909 ∠ −157.5°
Case 7A0.294 ∠ 26.2°7.085 ∠ −128.0°0.391 ∠ 48.4°6.849 ∠ 52.4°
Case 7B0.642 ∠ −22.0°6.507 ∠ 114.9°0.391 ∠ 48.4°6.849 ∠ −67.6°
Case 7C0.292 ∠ −70.0°6.421 ∠ −10.6°0.391 ∠ 48.4°6.849 ∠ 172.4°
Case 8A0.284 ∠ −2.0°0.354 ∠ −132.7°6.298 ∠ −127.6°6.849 ∠ 52.4°
Case 8B0.462 ∠ −21.9°0.354 ∠ −132.7°6.961 ∠ 109.4°6.849 ∠ −67.6°
Case 8C0.283 ∠ −41.7°0.354 ∠ −132.7°7.041 ∠ −5.0°6.849 ∠ 172.4°
Table 7. Inducing and shielding side analysis in DCL.
Table 7. Inducing and shielding side analysis in DCL.
Case StudyInducing Side
VU + VL (V)
Shielding Side
VG + VN (V)
Induced Voltage
VP (V)
Case 40.037 ∠ 58.3°0.322 ∠ −21.9°0.330 ∠ −15.5°
Case 513.446 ∠ −7.7°13.644 ∠ 173.8°0.401 ∠ −126.1°
Case 67.785 ∠ 22.4°7.813 ∠ −154.6°0.406 ∠ −70.0°
Case 7A6.695 ∠ −127.8°7.114 ∠ 51.4°0.430 ∠ 38.5°
Case 7B6.672 ∠ 111.8°7.312 ∠ −64.0°0.819 ∠ −27.5°
Case 7C6.631 ∠ −7.7°6.719 ∠ 174.6°0.284 ∠ −114.6°
Case 8A6.651 ∠ −127.9°7.017 ∠ 50.5°0.412 ∠ 24.2°
Case 8B6.802 ∠ 112.1°7.179 ∠ −64.9°0.525 ∠ −22.3°
Case 8C6.831 ∠ −7.3°6.616 ∠ 173.8°0.253 ∠ −38.9°
Table 8. Comparison of Case 2C and Case 5.
Table 8. Comparison of Case 2C and Case 5.
Case StudyInducing SideShielding Side
VG + VN (V)
Induced Voltage
VP (V)
Case 2C6.976 ∠ −4.71°0.198 ∠ −57.54° + 6.783 ∠ 172.7°0.353 ∠ 20.23°
Case 513.446 ∠ −7.7°0.343 ∠ −88.0° + 13.697 ∠ 172.4°0.401 ∠ −126.1°
Table 9. Induced voltage according to the separation distance in SCL.
Table 9. Induced voltage according to the separation distance in SCL.
Case
Study
Separation Distance (VEMTP (V))
−90
(m)
−60
(m)
−30
(m)
0
(m)
+30
(m)
+60
(m)
+90
(m)
Case 10.1110.1410.2161.5400.2870.1870.142
Case 2A0.0990.1360.2360.3320.5330.3360.250
Case 2B0.2710.3260.4480.5650.5520.3870.309
Case 2C0.2070.2750.4340.1770.3410.2130.156
Case 30.2090.2680.4070.2430.4110.2710.207
Table 10. Inducing and shielding side analysis in Case 2A of Table 9.
Table 10. Inducing and shielding side analysis in Case 2A of Table 9.
Separation DistanceInducing Side Voltage
(V)
Shielding Side Voltage
(V)
Induced Voltage
(V)
−90 (m)4.164 ∠ −135.51°4.118 ∠ 43.26°0.099 ∠ −73.53°
−60 (m)4.862 ∠ −132.86°4.768 ∠ 45.99°0.136 ∠ −87.79°
−30 (m)6.037 ∠ −129.93°5.827 ∠ 49.09°0.236 ∠ −104.45°
0 (m)7.086 ∠ −127.30°7.411 ∠ 52.05°0.332 ∠ 38.28°
+30 (m)5.301 ∠ −129.94°5.827 ∠ 49.09°0.533 ∠ 39.34°
+60 (m)4.438 ∠ −133.13°4.768 ∠ 45.99°0.336 ∠ 34.56°
+90 (m)3.874 ∠ −135.94°4.118 ∠ 43.26°0.250 ∠ 30.81°
Table 11. Induced voltage according to the separation distance in DCL.
Table 11. Induced voltage according to the separation distance in DCL.
Case StudySeparation Distance (VEMTP (V))
−90
(m)
−60
(m)
−30
(m)
0
(m)
+30
(m)
+60
(m)
+90
(m)
Case 40.1970.2280.2750.3110.2800.2280.197
Case 50.2020.2370.3030.5420.2890.2350.202
Case 60.2470.2860.3510.4390.3470.2860.247
Case 7A0.1990.2390.3230.5500.2550.1970.170
Case 7B0.3830.4470.5620.8640.6050.4700.399
Case 7C0.1790.2090.2540.3230.2340.1970.171
Case 8A0.1710.1960.2340.3810.2900.2190.184
Case 8B0.2890.3390.4310.5770.3960.3170.273
Case 8C0.1580.1780.2010.1810.2650.2160.184
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kim, H.-S.; Min, H.-Y.; Chase, J.G.; Kim, C.-H. Analysis of Induced Voltage on Pipeline Located Close to Parallel Distribution System. Energies 2021, 14, 8536. https://doi.org/10.3390/en14248536

AMA Style

Kim H-S, Min H-Y, Chase JG, Kim C-H. Analysis of Induced Voltage on Pipeline Located Close to Parallel Distribution System. Energies. 2021; 14(24):8536. https://doi.org/10.3390/en14248536

Chicago/Turabian Style

Kim, Hyoun-Su, Hae-Yeol Min, J. Geoffrey Chase, and Chul-Hwan Kim. 2021. "Analysis of Induced Voltage on Pipeline Located Close to Parallel Distribution System" Energies 14, no. 24: 8536. https://doi.org/10.3390/en14248536

APA Style

Kim, H. -S., Min, H. -Y., Chase, J. G., & Kim, C. -H. (2021). Analysis of Induced Voltage on Pipeline Located Close to Parallel Distribution System. Energies, 14(24), 8536. https://doi.org/10.3390/en14248536

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop