Open access peer-reviewed chapter

Determination of Optimal Transformation Ratios of Power System Transformers in Conditions of Incomplete Information Regarding the Values of Diagnostic Parameters

By Lezhniuk Petro Demianovych, Rubanenko Oleksandr Evgeniovich and Rubanenko Olena Oleksandrivna

Submitted: October 8th 2018Reviewed: February 5th 2019Published: December 9th 2019

DOI: 10.5772/intechopen.84959

Downloaded: 66

Abstract

On the base of damage rate analysis of power transformers and methods of electrical energy system (EES) modes control the necessity of using the results of on-line diagnostics of LTC transformers not only for determinations of the expending of further operation or equipment repair but also for calculation of optimal transformation coefficients (with account of the suggested RRCT) for their application in the process of modes control has been proved. Improved method of determination of control action, realized by the LTC transformers by means of comparative analysis of the results calculation of EES modes with quasi resistances of the circuit branches. Such peculiarity of the suggested method of determination of control actions by LTC transformers, as the account of RRCT, in the process of EES mode control provides such advantages as reduction of the damage rate of the equipment, reduction of active power losses in EES. Due to the peculiarities of the method of determination of control actions by LTC transformers, with the account of their technical state, perspectives of the development and introduction in EES modern microprocessor-based systems of automatic control of transformers LTC open.

Keywords

  • on-line diagnostics
  • control
  • normal modes
  • active power losses
  • similarity theory
  • neuro-fuzzy modeling
  • basic similarity criteria
  • membership functions
  • uncertainty

1. Introduction

Characteristic feature of present day situation are the attempts of utility companies to increase energy efficiency in conditions of continuing aging of high voltage equipment.

Practice of large-scale introduction of the intelligent support of decision making of solution processes proves their efficiency. One of the directions of efficiency increase of electric energy transportation is improvement of methods and means of active power losses reduction on conditions of maintaining reliable operation of high voltage equipment, including outdated equipment.

The set of electric energy system (EES) states and the processes of transition from one state into another her is EES mode (further-mode), characterized by the parameters, for instance, electrical voltages an substations loads, currents in transmission lines transformation ratios of the transformers, etc., normal operation mode 07 of EES modes. Control by the system of on-line dispatching control (LDC). Modern LDS technologies, for instance, provided by smart grids concept, are aimed at improvement of its information support. This enables the optimal implement more efficiently energy-saving technologies in electric systems, when out of date high-voltage equipment is using.

Means of similarity theory, in particular the criterion method (KM), can effectively solve and analyze optimization problems [1]. Criteria-based method can be defined as a set of techniques and principles, according to which the analysis, comparison and interpretation of baseline data to provide scientific and practical conclusions. The ultimate goal of studies using the criterial method is to reveal regularities, which under certain conditions can be represented as the law control [2]. The main purpose of KM is to find the variant of the process or the object. Most often, in the further analysis, the parameters are used as reference. According to his idea, KM is close to the geometrical programming [3]. This use of duality of optimization problems is the replacement of the direct problem for the corresponding dual. The main difference is that the basis for geometric programming is inequality between geometric and arithmetic averages, and the background of the KM—matrix properties of dimensions or indicators [4]. This is the meaning of dual variables. In geometric programming are weighting factors, in KM—similarity criteria. That is, the result of solving the tasks of the KM is values of criteria of similarity or, in other words, the optimum ratio of the individual parameters, and not they themselves. This specific feature characteristic only of KM, determines its scope.

2. Analysis of the literature data and problem set up

In [5] the technique of voltage drop decrease in separate parts of distribution electric grids is, suggested, but it does not take into account technical state of regulating devices [6]. In [7] on the example of Indian electric grids the statistics of the increase of power, transmitted in electric grid is considered.

It is stated that in order to improve the grid reliability and efficiency of energy transmission it is necessary to use power transformers, equipped with LTC and automatic or automated control systems for such LTC control. It enables to control power flows in EES by means of LTC so that the parameters of electric grids modes were within the limits of normal values of equipment (transmission lines, switching devices, transformers, etc.) parameter.

It is provided by usage of FACTS technologies. For instance, in [8], the possibilities of using phase-shift transformer for power flows change in electric energy system to reduce power losses in the process of energy transmission in transmission lines of Slovak Republic are considered. In [9] high price of FACTs technologies usage in energy branch for the reduction of electrical power is proved. In [10] three variants of power flows control in electric grid, using three FACTs devices are considered with Thyristor Controlled Series Capacitor (TCSC), Static Synchronous Series Compensator (SSSC) and phase-shift transformer (PST), but in [11] any attention is paid to the state of equipment, used for modes control. In [12] the conclusion is made that prolongation of power transformers operation term for 20–30 years is more profitable than their replacement by new ones, and the quantity of power transformers in the USA, that have been in operation for more than 25 years (certificate resource – 25 years) is approximately 65%.

In [13] attention is paid to the system of continuous monitoring of technical state of power transformers, the given system is used at the transformers of joint-stock company “Magnitogorsk Metallurgical Complex”, HYDRAN analyzer, methods of localization and identification of faults, practical necessity of partial discharges control is underlined, but the results of diagnostics during modes control are not paid attention to.

In [14] it is noted that in local electric system in order to provide stable operation and indices of electric energy quality it is necessary to use modern control systems that take into consideration voltages in nodes and frequency and eliminate emergency deviations. At the same time, in [15] modeling of non-stationary critical operation modes of EES in the process of parameters change in wide limits by means of application of non-linear mathematical models attention is paid to. This enables to study the consequences of such modes, promptly take measures, aimed at their prevention or elimination. In pages [16] technical state of the equipment of these systems is not taken into account, this can lead to the damage of the equipment and undersupply of energy to the consumers.

Thus, the problem of development of the methods of diagnostics results account during control of EES modes is not solved.

