Elevated Neutral to Earth Voltages Due to Harmonics A T&D Update E. R. (Randy) Collins, PhD, PE Dept. of Electrical and Computer Engineering Clemson University Clemson, South Carolina Stray Voltage Panel Session 8 IEEE T&D Meeting - Chicago April, 8 1 Two Publications of Interest Analysis of Elevated Neutral-to-Earth Voltage in Distribution Systems with Harmonic Distortion to appear in 8, IEEE Trans. on Power Delivery, Paper TPWRD-65-6 by Randy Collins and Jian Jiang. Elevated Neutral-To-Earth Voltage In Distribution Systems With Harmonic Distortion, PhD Dissertation, Clemson University, by Jian Jiang, 7. 1
Three Phase Distribution Feeder with Single Phase Lateral 3 Earth Equipotentials surround ground rods
What do we mean by NEV and how does this relate to stray voltage? 5 Neutral-to-earth voltage is the neutral conductor voltage referred to remote earth. The NEV can be much higher than the stray voltage that can be contacted by humans or animals, so the results must be tempered with this understanding. Knowledge of NEV under different system conditions is extremely important in analyzing and mitigating stray voltage problems arising from NEV. 6 Common Sources and Conventional Resolutions for Elevated NEV and Stray Voltage Sources Power system grounding Load unbalance Transformer connections Neutral conductor impedance Alleviation methods Balance the loads Three-phase single phase laterals Increase the neutral conductor size Improve grounding connection Repair bad neutral connections and splices 3
Trends of Consumer Electronics and Their Impact on NEV 7 Proliferation of single phase nonlinear power supplies in electronic appliances Non-linear response of single-phase motors to voltage distortion Increased triplen harmonic distortion Elevated neutral conductor current 8 Our Research Objectives Develop means for systematic analysis of NEV in distribution power systems including harmonic distortion Develop procedures to predict NEV for planning and assessment, and to troubleshoot existing problems Assess the effects of NEV mitigation techniques in an environment with harmonic distortion
9 Steps of This Project Develop detailed power system models, including the nonlinear loads, for NEV analysis. Construct a multiphase harmonic load flow algorithm to calculate NEV. Test the algorithm accuracy by comparing with the actual field measurement. Examine the algorithm capabilities using a standard IEEE example system. Develop procedures for locating sources of elevated NEV and assess common methods in reducing NEV with nonlinear loads. Single-Phase Distribution Feeder Modeling 1 Phase Conductor Neutral Parallel Ground Current 5
6 11 Conventional Transmission Line Model (Carson s Line) ( ) 1 1.6935 / 1 8.65... _ ln 1.386... cos 3 1 8,, ln ln = = = = = = = G f S k k Q k P n a j i G Q D S j G P Z G ii Q Dii ii S j G ii P i r ii Z ρ θ π ω ω ω ω ω ω S Distance between conductors and the images D Distance between conductors = n a g nn g na g an g aa nn aa I I Z Z Z Z V V ' ' 1 Practical Distribution Feeder with Nonzero Grounding Resistance
13 Equivalent PI Model NEV 1 Z se [R] -1 [R] -1 NEV The series impedance is the same as in Carson s line Z aa Z an [ Z se ] = Z an Z nn The shunt admittance matrix is obtained from the grounding resistance 1 R = 1 R g [ ] 1 Multiple Segments Assembling I a a 1 Z a1 a Z a a 3 I a Z an1 Z na1 Z an Z na n 1 Z n1 n Z n n 3 I n Z ag1 Z ag I n R g Z ng1 R g Z ng R g I g I g g 1 g [ Y 1 ] [ Y ] g 3 7
15 Equivalent PI Model Apply Kron reduction to Ybus to eliminate the nodes in the middle of the feeder Derive the equivalent pi model from the resultant matrix Y new bus Z se = [ A] [ B] [ B] [ A] [ ] 1 = B = [ A B] Y = [ A B] Y sh sh 16 Step of This Project Linear modeling Nonlinear modeling Multiphase harmonic load flow Examine algorithm capability using IEEE example system Test algorithm accuracy by comparing with actual field measurements Assess three-phasing method with harmonic distortion present 8
Typical Non-Linear Loads in Distribution Systems 17 Three phase AC/DC converter Single phase load with uncontrolled rectifier front end Single phase induction motor Fluorescent lighting 8.9 1. Amps -1.. 8.3 1.5 16.6-8.9 18 Single Phase Rectifier ω di dθ s V th = RT is LT dv dθ o i s = ωc V R o eq V o d Vo dθ d is dθ dv dθ o ( α1 α ) ( α1α α α 3) Vo = α α 3Vth di dθ s ( α α ) ( α α α α ) i 1 1 3 s = α α V th dvth α dθ 9
Typical Waveforms of a Single Phase Rectifier 19 6 V o I s Vol t s/ Amps - - -V th V th θ1 θ - 6 5 1 15 5 3 35 Degr ee Step of This Project Linear modeling Nonlinear modeling Multiphase harmonic load flow Test algorithm accuracy by comparing with actual field measurement Examine algorithm capability using IEEE example system Assess three phasing method with harmonic distortion 1
Multiphase Harmonic Load Flow for NEV Analysis 1 The nature of NEV problems demands a multiphase load flow solution NEV problems are often more prevalent in radial distribution systems Simple load flow program with robust convergence is desired A Typical Radial Distribution Topology Root Node 1 3 5 Layer 1 6 7 8 9 1 11 Layer 1 13 1 15 16 17 18 19 Layer 3 1 3 5 6 7 8 Layer 9 3 31 3 33 3 35 Layer 5 11
Multiphase Load Flow Algorithm - Iterative Solution 3 Calculate the series impedance matrix for each branch and sum up the shunt admittance at every node except for the root node Assume the node voltages for initialization Find the nodal injection current using the node voltage in the previous iteration and determine the division between the neutral conductor and earth currents Back sweep the system to find the branch current from the last layer by applying KCL Forward sweep the system to find the voltage drop for each branch from the first layer and update the nodal voltage Check convergence and start a new iteration if necessary Steps of This Project Linear modeling Nonlinear modeling Multiphase harmonic load flow Test algorithm accuracy by comparing with actual field measurement Examine algorithm capability using IEEE example system Assess three phasing method with harmonic distortion 1
Algorithm Verification Using Actual System Measurement 5 Δ /YG Three phase feeder with a single phase lateral (all conductor sizes are 336 kcmil) Transformer impedance: 7% Assume uniformly distributed loads after the measuring point and ignore any load before it 6 Actual Measurement Phase Neutral Meter ratios: Voltage (5 V/div), Current (1 A/div) 13
7 Synchronized Measurement Waveforms 15 Voltage (kv) and Current (A) 1 5-5 Phase voltage Phase current Neutral plus messenger current -1-15 5 1 15 5 3 35 Degree 8 Synchronized Measurement Waveform 15 Voltage (kv) and Current (A) 1 5-5 Phase voltage Phase current Neutral plus messenger current -1-15 5 1 15 5 3 35 Degree 1
9 NEV Simulation Δ /YG The measured current is injected half way from the measuring point due to the uniform load distribution The calculated neutral current in the single phase lateral is compared with the measurement to check the modeling and algorithm accuracy Waveform Comparison Using Simulation Results 3 15 Voltage (kv) and Current (A) 1 5-5 Phase voltage Phase current Calculated neutral current Measured neutral current -1-15 5 1 15 5 3 35 Degree Close match-up between measurement and calculation 15
Comparison of Measurement and Simulation Results 31 15 Voltage (kv) and Current (A) 1 5-5 Phase voltage Phase current Calculated neutral current Measured neutral current -1-15 5 1 15 5 3 35 Degree Close match between measurement and calculation 3 Step of This Project Linear modeling Nonlinear modeling Multiphase harmonic load flow Test algorithm accuracy by comparing with actual field measurement Examine algorithm capability using IEEE example system Assess three phasing method with harmonic distortion 16
Test on IEEE Example Distribution System 33 Load flow algorithms test system Relative short feeder High loading level Various line configurations Different load connections Black Three phase Green Two phase Red One phase 3 Linear Loads on The System Load Type Bus No. Phase 1 kw kvar Phase kw kvar Phase 3 kw kvar Y 8 16 11 1 9 1 9 Y 17 15 Δ 5 3 13 Δ 3 385 385 385 Δ 3 17 151 Y 9 17 8 17
NEV Test with Linear Loads at Different Earth Resistivities 35 Neutral to Earth Voltage (rms Volts) 1 1 1 8 6 1 3 5 6 7 8 9 1 Ω-m 1 Ω-m 1 Ω-m Bus Number The Effect of Non-Linear Loads on NEV 36 Root Single-phase rectifier: V, 3kW with a μf smoothing capacitor 3 5 6 7 8 9 1 Nonlinear load bus Layer 1 Layer Layer 3 Layer Case 1: A single rectifier on a single-phase line at bus 3 Case : Three balanced rectifiers on a three-phase line at bus 3 Case 3: Two identical rectifiers on two phases of a three-phase line at bus 3 to form an unbalanced load 18
NEV RMS Values with Different Nonlinear Load Combinations 37 Neutral to Earth Voltage (rms volts) 35 3 5 15 1 5 1 3 5 6 7 8 9 Linear load only Case 1 Case Case 3 1 Ω-m ground resistance Bus Number Case : Three identical single-phase rectifiers at bus 3 38 Step of This Project Linear modeling Nonlinear modeling Multiphase harmonic load flow Test algorithm accuracy by comparing with actual field measurement Examine algorithm capability using IEEE example system Assess three phasing method with harmonic distortion 19
39 Three Phasing Single Phase Lateral Δ /YG A common practice for distribution system upgrade Simulation setup Three phasing the single phase lateral Evenly distribute the measured current among the three phases Inject the currents half way from the measuring point Detailed Results of Harmonic Components Harmonic Order NEV (V) Single Phase Current (A) I ph I n I E NEV (V) Three Phase Current (A) I ph I n I E 1 1.7 9.15 6.66.7.18 3.5.. 3.93..5.18.93.1.5.18 5.7.3.1.1..1.. 7.6.1.7.7.1.5.. 9.8.8.3..8.3.3. 11.1.3....1.. 13.19.....1.. 15.17....17.1.. TOT RMS 1.8 1.5
Assessing The Effect of Three Phasing with Nonlinear Loads 1 Δ /YG Same system configuration The linear loads are modeled in constant power mode with its value doubled: 1kW @.9 pf lagging Nonlinear load: single phase rectifier at V, 6kW with 6 μf smoothing capacitors Single Phase Lateral NEV Spectrum with nonlinear load Voltage (V) 6 5 3 1 At each bus, the NEV spectrum include the following components: Total RMS Fund. 3 rd 5 th 7 th 9 th 1 3 Bus number 1
NEV Spectrum with Nonlinear load After Three Phasing 3 Voltage (V) 35 3 5 15 1 5 At each bus, the NEV spectrum include the following component: Total RMS Fund. 3 rd 5 th 7 th 9 th 1 3 Bus number Detailed Results for Harmonic Components with Nonlinear Load Harmonic Order NEV (V) Single Phase Current (A) I ph I n I E NEV (V) Three Phase Current (A) I ph I n I E 1 38.1 7.. 8.9.55 9.69.1.1 3 3.56 1.5 8.86 6. 36.3 5. 9.86 7.19 5 16.68 5.8 3.7 3..5..6. 7 5.58 1.67.8.91.13.6.3. 9 5.9 1.58.7.87 6.35.56.78.93 11.93.71.3.39.5.31.1.1 13 3.18.7.3..7.5.1.1 15 1.83.39.18.1..16..6 RMS TOT 57. 37
5 Conclusions Carson s line model is adapted for NEV analysis A multiphase harmonic load flow algorithm is developed for NEV analysis in distribution systems The simple current injection method used in harmonic analysis achieves similar accuracy as the detailed nonlinear device model It is possible to approximately predict the NEV profile and help locate potential stray voltage problems Triplen harmonics elevates the NEV due to their additive nature in the neutral Common methods for stray voltage mitigation, e.g. three phasing, may not be adequate when there is high harmonic distortion in the system 6 Future Research Opportunities More field measurements to validate model. Improve models to incorporate more complex situations Improve load flow program to handle large scale systems Better understand the human sensitivity to harmonic frequency voltages 3