Optimized Economic Loading ofDistributionTransformersUsing Minimum Energy Loss Computing

)is paper presents an approach to achieve the Economic Loading of Distribution Transformer (ELDT) based on minimizing the active energy loss.)e effects of the transformer on-load and no-load losses, load factor (LF), and load loss factor (LSF) on the total energy losses are considered.)e adopted technique in this paper consists of two phases, where ELDTis determined analytically in the first phase, and the load profile adequated (LSF) is statistically obtained in the second one.)e results declare that the proposed technique is suitable for the shifts system mode of operation such as in industrial plants. Moreover, this paper investigates the effect of the total active and reactive power and energy losses on ELDT. Finally, numerical examples with software analyses are performed as a valuable tool, which supports the decision-makers to decide trustfully the size of the transformer and its capacity (kVA) during the design stage, as well as to determine the economic loading during the operation based on the effective factors, that is, total power, energy losses, and the load profile.


Introduction
e key power system elements that affect the cost are the generators, transmission lines, and transformers; consequently, the maximum power system benefit is achieved when these elements operate at their maximum efficiency and minimum operational cost. Indeed, the appropriate design, selection, and operational conditions for the power system elements play a crucial factor in reducing the overall operational cost and hence improve the overall system efficiency and availability. e proposed work in this paper focuses on the power and energy losses of the transformer and its effect on the ELDT when waving the other system elements (generators and transmission lines). e transformer is regarded as static equipment where most of the generated power passes through it at different stages. Moreover, the transformer is usually designed for a lifetime of 20 years or more. Apart from small-size transformers, the power and distribution transformer efficiencies usually vary from 95% to 98.5%, taking into consideration the fact that these percentages are even higher for large-size transformers. Moreover, the vast number of transformers installed in a single power system makes the transformer design and its load be significant issues on the overall system cost, efficiency, and availability. Hence, the improper selection of transformers results in additional expenses added to the total investment value, extra cost to the power in terms of power and energy losses, and the extra availability/reliability of the overall system. Based on the above, it is essential to optimize the transformer's loading versus the technical losses, which might have been placed during the power system operation period.
Different approaches were proposed to assess the technical and economic transformers' cost and design issues, such as in [1], where the reduction of overall losses in the distribution transformer was discussed. e impact of the transformer overexcitation on the loading conditions was investigated in [2]. e effect of the transformer location on cost minimization was introduced by [3], where [4] showed how to minimize the iron losses. e research work in [5] classified the losses of transformers into three groups of losses, which are tank losses, core joints losses, and stray losses. e determination of the no-load losses magnitude by using finite element analysis was presented in [6]. e authors of [7] proposed a mathematical model to determine the uncertainty of the transformer on-load and no-load losses. e usage of load loss factor is utilized in [8] to find the amount of the transformer power and energy losses. e reduction of the core cost which resulted from newly developed materials was given in [9]. In [10], the relationship between the no-load current and the level of the total harmonic distortion was investigated, wherein, in [11], the impact of the harmonic level over the nonlinear load was investigated in good detail. e work in [12] proposed a thermal model to monitor the operation and identify the hotspot points.
e oil breakdown strength based on nanofluids was discussed in [13]. e developed methods to determine the transformer's life span considering the annual load factor were presented in [14,15]. In addition, the economics of the transformer and the evaluation of the total owning cost have been thoroughly discussed and analyzed in [3,[16][17][18].
In this paper, the ELDT is achieved for the transformer's capacity, based on the effective factors, parameters, and nature of the working mode during the transformer's lifetime operation. e rest of the paper is arranged as follows. In Section 2, the formulation of the problem is presented, supported by the required mathematical equations, and the transformer losses calculations are presented in Section 3. In Section 4, numerical examples are shown for the proposed technique to simplify the understanding of implementing the technique. Finally, the results, conclusions, and future work proposals are presented in Section 5.

Problem Formulation
e optimum transformer loading is subject to all power and energy losses throughout the operational time, where one year is usually considered enough for the propose of the economic evaluation. Moreover, the objective function can be mathematically formulated as follows: Find It is subject to And where S L(econ.) is the economic loading of the transformer (kVA), ∆E LL is the on-load energy loss of the transformer (kWh), E NLL is the no-load energy loss of the transformer (kWh), S n is the nominal capacity of the transformer (kVA), ∆E represents overall average active energy losses (kWh), and S D is the designed capacity of the transformer (kVA).
In the other sections of this work, the authors look out to optimize the loading of the transformer by finding ELDT, which leads to finding the value of S L(econ.) in (1) in (kVA); this obtained value shall be less than the capacity of the transformer, which was determined during the design stage S D in (2) to avoid the overloading of the transformer (i.e., S D ≥ S L (econ.) ). Usually, the designed capacity (kVA) of the transformer is not matched with the standard unit capacity of the transformers; therefore, the designer should select a transformer with a unit size capacity that is nearest to greater standard capacity (S n ≥ S D ). Consider the total energy losses (∆E) (kWh) in (3), taking into account the future expansion in operation (if any) and the increase in the demand on power as illustrated in numerical examples in Section 4.