It is known that operation control (RTOC) in Ukraine is carried by a man. Overloading of this person with a great volume of diagnostic parameters data, especially in conditions of limited time for decision-making, leads to their actions. In the process of modes on-line control, especially post-occident modes. It is expedient to assess the state of equipment by generalized indices, for instance, by residual resources coefficient of the transformers (RRCT). The development of the method of on-line diagnostics of the transformers and the account of RRCT in the process of ESS modes control for minimization of total losses of active power are not considered in literature sources and is the subject of authors study.

3. Objectives and tasks of the research

The objective of the research is the development of the method of diagnostics of the transformers with LTC and account of PRCT values in the process of EES modes control for minimization of active power losses. To realize this objective the following problems are to be solved:

  • substantiate the expediency of applying the results of diagnostics of LTC- transformers in the process of optimal control of EES modes;

  • develop neuro-fuzzy model of residual resource coefficient of the transformers (RRCT);

  • develop the method of RRCT values and power transformers with LTC state account in the process of EES modes control.

4. Materials and methods of transformers diagnostic study

Automation of the process of power flow control may be provided by means of centralized remotely controlled alternative usage of switching devices (LTC) of the transformers. Under such conditions, there appears the possibility of the analysis of control actions of separate LTC on mode parameters of EES by means of the feedback. This approach improves the operation quality of adaptive control automatic systems of the LTC position control. For this purpose, at considerable changes of load schedule it is necessary to perform ranking of the transformers with LTC by the quality of their impact on maintaining parameters of the modes.

Realization of measures, aimed at reduction of power losses is limited by the possibilities of the equipment involved in the provision of mode; namely, by its technical state. It is known, that the damage of high voltage equipment during mode control (for instance, power transformers) leads to losses, which considerably exceed the cost of electric energy, saved as a result of losses decrease. Failure rate of the outdated high voltage equipment (power transformers, shunting reactors, instrument current and voltage transformers, switches, etc.) increases, when such equipment has been in operation for more than 25 years [17]. Taking into consideration the fact that the control of EES modes is accompanied by the operation of switching devices, regulation devices of transformers, emergence of switching surges, ferro-resonances, currents increase in power and instrument transformers, transmission lines, etc., then the control of modes must be realized, taking into consideration their technical state [18] and possible expenses for their replacement or repair.

Thus, it is necessary to know current state of high voltage electric equipment of EES, which is in operation during modes control.

5. Determination of current technical state of power transformers

We will consider the method of determination of power transformers current state and RRCT values in the process of EES modes control on the example of power high voltage transformers, which have on-load-tap changing device.

We suggest to evaluate technical state of power transformer by means of the analysis of the value of its residual resource coefficient. Power transformer residual resource coefficient has the dimensionality in relative units and can change in the process of operation in the range from one (the best technical state) to zero (the worst technical state, when the transformer must be removed out of service for inspection, repair, replacement, etc.).

Then, we will consider the example of residual resource coefficient determination of the transformer АТDCTN 125000–330/110. First we will study the statistics of failure rate of such transformers. Table 1 contains the example of possible reasons and amount of transformers removal out of service, that is close to data, published in studies [19].

Transformer elementDesignationParameter nameUnits%
WindingsZkWinding deformation81.6
t0Deterioration of contact joints state102
Pi.pIdle power that characterizes of the magnetic quality153
InsulationRinRContamination of isolation6513.4
WHumidification of the isolation4810
BushingskbushDefects of bushings7415.2
OilCADGcContent of dissolved gases7114.6
PCAHigh moisture content and deviations of other parameters of the oil439
CADGdDischarges in oil6413.2
LTCkdef.LTCLTC defects459.3
Cooling systemImotor or ImtThe current of oil pump drive motor142.9
tcooloCoolers temperature163.3
TankktankTank leakage122.5
Total485100

Table 1.

Reasons of removing out of service power transformers.

In Table 1 such symbols are used: Zk is the resistance of the transformer windings (during measurements in short – circuit mode); t° is the temperature of contact points (for instance, bushing of the bus duct or with winding lead); Pi.p. idle mode power, that characterizes the quality of magnetic circuit; Rin is the resistance of the insulation for revealing the contamination and aging of solid and liquid insulation (also it is necessary to determine the capacity and dielectric loss tangent, also it is desirable to determine the degree of polymerization); W humidification of the isolation; kresid.res.bush or kbush is the residual resource coefficient of the bushings; CADGC is residual resource coefficient of the transformer by the results of chromatographic analysis of dissolved gas in the transformer oil of the tank and LTC (ethylene, ethane, methane) of the transformer, that characterizes oil contamination by the gases, dissolved in it and among them acetylene and hydrogen (for revealing of discharges); PCA residual resource coefficient of the transformer by the results of physical–chemical analyses of transformer oil from transformer tank, contactor and LTC tap changer; CADGd is residual resource coefficient of the transformer by the results of chromatographic analysis of the dissolved hydrogen and acetylene in the transformer oil of the tank and LTC of the transformer tap changer in order to reveal the discharge; kdef.LTC or kLTC coefficient of the transformer LTC residual resource; Im is the current of electric motors oil pumps and fans of cooling system; t°cool is coolers temperature; k tank is the residual resource coefficient of the transformer tank, determined by the availability (takes the value “0”) or absent of oil leakage (takes the value 1).

From Table 1 shows that the transformers are often displayed in repairs due to moisture and oil contamination, insulation and high-voltage inputs defects.

The task of creating a mathematical model complicated with incomplete initial data as part of the parameters known at the time of payment, such as the reasons for the need for additional studies. To establish reciprocal links diagnostic parameters very constructive simulation technology is unclear. This simulation allows to obtain more reliable results compared to the results of existing diagnostic systems.

