Abstract:
A fault detector for operating process of electric arc furnace and method thereof are disclosed, which belong to the technical field of fault detection. The fault detector includes the temperature signal acquisition equipment, the current signal acquisition equipment and a host computer. The multimode fault monitoring and diagnosis method comprises the following steps: acquiring and standardizing the data; establishing the preliminary monitoring model for the operating process to obtain common subsets of M operating modes and typical subsets of every operating mode; calculating the T 2  statistics and the SPE statistics, and monitoring and diagnosing fault in the operating process. The present invention has the advantages that the colorimetric temperature measurement can improve the calculation accuracy, different equipment becomes compatible, and the fault detector is suitable for operating in industrial production process with a variety of steady modes and can diagnose faults in a certain operating mode.

Description:
BACKGROUND OF THE INVENTION 
       [0001]    1. Field of the Invention 
         [0002]    The present invention relates to a fault detector for operating process of an electric arc furnace and a method thereof, which belongs to the technical field of fault detection. 
         [0003]    2. The Prior Arts 
         [0004]    Electric arc furnace operates by three-phrase AC power, and is used to melt ore and metal with the electric arc between the electrodes and the molten materials. In the operating process, the electric arc furnace produces lots of waste gas which is exhausted by raising and lowering the electrodes. If the electrodes are not adjusted in time, a mass of gas can accumulate in the electric arc furnace, which can result in voids in the products. If the backlog gas is excessive, the molten material in the furnace may spray. Moreover, the electrodes position appears excursion during adjustment, and distribution of the materials in the electric arc furnace is not homogeneous, so the temperature around electrodes is high, however the temperature around other electrodes is low. The imbalance distribution of heat would influence the quality of products. More serious, if the imbalance distribution of heat exists for a long time, in the higher temperature portion of the electric arc furnace, the temperature rises sharply. And the furnace could be even worn out without effective detection and prevention. 
         [0005]    Fault detection and diagnosis have been the key means to ensure the safety of equipment and personnel and the quality of products. In the past several decades, Principal Component Analysis (PCA) and Partial Least Square (PLS) have been used widely in industrial production and got lots of satisfying results. Multivariate Statistics Process Control (MSPC) has also been studied deeply. However, there is often more than one operating mode in the electric arc furnace, so fault detection of the multiple modes of production and multimode productive process becomes more complicated. There is often a lot of problems in the multivariate statistical process control (MSPC) applied to the multimode production process. MSPM may sometimes make some mistakes, even if the production process is working in another steady-state operation mode. The recursive PCA and self-adaptive PLS arise over these years. Although these methods can also be used for online fault detection, it still lack the ability to deal with the multimode operating process, and the detection results are unsatisfactory. 
         [0006]    Moreover, because when the arc forms the energy is concentrated and the temperature of the arc zone near the electrodes could be above 3000° C., it is hard to measure. Temperature can be got by observing the temperature of the furnace wall in the traditional production process. The temperature of the furnace wall and that of the molten materials inside the electric arc furnace are usually different because there are specific heat capacity and thermal resistivity. But the temperature in the electric arc furnace is hard to measure. The traditional thermocouple measures temperature by contact, but the thermocouple will melt at high temperature. Therefore, the costs of temperature measurement increase greatly with the raising of the temperature. In this case, some engineers and technicians try to invent advanced thermocouples which can be used at higher temperature by changing the materials. Although the range of temperature which can be measured is widen, due to the restriction of the contact temperature measurement technology, the equipment is worn out seriously. Therefore, the durability of the equipment is affected, and more impurities are generated because the materials are molten in the electric arc furnace, which affects the quality of the products. 
       SUMMARY OF THE INVENTION 
       [0007]    In consideration of the shortage of the prior art, the present invention provides a detector of the production process with an electric arc furnace and a method thereof. 
         [0008]    The detector comprises the temperature signal acquisition equipment, the current signal acquisition equipment and a host computer. The current signal acquisition equipment comprises an AC transformer and a Programmable Logic Controller (PLC). The output terminal of the AC transformer is connected to the input terminal of the PLC, and the output terminal of PLC is connected to the host computer. The temperature signal acquisition equipment comprises an imaging device, a color sensor, two transceiver integrated circuits and a serial communication circuit. The imaging device comprises an objective lens and an eyepiece. The transceiver integrated circuit comprises a microcontroller and a RF transceiver. One RF transceiver of the two transceiver integrated circuit works as transmitter, and the other works as an acceptor. The input terminal of the color sensor is connected to the imaging device. In the first group transceiver integrated circuit the output terminal of the color sensor is connected to the input terminal of the microcontroller, and the output terminal of the microcontroller is connected to the input terminal of the RF transmitter. The output terminal of the RF transmitter is connected to the input terminal of the RF receiver in the second group transceiver integrated circuit by wireless network. The output terminal of the RF receiver is connected to the input terminal of the microcontroller in the second group of transceiver integrated circuit. The output terminal of the microcontroller is connected to the input terminal of the serial communication circuit. The output terminal of the serial communication circuit is connected to the host computer. The objective lens of the imaging device aims the measured target. The photosensitive sensor is fixed in the eyepiece at one end, so that the measured target radiation light can be projected to the photosensitive parts of the color sensor by the imaging device. The color sensor is connected to the microcontroller. When collecting the signal, the microcontroller controls the RF transmitter circuit to transmit wireless RF signals. The RF receiver circuit controls the serial communication circuit to send data to the host computer under the coordination of the microcontroller after receiving signals. Two identical transceiver integrated circuits are used to transmit and receive wireless signals and to play a different role in the control of the microcontroller. The configuration software of the host computer has the function of real-time display, history enquiry, data storage and alarm processing, etc. 
         [0009]    Field data are transmitted to the host computer by two acquisition equipments. The host computer calculates and analyzes the field data by the method in the present invention to detect faults. 
         [0010]    The function of each component of the fault detector is as follows: 
         [0011]    (1) Current Signal Acquisition equipment: The current signal acquisition equipment comprises an AC transformer and a Programmable Logic Controller (PLC). The electric arc furnace has the characteristics of low voltage and high current, so the voltage of the wire between the transformer and the electric arc furnace is not high, but the current is very high. When the voltage is 100 V in the operating process, the current can be up to 10000 A, which is hard to measure. Therefore, the AC transformer is used to transform the high current in the electrodes to the normal signal of 0˜5 V. The AC transformer is connected to the input terminal of the PLC. After the PLC completes the analogue-to-digital conversion, current signals in the three electrodes are acquired by the PLC simultaneously, and then, the current signals are transmitted to the host computer by the serial ports of the PLC. 
         [0012]    (2) Temperature Signal Acquisition equipment: The temperature signal acquisition equipment comprises an imaging device, a color sensor, a microcontroller, transceiver integrated circuit and a serial communication circuit. The electromagnetic spectrum radiated by the high-temperature molten materials contains electromagnetic waves with all the wavelength frequency from infrared rays to ultraviolet rays. The electromagnetic waves with different wavelengths have different radiation intensity in different temperatures, which also contains the visible part. Two kinds of visible light with different wavelengths, such as red light and green light, have different ratios of the radiation intensity in different temperatures. Therefore, temperature can be measured by indirectly measuring the ratio of the radiation intensity of two kinds of visible light with different wavelengths. The temperature signal acquisition equipment is used to make the measured target radiation light projected to the photosensitive parts of the color sensor by the imaging device. The color sensor measures the radiation intensity of two kinds of visible light such as red light and green light with the help of the microcontroller. After the data are transmitted from the color sensor to the microcontroller, the transmitter circuit of the microcontroller transmits the data in the form of wireless signals. After the signals are received by the receiver circuit, the receiver circuit is controlled by the microcontroller to transmit the signals to the host computer through the serial communication circuit. The program to control the color sensor and to control the transmission of wireless data is written in the microcontroller at the sending end. The program to control serial communication and to receive wireless data is written in the microcontroller at the receiving end. 
         [0013]    (3) Host Computer: After the data are transmitted to the host computer, firstly, the monitoring software transforms the data acquired by the color sensor to temperature data by the method of the colorimetric temperature measurement, and then, the current data and the temperature data are written into the database. At the same time, the current data and the temperature data as well as the history temperature curves thereof are displayed on the monitor interface. The host computer also detects faults of the operating process by analyzing the history data. 
         [0014]    The host computer needs to calculate the data to obtain the related temperature after the host computer receives the signals of radiation intensity of the two kinds of visible light such as red light and green light, and the method of the colorimetric temperature measurement is as follows: 
         [0015]    Based on Planck&#39;s law of radiation: 
         [0000]        M (μ, T )=ε(μ, T ) c   1 ·μ −5 /( e   c     2     /(μ·T) −1)  (1)
 
         [0016]    In the equation (1), μ is the wavelength of the electromagnetic wave, M(μ,T) is the spectral radiant emittance, ε(μ,T) is the gray body emissivity, T is degree Kelvin, c 1  is the first thermal radiation constant, c 2  is the second thermal radiation constant, and e is Natural logarithm. 
         [0017]    When c 2 /(μ·T)&gt;&gt;1, Planck&#39;s law can be superseded by Wien&#39;s law, then 
         [0000]        M (μ, T )=ε(μ, T ) c   1 ·μ −5   ·e   −c     2     /(μ·T)   (2)
 
         [0018]    When two wavelengths μ 1 , μ 2  are given, the corresponding spectral radiant emittances are: 
         [0000]        M (μ 1   ,T )=ε(μ 1   T ) c   1 ·μ 1   −5   ·e   −c     2     /(μ     1     ·T)   (3)
 
         [0000]        M (μ 2   ,T )=ε(μ 2   T ) c   1 ·μ 2   −5   ·e   −c     2     /(μ     2     ·T)   (4)
 
         [0019]    The temperature measurement equation can be obtained by comparing the equations (3) and (4) above: 
         [0000]    
       
         
           
             
               
                 
                   T 
                   = 
                   
                     
                       
                         c 
                         2 
                       
                       · 
                       
                         ( 
                         
                           
                             1 
                             / 
                             
                               μ 
                               2 
                             
                           
                           - 
                           
                             1 
                             / 
                             
                               μ 
                               1 
                             
                           
                         
                         ) 
                       
                     
                     
                       
                         ln 
                          
                         
                           ( 
                           
                             
                               M 
                                
                               
                                 ( 
                                 
                                   
                                     μ 
                                     1 
                                   
                                   , 
                                   T 
                                 
                                 ) 
                               
                             
                             / 
                             
                               M 
                                
                               
                                 ( 
                                 
                                   
                                     μ 
                                     2 
                                   
                                   , 
                                   T 
                                 
                                 ) 
                               
                             
                           
                           ) 
                         
                       
                       - 
                       
                         ln 
                          
                         
                           ( 
                           
                             
                               ɛ 
                               1 
                             
                             / 
                             
                               ɛ 
                               2 
                             
                           
                           ) 
                         
                       
                       - 
                       
                         5 
                          
                         
                             
                         
                          
                         
                           ln 
                            
                           
                             ( 
                             
                               
                                 μ 
                                 2 
                               
                               / 
                               
                                 μ 
                                 1 
                               
                             
                             ) 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   5 
                   ) 
                 
               
             
           
         
       
     
         [0020]    The spectral radiant emittances of different wavelengths gathered by red filter and green filter are used as the input of equation (5). Because the wavelengths of two colors are approximate, the ratio of the gray body emissivity is approximately 1, so that the equation can be simplified as: 
         [0000]    
       
         
           
             
               
                 
                   T 
                   = 
                   
                     
                       
                         c 
                         2 
                       
                       · 
                       
                         ( 
                         
                           
                             1 
                             / 
                             
                               μ 
                               2 
                             
                           
                           - 
                           
                             1 
                             / 
                             
                               μ 
                               1 
                             
                           
                         
                         ) 
                       
                     
                     
                       
                         ln 
                          
                         
                           ( 
                           
                             
                               M 
                                
                               
                                 ( 
                                 
                                   
                                     μ 
                                     1 
                                   
                                   , 
                                   T 
                                 
                                 ) 
                               
                             
                             / 
                             
                               M 
                                
                               
                                 ( 
                                 
                                   
                                     μ 
                                     2 
                                   
                                   , 
                                   T 
                                 
                                 ) 
                               
                             
                           
                           ) 
                         
                       
                       - 
                       
                         5 
                          
                         
                             
                         
                          
                         
                           ln 
                            
                           
                             ( 
                             
                               
                                 μ 
                                 2 
                               
                               / 
                               
                                 μ 
                                 1 
                               
                             
                             ) 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   6 
                   ) 
                 
               
             
           
         
       
     
         [0021]    The electric arc furnace has multiple modes in the operating process. There always exits a stable working state in every mode. The current data and the temperature data will change when the work mode switches, leading to omission and misstatement of fault diagnosis. The present invention provides a multimode fault diagnosis method which can intelligently recognize the current operating mode when the operating mode changes and diagnose the fault for this mode. 
         [0022]    The fault detection device of the electric arc furnace is used to take multimode fault monitoring and diagnosis, and it comprises the following steps: 
         [0023]    Step 1: Acquiring and standardizing data. M operating modes are supposed in the operating process of the electric arc furnace. Three-phase current I A , I B  and I C  in the M operating modes and temperature t in the electric arc furnace are acquired when the electric arc furnace operates steadily in different modes. Then, the initial dataset {dot over (X)} m =[x 1,m   T , x 2,m   T , . . . x l,m   T  . . . , x N     m     T ] T (N m ×J) (m=1, 2, . . . , M) is obtained, where the subscript m is the m th  operating mode, N m  is the number of the data sample in the m th  operating mode, and J is the number of the variables which is defined as 4 in the present invention. To explain the calculating method clearly, J is still used. The initial dataset is standardized by the methods of calculating the average and standard deviation, and then, the standard dataset X m =(N m ×J) (m=1, 2, . . . , M) in the normal operation state can be obtained. 
         [0024]    Step 2: Decompose the data and establish the preliminary monitoring model for the operating process of the electric arc furnace. The M datasets are obtained in Step 1 by standardized. And then build the preliminary monitoring model for the operating process of the electric arc furnace, that calculate the global main factor matrix P g  under the M different operating modes of the monitoring model to decompose the datasets. The datasets is obtained in the M different operating modes of the electric arc furnace. Then close common subsets in these M operating modes and special subsets of every operating mode are obtained. 
         [0025]    To establish the preliminary monitoring model for the operating process of the electric arc furnace, the modeling has the following two steps of: 
         [0026]    1. Calculating the global main factor matrix P g  for the operating process of the electric arc furnace. 
         [0027]    The datasets in different operating modes of the electric arc furnace have the same variable correlation. The linear combination of the global main factor can represent other samples and can even replace all the other samples. The global main factors can also represent the main characteristic of the initial data, the differences of data distribution and the change information of themselves. Therefore, the global main factors for the operating process of the electric arc furnace can be used as the similarity evaluation index of datasets for different operating modes. 
         [0028]    In the operating process of the electric arc furnace, the main factor of dataset of every operating mode are defined as p m  (m=1, 2, . . . , M), where the subscript m represents the number of different operating modes. The relation between the main factor p m  of dataset of every operating mode and this dataset is: 
         [0000]    
       
         
           
             
               
                 
                   
                     p 
                     m 
                   
                   = 
                   
                     
                       
                         ∑ 
                         
                           n 
                           = 
                           1 
                         
                         
                           N 
                           m 
                         
                       
                        
                       
                           
                       
                        
                       
                         
                           a 
                           
                             n 
                             , 
                             m 
                           
                         
                          
                         
                           x 
                           
                             n 
                             , 
                             m 
                           
                         
                       
                     
                     = 
                     
                       
                         X 
                         m 
                         T 
                       
                        
                       
                         a 
                         m 
                       
                     
                   
                 
               
               
                 
                   ( 
                   7 
                   ) 
                 
               
             
           
         
       
     
         [0029]    α m =[a 1,m , a 2,m , . . . , a n,m ] T  is the linear combination coefficient, where the letter n represents the number of rows of α m , i.e., p m  can be represented as the linear combination of the N m  vectors of X m . Then, the global main factor for the operating process of the electric arc furnace p g  is introduced, and make it approximates to the main factor of the dataset obtained in each operating mode of the electric arc furnace. Global main factor p g  is used to describe the main factor of the dataset of each operating mode. The method to obtain p g  is as follows: 
         [0030]    To make p g  sufficiently approximate to p m,j  which is the main factor of the subset for each operating mode of the electric arc furnace, (p g   T X m   T α m ) 2  is maximized, and R 2 =(p g   T X 1   T α 1 ) 2 +(p g   T X 2   T α 2 ) 2 + . . . +(p s   T X M   T α M ) 2  is introduced. Then, the problem changes to calculating the maximum of R 2 . At the same time, the inequality 
         [0000]    
       
         
           
             
               ɛ 
               ≤ 
               
                 
                   
                     ( 
                     
                       
                         p 
                         g 
                         T 
                       
                        
                       
                         X 
                         1 
                         T 
                       
                        
                       
                         a 
                         1 
                       
                     
                     ) 
                   
                   2 
                 
                 
                   
                     ( 
                     
                       
                         p 
                         g 
                         T 
                       
                        
                       
                         X 
                         2 
                         T 
                       
                        
                       
                         a 
                         2 
                       
                     
                     ) 
                   
                   2 
                 
               
             
             = 
             
               
                 
                   
                     ( 
                     
                       
                         p 
                         g 
                         T 
                       
                        
                       
                         X 
                         2 
                         T 
                       
                        
                       
                         a 
                         2 
                       
                     
                     ) 
                   
                   2 
                 
                 
                   
                     ( 
                     
                       
                         p 
                         g 
                         T 
                       
                        
                       
                         X 
                         3 
                         T 
                       
                        
                       
                         a 
                         3 
                       
                     
                     ) 
                   
                   2 
                 
               
               ≤ 
               
                 1 
                 ɛ 
               
             
           
         
       
     
         [0000]    is introduced, where ε is a constraint parameter and satisfies the condition of 0.8≦ε&lt;1. To make p g  sufficiently approximate to p m,j , s is defined as 0.9 in the present invention. The problem changes to solve the following equation: 
         [0000]    
       
         
           
             
               
                 
                   
                     max 
                      
                     
                         
                     
                      
                     
                       R 
                       2 
                     
                   
                   = 
                   
                     
                       max 
                        
                       
                         ( 
                         
                           
                             ∑ 
                             
                               m 
                               = 
                               1 
                             
                             M 
                           
                            
                           
                               
                           
                            
                           
                             
                               ( 
                               
                                 
                                   p 
                                   g 
                                   T 
                                 
                                  
                                 
                                   X 
                                   m 
                                   T 
                                 
                                  
                                 
                                   a 
                                   m 
                                 
                               
                               ) 
                             
                             2 
                           
                         
                         ) 
                       
                     
                     . 
                   
                 
               
               
                 
                   ( 
                   8 
                   ) 
                 
               
             
           
         
       
     
         [0000]    Satisfy constraints are as follows: 
         [0000]    
       
         
           
             
               
                 
                   s 
                   . 
                   t 
                   . 
                   
                     { 
                     
                       
                         
                           
                             
                               
                                 
                                   ɛ 
                                    
                                   
                                     ( 
                                     
                                       
                                         p 
                                         g 
                                         T 
                                       
                                        
                                       
                                         X 
                                         
                                           m 
                                           - 
                                           1 
                                         
                                         T 
                                       
                                        
                                       
                                         a 
                                         
                                           m 
                                           - 
                                           1 
                                         
                                       
                                     
                                     ) 
                                   
                                 
                                 2 
                               
                               - 
                               
                                 
                                   ( 
                                   
                                     
                                       p 
                                       g 
                                       T 
                                     
                                      
                                     
                                       X 
                                       m 
                                       T 
                                     
                                      
                                     
                                       a 
                                       m 
                                     
                                   
                                   ) 
                                 
                                 2 
                               
                             
                             ≤ 
                             0 
                           
                         
                       
                       
                         
                           
                             
                               
                                 
                                   ɛ 
                                    
                                   
                                     ( 
                                     
                                       
                                         p 
                                         g 
                                         T 
                                       
                                        
                                       
                                         X 
                                         m 
                                         T 
                                       
                                        
                                       
                                         a 
                                         m 
                                       
                                     
                                     ) 
                                   
                                 
                                 2 
                               
                               - 
                               
                                 
                                   ( 
                                   
                                     
                                       p 
                                       g 
                                       T 
                                     
                                      
                                     
                                       X 
                                       
                                         m 
                                         - 
                                         1 
                                       
                                       T 
                                     
                                      
                                     
                                       a 
                                       
                                         m 
                                         - 
                                         1 
                                       
                                     
                                   
                                   ) 
                                 
                                 2 
                               
                             
                             ≤ 
                             0 
                           
                         
                       
                       
                         
                           
                             
                               
                                 p 
                                 g 
                                 T 
                               
                                
                               
                                 p 
                                 g 
                               
                             
                             = 
                             1 
                           
                         
                       
                       
                         
                           
                             
                               
                                 a 
                                 m 
                                 T 
                               
                                
                               
                                 a 
                                 m 
                               
                             
                             = 
                             1. 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   9 
                   ) 
                 
               
             
           
         
       
     
         [0031]    The present invention deals with fault diagnosis for the different operating modes of the electric arc furnace. Different mode is reflected in the different raw materials. Raw materials will be adjusted due to production needs. When raw material is different, the current for normal production is also different, but it still operates in a stable state. To keep the generality, suppose there are two kinds of production modes corresponding to the two different kinds of raw materials A and B. Then, mode A and mode B respectively represent the operating modes of the electric arc furnace, of which the raw material is material A and material B respectively. That is to say M is equal to 2. To achieve the above objective, an initial objective function is introduced, and the maximum of the objective function is calculated as follows: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       F 
                        
                       
                         ( 
                         
                           
                             p 
                             g 
                           
                           , 
                           a 
                           , 
                           λ 
                         
                         ) 
                       
                     
                     = 
                     
                       
                         
                           ( 
                           
                             
                               p 
                               g 
                               T 
                             
                              
                             
                               X 
                               A 
                               T 
                             
                              
                             
                               a 
                               A 
                             
                           
                           ) 
                         
                         2 
                       
                       + 
                       
                         
                           ( 
                           
                             
                               p 
                               g 
                               T 
                             
                              
                             
                               X 
                               B 
                               T 
                             
                              
                             
                               a 
                               B 
                             
                           
                           ) 
                         
                         2 
                       
                       - 
                       
                         
                           λ 
                           g 
                         
                          
                         
                           ( 
                           
                             
                               
                                 p 
                                 g 
                                 T 
                               
                                
                               
                                 p 
                                 g 
                               
                             
                             - 
                             1 
                           
                           ) 
                         
                       
                       - 
                       
                         
                           λ 
                           A 
                         
                          
                         
                           ( 
                           
                             
                               
                                 a 
                                 A 
                                 T 
                               
                                
                               
                                 a 
                                 A 
                               
                             
                             - 
                             1 
                           
                           ) 
                         
                       
                       - 
                       
                         
                           λ 
                           B 
                         
                          
                         
                           ( 
                           
                             
                               
                                 a 
                                 B 
                                 T 
                               
                                
                               
                                 a 
                                 B 
                               
                             
                             - 
                             1 
                           
                           ) 
                         
                       
                       - 
                       
                         
                           λ 
                           1 
                         
                          
                         
                           ( 
                           
                             
                               
                                 ɛ 
                                  
                                 
                                   ( 
                                   
                                     
                                       p 
                                       g 
                                       T 
                                     
                                      
                                     
                                       X 
                                       A 
                                       T 
                                     
                                      
                                     
                                       a 
                                       A 
                                     
                                   
                                   ) 
                                 
                               
                               2 
                             
                             - 
                             
                               
                                 ( 
                                 
                                   
                                     p 
                                     g 
                                     T 
                                   
                                    
                                   
                                     X 
                                     B 
                                     T 
                                   
                                    
                                   
                                     a 
                                     B 
                                   
                                 
                                 ) 
                               
                               2 
                             
                           
                           ) 
                         
                       
                       - 
                       
                         
                           λ 
                           2 
                         
                          
                         
                           ( 
                           
                             
                               
                                 ɛ 
                                  
                                 
                                   ( 
                                   
                                     
                                       p 
                                       g 
                                       T 
                                     
                                      
                                     
                                       X 
                                       B 
                                       T 
                                     
                                      
                                     
                                       a 
                                       B 
                                     
                                   
                                   ) 
                                 
                               
                               2 
                             
                             - 
                             
                               
                                 ( 
                                 
                                   
                                     p 
                                     g 
                                     T 
                                   
                                    
                                   
                                     X 
                                     A 
                                     T 
                                   
                                    
                                   
                                     a 
                                     A 
                                   
                                 
                                 ) 
                               
                               2 
                             
                           
                           ) 
                         
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   10 
                   ) 
                 
               
             
           
         
       
     
         [0000]    Where λ g , λ A , λ B , λ 1  and λ 2  are constants. The partial derivative of p g , α A  and α B  for F(p g , α, λ) in above equation is solved to obtain the following equations respectively: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         ( 
                         
                           1 
                           - 
                           
                             
                               λ 
                               1 
                             
                              
                             ɛ 
                           
                           + 
                           
                             λ 
                             2 
                           
                         
                         ) 
                       
                        
                       
                         X 
                         A 
                         T 
                       
                        
                       
                         X 
                         A 
                       
                        
                       
                         p 
                         g 
                       
                     
                     + 
                     
                       
                         ( 
                         
                           1 
                           - 
                           
                             
                               λ 
                               2 
                             
                              
                             ɛ 
                           
                           + 
                           
                             λ 
                             1 
                           
                         
                         ) 
                       
                        
                       
                         X 
                         B 
                         T 
                       
                        
                       
                         X 
                         B 
                       
                        
                       
                         p 
                         g 
                       
                     
                   
                   = 
                   
                     
                       λ 
                       g 
                     
                      
                     
                       p 
                       g 
                     
                   
                 
               
               
                 
                   ( 
                   11 
                   ) 
                 
               
             
             
               
                 
                   
                     a 
                     A 
                   
                   = 
                   
                     
                       
                         
                           
                             ( 
                             
                               1 
                               - 
                               
                                 
                                   λ 
                                   1 
                                 
                                  
                                 ɛ 
                               
                               + 
                               
                                 λ 
                                 2 
                               
                             
                             ) 
                           
                           
                             λ 
                             A 
                           
                         
                       
                        
                       
                         X 
                         A 
                       
                        
                       
                         p 
                         g 
                       
                     
                     = 
                     
                       
                         
                           1 
                           
                             
                               
                                 p 
                                 g 
                                 T 
                               
                                
                               
                                 ( 
                                 
                                   
                                     X 
                                     A 
                                     T 
                                   
                                    
                                   
                                     X 
                                     A 
                                   
                                 
                                 ) 
                               
                             
                              
                             
                               p 
                               g 
                             
                           
                         
                       
                        
                       
                         X 
                         A 
                       
                        
                       
                         p 
                         g 
                       
                     
                   
                 
               
               
                 
                   ( 
                   12 
                   ) 
                 
               
             
             
               
                 
                   
                     a 
                     B 
                   
                   = 
                   
                     
                       
                         
                           
                             ( 
                             
                               1 
                               - 
                               
                                 
                                   λ 
                                   2 
                                 
                                  
                                 ɛ 
                               
                               + 
                               
                                 λ 
                                 1 
                               
                             
                             ) 
                           
                           
                             λ 
                             B 
                           
                         
                       
                        
                       
                         X 
                         B 
                       
                        
                       
                         p 
                         g 
                       
                     
                     = 
                     
                       
                         
                           1 
                           
                             
                               
                                 p 
                                 g 
                                 T 
                               
                                
                               
                                 ( 
                                 
                                   
                                     X 
                                     B 
                                     T 
                                   
                                    
                                   
                                     X 
                                     B 
                                   
                                 
                                 ) 
                               
                             
                              
                             
                               p 
                               g 
                             
                           
                         
                       
                        
                       
                         X 
                         B 
                       
                        
                       
                         p 
                         g 
                       
                     
                   
                 
               
               
                 
                   ( 
                   13 
                   ) 
                 
               
             
           
         
       
     
         [0000]    λ A  and λ B  can be calculated from equations (11), (12) and (13): 
         [0000]      λ A =(1−λ 1 ε+λ 2 ) p   g   T ( X   A   T   X   A ) p   g   (14)
 
         [0000]      λ B =(1−λ 2 ε+λ 1 ) p   g   T ( X   B   T   X   B ) p   g   (15)
 
         [0032]    The equation (11) is analyzed to calculate p g . Based on the previously mentioned inequality, p g  should meet the following conditions: 
         [0000]    
       
         
           
               
             
               { 
               
                 
                   
                     
                       
                         
                           
                             λ 
                             1 
                           
                            
                           
                             ( 
                             
                               
                                 
                                   ɛ 
                                    
                                   
                                     ( 
                                     
                                       
                                         p 
                                         g 
                                         T 
                                       
                                        
                                       
                                         X 
                                         A 
                                         T 
                                       
                                        
                                       
                                         a 
                                         A 
                                       
                                     
                                     ) 
                                   
                                 
                                 2 
                               
                               - 
                               
                                 
                                   ( 
                                   
                                     
                                       p 
                                       g 
                                       T 
                                     
                                      
                                     
                                       X 
                                       B 
                                       T 
                                     
                                      
                                     
                                       a 
                                       B 
                                     
                                   
                                   ) 
                                 
                                 2 
                               
                             
                             ) 
                           
                         
                         = 
                         0 
                       
                       , 
                       
                         ( 
                         16 
                         ) 
                       
                     
                   
                 
                 
                   
                     
                       
                         
                           
                             λ 
                             2 
                           
                            
                           
                             ( 
                             
                               
                                 
                                   ɛ 
                                    
                                   
                                     ( 
                                     
                                       
                                         p 
                                         g 
                                         T 
                                       
                                        
                                       
                                         X 
                                         B 
                                         T 
                                       
                                        
                                       
                                         a 
                                         B 
                                       
                                     
                                     ) 
                                   
                                 
                                 2 
                               
                               - 
                               
                                 
                                   ( 
                                   
                                     
                                       p 
                                       g 
                                       T 
                                     
                                      
                                     
                                       X 
                                       A 
                                       T 
                                     
                                      
                                     
                                       a 
                                       A 
                                     
                                   
                                   ) 
                                 
                                 2 
                               
                             
                             ) 
                           
                         
                         = 
                         0 
                       
                       , 
                       
                         ( 
                         17 
                         ) 
                       
                     
                   
                 
                 
                   
                     
                       
                         λ 
                         1 
                       
                       , 
                       
                         
                           λ 
                           2 
                         
                         ≥ 
                         0. 
                       
                     
                   
                 
               
             
           
         
       
     
         [0033]    Equations (12) and (13) are submitted into equations (16) and (17) to obtain the result as follows: 
         [0000]    
       
         
           
               
             
               { 
               
                 
                   
                     
                       
                         
                           
                             λ 
                             1 
                           
                            
                           
                             ( 
                             
                               
                                 ɛ 
                                  
                                 
                                     
                                 
                                  
                                 
                                   
                                     p 
                                     g 
                                     T 
                                   
                                    
                                   
                                     ( 
                                     
                                       
                                         X 
                                         A 
                                         T 
                                       
                                        
                                       
                                         X 
                                         A 
                                       
                                     
                                     ) 
                                   
                                 
                                  
                                 
                                   p 
                                   g 
                                 
                               
                               - 
                               
                                 
                                   
                                     p 
                                     g 
                                     T 
                                   
                                    
                                   
                                     ( 
                                     
                                       
                                         X 
                                         B 
                                         T 
                                       
                                        
                                       
                                         X 
                                         B 
                                       
                                     
                                     ) 
                                   
                                 
                                  
                                 
                                   p 
                                   g 
                                 
                               
                             
                             ) 
                           
                         
                         = 
                         0 
                       
                       , 
                       
                         ( 
                         18 
                         ) 
                       
                     
                   
                 
                 
                   
                     
                       
                         
                           
                             λ 
                             2 
                           
                            
                           
                             ( 
                             
                               
                                 ɛ 
                                  
                                 
                                     
                                 
                                  
                                 
                                   
                                     p 
                                     g 
                                     T 
                                   
                                    
                                   
                                     ( 
                                     
                                       
                                         X 
                                         B 
                                         T 
                                       
                                        
                                       
                                         X 
                                         B 
                                       
                                     
                                     ) 
                                   
                                 
                                  
                                 
                                   p 
                                   g 
                                 
                               
                               - 
                               
                                 
                                   
                                     p 
                                     g 
                                     T 
                                   
                                    
                                   
                                     ( 
                                     
                                       
                                         X 
                                         A 
                                         T 
                                       
                                        
                                       
                                         X 
                                         A 
                                       
                                     
                                     ) 
                                   
                                 
                                  
                                 
                                   p 
                                   g 
                                 
                               
                             
                             ) 
                           
                         
                         = 
                         0 
                       
                       , 
                       
                         ( 
                         19 
                         ) 
                       
                     
                   
                 
                 
                   
                     
                       
                         λ 
                         1 
                       
                       , 
                       
                         
                           λ 
                           2 
                         
                         ≥ 
                         0. 
                       
                     
                   
                 
               
             
           
         
       
     
         [0034]    In above equation, in order to calculate p g , the values of λ 1  and λ 2  are analyzed as follows: 
         [0035]    (1) If λ 1 ≠0 and λ 2 ≠0, the equations (18) and (19) are contradictory. Therefore, the equation has no solution in this case. 
         [0036]    (2) If λ 1 =0 and λ 2 =0, the equations (18) and (19) are constantly equal to 0, i.e., equations (18) and (19) lose binding effect, and equation (11) can be rewritten as follows: 
         [0000]        X   A   T   X   A   p   g   +X   B   T   X   B   p   g =λ g   p   g   (20)
 
         [0037]    From equation (20), calculating p g  changes to calculating the eigenvector of (X A   T X A +X B   T X B ), by which the maximum of (p g   T X A   T α A ) 2 +(p g   T X B   T α B ) 2  can be obtained. However, because equations (18) and (19) lose binding effect, p g  may not sufficiently approximate to the main factor of the subset for every mode. Therefore, the solution p g  cannot meet requirements. 
         [0038]    (3) If λ 1 =0 and λ 2 ≠0, from equations (18) and (19), the following equation is obtained: 
         [0000]      ε p   g   T   X   B   T   X   B   p   g   −p   g   T   X   A   T   X   A   p   g =0  (21)
 
         [0039]    Where the parameters should meet the condition of 0.8≦ε&lt;1. Here, s is defined as 0.9. Equation (21) has only one unknown vector p g , and some groups of solution obtained by solving the equation (21) are as the solution of p g , which is defined as p g1 . 
         [0040]    (4) If λ 1 ≠0 and λ 2 =0, the same method is used, and the following equation is obtained: 
         [0000]      ε p   g   T   X   A   T   X   A   p   g   −p   g   T   X   B   T   X   B   p   g =0  (22)
 
         [0000]    The analysis method similar to the method in step (3) is used to obtain the solution of another group of solution of p g , which is defined as p g2 . 
         [0041]    P g1  and p g2  can be the solutions of p g  which meet the conditions. P g1  and p g2  are combined together to constitute a matrix p g =[p g1 , p g2 ] of the R row, where the sum of the rows of P g1  and p g2  is R. Then, the global main factor matrix P g  (J×R) for the operating process of the electric arc furnace is obtained as the solution which meets the assuming conditions. 
         [0042]    2. Decomposed the dataset in the two working modes of the electric arc furnace working process. 
         [0043]    The global main factor matrix P g  (J×R) of the electric arc furnace operating process is obtained by the modeling above. The column vectors of P g  should meet the linear conditions p g =X A   T α 1  and p g =X B   T α 2 , where α 1  and α 2  are the linear combination coefficients, p g   T p g =1, α 1   T α 1 =1, and α 2   T α 2 =1. Then, the equations above can be transformed to α 1 =X A P g  and α 2 =X B P g . R groups of linear combination coefficients α 1  and α 2  are made respectively to constitute the score matrices T A   C  and T B   C  of the two common subsets, and the following equation can be obtained: 
         [0000]        T   A   C   =X   A   P   g   ,T   B   C   =X   B   P   g   (23)
 
         [0044]    By the following calculation, the dataset X A  for mode A of the electric arc furnace can be decomposed into the common dataset X A   C  containing common information and the special dataset X A   S . 
         [0000]    
       
         
           
             
               
                 
                   { 
                   
                     
                       
                         
                           
                             X 
                             A 
                           
                           = 
                           
                             
                               X 
                               A 
                               C 
                             
                             + 
                             
                               X 
                               A 
                               S 
                             
                           
                         
                       
                     
                     
                       
                         
                           
                             X 
                             A 
                             C 
                           
                           = 
                           
                             
                               
                                 T 
                                 A 
                                 C 
                               
                                
                               
                                 P 
                                 g 
                                 T 
                               
                             
                             = 
                             
                               
                                 X 
                                 A 
                               
                                
                               
                                 P 
                                 g 
                               
                                
                               
                                 P 
                                 g 
                                 T 
                               
                             
                           
                         
                       
                     
                     
                       
                         
                           
                             X 
                             A 
                             S 
                           
                           = 
                           
                             
                               
                                 X 
                                 A 
                               
                               - 
                               
                                 X 
                                 A 
                                 C 
                               
                             
                             = 
                             
                               
                                 
                                   X 
                                   A 
                                   T 
                                 
                                  
                                 
                                   ( 
                                   
                                     I 
                                     - 
                                     
                                       
                                         P 
                                         g 
                                       
                                        
                                       
                                         P 
                                         g 
                                         T 
                                       
                                     
                                   
                                   ) 
                                 
                               
                               . 
                             
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   24 
                   ) 
                 
               
             
           
         
       
     
         [0045]    The dataset X B  for mode B of the electric arc furnace can be decomposed into the common dataset X B   C  and special dataset X B   S  in the same method. 
         [0000]    
       
         
           
             
               
                 
                   { 
                   
                     
                       
                         
                           
                             X 
                             B 
                           
                           = 
                           
                             
                               X 
                               B 
                               C 
                             
                             + 
                             
                               X 
                               B 
                               S 
                             
                           
                         
                       
                     
                     
                       
                         
                           
                             X 
                             B 
                             C 
                           
                           = 
                           
                             
                               
                                 T 
                                 B 
                                 C 
                               
                                
                               
                                 P 
                                 g 
                                 T 
                               
                             
                             = 
                             
                               
                                 X 
                                 B 
                               
                                
                               
                                 P 
                                 g 
                               
                                
                               
                                 P 
                                 g 
                                 T 
                               
                             
                           
                         
                       
                     
                     
                       
                         
                           
                             X 
                             B 
                             S 
                           
                           = 
                           
                             
                               
                                 X 
                                 B 
                               
                               - 
                               
                                 X 
                                 B 
                                 C 
                               
                             
                             = 
                             
                               
                                 
                                   X 
                                   B 
                                   T 
                                 
                                  
                                 
                                   ( 
                                   
                                     I 
                                     - 
                                     
                                       
                                         P 
                                         g 
                                       
                                        
                                       
                                         P 
                                         g 
                                         T 
                                       
                                     
                                   
                                   ) 
                                 
                               
                               . 
                             
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   25 
                   ) 
                 
               
             
           
         
       
     
         [0046]    After the decomposition above is completed, the preliminary monitoring model of the electric arc furnace operating process is established. 
         [0047]    Step 3: Monitoring and diagnosing the fault of the electric arc furnace operating process; the preliminary monitoring model of the electric arc furnace is used to calculate the T 2  statistics and the SPE statistics. 
         [0048]    The score matrices T A   S  and T B   S  of special parts of two modes datasets and the residual error matrices E A   S  and E B   S  of the special parts can be obtained combination of Principal Component Analysis (PCA) algorithm. 
         [0000]    
       
         
           
             
               
                 
                   { 
                   
                     
                       
                         
                           
                             T 
                             A 
                             S 
                           
                           = 
                           
                             
                               
                                 X 
                                 A 
                                 S 
                               
                                
                               
                                 P 
                                 A 
                                 S 
                               
                             
                             = 
                             
                               
                                 
                                   X 
                                   A 
                                 
                                  
                                 
                                   ( 
                                   
                                     I 
                                     - 
                                     
                                       
                                         P 
                                         g 
                                       
                                        
                                       
                                         P 
                                         g 
                                         T 
                                       
                                     
                                   
                                   ) 
                                 
                               
                                
                               
                                 P 
                                 A 
                                 S 
                               
                             
                           
                         
                       
                     
                     
                       
                         
                           
                             
                               X 
                               ^ 
                             
                             A 
                             S 
                           
                           = 
                           
                             
                               
                                 T 
                                 A 
                                 S 
                               
                                
                               
                                 ( 
                                 
                                   P 
                                   A 
                                   S 
                                 
                                 ) 
                               
                             
                             T 
                           
                         
                       
                     
                     
                       
                         
                           
                             E 
                             A 
                             S 
                           
                           = 
                           
                             
                               X 
                               A 
                               S 
                             
                             - 
                             
                               
                                 X 
                                 ^ 
                               
                               A 
                               S 
                             
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   26 
                   ) 
                 
               
             
             
               
                 
                   { 
                   
                     
                       
                         
                           
                             T 
                             B 
                             S 
                           
                           = 
                           
                             
                               
                                 X 
                                 B 
                                 S 
                               
                                
                               
                                 P 
                                 B 
                                 S 
                               
                             
                             = 
                             
                               
                                 
                                   X 
                                   B 
                                 
                                  
                                 
                                   ( 
                                   
                                     I 
                                     - 
                                     
                                       
                                         P 
                                         g 
                                       
                                        
                                       
                                         P 
                                         g 
                                         T 
                                       
                                     
                                   
                                   ) 
                                 
                               
                                
                               
                                 P 
                                 B 
                                 S 
                               
                             
                           
                         
                       
                     
                     
                       
                         
                           
                             
                               X 
                               ^ 
                             
                             B 
                             S 
                           
                           = 
                           
                             
                               
                                 T 
                                 B 
                                 S 
                               
                                
                               
                                 ( 
                                 
                                   P 
                                   B 
                                   S 
                                 
                                 ) 
                               
                             
                             T 
                           
                         
                       
                     
                     
                       
                         
                           
                             E 
                             B 
                             S 
                           
                           = 
                           
                             
                               X 
                               B 
                               S 
                             
                             - 
                             
                               
                                 
                                   X 
                                   ^ 
                                 
                                 B 
                                 S 
                               
                               . 
                             
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   27 
                   ) 
                 
               
             
           
         
       
     
         [0049]    Where {circumflex over (X)} A   S  and {circumflex over (X)} B   S  are respectively the estimated matrices of X A   S  and X B   S , P A   S (J×R A   S ) and P B   S (J×R B   S ) are respectively the PCA loading matrices of two special subsets, R A   S  and R B   S  are respectively the number of principle components for each subset, and E A   S  and E B   S  are respectively the residual errors of the special subsets and also the residual errors of the final model. The loading matrices P A   S  and P B   S  of special subsets for mode A and mode B of the electric arc furnace are calculated as follows: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       1 
                       
                         
                           N 
                           A 
                         
                         - 
                         1 
                       
                     
                      
                     
                       ( 
                       
                         
                           
                             X 
                             A 
                             S 
                           
                            
                           
                             ( 
                             
                               X 
                               A 
                               S 
                             
                             ) 
                           
                         
                         T 
                       
                       ) 
                     
                      
                     
                       P 
                       A 
                       S 
                     
                   
                   = 
                   
                     
                       γ 
                       A 
                     
                      
                     
                       P 
                       A 
                       S 
                     
                   
                 
               
               
                 
                   ( 
                   28 
                   ) 
                 
               
             
             
               
                 
                   
                     
                       1 
                       
                         
                           N 
                           B 
                         
                         - 
                         1 
                       
                     
                      
                     
                       ( 
                       
                         
                           
                             X 
                             B 
                             S 
                           
                            
                           
                             ( 
                             
                               X 
                               B 
                               S 
                             
                             ) 
                           
                         
                         T 
                       
                       ) 
                     
                      
                     
                       P 
                       B 
                       S 
                     
                   
                   = 
                   
                     
                       γ 
                       B 
                     
                      
                     
                       
                         P 
                         B 
                         S 
                       
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   29 
                   ) 
                 
               
             
           
         
       
     
         [0050]    Where 
         [0000]    
       
         
           
             
               1 
               
                 
                   N 
                   A 
                 
                 - 
                 1 
               
             
              
             
               ( 
               
                 
                   
                     X 
                     A 
                     S 
                   
                    
                   
                     ( 
                     
                       X 
                       A 
                       S 
                     
                     ) 
                   
                 
                 T 
               
               ) 
             
              
             
                 
             
              
             and 
              
             
                 
             
              
             
               1 
               
                 
                   N 
                   B 
                 
                 - 
                 1 
               
             
              
             
               ( 
               
                 
                   
                     X 
                     B 
                     S 
                   
                    
                   
                     ( 
                     
                       X 
                       B 
                       S 
                     
                     ) 
                   
                 
                 T 
               
               ) 
             
           
         
       
     
         [0000]    are the covariance matrices of datasets of mode A and mode B of the electric arc furnace, γ A  and γ B  are respectively the eigenvalues of the covariance matrices, and P A   S εR J×1  and P B   S εR J×1  are respectively consist of the eigenvectors of the covariance matrices of mode A and mode B of the electric arc furnace. 
         [0051]    The datasets for mode A and mode B are decomposed into common datasets and special datasets, that is to say, X A   C  and X A   S  obtained from equation (24), and X B   C  and X B   S  obtained from equation (25). Then, the diagonal matrix which consists of these four matrix eigenvalues is calculated as follows: 
         [0052]    The eigenvalue of X A   C  is supposed to be φ, so the characteristic polynomial of X A   C  is |φ·E−X A   C | where E is the unit matrix. Make the characteristic polynomial equal to zero, that is to say, |φ·E−X A   C |=0. Then, the characteristic equation G(φ)=0 whose variable is φ can be obtained. By solving this characteristic equation, the eigenvalue of X A   C  can be obtained. The eigenvalue of X A   C  is used as the diagonal elements to obtain the diagonal matrix Λ A   C  and diagonal matrixes Λ A   S , Λ B   C  and Λ B   S  which consist of the eigenvalues of another three matrixes X A   S , X B   C  and X B   S  can be obtained in the same method. 
         [0053]    Fault detection analysis is performed in the common datasets and the special datasets of mode A and mode B respectively. In statistical process monitoring, based on the principle of hypothesis testing to determine whether the failure occurred in the process, and the Hotelling (T 2 ) statistics and the squared prediction error (SPE) statistics are used for fault detection, and faults can be detected by plotting control charts. 
         [0054]    X new (J×1) is the new observation data of the electric arc furnace operating process. By projecting X new  to the global main factor P g , the common datasets score vector of the new sample data for two modes t new   C =P g   T x new  can be obtained. At the same time, the T 2  statistics T C,A   2  and T C,B   2  of common parts of mode A and mode B can be obtained by the following equations (30) and (31). The T 2  statistics of the common parts of the two modes are calculated as follows. 
         [0000]        T   C,A   2 =( t   new   C ) T (Λ A   C ) −1   t   new   C   (30)
 
         [0000]        T   C,B   2 =( t   new   C ) T (Λ B   C ) −1   t   new   C   (31)
 
         [0055]    The scores of the special subsets of the electric arc furnace operating process, to t A,new   S =(P A   S ) T x new  and t B,new   S =(P B   S ) T x new , are calculated in the similar method. Meanwhile, equations (32) and (33) below are used to calculate the T 2  statistics of the special subsets of mode A and mode B. 
         [0000]        T   S,A   2 =( t   A,new   S ) T (Λ A   S ) −1   t   A,new   S   (32)
 
         [0000]        T   S,B   2 =( t   B,new   S ) T (Λ B   S ) −1   t   B,new   S   (33)
 
         [0000]    Λ A   S  and Λ B   S  respectively represent the diagonal matrixes which consist of the eigenvalues of the special subsets. 
         [0056]    To calculate the final residual errors e A,new   S  and e B,new   S , {circumflex over (x)} A,new   S  and {circumflex over (x)} B,new   S  are supposed to be the estimated values of x new   S  in mode A and mode B. {circumflex over (x)} A,new   S , {circumflex over (x)} B,new   S  and x new   S  are submitted into equations (26) and (27) to obtain following equations: 
         [0000]        e   A,new   S   =x   new   S   −{circumflex over (x)}   A,new   S =( I−[P   g   ,P   A   S   ][P   g   ,P   A   S ] T ) x   new   (34)
 
         [0000]        e   B,new   S   =x   new   S   −{circumflex over (x)}   B,new   S =( I−[P   g   ,P   B   S   ][P   g   ,P   B   S ] T ) x   new   (35)
 
         [0057]    The SPE statistics SPE A,new   S  and SPE B,new   S  of mode A and mode B are calculated by the residual errors as follows: 
         [0000]        SPE   A,new   S =( e   A,new   S ) T   e   A,new   S   (36)
 
         [0000]        SPE   B,new   S =( e   B,new   S ) T   e   B,new   S   (37)
 
         [0058]    The control limit of the Hotelling T 2  statistics is calculated as follows: 
         [0000]    
       
         
           
             
               
                 
                   
                     T 
                     
                       lim 
                       , 
                       A 
                     
                     2 
                   
                   = 
                   
                     
                       
                         p 
                         · 
                         
                           ( 
                           
                             
                               n 
                               A 
                             
                             - 
                             1 
                           
                           ) 
                         
                       
                       
                         
                           n 
                           A 
                         
                         - 
                         p 
                       
                     
                      
                     
                       F 
                        
                       
                         ( 
                         
                           p 
                           , 
                           
                             
                               n 
                               A 
                             
                             - 
                             1 
                           
                           , 
                           α 
                         
                         ) 
                       
                     
                   
                 
               
               
                 
                   ( 
                   38 
                   ) 
                 
               
             
             
               
                 
                   
                     T 
                     
                       lim 
                       , 
                       B 
                     
                     2 
                   
                   = 
                   
                     
                       
                         p 
                         · 
                         
                           ( 
                           
                             
                               n 
                               B 
                             
                             - 
                             1 
                           
                           ) 
                         
                       
                       
                         
                           n 
                           B 
                         
                         - 
                         p 
                       
                     
                      
                     
                       
                         F 
                          
                         
                           ( 
                           
                             p 
                             , 
                             
                               
                                 n 
                                 B 
                               
                               - 
                               1 
                             
                             , 
                             α 
                           
                           ) 
                         
                       
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   39 
                   ) 
                 
               
             
           
         
       
     
         [0059]    Where F(p, n A −1, α) represents F distribution, of which the confidence limit is α, and degrees of freedom are respectively P and n A −1; F(p, n B −1,α) represents F distribution, of which the confidence limit is α, and degrees of freedom are respectively p and n B −1. 
         [0060]    In the residual subspace, the control limit of the SPE statistics SPE A  and SPE B  of mode A and mode B can be calculated by the equations as follows: 
         [0000]        SPE   A   ˜g   A χ h     A     2   (40)
 
         [0000]        SPE   B   ˜g   B χ h     B     2   (41)
 
         [0061]    Where χ h     A     2  represents χ 2  distribution, of which the degree of freedom is h A . g A =b A /2a A , h A =2(a A ) 2 /b A , a A  and b A  are respectively the SPE estimated mean and variance of mode A in normal operation. χ h     B     2  represents χ 2  distribution, of which the degree of freedom is h B . g B =b B /2a B , h B =2(a B ) 2 /b B , a B  and b B  are respectively the SPE estimated mean and variance of mode B in normal operation. 
         [0062]    The T 2  statistics and the SPE statistics are used for the fault monitoring. By comparing the control charts plotted respectively in mode A and mode B, the current operating mode of the electric arc furnace can be determined, and online fault detection can be performed in the corresponding operating mode. When the T 2  statistics and the SPE statistics of the corresponding operating mode go beyond the control limits, some fault may occur, otherwise, the whole operating process is normal. 
         [0063]    The advantage of present invention is: This device uses the radiation thermometry method can conduct long-range non-contact temperature measurement; and because the radiation intensity of two different kinds of visible light can be acquired, the colorimetric temperature measurement can be used to improve the calculation accuracy; the wireless communication protocol can form the wireless sensor network to make field wiring convenient, and the receiving and responding mechanism is used to ensure the accuracy of data transmission, and when communication failures occur due to the special situation, this mechanism also can intelligently determine the cause of the error and attempts to restore communication; the microcontroller of the receiving end realizes the same communication protocol as the PLC and makes different equipment compatible; this invention is suitable for working in a variety of different stable modes of industrial production process and it can identify the current operating mode of the electric arc furnace intelligently, and can diagnose faults for this work mode; this device has a good human-computer interface which can display temperature history curves to show the temperature change trend intuitively in the operating process, and it has the alarm function as well. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0064]    The present invention can be understood in more detail by reading the subsequent detailed description in conjunction with the examples and references made to the accompanying drawings, wherein: 
           [0065]      FIG. 1  is the overall structural diagram of the fault detector in the present invention; 
           [0066]      FIG. 2  is the temperature signal acquisition equipment structural diagram in the present invention; 
           [0067]      FIG. 3  is the sender controller program flowchart in the present invention; 
           [0068]      FIG. 4  is the receiver controller program flowchart in the present invention; 
           [0069]      FIG. 5  is the transceiver integrated circuit schematic in the present invention; 
           [0070]      FIG. 6  is the serial communication circuit schematic in the present invention; 
           [0071]      FIG. 7  is the electric arc furnace operating process schematic in the present invention; 
           [0072]      FIG. 8  is the control chart of the T 2  statistics of the common part of the dataset in normal operation; 
           [0073]      FIG. 9  is the control chart of the T 2  statistics of the special part of the dataset in normal operation; 
           [0074]      FIG. 10  is the control chart of the SPE statistics of the complete dataset in normal operation; 
           [0075]      FIG. 11  is the control chart of the T 2  statistics of the common part of the new observation dataset; 
           [0076]      FIG. 12  is the control chart of the T2 statistics of the special part of the new observation dataset; 
           [0077]      FIG. 13  is the control chart of the SPE statistics of the complete new observation dataset; and 
           [0078]      FIG. 14  is the flow chart of the fault diagnosis method. 
       
    
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
       [0079]    In the embodiment, the type of the controller is CC2430. The type of the RF transceiver is CC2591. The AC transformer is a LYM-0.5 bus type low voltage current transformer. The type of the PLC is Siemens S7-200. The imaging device is a telescope. The type of the color sensor is TCS320. 
         [0080]    As shown in  FIG. 1 , the fault detector comprises the temperature signal acquisition equipment, the current signal acquisition equipment and a host computer. The current signal acquisition equipment comprises an AC transformer and a PLC. The AC transformer is respectively connected with the three electrodes of the electric arc furnace in series. The large current signals which go through the electrodes is transformed into standard signals of 0 to 5 V, and then, the standard signals are transmitted to the analog input terminal of the PLC. After analogue-to-digital conversion by the PLC, the digital signals are transmitted to the host computer through the serial port of PLC. As shown in  FIG. 2 , the temperature signal acquisition equipment comprises an imaging device, a color sensor, two groups of transceiver integrated circuits and a serial communication circuit. The imaging device comprises an objective lens and an eyepiece; the transceiver integrated circuit comprises a microcontroller and a RF transceiver. Of two groups of transceiver integrated circuits, one RF transceiver works as a transmitter, and another RF transceiver works as a receiver. The input terminal of the color sensor is connected to the imaging device. The output terminal of the color sensor is connected to the input terminal of microcontroller of the first group transceiver integrated circuit. The transceiver integrated circuit is shown in  FIG. 5 . The output terminals of the microcontroller RF_N, TXRX_SWITCH, P1 — 1 and P1 — 7 are respectively connected to the input terminals of the RF transmitter RF_N, RF_P, HGM, EN and PAEN. The output terminal of the RF transmitter is connected to the input terminal of the RF receiver of the second group transceiver integrated circuit through wireless network. The output terminals of the RF receiver RF_N, RF_P, HGM, EN and PAEN are respectively connected to the input terminals of the microcontroller RF_N, TXRX_SWITCH, P1 — 1 and P1 — 7. The P0 — 2 terminal of the microcontroller is connected to the P0-2/RX terminal of the serial communication circuit. The P0 — 3 terminal of the microcontroller is connected to the P0-3/TX terminal of the serial communication circuit. The TIOUT and RIOUT terminals of the serial communication circuit are connected to the host computer. The serial communication circuit is shown in  FIG. 6 . The objective lens of the imaging device aims the measured target. The photosensitive sensor is fixed in the eyepiece at one end, so that the measured target radiation light can be projected to the photosensitive parts of the color sensor by the imaging device. The color sensor is connected to the microcontroller. When collecting the signal, the microcontroller controls the RF transmitter circuit to transmit wireless RF signals. The RF receiver circuit controls the serial communication circuit to send data to the host computer under the coordination of the microcontroller after receiving signals. Two identical transceiver integrated circuits are used to transmit and receive wireless signals and to play a different role in the control of the microcontroller. The configuration software of the host computer has the function of real-time display, history enquiry, data storage and alarm processing, etc. 
         [0081]    In the electric arc furnace production site, the electromagnetic environment is complex. Therefore, the fault detector in the present invention uses the wireless network protocol based on IEEE 802.15.4 specification. The protocol has the functions of automatic search network, automatic recovery after errors and data transmission correction. The protocol is also used in the input terminal and the output terminal of the microcontroller, which ensures the accuracy of data transmission. The detail is shown in  FIG. 3 , the sender controller program flowchart. Acquired data are encapsulated as a frame by the microcontroller, including loading the local address, the target address and loading request confirmation flag, adding the frame length and calculating the CRC cycle test code. The transmitter circuit is provided by Texas Instruments Inc. and has the function of demodulation. The transmitter circuit contains RF SOC (system on Chip) chip CC2430 of microcontroller core 8051 and can support the IEEE 802.15.4 specification at the physical layer to ensure that data can be transmitted correctly. 
         [0082]    The communication protocol between S7-200 PLC and the configuration software is obtained by analyzing the communication process between S7-200 PLC and the configuration software, and the C language is used to realize the protocol on the input terminal of the microcontroller, and realizes the goal of device compatibility finally. The specific implementation method is shown in  FIG. 4 . The data are packaged according to the format of the communication between the PLC and the host computer before the data are transmitted, including loading the local address, write data command, registers type, offset address. Calculate the number of data, load data and load calibration bytes. After packaging, the data are transmitted to the host computer by the serial port. 
         [0083]    Operating process of the electric arc furnace: 
         [0084]    The schematic diagram of the operating process of the electric arc furnace is shown in  FIG. 7 . Number  1  is the current transformer. Number  2 ,  3  and  4  are three electrodes. Number  5  is the furnace wall. Number  6  is the molten pool. Number  7  is the discharged gas in the production process. Number  8  is the fixed pulley. 
         [0085]    The imaging device, the color sensor and the transceiver integrated circuit are fixed in the supporting structure with the fixed pulley, and the objective aims at the molten pool. The transceiver integrated circuit, the serial communication circuit and the host computer are placed in the workbench. After start of the production, two kinds of red and green visible light signals radiated from the molten pool are acquired, and then, the signals are transmitted to the receiving circuit through the transceiver integrated circuits. Under the control of the microcontroller, the signals are transmitted to the host computer. The PLC and the current transformer are linked together to acquire current data and transmit the data to the host computer. 
         [0086]    The electric arc furnace produces a lot of waste gas in the operating process. If too much gas is accumulated, the molten materials in the electric arc furnace may erupt and even cause failure of material spray. The raw materials of the electric arc furnace operating process include magnesium powder and magnesium pieces. Some of the data acquired for two modes are shown as follows, where mode A is corresponding to the magnesium powder and mode B is corresponding to the magnesium pieces. 
         [0000]    
       
         
               
             
               
               
               
             
               
               
               
               
               
             
               
               
               
               
               
             
           
               
                 TABLE 1 
               
             
             
               
                   
               
               
                 Some data of the modeling data in mode A 
               
             
          
           
               
                   
                 Variable 
                   
               
             
          
           
               
                   
                 Current I A   
                 Current I B   
                 Current I C   
                 Temperature t 
               
               
                 No. 
                 (Ampere) 
                 (Ampere) 
                 (Ampere) 
                 (Centigrade) 
               
               
                   
               
             
          
           
               
                 1 
                 8264.47 
                 10434.55 
                 8824.69 
                 2853.2 
               
               
                 2 
                 8043.67 
                 9955.90 
                 8534.97 
                 2860.4 
               
               
                 3 
                 8111.64 
                 9807.76 
                 8392.92 
                 2858.5 
               
               
                 4 
                 8199.31 
                 10056.70 
                 8415.42 
                 2861.1 
               
               
                 5 
                 8647.01 
                 10402.67 
                 8082.11 
                 2866.7 
               
               
                 6 
                 8962.98 
                 10682.08 
                 7913.34 
                 2855.4 
               
               
                 7 
                 8880.94 
                 10739.74 
                 8495.12 
                 2864.3 
               
               
                 8 
                 8625.92 
                 10835.38 
                 8957.83 
                 2860.9 
               
               
                 9 
                 8551.85 
                 10704.11 
                 8975.64 
                 2850.5 
               
               
                 10 
                 8659.20 
                 10614.10 
                 8686.86 
                 2855.7 
               
               
                   
               
             
          
         
       
     
         [0000]    
       
         
               
             
               
               
               
             
               
               
               
               
               
             
               
               
               
               
               
             
           
               
                 TABLE 2 
               
             
             
               
                   
               
               
                 Some data of the modeling data in mode B 
               
             
          
           
               
                   
                 Variable 
                   
               
             
          
           
               
                   
                 Current I A   
                 Current I B   
                 Current I C   
                 Temperature t 
               
               
                 No. 
                 (Ampere) 
                 (Ampere) 
                 (Ampere) 
                 (Centigrade) 
               
               
                   
               
             
          
           
               
                 1 
                 11353.87 
                 9858.86 
                 9616.025 
                 2816.3 
               
               
                 2 
                 11199.16 
                 9947.46 
                 9539.61 
                 2820.5 
               
               
                 3 
                 11512.79 
                 10553.63 
                 11540.92 
                 2826.8 
               
               
                 4 
                 11335.58 
                 10410.64 
                 10027.16 
                 2830.1 
               
               
                 5 
                 11290.58 
                 10220.31 
                 9491.79 
                 2837.2 
               
               
                 6 
                 11475.76 
                 10212.81 
                 10669.42 
                 2835.9 
               
               
                 7 
                 11344.49 
                 9866.36 
                 9772.60 
                 2840.4 
               
               
                 8 
                 11200.10 
                 9642.74 
                 9315.05 
                 2833.3 
               
               
                 9 
                 11315.89 
                 9970.90 
                 9643.68 
                 2844.7 
               
               
                 10 
                 11502.95 
                 9994.81 
                 10036.07 
                 2839.3 
               
               
                   
               
             
          
         
       
     
         [0087]    By identifying the operating mode of the electric arc furnace, targeted fault detection is performed, and detecting results are analyzed. 
         [0088]    The electric arc furnaces running fault detection device is used for multimode fault monitor and diagnosis. The multimode fault monitoring and diagnosis contains the following steps of (see  FIG. 14 ): 
         [0089]    Step 1: Acquiring and standardizing data. 
         [0090]    The training data and the real-time data are acquired at the sampling interval of 1.2 seconds. The data in one of the normal operating modes are used to build mode A, the data in another normal operating mode is used to build mode B. In mode A, 349 sample points are selected, and in mode B, 403 sample points are selected, where the confidence limit is 99%. Each sample includes four variables, and the mean and the standard deviation is used to make data standard. 
         [0091]    Step 2: Decomposing the data. The common part and the special part of the two modes are decomposed, and build models of two parts. Through calculating the model global factor matrix P g  of the monitoring model for mode A and mode B, datasets obtained in mode A and mode B of the electric arc furnace are decomposed, and the common subsets with close relation and the special subsets of mode A and mode B are obtained. 
         [0092]    The sample datasets in mode A and mode B are standardized. The global factor matrix P g  is obtained through the methods mentioned above, and then, the datasets for mode A and the datasets for mode B are decomposed to obtain the common subsets and the special subsets. The PCA model is established for each subset. Loading vectors P A   S  and P B   S  of the special subsets for the two modes are calculated respectively. x new  is supposed to be an available new sample in the observation data, and then, x new  is projected by P g  to get the score vector of the common part. x new  is projected by P A   S  and P B   S  to obtain the score vector of the special part. Then, the T 2  statistics are obtained, and the SPE statistics are obtained by calculating the residual error. 
         [0093]    Step 3: Fault detection and diagnosis are used for the operating process of the electric arc furnace. 
         [0094]    Normal operation data for mode A and mode B are selected. There are 349 sample points and 3 variable datasets for mode A, and there are 403 sample points and 3 variable datasets for mode B. Models of mode A and B are established respectively. First of all, a group of new observation is tested. As shown in  FIG. 8  and  FIG. 9 , both the T 2  statistics of the common part and the T 2  statistics of the special part in mode A are not beyond the confidence area. However, the T 2  statistics of the special part in mode B is far beyond the control limit. As shown in  FIG. 10 , the SPE statistics of two modes are not beyond the control limit. Thus, it can determine the testing data is normal data. The comparison of the three figures shows that current data is much closer to mode A, and the current operating mode of the electric arc furnace can be determined in mode A. 
         [0095]    The analysis above shows that the electric arc furnace operates in mode A, so fault detection is performed in mode A. As shown in  FIGS. 11 ,  12  and  13 , the comparison of three figures shows that the T 2  statistics of the special part overrun most seriously, but the T 2  statistics of the common part overrun only in part of the areas, and most of the areas do not overrun. Therefore, it can determine a failure has occurred, and it concentrates in the special part. 
         [0096]    The example above illustrates the effectiveness of the fault detector for the operating process of the electric arc furnace and the method thereof, and fault monitoring for the operating process of the electric arc furnace is realized.