Transformer Losses Calculation
e current flows through the transformer windings, causing a power loss even if the secondary winding of the transformer is open. However, when the transformer is loaded, an additional power loss occurs.
e flow of the current in the transformer windings causes active and reactive power losses. e active and reactive power losses in the transformer are explained as follows.

Transformer Active Power Losses.
e active power losses (ΔP) caused by the current flow in the transformer are of two types, which are the iron losses or the no-load losses (ΔP Fe ) and the copper losses or the on-load losses(ΔP Cu ). e first active power losses are independent of the load current, which are called the iron losses or the no-load losses (ΔP Fe ); they are hysteresis losses and eddy current losses. e second type of active power loss is due to the flow of the load current in the transformer's windings, which is called the copper losses or load losses ( ΔP Cu ). e total power losses (ΔP ) then can be calculated from the sum of the power losses caused by the flow of active and reactive power which consist of two components. e first is the active power losses due to flow of active power (ΔP p ) and the second is active power losses due to the flow of reactive power (ΔP q ), expressed as follows: where ΔP represents the total active power losses (kW). ΔP p represents the active power losses due to the flow of active power (kW). ΔP q represents the active power losses due to the flow of reactive power (kW). e active power losses of the transformer with respect to the transformer loading can be expressed as follows: where ΔP Cu represents active power on-load losses (kW), ΔP Fe represents active power no-load losses (kW), and S L represents the loading of the transformer (kVA). e values of ΔP Fe and ΔP Cu are usually a part of the transformer nameplate information or the technical data sheet of the transformer provided by the manufacturers.
In addition, the flow of the reactive power is causing reactive power losses (ΔQ) related to the no-load current (ΔQ o ) and the on-load current (ΔQ L ), respectively. When the transformer is loaded with a power of S L of its nominal capacity S n , then ΔQ with respect to the transformer loading can be presented as follows: where ΔQ represents reactive power losses (kVAR), ΔQ L represents reactive power on-load losses (kVAR), and ΔQ o represents reactive power no-load losses (kVAR).
ΔQ o and ΔQ L can be calculated from the two following formulas, respectively: where I o% is the no-load current (%) and ΔV x% is the impedance voltage of the transformer (%).
Typically the value of current (I o% ) is about 2.20% to 8.5% for distribution transformer and can even vary between 0.1% and 10% of total load current in some special transformers [10]. e value of ∆V x% usually varies between 3.5% and 6.5% and is up to 10% or more in some particular case transformers [19]. e flow of the reactive power causes a loss in active power which can be represented as an equivalent coefficient (C e ). us the active power losses due to the flow of reactive power are expressed as follows: where C e is the equivalent coefficient converting the reactive power losses to active power losses (kW/kVAR). Substituting (5) and (6) into (4) with considering (9) yields the total active power losses that are shown in the following formula: e equivalent coefficient C e can be calculated from the expression in (11) according to [25] as follows: where P is the active load power (kW), Q is the reactive load power (kVAR), and R Tr is the transformer windings resistance (Ω). e transformer windings resistance R Tr according to [18] can be calculated as follows:

Transformer Active Energy
Losses. e product of the power losses (ΔP (t) ) integration over a time (T) gives the energy losses (ΔE), and if the time-varying power losses ΔP (t) over a time (T) are arranged in descending order from ΔP (max) to ΔP (min) in equal integrated areas, then the active energy losses can be expressed as follows: where τ is the equivalent load losses hours (h/year), T represents the year's hours � 8760 h, ΔP (t) represents the time-varying active power losses (kW), and ΔP (max) is maximum value of active power losses (kW). Substituting (10) into (13) yields the following formula: 3.3. e Transformer Power Economic Loading due to the Transferred Energy. If the time-varying loading power (P (t) ) over a time (T) is arranged in descending order from P (max) to P (min) in equal integrated areas, then the active energy consumed over a time period (T eq ) becomes where T eq represents equivalent working hours over a year (h). Based on the triangular power, the relation between the active power, apparent power of the transformer, and power factor (cos φ) is as follows: where cos φ is the power factor. Substituting (15) into (16) yields the following formula: where E is the transferred energy in one year (kWh). e relative energy losses for the transformed energy over the year (T � 8760 h) can be calculated by dividing transformer energy losses in (14) by the transformed energy in (17) which yields the following formula: Equation (18) represents the relative active energy losses of the transformer. When performing some mathematical manipulation, the final relative active energy loss is Mathematical Problems in Engineering e minimum function is obtained when differentiating (19) with respect to (S L ) and equalizing the result to zero as in (20), resulting in (21).
In this case, S L is the optimum loading that is equal to S L(econ.) ; accordingly, the transformer optimum loading is obtained by performing some mathematical arrangement as shown in the following equation: e second derivative of (19) is checked and it is found that (d 2 (ΔE /E))/(dS 2 l ) > 0, which means that the obtained result has a local minimum.
If the solution of equation (22) equals the standard unit size of the transformer, then the achieved solution is regarded as an optimum; however, if it does not equal the standard unit size (which is most probably), then the optimum solution will be obtained by selecting the nearest greater standard unit size; hence, the obtained solution is the economic loading.
In equation (22), the values of T eq and τ are still unknown. However, T eq is determined in the following clause (4) of this section.