In Table 1 under the term the controlled diagnostic parameter we mean the parameter deviation of which from the norm helped to remove the transformer out of service or was taken into account in the process of its removal out of service. In Table 1 the following diagnostic parameters are given: parameters, that characterize the state of the windings, insulation, bushings, oil, LTС, cooling systems, tank.

Having analyzed the data of Table 1 the scheme was created that shows whether dependent or independent is the impact of diagnostic parameters on the coefficient of total residual resource of the transformer (Figure 1).

Figure 1.

Structural diagram of the model of total residual resource coefficient of the transformer.

Figure 1 does not show mutual impact of one controlled diagnostic parameter on the other one; it is shown either in dependent or independent manner how these parameters influence the coefficient of total residual resource of power transformer (PT).

In Figure 1 over the parameter the percentage amount of revealed faulty transformers by the given parameter is shown, that is given in percent from the total amount of faulty transformers.

The blocks with parameters whose deviations from the norm substantiate the necessity of the output of the transformers for repair, are shown sequentially and are shown kres – residual resource coefficient of the transformer (RRCT). In parallel, blocks with parameters are also depicted. A large change in these parameters proves the necessity of outputting a power transformer (PT) for repair. PT is repaired in case of deviation from the norms of these parameters. This is due to the requirements for the reliability of the transformers. In each of the given blocks parallel can be allocated but it are not shown to simplify the calculation (for instance, currents of electric motors of oil pumps and fans).

In order to obtain the generalized parameter of the residual resource of the transformer, it is proposed from the known values of diagnostic parameters to pass to the corresponding values of residual resources coefficients (in relative units) by each diagnostic parameter. This will allow you to take into account the value of all diagnostic parameters and the impact of each of them.

These coefficients are defined in relative units by (1) and that is why they characterize total output of the transformers from the moment of their technical state control to transition to boundary state that is residual technical resource (12). Residual resource coefficient ki1by i1th diagnostic parameter:

ki1=xi1,limхі1,curхі1,limхі1,in,E1

where xi1,limis admissible limit normative value of i1thdiagnostic parameter; xi1,curis value of i1thdiagnostic parameter at the moment of control; xi1inis initial value of i1thdiagnostic parameter (at the moment of putting into operation of new equipment or after repair), i1 is number of diagnostic parameter.

We perform the reduction of the circuit by the following expressions. For serial part of the circuit (Figure 1) the coefficient of total residual resource is found by the expression:

ktot.resid.res.=τ=1νkτpτ,E2

where kτ is the coefficient of residual resource of PT by τth diagnostic parameter; τ is τth diagnostic parameter; ν is the amount of blocks in the serial part of the circuit of Figure 1, pτ probability of control parameters deviator from maximum permissible normalized value of this parameter is found by means of the expression (3):

pτ=yτm2,E3

where yτ is a number of controlled parameter deviations from admissible limiting normalized value of this parameter, which were revealed by means of τth diagnostic parameter control (τ for serial part of the circuit) from the total number of the revealed deviations of controlled parameters from admissible limiting normalized value; m2 is total quantity of the revealed deviations of controlled diagnostic parameter from their admissible limiting normalized values.

For parallel part of the circuit the coefficient of total residual resource is found by the expression (4)

ktot.resid.res.=1j=1m11kres,jpj,E4

where kres,j is the coefficient of residual resource of РТ by jth diagnostic parameter; j is number of jth diagnostic parameter; m1 is a quantity of blocks(parameters) in parallel part of the circuit that is reduced.

The coefficient of total residual resource of РТ is determined by the expression (5):

kres=kwind.kin.kbushkoilkLTCkcoolktank,E5

where kwind, kin, kbush, koil, kLIC, kcool, ktank are known at the moment of calculation values of the coefficient of residual resource: of the windings, of the insulation, of bushings, of the oil, of LTC, of system of the cooling, of tank of the transformer, by the elements of the transformer, correspondingly.

For the creation of mathematical model of residual resource coefficient of the transformer parameters were used, by each of these parameters the conclusion regarding the state of the transformer can be made. But none of these parameters completely characterizes technical state of the transformer, it only shows certain changes of technical state of power transformer.

Mathematical model of residual resource coefficient of the transformer was created by means of MatLab. Using this model, it is possible to edit the already created (5) probabilistic sample of teaching data. These data help to obtain analytical dependence of residual resource coefficient of the transformer on diagnostic parameters in the form of the polynomial. For seven input parameters of the model, that randomly changed from 0 to 1, the coefficient of total residual resource of the transformer (5) was determined, where input parameters of the model were reduced to relative units of their deviation from the norm.

By means of Anfis Editor using hybrid teaching algorithm and applying Sugeno algorithm of neuro-fuzzy conclusion neuro-fuzzy model of residual resource coefficient of the transformer (using subclusterization method) was obtained.

Figure 2 contains the copy of screen saver in Matlab environment where the structure of the obtained neural network is shown.

Figure 2.

Structure of Anfis-network of the transformer.

For each input variable of neuro-model four linguistic terms with Gaussian membership functions were used:

kres.і1=fxі1σі1cі1=exі1cі122σі12,E6

where δi1and ci1are numerical parameter; δi12in probability theory it is called dispersion of the distribution(14), and the second parameter ci1is mathematic expectation; i1 is input parameter of neuro-fuzzy model, that corresponds to diagnostic parameter (i1 = 1, 2, 3, 4, 5, 6, 7), xi is value of i1th input parameter of the model: x1 – kwind, x2 – kin, x3 – kbush, x4 – koil, x5 – kLTC, x6 – kcool, x7 – ktank.

These are such terms as: “normal” values of diagnostic parameter, “minor deviations” of diagnostic parameter value, “prefault” values of diagnostic parameter, “emergency” value of diagnostic parameter.

