390 likes | 579 Views
Title:. A New Hysteretic Reactor Model for Transformer Energization Applications. By : Afshin Rezaei-Zare & Reza Iravani University of Toronto. June 2011. Outline. Existing hysteresis models in EMT programs Drawbacks of the existing models New hysteretic reactor model
E N D
Title: A New Hysteretic Reactor Model for Transformer Energization Applications By: Afshin Rezaei-Zare & Reza Iravani University of Toronto June 2011
Outline • Existing hysteresis models in EMT programs • Drawbacks of the existing models • New hysteretic reactor model • Impact on Remnant Flux (Lab. Measurement) • De-energization / Re-energization • 33kV-VT Ferroresonance Lab. test results • Conclusions • Applications
Existing Hysteresis Models in EMT Programs • EMTP Type-96 • EMTP Type-92 (Current Hysteresis model of the EMTP-RV) • PSCAD/EMTDC Jiles-Atherton model(Not a reactor but incorporated in the CT model) • Proposed New Hysteretic Reactor
Type 96 model • Piecewise linear model • Originally developed by Talukdar and Bailey in 1976 and modified in 1982 by Frame and Mohan • Simple and computationally efficient • Minor loops are obtained by linearly decreasing the distance between the reversal point and the penultimate reversal point
Drawbacks of Type 96 Model • No stack is used to store the extrema of excitation which leads to open cycles • Similarity of minor loops to the major loop due to scaling approach used by the model (such a similarity is not valid in reality) • The existence of a saturation point is not verified experimentally • The model implemented in EMTP-V3 is pseudo nonlinear
Drawbacks of Type 96 Model – Cont’d • Noisy behavior and erroneous results (in some transients such as ferroresonance) due to switching the operating point between two adjacent branches of the piece-wise linear characteristic (Artificial switching & Numerical oscillations) Piece-wise linear model Smooth nonlinear model
Type 92 model • Developed in 1996 by Ontario Hydro • Based on hyperbolic functions • Instantaneous flux is separated in two components: i) hysteresis (irreversible) ii) saturation (reversible) • Current EMTP-RV model is based on this approach
Model Type 92 Hyperbolic functions in Type 92: instantaneous flux is used to find unsaturated flux which is then used to find instantaneous current Unsaturated flux vs. Current (to describe hysteresis) Saturated flux vs. Unsaturated flux (to describe saturation) slope slope
Drawbacks of Type 92 Model Inaccuracy (1) • Limited flexibility to fit to the hysteresis major loop (only based on one hyperbolic term) Raw data Fitted data
Drawbacks of Type 92 Model Inaccuracy (2) • Only upper part of the trajectory is used and the lower part is assumed to be symmetric to the upper part (while in reality the shapes of the two parts are independent) Current (A)
Jiles-Atherton Model • Physically correct model • decomposes the magnetization into “reversible anhysteretic” and “irreversible” components based on a weighted average: • Reversible part is based on Langevin function: • Irreversible part is based on the differential equation:
Drawbacks of the Jiles-Atherton Model • Limited flexibility to fit to the measurements due to the utilized Langevin function, and the model very few parameters (5 parameters) • In some cases, non-physical results as the input current changes the direction • Formations of minor loops and the major loop are dependent (changing the parameters changes both minor and major loop shapes) • In the PSCAD/EMTDC, it is not available as a reactor to build a desired general system for transient studies.(Only incorporated in a CT model)
New Hysteretic Reactor Model • A modified Preisach Model- a time-domain implementation with true-nonlinear solution within the EMTP-RV - Independent formation of minor loops from the major loop (consistent with the observed hysteresis loops of the magnetic materials) • Physically correct hysteresis model • Memory dependent model: past excitation extrema are stored in memory to form the magnetization trajectories. • Representing wiping-out property, (a well-known physical property of the ferromagnetic materials)
New Hysteretic Reactor Model Forms major loop Forms minor loop Same major loops – Different minor loops
Remnant Flux 80% 40% -50% 0%
Impact on Remnant Flux (Lab. Measurement) im Close-up Window
De-energization / Re-energization Auto-reclosure operations on a 12kA Fault Current
Remnant Flux subsequent to the second current interruption Remnant flux Different Minor loop shapes Different Remnant Flux (for the same switching scenario)
Impacts on CT Saturation and protection (Following the final reclosure on the permanent fault)
33kV-VT Ferroresonance Laboratory test results 63.5 kV (2.36pu) 52.4 kV (1.94pu) 38.4 kV (1.43pu) 30.7 kV (1.14pu) 11.8 kV (0.44pu) Source peak voltage Measured VT voltage
33kV-VT Ferroresonance Lab Test 63.5 kV (2.36pu) Voltage 52.4 kV (1.94pu) 218 W Power Loss 103 W 29 W
33kV-VT Ferroresonance Lab Test Model Type-92 Hysteresis Loop
33kV-VT Ferroresonance Lab Test New Reactor Hysteresis Loops Hysteresis loop at rated voltage Major loop
33kV-VT Ferroresonance Lab Test Bifurcation Points
33kV-VT Ferroresonance Lab TestCore Power Loss Power Loss Comparison
33kV-VT Ferroresonance Lab Test – Cont’d Hysteresis Loops Comparison New Reactor Measurement EMTP Type-96 EMTP-RV (Type 92)
33kV-VT Ferroresonance Lab Test Dynamic Inductance ( Slope of magnetization trajectories ) Before Ferroresonance (Normal conditions) Under Ferroresonance conditions
Capability of the models to represent the core dynamic behaviors Core Inductance change As the core is driven into ferroresonance with respect to normal operation Model Change direction Measurement New Reactor EMTP Type-96 EMTP-RV (Type-92) No change Single-valued saturation curve No change
Another Example – Comparison between two hysteresis models with the same major loop but different minor loop formations Bifurcation diagrams Model 1 Model 2
Conclusions New Model Features • The model is based on widely-verified and accepted Preisach model of hysteresis • Independent formation of major and minor loops • True nonlinear solution within the EMTP-RV • Can accurately represent the physical properties of the magnetic core materials • Can accurately represent the dynamic core behavior under electromagnetic transients
Applications For accurate EMTP studies on : • De-energizing/re-energizing of transformers • Ferroresonance phenomena in power and instrument transformers • Determination of the core remnant flux • Precise modeling of VTs, CTs, and CVTs for protection studies • Accurate modeling of electrical machines • Efficient design of control systems for power-electronic based drives by taking into account the machine nonlinearity and actual power loss
Important points • It is evident that a more detailed model needs more parameters. although, a model with simple implementation and with less required parameters is generally preferable, the accuracy of such models are limited. • For a sophisticated hysteresis model, “not needing minor loop data”, is a drawback not an advantage. Due to different behavior of minor loops (extensively verified by experiments), neglecting the minor loop parameters can result in completely different and unexpected results. • For the new reactor, if the minor loop data are not available, a set of pre-specified default values can be considered.