Equivalent Working Hours Determination.
Based on the formula published in [20][21][22] which determines the equivalent working hours per month (T (w/m)eq ) (also called maximum load utilization time), the expression of T (w/m)eq is shown in the following equation: T eq � N · T (w/m)eq .
e proof of equations (23) and (24) will be introduced in Appendix A at the end of this study. In equations (23) and (24), A p(w) is the sum of the consumed active energies (day and night) during the time of operation per month (kWh/N), A p(aw) is the consumed active energy after the time of operation per month (kWh/N), D (w/y) is the number of working days over the year (excluding holidays, shutdowns, weekends, etc.), n (h/m) is the number of hours per month (i.e., year's hours/ 12 months � 730 h), n s is the number of working shifts during the time of operation in one day, T eq represents the equivalent working hours (h), T (w/m)eq represents the equivalent working time over the month in hours (h/month), and N is e number of months (N � 1, 2, . . ., 12, where N � 12 for one year). e equivalent working hours (T eq ) in (24) can be calculated by multiplying T (w/m)eq in (23) with the number of the months, which is normally one year (12 months). (13), where the energy loss is determined, and when the timevarying power losses values (ΔP (t) ) are arranged in descending order from ΔP (max) to ΔP (min) as depicted in Figure 1, the areas under the curve of the dotted and straight line are equal.

Maximum Loss Time Determination. Based on
Consider the relation between LSF and the year's hours (T � 8760 h) to obtain τ as in the following equation: where LSF is load loss factor (unit less). When multiplying the LSF factor by the peak power losses (ΔP (max) ) for a period (T), the overall average energy lost (∆E) was given in [23]. Substituting (25) into (13) yields the following formula: e relationship between LF and LSF is given in [24], where the classical computational formula is presented as follows: e value of k in (27) could be any value between zero and one. For instance, k � 0.2 in the UK, 0.3 in the USA, and 0.33 in PL. In this paper and for the proceeding calculations, the value of k is assumed as 0.33.
Many authors in many literature works are using k � 0.30 or k � 0.33, which is appropriate to the urban power grid. In fact, the value of k in (27) might have different values subject to the loads' profile curve. Hence, it is possible to obtain different mathematical equations by using the fitting curve techniques to represent the relation between LSF and LF for the analyzed period.
LF in equation (27) is defined as the ratio of the average load (P (avg) ) supplied during the designated period (T) to the maximum load (P (max) ) occurring in the same period (T) [25]. Consequently, LF can be expressed as follows: where P (avg) is the average load power (kW).

Numerical Example
e proposed method is implemented in a numerical example, where the calculations are as follows.

Load Input Data.
It is intended to build an industrial plant with an expected maximum power of 584 kW to supply the processes requirement and associated buildings, with an expected power factor equaling 0.80. e load is expected to increase in the future by 20%. e plant could initially work five days a week in one shift or two shifts or work continuously in three shifts and one shift off.
To find the ELDT, the calculation in this example is conducted to one shift of work: 8 hours a day, five days a week, excluding holidays. To maintain the process outside the working hours' time some of equipment and lighting remain ON (essential load). e amount of energy required to maintain the process outside the time of operation is equal to 20% of the energy consumed during the normal time of operation.
However, at the end of the example, the results of twoand three-shift modes of work are presented.

Solution:
Calculate: (i) T eq represents the equivalent working hours (h). e consumed energy (E) over the year is estimated by considering the essential load to maintain the plant in the running state after the working hours. It is calculated using the following formula: ( A p (aw) /A p(w) ) � 20%. T eq is calculated for one shift (8 hours a day, five days a week, excluding holidays, etc.), as in (23).
And, based on equation (24), (ii) E is the active energy consumed over the year. From equation (15), E is obtained: (iii) S L is the apparent power of the load. From equation (16), we obtain (iv) LF is the load factor. Equation (28) (vi) τ represents the equivalent load losses' hours (h/ year). By using equation (25), we obtain (v) Q is reactive power load. From the triangular power, for the relation between active, reactive, and apparent power, the reactive power equals (vi) S n is the nominal capacity of the transformer (kVA). Considering the constraints in (2), the selection of the transformer capacity shall be equal to or more than the designed value; also, taking into account the future expansion in the load by 20%, the capacity of the transformer (S n ) is Round to the nearest standard available size, select a transformer with a capacity equal to S n � 1000 kVA.
Based on (24), (25), (27), and (28), the load data for the one-shift mode of work is calculated. e result of the calculations is summarized in Table 1.

Transformers Technical Data.
e technical data for the transformer (data belonging to transformer type (THUE23/ 1000) manufactured by ALSTHOM) which is provided in the manufacturer's technical data sheet and the calculated data are given in Table 2. e given data of the transformer and the calculated data are summarized in Table 2, where ΔQ o and ∆Q L are calculated based on (7) and (8), respectively, as follows: ΔQ o � 2.12 100 × 1000 � 21.20 kVAR, × 1000 � 50.00 kVAR.
Moreover, R Tr and C e can be calculated based on (11) and (12), respectively, as follows:

Transformer Power Economic
Loading. Based on the given data, the ELDT can be calculated from equation (22) for one working condition for one shift.
However, the ELDT for two shifts and that for three shifts are similar; the difference is only in calculating the time of T eq and τ. e results for the work in shifts system are presented in Table 3.
However, when neglecting the effect of T eq time and τ in (22), the ELDT reduces to 511 [kVA] for this transformer.
A lot of exercises can be performed by the recently devolved software.

Results, Conclusions, and Future Work
is paper presents a detailed approach for studying and calculating the ELDT. Different schedules of working hours for different transformers input technical data were taken into account. e evaluation shows that the ELDT is accruing when the objective is to minimize the active energy losses rather than active power losses, which was justified by the product of the equivalent working hours to the equivalent load losses, which is always more than or equal to one; the results of calculations are presented in Table 3. e submitted formula for the ELDT shows the effect of the on-load losses (the copper losses) and the no-load losses (the core losses) on the result of the transformer economic loading. e optimality is achieved when the load losses are equal to no-load losses; however, this solution is very costly and sometimes practically is impossible; hence, a compromise is needed. e explained numerical examples show the impact of the technical parameters (in terms of load and no-load losses) and operational time on the ELDT, where the ELDT could reach a value equal to or less than the transformer rated capacity. e more the relative value of the on-load losses to the no-load losses is reduced, the more enhancement of the ELDT is achieved, the utilization end is enhanced, and the operational cost of the transformer is reduced.
In addition, the software has been devolved to find the ELDT faster and easier, which is very beneficial, especially when having two or more transformers with different capacity and technical data working in parallel to feed a variable load.
is paper, along with the numerical examples, could be a base (platform) for other researches or future work, such as the effect of energy loss cost in terms of money on ELDT and the effect of load changing and power factor on the ELDT.     where (365: D (w/y) ) is the number of working days/ year (excluding holidays, etc., as above); and the value (A p(aw) /A p(w) ) represents the ratio between the active energies consumed (day and night) after and during the time of operation, multiplied by the number of hours per day (24(h/d)) and divided by the number of months per year (12(m/y)) to obtain the value of this part in hours per month. (iii) When the plant works in one or two shifts, that is, n s � 1 or n s � 2, respectively, 5 days/week (excluding the public holidays and the time for maintenance and shutdown), such as in most public institutions, all the three parts of equation (A.1) will not equal zero as illustrated in the numerical example.

Conflicts of Interest
e authors declare that they have no conflicts of interest.
Mathematical Problems in Engineering 7