For determining the value of total residual coefficient neuro-fuzzy non-linear autoregressive model of the total residual resource coefficient of the transformer is used. This model establishes neuro-fuzzy non-linear transformation between the values of residual resource coefficients by diagnostic parameters and total residual resource coefficient of the transformer (7):

ktot.resid.res=Fkwind.kin.kbushkoilkLTCkcoolktank,E7

where F is neuro-fuzzy functional transformation.

For determination of the value of total residual resource coefficient of the transformer, we use Takagi-Sugeno model of logic conclusion.

Mathematical model of total residual resource coefficient is the system of logic equations (8).

IFkwind."normal"ANDkin."normal"ANDkвush."normal"ANDkoil"normal"ANDkLTC"normal"ANDkсоol."normal"ANDktank"normal"THENktot.resid.res=a11kwind.+a12kin.+a13kbush.+a14koil+a15kLTC+a16kсоol.+a17ktank+c1IFkwind"minordeviationsANDkin"minordeviation"ANDkвв"minordeviation"ANDkoil"minordeviation"ANDkLTC"minordeviation"ANDkсоol"minordeviation"ANDktank"minordeviation"THENktot.resid.res=a21kwind.+a22kin.+a23kbush.+a24koil+a25kLTC+a26kсоol.+a27ktank+c2IFkwind"prefault"ANDkіn"prefault"ANDkвв"prefault"ANDkoil"prefault"ANDkLTC"prefault"ANDkсоol"prefault"ANDktank"prefault"THENktot.resid.res=a31kwind.+a32kin.+a33kbush.+a34koil+a35kLTC+a36kсоol.+a37ktank+c3IFkwind"emergency"ANDkіn"emergency"ANDkвв"emergency"ANDkoil"emergency"ANDkLTC"emergency"ANDkсоol"emergency"ANDktank"emergency"THENktot.resid.res=a41kwind.+a42kin.+a43kbush.+a44koil+a45kLTC+a46kсоol.+a47ktank+c4E8

Output of the model total ktot.resid.res. is found as weighted sum of conclusions (8) of rules base written in the form the system of logic equations:

ktot.resid.res.=j2=1m3wj2aj21kwind.+aj22kіn.+aj23kвush.+aj24koil++aj25kLTC+aj26kсоol.+aj27ktank+cj2,E9

where 0wj21is the degree of execution (weight) of the j2th rule that is determined by the correspondence of real changes of diagnostic parameters of the transformer.

ANFIS is the simplest network of direct propagation that contains adaptive nodes, using the teaching rules the parameters of these nodes are arranged to minimize the error between the real output of the model ktot.resid.mod. and real total residual resource coefficient ktot.resid.res of the transformer

δ=1N1k3=0N11ktot.resid.res.modk3ktot.resid.res.k32min,E10

where N is number of rows in teaching sample; k3 is the number of the row in teaching sample, starting from the row with consecutive number “0”.

Taking into account the iterative computation experiments carried out the vector of membership functions parameters is determined in Table 2.

ParametersInput parameters of the modelName of the termNumber of the ruleParameters of membership function
σС
Winding stateKwind.Normal10.38250.7944
Minor deviation20.4790.5197
Prefault30.49030.5668
Emergency40.40.1697
Insulation statekіn.Normal10.36530.8698
Minor deviation20.46420.6104
Prefault30.51020.5267
Emergency40.39490.1742
State of BBkbush.Normal10.32020.9221
Minor deviation20.34190.7649
Prefault30.49140.5376
Emergency40.40320.1925
State of oilKoilNormal10.43690.9273
Minor deviation20.34040.9674
Prefault30.4120.599
Emergency40.40310.2057
State LTCKLTCNormal10.39840.973
Minor deviation20.33160.963
Prefault30.44680.5881
Emergency40.44280.2349
State of cooling systemkсоol.Normal10.34391153
Minor deviation20.35070.9706
Prefault30.4370.597
Emergency40.42630.2397
State of tankKtankNormal10.34540.9506
Minor deviation20.38011017
Prefault30.45820.6273
Emergency20.54510.564

Table 2.

Parameters of membership function.

It is seen from Figure 2 that in the process of formation of the structure of neuro fuzzy model of the transformer seven inputs and one output of this model were set. Each of seven inputs has four terms. That is, each set of possible values of input parameters of the model is conventionally divided into four subsets: “normal” values of input parameter, “miner deviations” of the values of input parameter, and “prefault” values of input parameter, “emergency” values of input parameter. Membership degree of each value of input parameter to corresponding set of values is determined by Gaussian membership function. The model is intended for determining the numerical value of total residual resource coefficient of the transformer, that is why it has one output. This numerical value is found by means of solution of linear equation, that describes the dependence of the coefficient of total residual resource of the transformer on input parameters.

The obtained neuro-fuzzy model allows to determine the value of total residual resource coefficient of the transformer depending on the values of input parameters residual resources coefficients by each of controlled diagnostic parameters. The error of PPCT mathematical model changes from +0,004 relative units, if PPCT equals 0, to −0,032, when PPCT equals 1.

Taking into account the data of the Table 1 and Table 2 and (9) we obtain mathematical model of the coefficient of total residual resource in the form:

IFkwind."normal"ANDkіn"normal"ANDkвв"normal"ANDkoil"normal"ANDkLTC"normal"ANDkсоol"normal"ANDktank"normal"THENktot.resid.res=0,6166kwind+0,4125kіn+0,4618kвв+1,83koil+1,804kLTC++0,0462kсоol.+1,96ktank5,377IFktank"minor deviation"ANDkin."minor deviation"ANDkвush"minor deviation"ANDkoil."minor deviation"ANDkLTC"minor deviation"ANDkсоol."minor deviation"ANDktank"minor deviation"THENktot.resid.res=0,0393kwind.+0,2609×kin.+0,1086×kвush.0,37koil0,1459kLTC0,02387kсоol.0,05863ktank+0,1288IFkwindprefaultANDkіnprefaultANDkвush.prefaultANDkoilprefaultANDkLTCprefaltANDkсоol.prefaultANDktankprefaultTHENktot.resid.res=0,2165kwind..0,3714kin.0,4678×kвush.0,514koil0,882kLTC0,5302kсоol.1,406ktank+3,88IFkwind."emergency"ANDkіn"emergency"ANDkвв"emergency"ANDkoil"emergency"ANDkLTC"emergency"ANDkсоol."emergency"ANDktank"emergency"THENktot.resid.res=0,03166kwind.0,06144kіn.0,387kвush.+0,06koil+0,3199kLTC0,026kсоol.0,006ktank+0,003E11

The obtained neuro-fuzzy model allows to determine the value of total residual resource coefficient of the transformer depending on the values of input parameters residual resources coefficients by each of controlled diagnostic parameters. The error of PPCT mathematical model changes from +0.004 relative units, if PPCT equals 0, to −0.032, when PPCT equals 1.

Despite the complexity of dependences, mathematical model of residual resource coefficient of the transformer (11) may be used for programming neuro-fuzzy controller in order to create the device for on-line determination of transformer state by means of analysis of residual resource coefficient of the transformer value.

6. Account of the forecast current value of residual resource of the transformers in the process of control of EES modes

It is known that in the process of operation, energy enterprise plans to remove out of service the equipment in the overhaul, cost of is forecast. Removal of the transformer into overhaul in a planned number of years (TWF) of trouble-free operation (12 years) provides certain list of works and their expected cost Boh pl. For instance, for 330/110 kV transformers of 125–250 MVA power the cost (B) of such repair is 770–11,550 $. We propose to assume that removal out of service the transformers into current repair requires unscheduled cost.

The cost of repair may increase by the value ΔВ1, replacement of damaged blocks of the transformer and additional work, connected with the replacement. These costs are not provided in case of “typical” planned overhaul

ΔB1=i=1nBieγikres.iβi,E12

where Bi is the cost of replacement of ith damaged block of the transformer and additional work, connected with this replacement, n is a number of damaged blocks that require unscheduled replacement; kres.iβiis residual resource coefficient of ith block that requires unscheduled replacement; γ and β are coefficients, that characterize the impact of residual resource coefficient on the expected cost of unscheduled repair or replacement of ith block of the transformer (is determined by means of processing of statistic data).

Repair cost may increase by the cost of ΔB2 (as compared with expected) in case of enlarged current (instead of planned overhaul) repair of the transformer, that did not operate for planned number of years:

ΔВ2j=1eαjTj1ВOH,E13

where j is a number of the transformer, Tj is time, the ith transformer functioned after putting into operation or after the last overhaul (enlarged current) repaint to the moment of mode control, λ is the coefficient, that characterizes the intensity of cost growth ∆B2 that depends on the construction of the transformer, conductions and operation mode (is determined experimentally), BOH is the cost of transformer overhaul.

It should be noted that removal the transformer out of service takes place not only as a result of relay protection, emergencies control automation operation but also by a person responsible for safety operation by the results of control of diagnostic parameter, values of which sometimes only approaches to limiting values.

Within the context of creation of modern Smart Grids and to provide safe, reliable, quality and economic efficient operation of EES it is necessary to perform the control over active power overflow to realize by means of the transformer, performing reliable and information archons on the mode. That is why, we suggest to take into account the coefficient of regulating transformer limitation:

kwind,j=1kres,jВcq,j,E14

where Bcq is the coefficient of repair cost value growth of the jth transformer.

Bcq,j=ΔB1,j+ΔB2,jΔB1,j+Bpl.,j.E15

As the example, we will consider 23 nodes 230/138 kV test circuit (Figure 3). In branches 11–9, 12–14, 12–9, 11–4, and 3–7 transformers АТDCTN-63000/230/138, АТDCTN-100000/230/138 and АТDCTN-125000/230/138 are installed. Initial node loads, complex transformation ratios and corresponding transformers LTC positions (number of taps) are given in Tables 3 and 4.

Figure 3.

Scheme of grid 230/138 kV for 23 nodes.

BranchesR (Ohm)X (Ohm)kactivekreactive
№ of beginning№ of the end
11100.6270.64870
1290.379.280.64980
1190.3130.64790
730.2111.530.650
120.495126.47110
1310.39840.22110
1541.51616.09210
2462.46424.12910
2694.64936.56510
399.88220.96210
4951.03819.74910
5104.34216.81610
898.8215.12410
8102.0672.14510
11235.20722.79310
111428.56622.11210
12235.20725.1810
121365.59651.10110
231358.71945.75910
14162.64520.57810
15162.3388.40410
161717.45713.70110
1619311711.20610
17180.952276.17610
172271.41555.70410
212246.02335.86610
7155.6824.86510
21180.8736.85110
21151.66612.9610
19201.34910.47410
20130.7415.71310
12100.31140.65240
10626.47111.52210

Table 3.

Information of circuit branches.

№ of the nodeU (kV)Phase (grad)Pload (MW)Qload (MVAr)Pgen (MW)Qgen (MW)
1136.34−5.731082218230
2135.54−6.27187.1576172.930
3141.22−0.94176.436.2600
4136−8.2741500
5137−8.1168.1613.4400
6136.36−10.3129.225.2700
7218.842.9919.41.9400
8142.57−7.72169.2934.65131.55131.75
9141.55−5.62756600
10141.64−7.29191.14900
11222.36−3.12401000
12219.92−4.1254.8817.6400
13236.157.8900495150
14228.341.56184.337.050101.39
15233.2211.21304.3261.44235.251.32
16233.45101002018580
17237.7915.2835.6413.8600
18241.516.92323.0165.96417.1176.19
19231.477.6177.3836.2600
20233.817.311282600
21241.517.8300425.7146
22241.525.9500420−3.69
23234.6026554417.85281.37

