It is justified that during geomagnetic storms the high voltage synchronous engines are being impacted by high current harmonics of even sequences powered by power transformer due to geo-induced high voltage currents flowed through the windings. Equivalent circuits of step down substation and HV synchronous motors are made for making it possible to consider a saturation of power transformer magnetic system and higher current harmonics availability in stator windings. Analytic expressions for higher current harmonics and extra capacity losses calculation in stator windings are received, as well as the calculation of induction torques allowing to denote a rate of geomagnetic processes impact on synchronous engine operation at various step down substation parameters.
During geomagnetic storm periods (GMS), the geomagnetically induced currents (GIC) are flowing into highvoltage power grids and high-voltage power transformer windings with dead-earthed neutral, which in turns create asymmetric saturation of magnetic system and multiple increase of nonsinusoidal magnetizing current [1,2]. As a result, power transformer becomes a powerful source of high order current harmonics, which impact on all power consumers of power distribution circuit [3,4]. Electrical motor load is the most sensitive to supply voltage quality degradation, especially high-voltage synchronous motors (SM) with direct connection to power distribution circuit 6 (10) kV. High order currents in stator SM windings result into additional power loss and additional electromagnetic torque constituents appearing.
The most dangerous thing for SM is additional power losses that cause insulation temperature increase, particularly in slot-embedded part of stator windings in the period of GMS. This is a factor of gradual accumulating of defects and accelerated slot-embedded insulation wear, damaging of which could cause up to 80% cases of highvoltage SM stator insulation failures. Nevertheless, the questions of GIC’s impact on work conditions of SM high-voltage are not reflected in science literature. In this relation it is up-to-date and practically significant task to evaluate a level of high order current in stator SM windings during GMS period.
Assigned task is solved based on typical step down substation of industrial enterprise, main consumers of which are SMs of turbine-driven air compressor connected to sections I c - II с of switch-yard bus SY-6kV (
Equivalent circuit has the following designations accepted:
—inductive and active resistances, and also capacitive conductance of phase wire PL-110 on n-harmonic frequency;
—inductive and active resistances of phase winding of high-voltage (HV) T1 power transformer on n-harmonic frequency;
—inductive and active resistances on n-harmonic frequency of T1 power transformer LV1 phase winding and cable line CL-6 feeding SM1 from section I c ;
—inductive and active resistances on n-harmonic frequency of T1 power transformer LV2 phase winding and cable line CL-6 feeding SM1 from section II c ;
—capacitive conductances on n-harmonic frequency of cable lines CL-6 feeding SM1 and SM2 correspondingly;
—n-harmonic of Т1 power transformer’s magnetizing current.
Magnetization current’s harmonics parameters for equivalent circuit
As it is shown on
Due to magnetic and rotor electric asymmetry, conditioned by excitating winding, SM equivalent circuits on longitudinal and transversal axes are different as it is shown on
High order harmonic currents are starting to flow when saturating of T1 power transformer under GIC impact in SM stator windings except basic current constituent with 50 Hz frequency. If saturation isn’t taken into account, SM is allowed to be considered as line system for which superposition principle being applicable. In this case analysis is hold independently for basic constituent and for each high order harmonic. Among other things, the consumed current and SM electromagnetic torque are defined by amount of constituent currents and torques from each harmonic. SM equivalent circuit parameters:
—active and inductive resistances of stator winding dissipation on n-harmonic frequency;
—mutual induction resistance between stator and rotor windings on longitudinal and transversal axes on n-harmonic frequency;
—active and inductive resistances of excitating winding dissipation on n-harmonic frequency;
—inductive dissipation resistances of longitudinal and transversal damping contours on n-harmonic frequency;
—active resistances of longitudinal and transversal damping contours on n-harmonic frequency.
The method stated in [
1) An increase of all inductive resistances on n times comparing to values on basic harmonic frequency;
2) n-harmonic’s slipping (s(n)) of magnetomotive force
(mmf) of direct and reverse sequence created by n-harmonic of stator winding phase currents, value of which is defined when rotor synchronous motion saving by following expression
where s—rotor slipping relatively to rotating magnet field of basic harmonic (during synchronous rotor motion at s = 0). Upper sign corresponds to stator winding current harmonics of direct sequence; lower sign corresponds to stator winding current harmonics of reverse sequence: an increase of active stator and rotor resistances due to skin-effect on higher frequencies. Dependence of active SM contours resistances on frequency is being approximated by empirical expression [
where—contour resistance on basic harmonic frequency.
In order to evaluate the individual harmonics impact rate, it is advisable to determine an average equivalent and complex SM resistance on two axes [
where—complex conductances of mutual induction between stator winding and rotor contours on longitudinal and transversal axes on n-harmonic frequency;—complex excitating winding conductance on n-harmonic frequency;—complex longitudinal and transversal damping contours conductances on n-harmonic frequency correspondingly.
Expression (3) can be folded and written in a from that is more convenient and compact for its further use
where
—are average equivalent active and reactive SM resistances on two axes;—summary active resistances of parallel branches of SM equivalent circuit on n-harmonic frequency on longitudinal and transversal axes correspondingly;—summary reactive resistances of SM equivalent circuit’s parallel branches on longitudinal and transversal axes correspondingly.
On the basis of (4) the SM equivalent circuit for high order harmonics (see
In order to define high order harmonic currents that flow into stator SM windings during GIC power transformer magnetic biasing, it is required to complement complex resistance of SM equivalent circuit to step down substation equivalent circuit showed on
It is advisable to have current source preliminarily being transformed into equivalent EMF source, the internal complex resistance of which is defined by the following expressions
As follows from the integration, the estimated equivalent circuit (see
From the condition
resonance frequency of internal complex resistance wp is defined
where w(1)—is angular frequency of basic harmonic;—wave resistance PL-110 kV.
It should be noted that n = 20 harmonic is resonance for stepdown substation showed on
Using equivalent circuit on
where KHV-LV—complex division ratio of magnetizing current’s n-harmonic between HV and LV sides; KLV1-LV2 —complex division ratio of magnetizing current’s n-harmonic between split LV windings; KD—complex division ratio of magnetizing current’s n-harmonic between SM and longitudinal capacitive conductance of power cable.
Expression (8) graphically shows which part of n-harmonic of T1 power transformer magnetizing current flows in SM stator winding and how it depends on step down substation parameters. Values of complex division ratios are defined by following expressions
where —complex SM resistances on n-harmonic frequency taking power cable parameters into account.
As per (8) and (9) it is possible to define high order harmonics in SM stator windings connected to bus sections I с (ISM(n)I) and II с (ISM(n)II), during GIC impacting on T1 power transformer from power grid side. Quantitative estimations made for synchronous turbomotors of STM series with 10 MW capacity allow to admit that majority level of high order harmonics in stator windings do not exceed (3 ÷ 7)% from same-name level harmonics of T1 power transformer magnetizing current. It’s only for 19th, 20th and 22nd harmonics that relative level is significantly higher and reaches the values of 19%, 48% and 13.2% correspondingly. Although in same-name units the indicated harmonics are not exceeding 2nd current harmonic level in SM stator windings.
The most serious effects for SM could arise providing that bC1 » 0 and high order harmonics of magnetizing current will “flow down” on T1 power transformer LV side, and SM power is carried out only from one bus section.
Finally magnetizing current’s high order harmonics of T1 power transformer will practically be flowing through SM stator windings, i.e. ISM(n) @ I(n) In this case additional capacity losses that occurred in SM due to high order current’s harmonics, at values set in
At ISM(n) @ Iμ(n) a necessity is also arisen to take into account magnetizing current of asynchronous electromagnetic torque constituents created by high order harmonics.
Permanent asynchronous moment on SM shaft, conditioned by interaction of same magnetic flux harmonics in air gap and rotor current, is made up of two constituents [
where—is asynchronous torque, conditioned by excitating winding reaction and longitudinal damping contour;—asynchronous torque conditioned by longitudinal damping contour.
In order to use the last formula it is necessary to know the n-harmonic currents values in separate SM contours. It is notably inconvenient for practical calculations carrying out. In this regard it is advisable to refer to simplified SM equivalent circuit with equivalent complex resistance. For such a scheme it is possible to write down
where р—is poles pairs number; m—number of SM stator windings phases.
Taking slipping into account regarding (1) and SM n-harmonic current values as per (8) and (9) for synchronous rotor motion (s = 0) the expression (12) is reduced to form
For instance, for the most significant 2nd harmonic of T1 power transformer magnetizing current and unfavorable conditions combination, when КHV-LV = 1, КLV1-LV2 = 1 and КD = 1, the relative value of asynchronous torque is equal to approx. This value of asynchronous torque is capable to affect to rotor mechanical movement and disturb stable SM operation in case of nominal load on shaft under low power voltage conditions.
For the first time the existence of indirect impact of geoinduced currents on high voltage synchronous engines operation within geomagnetic storms is justified. Danger level of geomagnetic activity against synchronous engines is defined by higher current harmonics level of even sequences occurred in stator windings which is significantly lower in the case when step down substation power supply is delivered from electrical transmission line having cable plants.