Abstract:
The present application discloses a method to optimize operation of a transformer cooling system, the corresponding cooling system, and a method to determine the capacity of Variable Frequency Drives (VFD) that are used in the transformer cooling system. The method comprises: preprocessing the initial data input by user; collecting the on-line data, and calculating the optimized control command to meet the requirement of the transformer loss, top-oll temperature variation and noise; and executing the control actions by controlling a controllable switch and/or sending a control command to a VFD. Compared with the existing prior arts, the proposed solutions are much more intuitive and practical in the field of the cooling system.

Description:
FIELD OF THE INVENTION 
       [0001]    This invention relates to the cooling technical field, and more particularly to a method to optimize operation of a transformer cooling system, the corresponding transformer cooling system, and a method to determine the capacity of Variable Frequency Drives (VFD) that are used in the said transformer cooling system. 
       BACKGROUND OF THE INVENTION 
       [0002]    Transformer is one of the most critical components of a substation, whose safety, reliability and efficiency are of high importance to the overall power grid. For each transformer, especially power transformers with voltage level 110 kV and above, a dedicated cooling system consisting of multiple motor-fan units is required to keep the winding temperature within an acceptable range. The operation of the transformer is therefore closely related to 1) how the cooling system is designed and 2) how the cooling system is operated. 
         [0003]    As to the cooling system design, it is common understanding that variable speed operation of these cooling fans can achieve higher efficiency compared with fixed speed operation. Therefore transformer cooling systems tend to install VFDs for motor-fan units to ensure high efficiency operation, the system architectures are shown in  FIG. 1A  and  FIG. 1B . However, these two types of architecture have their own disadvantages. The first architecture as shown in  FIG. 1A  requires high capital investment because it installs VFD for each motor-fan chain; plus if the motor-fan chain is mostly working at rated speed, VFD solution might lower the efficiency due to its own power losses. The second architecture as shown in  FIG. 1B  can relatively reduce the capital investment because it uses one big VFD to drive a plurality of motor-fan chains jointly at the same operation point. But the disadvantages are also obvious: Firstly, each motor-fan chain has low efficiency when the VFD utilized capacity is relatively low; secondly, there are different ways for load distribution among different VFD-fed motor-fan chains to meet the same total output requirement. It is not always true to distribute the load evenly among individual chains in order to have optimal system efficiency. 
         [0004]    As to cooling system operation, the core is how to control the winding temperature. Normally, lower winding temperature leads to the lower copper loss of winding. However, the power consumption of the cooling system will be higher at the same time, meaning that the overall efficiency, considering both transformer winding and the cooling system itself, might be less optimal. 
         [0005]    Besides efficiency, the variation of the winding temperature is also one key factor which will affect the lifecycle of the transformer. The more frequency the temperature varies, the faster the transformer aging will be. It could be so that the efficiency of the transformer is optimized, however at a cost of shortened transformer lifetime. 
         [0006]    For transformer operated at urban area, noise level is also one important criterion to consider in order to reduce the impact on the neighbouring residents especially at night. Currently, few solution is available to control the cooling system to tackle the noise problem. 
         [0007]    To overcome above shortcomings, the person skilled in the art aims to solve two problems as follows. 
         [0008]    1) How to design the cooling system to realize speed regulation for the motor-fan loads selectively with less capital investment on VFDs. 
         [0009]    2) How to improve the operation efficiency of transformer by cooling control considering the transformer copper loss, the motor-fans power consumption and the speed regulation of VFD. 
         [0010]    3) How to control the winding temperature as well as its variation in order to extend the lifecycle of the transformer and meanwhile achieve the best overall system efficiency. 
         [0011]    4) How to operate the cooling system to optimize not only the efficiency and lifecycle, but also minimize the noise level so as to reduce the negative impact on the surrounding environment. 
       SUMMARY OF THE INVENTION 
       [0012]    The objects of the present invention are achieved by a method to optimize operation of a transformer cooling system, the corresponding cooling system, and a method to determine capacity of the VFDs that are used in the said transformer cooling system, in order to improve the operation efficiency of the whole transformer with limited capital investment on cooling system hardware upgrade, and meanwhile to extend the transformer lifecycle and lower the noise level of the transformer system. 
         [0013]    According to one aspect of the invention, said method to optimize the operation of the transformer cooling system, comprises the following steps: preprocessing the initial data input by user; collecting the on-line data, and calculating the cooling capacity required to meet the requirements of transformer loss; and executing the control actions by controlling a controllable switch and/or sending a control command to a VFD. 
         [0014]    According to a preferred embodiment of the present invention, said calculating the optimized control command step further considers the requirement of the top-oil temperature variation and/or the noise level. 
         [0015]    According to a preferred embodiment of the present invention, said calculating the optimized control command step further considers the requirements according to the weighting factors of the transformer loss, the top-oil temperature variation and the noise level, which are capable of pre-defining by the user. 
         [0016]    According to a preferred embodiment of the present invention, said preprocessing step comprises the following steps: collecting parameters of the transformer type, the transformer ratio, and the ratio of load losses at rated current to no-load losses; collecting parameters of the transformer thermal model; collecting parameters of the tap changer mid position, the step voltage and the present tap changer position; collecting parameters of the cooler type, the fan number and the power of the radiator; and collecting the relationship curve between the fan noise and the fan capacity. 
         [0017]    According to a preferred embodiment of the present invention, said preprocessing step further includes the following steps: calculating the transformer copper loss; calculating the winding temperature; calculating the load current of different sides; and calculating the power consumption of cooling system. 
         [0018]    According to a preferred embodiment of the present invention, said on-line data includes: the load current, the temperatures and the status of the cooler; and said calculating step comprises the following steps: calculating the cooling capacity required to meet said requirement; calculating the number of fans including the fan driven by the VFD; comparing the fans required with the existing fans in operation; and leading to different possible operation solutions in accordance with the comparison. 
         [0019]    According to a preferred embodiment of the present invention, the actual transformer loss P K ′ under specific load level for three-winding transformer is calculated by the following equation: 
         [0000]    
       
         
           
             
               P 
               k 
               ′ 
             
             = 
             
               
                 
                   1 
                   + 
                   
                     α 
                      
                     
                       
                         
                             
                         
                          
                         
                           θ 
                           w 
                         
                       
                       _ 
                     
                   
                 
                 
                   1 
                   + 
                   
                     75 
                      
                     
                         
                     
                      
                     α 
                   
                 
               
                
               
                 ( 
                 
                   
                     
                       β 
                       1 
                       2 
                     
                      
                     
                       P 
                       
                         k 
                          
                         
                             
                         
                          
                         2 
                          
                         
                             
                         
                          
                         N 
                       
                     
                   
                   + 
                   
                     
                       β 
                       2 
                       2 
                     
                      
                     
                       P 
                       
                         k 
                          
                         
                             
                         
                          
                         2 
                          
                         
                             
                         
                          
                         N 
                       
                     
                   
                   + 
                   
                     
                       β 
                       3 
                       2 
                     
                      
                     
                       P 
                       
                         k 
                          
                         
                             
                         
                          
                         3 
                          
                         
                             
                         
                          
                         N 
                       
                     
                   
                 
                 ) 
               
             
           
         
       
     
         [0000]    Wherein,  θ   w  is the average winding temperature; a is temperature factor; β 1 , β 2 , β 3  are load factors; P  k1N , P k2N , P k3N  are the winding losses at rated current. 
         [0020]    According to a preferred embodiment of the present invention, said top-oil temperature variation Dθ 0  over time dt is calculated by the following equation: 
         [0000]    
       
         
           
             
               D 
                
               
                   
               
                
               
                 θ 
                 o 
               
             
             = 
             
               
                 { 
                 
                   
                     
                       
                         
                           [ 
                           
                             
                               1 
                               + 
                               
                                 RK 
                                 2 
                               
                             
                             
                               1 
                               + 
                               R 
                             
                           
                           ] 
                         
                         x 
                       
                       · 
                       Δ 
                     
                      
                     
                         
                     
                      
                     
                       
                         θ 
                         or 
                       
                       · 
                       
                         100 
                         
                           X 
                           cor 
                         
                       
                     
                   
                   - 
                   
                     ( 
                     
                       
                         θ 
                         oi 
                       
                       - 
                       
                         θ 
                         a 
                       
                     
                     ) 
                   
                 
                 } 
               
               · 
               
                 dt 
                 
                   τ 
                   o 
                 
               
             
           
         
       
     
         [0000]    Wherein, Δθ or  is top-oil temperature rise in the steady state at rated losses (K); R is ratio of load losses at rated current to no-load losses; K is load factor; τ 0  is average oil time constant; θ oi  is the top-oil temperature at prior time; θ a  is the ambient temperature; X cor  is the rate of cooling in operation. 
         [0021]    According to a preferred embodiment of the present invention, the total noise from the transformer and the fan Lp t  is calculated by the following equation: 
         [0000]    
       
         
           
             
               Lp 
               t 
             
             = 
             
               { 
               
                 
                   
                     
                       
                         Lp 
                         
                           N 
                            
                           
                               
                           
                            
                           1 
                         
                       
                       , 
                     
                   
                   
                     
                       
                         Lp 
                         fan 
                       
                       = 
                       0 
                     
                   
                 
                 
                   
                     
                       
                         
                           Lp 
                           
                             N 
                              
                             
                                 
                             
                              
                             1 
                           
                         
                         + 
                         
                           10 
                            
                           
                               
                           
                            
                           
                             lg 
                              
                             
                               [ 
                               
                                 1 
                                 + 
                                 
                                   10 
                                   
                                     - 
                                     
                                       
                                         
                                           Lp 
                                           
                                             N 
                                              
                                             
                                                 
                                             
                                              
                                             1 
                                           
                                         
                                         - 
                                         
                                           Lp 
                                           fan 
                                         
                                       
                                       10 
                                     
                                   
                                 
                               
                               ] 
                             
                           
                         
                       
                       , 
                     
                   
                   
                     
                       
                         Lp 
                         
                           N 
                            
                           
                               
                           
                            
                           1 
                         
                       
                       &gt; 
                       
                         Lp 
                         fan 
                       
                     
                   
                 
                 
                   
                     
                       
                         
                           Lp 
                           fan 
                         
                         + 
                         
                           10 
                            
                           
                               
                           
                            
                           
                             lg 
                              
                             
                               [ 
                               
                                 1 
                                 + 
                                 
                                   10 
                                   
                                     - 
                                     
                                       
                                         
                                           Lp 
                                           fan 
                                         
                                         - 
                                         
                                           Lp 
                                           
                                             N 
                                              
                                             
                                                 
                                             
                                              
                                             1 
                                           
                                         
                                       
                                       10 
                                     
                                   
                                 
                               
                               ] 
                             
                           
                         
                       
                       , 
                     
                   
                   
                     
                       
                         Lp 
                         fan 
                       
                       &gt; 
                       
                         Lp 
                         
                           N 
                            
                           
                               
                           
                            
                           1 
                         
                       
                     
                   
                 
               
             
           
         
       
     
         [0000]    Wherein, Lp fan  is the fan noise; Lp N1  is the transformer noise. 
         [0022]    According to a preferred embodiment of the present invention, said different possible operation solutions comprises: switching on the integer fans with lower utilization rate and driving the rest fans by VFD with calculated frequency; switching off the integer number of fans with higher utilization rate and driving the rest fans by the VFD with calculated frequency; or changing the fan driven by the VFD with calculated frequency. 
         [0023]    According to a preferred embodiment of the present invention, said control actions includes: the start or stop of the fans; controllable switch operation; or VFD frequency regulation. 
         [0024]    According to another aspect of the invention, said method to determine capacity of the VFDs used in the said transformer cooling system, comprises the following steps: inputting parameters and the objectives of the transformer loss, the top-oil temperature variation and the noise; calculating the Net Present Value (NPV) curve versus of the VFD capacity which shows the relationship between the saved energy loss and the VFD cost; calculating the VFD capacity limit for the pre-defined top-oil temperature variation; calculating the VFD capacity limit for the pre-defined noise level; determining the VFD capacity which has highest NPV, meanwhile within the limits to fulfil both top-oil temperature variation and noise level requirements. 
         [0025]    According to a preferred embodiment of the present invention, said highest NPV is determined with the following steps: calculating saved energy loss of the cooling system due to the VFD; calculating the capital cost of the VFD; evaluating the NPV of the VFD considering both benefit and cost; and selecting the VFD capacity with the highest NPV. 
         [0026]    According to another aspect of the invention, said transformer cooling system, comprises a central controller, a transformer and a plurality of fans to cool down said transformer. Said transformer cooling system further comprises a shared VFD bus fed by VFD and an AC bus fed by AC power source, both of which being controlled by said central controller. Said shared VFD bus is shared by two or more motor-fan chains and selectively driving one, two or more said motor-fan chains. 
         [0027]    According to a preferred embodiment of the present invention, each of said motor-fan chain connects to a controllable switch, which switches said motor-fan chain among connecting to said AC bus, connecting to said shared VFD bus, and disconnecting from power supplies. 
         [0028]    Compared with the existing prior arts, the solution of the present invention saves the capital investment to upgrade cooling system hardware for transformer cooling system operation optimization. Another benefit of the present invention is that it can optimize the real-time operation efficiency of transformer by coordinating the transformer copper loss, cooling system power consumption, and VFD settings for individual motor-fan chain, meanwhile realize transformer lifecycle extension and noise level limitation. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0029]    The subject matter of the invention will be explained in more details in the following description with reference to preferred exemplary embodiments which are illustrated in the drawings, in which: 
           [0030]      FIGS. 1A and 1B  show an electrification scheme of the conventional transformer cooling system; in which  FIG. 1A  illustrates the structure of respectively installing VFD for each motor-fan chain, and  FIG. 1B  illustrates the structure of a plurality of motor-fan chains jointly driven by one VFD; 
           [0031]      FIG. 2  shows an electrification scheme of the transformer cooling system according to an embodiment of the present invention; 
           [0032]      FIG. 3  is the overall flow-chart for VFD capacity determination according to an embodiment of the present invention; 
           [0033]      FIG. 4  is the flow-chart for net present value calculation due to transformer efficiency improvement by installing different capacity of VFD in the cooling system according to an embodiment of the present invention; 
           [0034]      FIG. 5  is the main flow-chart for operation optimization of transformer cooling system according to an embodiment of the present invention; 
           [0035]      FIG. 6  illustrates a flow chart of parameters preprocessing procedures according to an embodiment of the present invention; 
           [0036]      FIG. 7  illustrates a flow chart of control command determination according to an embodiment of the present invention; 
           [0037]      FIG. 8  illustrates a flow chart of control command execution according to an embodiment of the present invention. 
       
    
    
     DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
       [0038]    Exemplary embodiments of the present invention are described in conjunction with the accompanying drawings hereinafter. For the sake of clarity and conciseness, not all the features of actual implementations are described in the specification. 
         [0039]    According to the first preferred embodiment, the electrical system design of the transformer cooling system is shown in  FIG. 2 , which consists of two power supply schemes for motor-fan loads, including an AC line supply and a VFD supply (e.g. VFD 1  in  FIG. 2 ). 
         [0040]    As shown in  FIG. 2 , one or more motor-fan chains can be connected to the VFD bus, the AC bus or disconnected from power supplies respectively through the controllable switches. That means, the motor-fan chains can only have one out of three statuses at one time: connecting to AC line, connecting to VFD, or disconnecting from power supplies. 
         [0041]    By coordinating the VFD and controllable switches, the start-up process of motor-fan loads can be optimized. As shown in  FIG. 2 , a motor-fan load can be switched to VFD for soft start. After completing the start-up process, it can be switched back to the AC line if it is operated at the rated output. In order to optimize the operation, the status information of VFD and controllable switches are all transmitted to a central controller. Besides these, the central controller also gets access to the real-time transformer load data, oil temperature and ambient temperature. With all these data, the controller performs the efficiency optimization calculation, top oil and its variation calculation, and noise level calculation of the whole transformer. After that, it will send out the control command to controllable devices, e.g. controllable switches for gross temperature regulation, and VFD for fine temperature regulation. 
         [0042]    According to the second preferred embodiment, the size of the VFD can be determined by techno-economic analysis to ensure best cost-effectiveness of the given type of transformer. The higher the VFD capacity is, the more accurate the temperature control will be, which can contribute to overall operation performance improvement. However, the cost of the VFD will also increase which will affect the business case. Meanwhile, different type of transformers have different cooling capacity requirement. The sizing of VFD should also take this into account.  FIG. 3  shows the overall procedures for VFD capacity determination. Firstly, the parameters and the operation objectives, e.g. transformer loss, top-oil temperature variation and expected noise level will be input by the users; secondly, the NPV curve which shows the relationship between transformer loss and VFD capacity will be calculated; thirdly, the VFD capacity limitations to achieve the predetermined top-oil temperature variation and noise level requirements will be calculated; fourthly, the VFD capacity can be determined which has the highest NPV for transformer loss reduction, and meanwhile can fulfill the lifecycle and noise level requirement. 
         [0043]      FIG. 4  illustrates how to calculate the NPV curve versus VFD capacity through transformer system efficiency improvement. In  FIG. 4 , P VFD represents the rated capacity of the VFD; P VFD0  and AP VFD  represent the initial capacity and incremental capacity of VFD used for iteration By calculating the save energy loss through VFD, and the corresponding capital investment of VFD, the net present value curve can be obtained versus different VFD capacities. 
         [0044]    According to another preferred embodiment, the central controller performs the optimization calculation in real-time. The flowchart is shown in  FIG. 5 . Whenever the optimization result changes, the central controller will update the control commands for VFD and/or controllable switches respectively. 
         [0045]    Step 1: the first step of the flowchart is to preprocess the initial data input by user. The detailed information is shown in  FIG. 6 , where totally five groups of data will be collected as follows:
       1) The transformer type, ratio, and ratio of load losses at rated current to no-load losses. The method uses them to calculate the copper loss.   2) Winding exponent, oil exponent, hot-spot to top-oil gradient, hot-spot factor, ambient temperature, average oil time constant, winding time constant, hot-spot-to-top-oil gradient at start, hot-spot-to-top-oil gradient at the rated current, top-oil temperature rise in steady state at rated losses, top-oil temperature rise at start, the load permissible in % of nameplate rating when all fans inoperative. The method uses them to calculate the hot-spot temperature which can be regarded as the winding temperature.   3) Tap changer mid position, step voltage, present tap changer position. The method uses them to calculate the load current of different sides.   4) Cooler type, fan number, the power of radiator. The method uses them to calculate the power consumption of cooling system.   5) Relationship curve between fan noise and fan capacity.       
 
         [0051]    After the preprocessing, all information except real-time data will be ready for calculation. 
         [0052]    Step 2: the second step, the central controller collects the load current, temperatures and the status of cooler. And then calculate the cooling capacity which can meet the requirements of transformer loss, top-oil temperature variation and/or transformer noise requirements. The detailed procedures for calculating winding loss, oil-temperature variation and noise are described from Section A to Section C; and the method to combine this three dimensional control objectives together using weighting factors are described in Section D. 
         [0053]    After the optimal cooling capacity is obtained by the central controller, the control strategy will lead to three possible operation solutions as shown in  FIG. 7 : if the number of fans required is greater, less than or equal to the number of existing fans in operation. 
         [0054]    If the number of fans required is nf_next, the number of existing fans is nf_prior, then 
         [0000]        n   fΔ=fix(   n   f   _   next )−fix( n   f   _ prior);
 
         [0000]        n   VFD   =n   f   _   next −( n   f   _   prior+nfΔ );
 
         [0055]    If n fΔ &gt;0, switch on the corresponding number of fans; otherwise, switch off the corresponding number of fans. And the rest fans driven by VFD should change n VFD . 
         [0056]    When to increase or decrease percentage of transformer cooling, the central controller calculates the number of motor-fan chains needed, it is assumed that the number of motor-fan chains in operation is m 1 .n 1 , the number calculated is m 2 .n 2 , where m i  is the integer number and n i  is the percentage of cooling capacity which will achieved by VFD. The central controller gets the integer number of motor-fan chains by m 2 −m 1 . The speed regulation of VFD can be calculated by n 2 . The priority of motor-fan chains depend on the utilization time. The central controller prioritizes the motor-fan chains according to the utilization time. Then, the central controller selects to start the motor-fan chain with lower utilization time, and selects to stop the motor-fan chain with higher utilization time. 
         [0057]    A. Basic Mathematics for Transformer Loss Calculation 
         [0058]    For three-winding transformer, the actual winding loss under specific load level is 
         [0000]    
       
         
           
             
               
                 
                   
                     P 
                     k 
                     ′ 
                   
                   = 
                   
                     
                       
                         1 
                         + 
                         
                           α 
                            
                           
                             
                               
                                   
                               
                                
                               
                                 θ 
                                 w 
                               
                             
                             _ 
                           
                         
                       
                       
                         1 
                         + 
                         
                           75 
                            
                           
                               
                           
                            
                           α 
                         
                       
                     
                      
                     
                       ( 
                       
                         
                           
                             β 
                             1 
                             2 
                           
                            
                           
                             P 
                             
                               k 
                                
                               
                                   
                               
                                
                               1 
                                
                               
                                   
                               
                                
                               N 
                             
                           
                         
                         + 
                         
                           
                             β 
                             2 
                             2 
                           
                            
                           
                             P 
                             
                               k 
                                
                               
                                   
                               
                                
                               2 
                                
                               
                                   
                               
                                
                               N 
                             
                           
                         
                         + 
                         
                           
                             β 
                             3 
                             2 
                           
                            
                           
                             P 
                             
                               k 
                                
                               
                                   
                               
                                
                               3 
                                
                               
                                   
                               
                                
                               N 
                             
                           
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   ( 
                   1 
                   ) 
                 
               
             
           
         
       
     
         [0059]    Where, 
         [0060]      θ   w : the average winding temperature; 
         [0061]    α: temperature factor; 
         [0062]    β 1 , β 2 , β 3 : load factor; 
         [0063]    P k1N , P k2N , P k3N : the winding loss at rated current; 
         [0064]    Assume n f  equals to the total required cooling power divided by rated cooling power of each motor-fan chain P f , which consists of two parts: n r , which is the integer part, and n v , which is the decimal part. 
         [0065]    Assume n r  is contributed by fans operated at rated speed; and n v  is contributed by fans controlled by VFD operated at partial speed. The total power demand can be expressed as (2), where η is the efficiency of the VFD. 
         [0000]        P   fans   =n   r   ×P   f   +n   v   ×P   f /η  (2)
 
         [0066]    If all fans are at the same speed and all driven by VFDs, we have 
         [0000]        P   fans   =n   f   ×P   f /η  (3)
 
         [0067]    The transformer loss can be calculated as formula (4) 
         [0000]        f   1   =P   t   =P   k   ′+P   fans   +C    (4)
 
         [0068]    Where, C is constant the power consumption of other parts. 
         [0069]    B. Basic Mathematics for Transformer Top-Oil Temperature Calculation 
         [0070]    The top-oil temperature variation over time dt is calculated by equation (5), 
         [0000]    
       
         
           
             
               
                 
                   
                     D 
                      
                     
                         
                     
                      
                     
                       θ 
                       o 
                     
                   
                   = 
                   
                     
                       { 
                       
                         
                           
                             
                               
                                 [ 
                                 
                                   
                                     1 
                                     + 
                                     
                                       RK 
                                       2 
                                     
                                   
                                   
                                     1 
                                     + 
                                     R 
                                   
                                 
                                 ] 
                               
                               x 
                             
                             · 
                             Δ 
                           
                            
                           
                               
                           
                            
                           
                             
                               θ 
                               or 
                             
                             · 
                             
                               100 
                               
                                 X 
                                 cor 
                               
                             
                           
                         
                         - 
                         
                           ( 
                           
                             
                               θ 
                               oi 
                             
                             - 
                             
                               θ 
                               a 
                             
                           
                           ) 
                         
                       
                       } 
                     
                     · 
                     
                       dt 
                       
                         τ 
                         o 
                       
                     
                   
                 
               
               
                 
                   ( 
                   5 
                   ) 
                 
               
             
           
         
       
     
         [0071]    Then the difference between the top-oil temperature and a given value is f 2 , 
         [0000]        f   2 =abs(θ oi   +Dθ   o −θ om )   (6)
 
         [0072]    Where, 
         [0073]    Δθ or : top-oil temperature rise in the steady state at rated losses (K); 
         [0074]    R: ratio of load losses at rated current to no-load losses; 
         [0075]    K: load factor; 
         [0076]    τ 0 : average oil time constant; 
         [0077]    θ oi : the top-oil temperature at prior time; 
         [0078]    θ a : the ambient temperature; 
         [0079]    θ om : the given value of top oil temperature; 
         [0080]    X cor : the rate of cooling in operation, which can be calculated by equation (7), where 
         [0081]    N is the rated current ratio of ONAN condition to ONAF condition; 
         [0000]    
       
         
           
             
               
                 
                   
                     X 
                     cor 
                   
                   = 
                   
                     
                       ( 
                       
                         N 
                         + 
                         
                           
                             X 
                             100 
                           
                            
                           
                             ( 
                             
                               1 
                               - 
                               N 
                             
                             ) 
                           
                         
                       
                       ) 
                     
                     × 
                     100 
                   
                 
               
               
                 
                   ( 
                   7 
                   ) 
                 
               
             
           
         
       
     
         [0082]    C. Basic Mathematics for Transformer Noise Level Calculation 
         [0083]    The transformer noise is Lp N1  at ON condition, and Lp N2  when all the fans are in operation at rated speed. The relationship between the noise Lp fan  caused by fans and the proportion of fans X is shown in equation (8): 
         [0000]        Lp   fan   =f ( X )   (8)
 
         [0084]    So when the proportion of fans in operation is X, the total noise from the transformer and the fan is: 
         [0000]    
       
         
           
             
               
                 
                   
                     f 
                     3 
                   
                   = 
                   
                     
                       Lp 
                       t 
                     
                     = 
                     
                       { 
                       
                         
                           
                             
                               
                                 Lp 
                                 
                                   N 
                                    
                                   
                                       
                                   
                                    
                                   1 
                                 
                               
                               , 
                             
                           
                           
                             
                               
                                 Lp 
                                 fan 
                               
                               = 
                               0 
                             
                           
                         
                         
                           
                             
                               
                                 
                                   Lp 
                                   
                                     N 
                                      
                                     
                                         
                                     
                                      
                                     1 
                                   
                                 
                                 + 
                                 
                                   10 
                                    
                                   
                                       
                                   
                                    
                                   
                                     lg 
                                      
                                     
                                       [ 
                                       
                                         1 
                                         + 
                                         
                                           10 
                                           
                                             - 
                                             
                                               
                                                 
                                                   Lp 
                                                   
                                                     N 
                                                      
                                                     
                                                         
                                                     
                                                      
                                                     1 
                                                   
                                                 
                                                 - 
                                                 
                                                   Lp 
                                                   fan 
                                                 
                                               
                                               10 
                                             
                                           
                                         
                                       
                                       ] 
                                     
                                   
                                 
                               
                               , 
                             
                           
                           
                             
                               
                                 Lp 
                                 
                                   N 
                                    
                                   
                                       
                                   
                                    
                                   1 
                                 
                               
                               &gt; 
                               
                                 Lp 
                                 fan 
                               
                             
                           
                         
                         
                           
                             
                               
                                 
                                   Lp 
                                   fan 
                                 
                                 + 
                                 
                                   10 
                                    
                                   
                                       
                                   
                                    
                                   
                                     lg 
                                      
                                     
                                       [ 
                                       
                                         1 
                                         + 
                                         
                                           10 
                                           
                                             - 
                                             
                                               
                                                 
                                                   Lp 
                                                   fan 
                                                 
                                                 - 
                                                 
                                                   Lp 
                                                   
                                                     N 
                                                      
                                                     
                                                         
                                                     
                                                      
                                                     1 
                                                   
                                                 
                                               
                                               10 
                                             
                                           
                                         
                                       
                                       ] 
                                     
                                   
                                 
                               
                               , 
                             
                           
                           
                             
                               
                                 Lp 
                                 fan 
                               
                               &gt; 
                               
                                 Lp 
                                 
                                   N 
                                    
                                   
                                       
                                   
                                    
                                   1 
                                 
                               
                             
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   9 
                   ) 
                 
               
             
           
         
       
     
         [0085]    Wherein, 
         [0086]    Lp fan : the fan noise; 
         [0087]    Lp N1 : the transformer noise. 
         [0088]    D. Objective Function with Weighting Factors 
         [0089]    When the cooling capacity varies, the variation of the loss f 1 , the top oil temperature f 2  and the noise f 3  are obviously different. In order to unify them, the maximum and minimum values of these three objectives f 1min , f 1max , f 2mim  f 2max , f 3min  and f 3max  are calculated at each moment and put into the objective function shown in (10). 
         [0090]    By using weighting factors w 1 , w 2 , w 3  for these three objectives, the objective function can be expressed as: 
         [0000]    
       
         
           
             
               
                 
                   
                     f 
                     obj 
                   
                   = 
                   
                     
                       
                         w 
                         1 
                       
                        
                       
                         
                           f 
                           1 
                         
                         
                           
                             f 
                             
                               1 
                                
                               
                                   
                               
                                
                               max 
                             
                           
                           - 
                           
                             f 
                             
                               1 
                                
                               
                                   
                               
                                
                               min 
                             
                           
                         
                       
                     
                     + 
                     
                       
                         w 
                         2 
                       
                        
                       
                         
                           f 
                           2 
                         
                         
                           
                             f 
                             
                               2 
                                
                               
                                   
                               
                                
                               max 
                             
                           
                           - 
                           
                             f 
                             
                               2 
                                
                               
                                   
                               
                                
                               min 
                             
                           
                         
                       
                     
                     + 
                     
                       
                         w 
                         3 
                       
                        
                       
                         
                           f 
                           3 
                         
                         
                           
                             f 
                             
                               3 
                                
                               
                                   
                               
                                
                               max 
                             
                           
                           - 
                           
                             f 
                             
                               3 
                                
                               
                                   
                               
                                
                               min 
                             
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   10 
                   ) 
                 
               
             
           
         
       
     
         [0000]    where, w 1 +w 2 +w 3 =1 
         [0091]    With formula (10), the optimal cooling capacity for all three objectives can be calculated. Also, each of objectives can be met individually when set its weight to 1, and set other weights to 0. 
         [0092]    Step 3: the third step, after the control commands calculation, the central controller will execute the results by controlling the switches directly or sending the control command to VFD, as shown in  FIG. 8 , where the control actions includes the start and stop of fans, controllable switch operation, and VFD frequency regulation. 
         [0093]    To start the fan, the central controller switches the motor-fan which does not need VFD directly to AC lines. For the motor-fan chain will be driven by VFD, the control center switches it to VFD, and sends the speed regulation reference to VFD. 
         [0094]    To stop the fan, the central controller directly switches the motor-fan chains off-line. 
         [0095]    The central controller repeats the Step 2 and Step 3in real-time. 
         [0096]    Advantages of the method and system according to this invention: 
         [0097]    This invention proposes a novel transformer cooling system and the corresponding operation method for optimal temperature control, which can improve the operation efficiency of the whole transformer with very limited capital investment on cooling system hardware upgrade, and meanwhile to extend the transformer lifecycle and lower the noise level of the transformer system. 
         [0098]    In this invention, the motor-fan loads of the cooling system will be controlled by one VFD selectively according to the temperature control requirement. For motor-fan loads needs to operate at rated power, they will connect to the AC bus directly. The temperature control will consider efficiency of the transformer windings and the cooling system together. Meanwhile, transformer top-oil temperature variation will be controlled in an coordinated way to extend the lifecycle. Furthermore, transformer noise level will be considered together in the cooling control in order to minimize the impact on the surrounding environment. With the proposed electrical design and the control method, the cooling system can be operated in an optimal way to achieve cost-effective efficiency improvement of the whole transformer. 
         [0099]    Though the present invention has been described on the basis of some preferred embodiments, those skilled in the art should appreciate that those embodiments should by no means limit the scope of the present invention. Without departing from the spirit and concept of the present invention, any variations and modifications to the embodiments should be within the apprehension of those with ordinary knowledge and skills in the art, and therefore fall in the scope of the present invention which is defined by the accompanied claims.