Table 4.

Information of the circuit nodes.

Knowing the circuit and normal node parameters we define transformation ratios.

ka.оpt=1diagReNk.bal.bZCeJUb1Ebal.aE16
kr.opt=diagImNk.bal.bZCeJUb1EbalrE17

where Nk.bal.b is the second matrix of branches connection in contour balanced transformation ratios; Z diagonal matrix of complex branches resistances; Ce is the matrix of currents distribution coefficients for economic mode of electric network (corresponds to minimal losses of electric energy); J vector-column of currents in nodes; Ub is the voltage of basic node; E*bal.a and E*bal.r are balancing is electric moving in relative units (EMF) in relative units (active and reactive components).

Taking in the consideration the discrete character of LTC switching’s, errors of instrument transforms, errors of data transmission channels and recommendations [20] we assume that non-sensitivity zone of active power losses may be considered as regulating actions on LTC of the transformer—to be 3% [21].

As initial conditions we assume that in accordance with load graph LTC of transformers 9–11, 9–12, 11–10, 12–10 have transformation ratios 0.6413 (14 tap), 0.6347 (14 tap), 0.6397 (14 tap), 06446 (14 tap).

It should be noted that further changes of operation mode were realized at admissible voltage deviations ±5%, from nominal voltage Unom.

Regulation of the transformer 7–3 is inexpedient on conditions of the usage of the given technique of determination of transformation ratios.

We define the losses of active and reactive power in the branches of the circuit at current transformation ratios (Table 5).

P-ters*№ of nodes
1234567891011121314151617181920212223
Pload, MW19.4437.4338.817.0213.6324.553.4927.0860.538.228.410.43029.4939.6208.5580.7535.4726.880053
Qload, MVAr3.9615.27.973.452.694.80.355.5414.529.82.13.3505.93843.3316.497.255.460010.8
ΔPƩ, МW4.49
ΔQƩ, МVar29.05
11–100.6437
№ tap14
12–100.6446
№ tap14
9–110.6413
№ tap14
9–120.6397
№ tap13
9–120.65
№ tap14

Table 5.

Parameters of the current mode.

Parameters.


ΔSΣbr=ΔPΣbr+jΔQΣbr=3j=1mΔSbr,j,E18

where ΔSbr=diagΔUbrÎbris vector-column of complete power losses in the branches of the circuit. Ebr is vector column of the current in branches, m is a number of the branch in the circuit, ΔUbr=MΣUnodeis vector-column of phase voltages in the nodes, ΔPbr, ΔQbr is vector-column of active and reactive power losses in the branches of the circuit(correspondingly).

We define transformation ratios (16–17) and position of LTC on condition of minimal amount of switchings (in order to maintain switching resources of LTC) to provide minimal losses of active power in branches of the circuit of Table 6.

P-ters№ of nodes
1234567891011121314151617181920212223
Pload, MW19.4437.4338.817.0213.6324.553.4927.0860.538.228.410.43029.4939.6208.5580.7535.4726.880053
Qload,. MVAr3.9615.27.973.452.694.80.355.5414.529.82.13.3505.93843.3316.497.255.460010.8
ΔPƩ,. МW4.42
ΔQƩ,. МVar28.69
11–100.6513
№ tap14
12–100.6542
№ tap15
9–110.6507
№ tap14
9–120.6521
№ tap14
9–120.65
№ tap14

Table 6.

Transformation ratios for the current mode.

As a result of realization of control actions, mode optimization power losses were reduced from ΔS1 = 4.49 + j29.05 (MVA) for mode (Table 5) to ΔS2 = 4.42 + j28.69 (MVA). Thus, the effect of realization of transformer LTC switchings is ΔS1-ΔS2 = 0.07 + j0.36 (MVA). The transition from the current to mode may be performed by switching the LTC of the transformer, installed in the branch 9–12 from 13 tap to 14 and carry out transformer regulation changing position of LTC from 14 tap to 15. We will consider the transition to another stage of daily load graph(load increase), its parameters are given in Table 7, and optimized transformation ratios and corresponding mode parameters - in Table 8.

P-ters№ of nodes
1234567891011121314151617181920212223
Pload, MW108187.1176.47468.16129.219.4169.227539.14054.880184.3304.310035.64323.0177.312800265
Qload, MVAr227636.261513.4425.271.9434.6566101017.64037.0561.442013.8665.9636.26260054
ΔPƩ, МW61.85
ΔQƩ, МVar412.73
11–100.6513
№ tap14
12–100.6542
№ tap15
9–110.6507
№ tap14
9–120.6521
№ tap14
9–120.65
№ tap14

Table 7.

Parameters of the mode after load change.

P-ters№ of nodes
1234567891011121314151617181920212223
Pload, MW108187.1176.47468.16129.219.4169.227539.14054.880184.3304.310035.64323.0177.312800265
Qload, MVAr227636.261513.4425.271.9434.6566101017.64037.0561.442013.8665.9636.26260054
ΔPƩ, МW61.80
ΔQƩ, МVar412.68
11–100,665
№ tap16
12–100.651
№ tap14
9–110.6753
№ tap17
9–120.659
№ tap15
7–30.65
№ tap14

Table 8.

Parameters of mode without taking into account technical state of the transformers.

As a result of performing control actions, mode optimization we succeeded in decreasing power losses from ΔS1 = 61.85 + j412.73 (MVA) for the mode (Table 7) to ΔS2 = 61.80 + j412.68 (MVA). Thus the effect of LTC transformer switchings is ΔS1-ΔS2 = 0.05 + j0.5 (MVA).

If as a result of determining the coefficient of regulating effect limitation for circuit transformers (Figure 3) the following values are obtained: kwind,9–11 = 0.85, kwind,12–9 = 0.4, kwind,11–10 = 0.3 а kwind,12–10 = 0,2, then expected quasi-decrease of losses, taking into account these coefficient will be defined.

Control actions are performed by the transformer, installed in the branch 9–11, namely, we change position of LTC with 14 tap on 15, in this case, the expected losses of active power are ΔP9–11 = 61.82 (MW). We find the decrease of active power losses ΔРΣ - ΔP9–11 = 61.85–61.82 = 0.03 (MW), however, having taken into account the coefficient of regulating effect limitation, losses decrease change δPquasi.9–11 = (ΔРΣ - ΔP9–11)∙kwind,9–11 = 0.0255 (MW). New quasi-losses ΔPquasi.9–11 = ΔP9–11 + δPquasi.9–11 = 61.82 + 0.0255 = 61.8455 (MW). Results of the calculation of other transformers are given in Table 9.

Transf.Ktr.cur.Ktr.opt.Kwind.jΔPtr.j (MW)ΔРΣ - ΔPtr.j (MW)δPquasi.j (MW)ΔPquasi j. (MW)
Ncur.Nopt.
9–110.65070.6590.8561.820.030.025561.8455
1415
12–90.65210.64190.4661.830.020.009261.8392
1413
11–100.65130.6650.3461.8350.0150.005161.8401
1416
12–100.65420.6510.2561.840.010.002561.8425
1514

Table 9.

Results of limiting effect coefficients calculation for circuit transformers.

We define mode parameters for the circuit with quasi-resistances from Table 11 and corrected transformation ratios from Table 12.

We find losses of active power in the branch, that contains the transformer, as function the element of vector-column of complete power losses in the branches of the circuit by the expression

ΔPα=ReΔSα,E19

where ΔSα = ΔUα∙Iα is element of vector-column of power losses in the branches, that contain transformers, ΔUα is kth element of vector-column of phase voltages drop in the branches, and Iα is the current of the branches with transformers couplings, α is the number of row, that correspond to the branch with transformer couplings in vector-column ΔSbr.

The value of quasi resistance in kth-branch:

Zα=ΔSαÎα2,E20

where α=k+β, where k – is the number of the row of the first branch, that contains the transformer, β is the coefficient of the change of consecutive number of branch, that contains the transformer, it changes in the range from 0 toψ1, ψis the number of branches, containing transformers.

Applying this algorithm, according to (20), quasi-resistances of the branches, containing transformers are found. The results of the calculations are given in Table 10. The aim of control is provider of minimum of all system active power losses that is determined by the expression

ParametersTransformer 9–11Transformer 12–9Transformer 11–10Transformer 12–10
Branch resistance, Ohm0.3 + j130.37 + j9.280.3 + j270.3 + j14
Quasi resistance of the branch, Ohm0.32 + j25.20.4 + j13.20.36 + j28.40.35 + j19.2

Table 10.

Quasi-resistances of transformers branches of the circuit.

P-ters№ of nodes
1234567891011121314151617181920212223
Pload, MW108187.1176.47468.16129.219.4169.227539.14054.880184.3304.310035.64323.0177.312800265
Qload, MVAr227636.261513.4425.271.9434.6566101017.64037.0561.442013.8665.9636.26260054
ΔPƩ, МW62.47
ΔQƩ, МVar422.16
11–100.665
№ tap16
12–100.651
№ tap14
9–110.659
№ tap15
9–120.6419
№ tap13
7–30.65
№ tap14

Table 11.

Parameters of normal mode after loads change, taking into account technical stale of transformers.

P-ters№ of nodes
1234567891011121314151617181920212223
Pload, MW108187.1176.47468.16129.219.4169.227539.14054.880184.3304.310035.64323.0177.312800265
Qload, MVAr227636.261513.4425.271.9434.6566101017.64037.0561.442013.8665.9636.26260054
ΔPƩ, МW62.46
ΔQƩ, МVar422.28
11–100.664
№ tap15
12–100.641
№ tap13
9–110.6753
№ tap17
9–120.6492
№ tap14
7–30.65
№ tap14

Table 12.

Parameters of mode after loads change, taking into account technical state of transformers and corrected transformation ratios.

ΔF=i=1nΔPimin.E21

If ΔFmin.=ΔPmin.is the minimum value of the efficiency function (active power losses), ΔFcur.=ΔPcur.current value of efficiency function (active power losses), n is the total number of branches in the circuits, kтr.min. value of transformation ratio at which calculated losses of active power are minimal, then the dependence of active power losses change (values of efficiency function ΔFcur.=ΔFcur.ΔFmin.) in relative units on the values of transformation ratios k=kcur.ktr.min.(Figure 4) for various transformers will be built.

Figure 4.

Charts of dependencies of changes in active power loss on the values of transformation ratios for small loads mode.

Thus, the transition from the current to the mode can be realized by switching the LTC of the transformer, installed in the branch 9–12 from 13 tap to 14tap and perform regulation of the transformer in branch 12–10, changing LTC position from tap 14 to 15 tap are shown in Figure 5.

Figure 5.

Charts of dependencies of changes in active power loss on the values of transformation ratios for large loads mode, without taking into account technical state of the transformers.

As a result of realization of control actions the mode will be reached by transformer switching of the branch 11–10 from 14 to 16 tap of LTC, transformer of the branch 12–10 from 15 tap to 14 tap, transformer of the branch 9–11 transformer from 14 to 15 tap and transformer of the branch 9–12 from 14 to 13 tap of LTC, respectively, how are shown in Figure 6.

Figure 6.

Transformation ratio for large loads mode, taking into account technical state of the transformers.

We see that due to consideration of technical state of transformers, their ranking occurred by the measure of impact on the reduction of active power losses. To reach mode now it is more expedient to use a transformer of 9–12 branch as it during one switching of LTC from 14 tap to 15 tap reduces most active power losses.

7. Discussion of the results of transformation ratios of EES transformers determination, taking into account the state of transformers

The analysis of the articles showed that neuro-fuzzy logic methods are used to solve various problems of operating electric power systems, such as improving power quality [18], classifying the faults of the electrical equipment based on sequence components [19], developing of the controller of the tariff, which based on neuro-fuzzy logic for the for distributing active power between a micro network and EPC for improving energy quality [20].

The error of RRCT determination by means of the developed mathematical neuro-fuzzy model, as compared with teaching sample and to the opinion of independent experts does not exceed the error of the devices, measuring diagnostic parameters. Such results are explained by complex usage of probability theory methods neuro-fuzzy modeling and modern software Matlab. These results also confirm the information provided in the article by Moudud Ahmed, Naruttam Kumar Roy. In their article [21], it is written that the use of automatic systems for adaptive control of electric power systems (EPS) based on neuro-fuzzy modeling and based on an inference system (ANFIS) is promising method. This improves EPS performance, for example, reduces power losses. Similar positive results of using neuro-fuzzy logic are described in the article by Priyanka Ray and A.K. Sinha [22]. This article says that the use of neuro-fuzzy logic has allowed the development of a hybrid control system that provides the maximum generated electrical power of hydro, wind and solar power plants even under incomplete data on current weather conditions and power consumption.

Also, in the works [23, 24, 25] of the authors H. Suna, R. M. Velasquez; J.W.M. Lara; Dong Ling; Yao-Yu Xu; Yu Liang; Yuan Li; Ning Liuand Quan and Jun Zhang were reviewed methods of intelligent diagnostics of transformers that use fuzzy logic and in the future can be applied to improve diagnostic systems and other power equipment.

Such feature of the suggested method for determining the control actions of LTC-transformers, as account of PPCT, in the process of ES modes control, provides such advantages as reduction of the equipment damage rate, decrease of active power losses in the EPS. Due to the of the peculiarities method of determination of control action of LTC-transformers, taking into account their technical state, perspectives of developments and introduction in EPS of modern microprocessor –based systems, automatic control of LTC of transformers become possible.

As compared with the known method of voltage drop control on the branches of EPS circuits, with the method of overloads decrease of transmission lines, at the expense of redistribution of power overflows in EPS, decrease of active power losses in the process of transportation by means of LTC- transformers, the suggested method allows to select, by means of account the suggested RRCT, the transformer for EPS mode control, that would simultaneously provide the reduction of power losses and is more reliable.

Usage of quasi-resistances of circuit branches, that unlike the transformers used, in the process of calculation of nominal resistances of the branches, take into account transformers state and possible losses of utility companies due to possible damages, allows to calculate EPS mode in rise of transformers transformation ratio change and by means comparison of calculated power losses select the most efficient transformer.

The suggested peculiarity of application the method of neuro-fuzzy modeling (usage in teaching sample the model of transformer resource instead of measured values of diagnostic parameters - calculated and partially corrected by independent experts of coefficients of residual resource) enables to take into account simultaneous impact on RRCT the results of both current and periodic control.

The drawback of the suggested mathematical neuro-fuzzy model of RRCT is necessity of large data base regarding coefficient of residual resource of diagnostic parameters CRRDP (Coefficient of residual resource of the diagnostic parameter) for specific transformers. Attempt to reduce database or use the model from other similar transformer results in the increase of model error. Limitation on the usage of RRCT model is the necessity of application only on one – investigated transformer. Therefore, we need models for each transformer. The method of determination of control actions by LTC transformers does not take into account voltage limitations in nodes and current limitations in the branches of the circuit.

Further development of the given research will be realized in the development of mathematical models of other types of high voltage equipment, involved in the process of EPS modes control, damage of which areas place (Figure 7).

Figure 7.

Damage of high voltage equipment in EPS: (a) 750 kV shunt reactor; (b) current transformer and 750 kV air circuit breaker; (c) 33 kV voltage transformer; (d) 750 kV SF6 circuit breaker; (e) air circuit breaker; (f) 110 kV SF6 circuit breaker.

Problems of the considered research development are caused by the necessity of long lasting experiments and observations over the processes of aging and development of high voltage equipment damage, processes of EPS modes parameters change not only on computer ad mathematical models of the equipment and EPS modes and on real equipment.

8. Conclusions

  1. Analysis of damager ate of power transformers and methods of the EPS modes control allows to state that it is a necessary to use the results of on-line diagnostics of LTC-transformers not only to determine the expediency of further operation or repair of the equipment and for calculation transformation values (with the account of the suggested RRCT) for their usage in the process of modes control.

  2. The model enables, by means of accounting, of both current and retrospective values of diagnostic parameters on RRCT and determine its current value. That is necessary for automatic and automated reliable and control of EPS modes.

  3. Improved method of determination of control actions by LTC- transformers, by means of comparative analysis of the results of EES modes with quasi resistances of circuit branches, enables to soled the transformer and calculate transformation ratio that provides minimal amount of LTC switching.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Lezhniuk Petro Demianovych, Rubanenko Oleksandr Evgeniovich and Rubanenko Olena Oleksandrivna (December 9th 2019). Determination of Optimal Transformation Ratios of Power System Transformers in Conditions of Incomplete Information Regarding the Values of Diagnostic Parameters, Fuzzy Logic, Constantin Volosencu, IntechOpen, DOI: 10.5772/intechopen.84959. Available from:

chapter statistics

66total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

The Fuzzy Logic Methodology for Evaluating the Causality of Factors in Organization Management

By Nazarov Dmitry Mikhailovich

Related Book

First chapter

Microassembly Using Water Drop

By Taksehi Mizuno

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us