Abstract:
The present invention is directed to a mathematical approach to detect faults by reconciling known data driven techniques with a physical understanding of the HVAC system and providing a direct linkage between model parameters and physical system quantities to arrive at classification rules that are easy to interpret, calibrate and implement. The fault modes of interest are low system refrigerant charge and air filter plugging. System data from standard sensors is analyzed under no-fault and full-fault conditions. The data is screened to uncover patterns though which the faults of interest manifest in sensor data and the patterns are analyzed and combined with available physical system information to develop an underlying principle that links failures to measured sensor responses. These principles are then translated into online algorithms for failure detection.

Description:
[0001]     This application claims priority to U.S. Provisional Application No. 60/605,080 which was filed on Aug. 27, 2004 and U.S. Provisional Application No. 60/635,523 which was filed on Dec. 13, 2004. 
     
    
     BACKGROUND OF THE INVENTION  
       [0002]     The present invention relates to heating, ventilation and cooling (HVAC) systems, and more particularly to fault diagnostics associated with early detection and isolation of failures in HVAC systems.  
         [0003]     HVAC systems often do not function as well as expected due to faults developed during routine operation. While these faults are indicative of a failure mode, many faults do not result in immediate system shut down or costly damages. However, most faults, if unnoticed for a long period of time, could adversely affect system performance, life, and lifecycle cost.  
         [0004]     While diagnostics refer to detection and isolation of faults, prognostic typically refers to predicting faults before they occur. In many applications, however, early detection and diagnostics may serve the same end as prognostics. This is the case when failure propagation happens at a reasonably slow pace. Small changes in system parameters typically do not have a substantial adverse effect initially. As such, accurate prediction of the time between detection of a fault, that is, a small change to one or more system parameters, to full system deterioration or shutdown is not crucial. For instance, detection of HVAC system refrigerant charge leakage and air filter plugging are examples of failure modes for which early detection of changes provides adequate information to take timely maintenance action.  
         [0005]     Approaches to diagnostics may be divided into two broad categories. One category deals with direct measurement of monitored quantities and another category combines sensing technologies with mathematical algorithms. The technical emphasis in these approaches is the development of dedicated sensors for measurement of crucial system parameters. While such approaches may be more accurate, they are typically costly as they involve adding dedicated hardware for each failure mode of interest. In the combined approach, algorithms play the major role since they allow inference about the health of the system from indirect measurements provided by the sensors. Because the addition of new sensors is more expensive and more difficult to manufacture, algorithms alone are incorporated to utilize available sensors that are configured for a control purpose.  
         [0006]     Design of failure detection and diagnostics algorithms have been subject of extensive research ranging from statistical approaches and reviews to techniques derived from artificial intelligence and reasoning, graph theory, and bond graphs. Several diagnostics techniques have been applied to the problem of chiller and HVAC fault isolation. Among known approaches, “black-box” or data-driven techniques (such as neural networks) have received considerable attention. Such approaches are well suited to domains where data is abundant but physical knowledge about the phenomenon is scarce. However, one problem with such approaches is that recalibrating the parameters of the black-box model typically requires extensive re-experimentation even if the system changes slightly as there is no direct linkage between model parameters and physical system quantities.  
         [0007]     As such, there is a desire for an analytical approach to detect faults which reconciles the results of known data driven techniques with a physical understanding of the HVAC system and provides a direct linkage between model parameters and physical system quantities to arrive at classification rules that are easy to interpret, calibrate and implement.  
       SUMMARY OF THE INVENTION  
       [0008]     The present invention is directed to an analytical approach to detect faults by reconciling known data driven techniques with a physical understanding of the HVAC system and providing a direct linkage between model parameters and physical system quantities to arrive at classification rules that are easy to interpret, calibrate and implement.  
         [0009]     The present invention focuses on two of the most common problems encountered a multi-modular split HVAC system, which are detecting low refrigerant conditions and air filter plugging. A method for refrigerant charge leak detection is disclosed that relies on a systematic technique for analysis of experimental data, extraction of fault signatures, formulation of fault detection principals, and development and implementation of diagnostic algorithms. A method of detecting air filter plugging is also disclosed that relies on reduced physics-based relationships in heat exchangers to estimate air mass flow through the heat exchangers.  
         [0010]     Both methods incorporate data filtering techniques to determine which portions of data carry the most information regarding the underlying failure, variable sub-selection based upon available sensors, calculation of a distance between faulty and normal data sets, and maximization of this distance with respect to filtering parameters and variable sub-selection. The sub-selected variables are then processed by classification techniques to generate easy to interpret and easy to implement classification rules.  
         [0011]     These and other features of the present invention can be best understood from the following specification and drawings, the following of which is a brief description. 
     
    
     BRIEF DESCRIPTION OF THE DRAWINGS  
       [0012]      FIG. 1  is schematic illustration of an example HVAC system according to the present invention;  
         [0013]      FIG. 2  is a graph illustrating how data filtering is used to zoom in on data depending on a transient response of a base signal;  
         [0014]      FIG. 3  is a flow chart detailing an example charge leakage calculation;  
         [0015]      FIG. 4  is a flow chart detailing an example air filter plugging calculation;  
         [0016]      FIG. 5  illustrates an algorithm response to induced low system refrigerant charge conditions for a first example MMS system;  
         [0017]      FIG. 6  illustrates an algorithm response to induced low system refrigerant charge conditions for a second example MMS system;  
         [0018]      FIG. 7  is a graph illustrating a receiver operating characteristic (ROC) for a low charge detection algorithm;  
         [0019]      FIG. 8  illustrates an algorithm response to induced air filter plugging conditions for a high wall MMS system;  
         [0020]      FIG. 9  illustrates an algorithm response to induced air filter plugging conditions for a 4-way MMS system; and  
         [0021]      FIG. 10  is a graph illustrating a ROC for an air filter plugging detection algorithm. 
     
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT  
       [0022]      FIG. 1  is a schematic illustration of an example HVAC system  10  according to the present invention. In this example, the HVAC system  10  is a duct-free heat pump system known as a multi-modular split system (MMS). The MMS  10  includes one outdoor unit  12  and two indoor units  14 A and  14 B, which operate in a cooling mode to provide cool air to an interior space during a warm or hot season and a heating mode to provide warm air to the interior space during a cool or cold season.  
         [0023]     The outdoor unit  12  includes a pair of parallel compressors  16 , which are variable speed, an outdoor expansion valve  18  to control a sub-cool in the cooling mode and to control a superheat in the heating mode, an outdoor heat exchanger  20 , which behaves as a condenser in the cooling mode and as an evaporator in the heating mode, and an outdoor fan  22 .  
         [0024]     Each of the two indoor units  14 A and  14 B includes an indoor expansion valve  24 , to control the sub-cool in the cooling mode and superheat in the heating mode, and an indoor heat exchanger  26 , which behaves as an evaporator in the cooling mode and as a condenser in the heating mode, and an indoor fan  28 .  
         [0025]     A 4-way valve  30  controls the mode of operation from the cooling mode to the heating mode and vice versa. The MMS  10  also includes a receiver tank  32  for storage of refrigerant charge that is operable to change the amount of refrigerant charge circulated depending on the conditions.  
         [0026]     The speed of the compressors  16  and the indoor fan  28  are adjusted in response to a deviation between a room temperature and a set point. The speed of the compressors  16  is further adjusted to match the total cooling or heating demand.  
         [0027]     Expansion valves are disposed throughout the MMS  10 . In the illustrated example, the expansion valves are pulse modulated valves  34  that are actuated via pulse modulation. The pulse modulated valves  34  are controlled by an actuation signal that adjusts an opening of the pulse modulated valves  34  to control a flow of refrigerant through the MMS  10 . Pulse modulated valves  34 A and  34 B are positioned in-line proximate to the indoor heat exchangers  26 . A pair of pulse modulated valves  35  is positioned in-line between a coil  36  and the receiver tank  32 .  
         [0028]     Multiple sensors are also disposed throughout the MMS  10 . In the illustrated example, the sensors include a plurality of refrigerant-side temperature sensors  38 , air-side temperature sensors  40 , and pressure sensors  42 .  
         [0029]     Refrigerant-side temperature sensors  38 A- 38 D are positioned near each end of each of the indoor heat exchangers  26 . Refrigerant-side temperature sensors  38 E and  38 F are positioned near one end of each of the compressors  16 . Refrigerant-side temperature sensor  38 G is positioned between the 4-way valve  30  and an accumulator  44 . Refrigerant-side temperature sensor  38 H is positioned between the pair of pulse modulated valves  35  and the receiver tank  34 . Refrigerant-side sensor  381  is positioned between the outdoor heat exchanger  20  and the coil  36 .  
         [0030]     Air-side temperature sensors  40 A and  40 B are positioned between the indoor fans  28  and the indoor heat exchangers  26  and air-side temperature sensor  40 C is positioned proximate to the outdoor fan  20 .  
         [0031]     Pressure sensor  42 A is positioned proximate to the accumulator  44  and pressure sensor  42 B is positioned between the compressors  16  and an oil separator  46 .  
         [0032]     The present invention focuses on developing algorithms to detect faults within the MMS  10  by utilizing existing system sensors and data. An analytical method of diagnosing at least one system fault utilizing existing system sensors is disclosed. The analytical method includes identifying sensors that are available within a given system, analyzing sensor data to determine which of the available sensors generate data indicative of a system fault based upon a maximum separation/minimum overlap between no-fault data and full-fault data for each available sensor, determining a fault relationship based on the analysis, comparing at least one measured system characteristic to the fault relationship and generating at least one fault code indicative of a failure mode when a system fault is identified. However, the system fault identified is not necessarily the system characteristic directly monitored by the sensor generating the data that is being analyzed. For example, the sensor may generate data associated with a pressure within the system, however, the failure code generated may be indicative of a low system refrigerant charge or an air filter plugging condition.  
         [0033]     The following discusses the development of both a low system refrigerant charge indicator and an air filter plugging indicator of the present invention.  
         [0000]     Development of a Low System Refrigerant Charge Indicator.  
         [0034]     A total MMS refrigerant charge is roughly proportional to a total volume filled with liquid refrigerant at any given point in time. Low system refrigerant charge occurs when the total system liquid volume drops. When the MMS  10  operates in the cooling mode, a total MMS refrigerant mass drops and the total MMS volume increases resulting in an overall increase in a volume of vapor refrigerant relative to a volume of liquid refrigerant. The overall increase in total MMS volume causes an increase in temperature of refrigerant exiting the indoor and outdoor heat exchangers  20  and  26  to a temperature above the boiling point of the refrigerant when the refrigerant is in a vapor state or a decrease in temperature of refrigerant exiting the indoor and outdoor heat exchangers  20  and  26  to a temperature below the boiling point of the refrigerant when the refrigerant is in a liquid state. These phenomena are known as superheat and sub-cool respectively. Under these conditions, the MMS  10  will tend to increase the vapor quality at those temperatures if saturated, thereby controlling the superheat.  
         [0035]     In most systems, the superheat is controlled by actuating an expansion valve. An increase in superheat requires a greater opening of the expansion valve to allow an increased flow of refrigerant to pass through the indoor and outdoor heat exchangers  20  and  26  to maintain the desired superheat. In the example illustrated, the expansion valve  34  is a pulse modulated valve (PMV) and actuation pf the PMV is controlled via pulse modulation. As such, an increase in superheat translates into a higher pulse actuation which in turn correlates to low system refrigerant charge. In addition, low system refrigerant charge is also associated with lower system suction pressures.  
         [0036]     While the present invention focuses on the cooling mode, the principles discussed are not limited to the cooling mode and may also extend to the heating mode.  
         [0037]     The present invention includes a fault detection principal for identifying low system refrigerant charge which may be summarized as follows:  
         [0038]     Low refrigerant charge indicator uses larger average PMV opening of the indoor units OR smaller values of suction pressure observed during the first few minutes (e.g. 10-20 minutes) after the compressor speed (rpm) exceeds a threshold (e.g. 40%-50% of its maximum rpm) excluding first 0-2 minutes of transient data.  
         [0039]     This fault detection principle is the result of a systematic approach for identification of fault principles and decision rules for low refrigerant charge detection. This approach applies several statistical classification techniques to pre-processed field trial data as described below. The pre-processing involves application of a data filtering process whose objective is to zoom into parts of the data that carry most information about the fault event of interest.  
         [0040]     The data filtering process decomposes a time-axis into time intervals during which “interesting” transient or steady state behavior occurs. Compressor speed (rpm), denoted by v(t) at time t, is used as the base-signal for filtering. The base signal is the signal whose temporal behavior is used to decompose the time axis as explained more closely below.  
         [0041]     The major step in the data filtering process is to decompose the time axis into a sequence of intervals F={I k }. This decomposition is a function of three filtering parameters κ≧0,Δ≧0,{overscore (v)}≧0 i.e. F=F(κ,Δ,{overscore (v)}) where {overscore (v)} is a threshold value, κ is a pre-determined time period during which an actual value measured by a sensor must remain above the threshold value {overscore (v)} before a fault can be generated, and Δ defines a start-up time period during which the data collected is not considered. F={I k } is constructed from another sequence {I′ k }. Each element, I k ′, is a closed time interval with the following three properties: (1) the length I′ k  is larger than or equal to κ, (2) v(t)≧{overscore (v)} for all t in I′ k , (3) I′ k  is the superset of all other overlapping closed intervals that possess the first two properties. Furthermore, I′ k ∩I′ j =Ø for k≈j. The sequence {I′ k }, k=1, 2, 3, . . . , is ordered in such a way that the supremum of all t belonging to I′ k , denoted by e(I′ k ), is smaller than the infimum of all t belonging to I′ k+1  denoted by b(I′ k+1 ). {I k } is constructed by letting b(I k )=b(I′ k ) while e(I k )=max(e(I′ k ),b(I k )+Δ).  
         [0042]     As illustrated in  FIG. 2 , the following example illustrates show how data filtering can be used to zoom in data depending on a transient response of the base signal: 
        1. HP(Δ,{overscore (v)})=F(0,Δ,{overscore (v)}) where Δ is small to medium, constructs a high-pass filter. Points A, B, and D are segmented out by HP(0,100%).     2. BP(κ,Δ,{overscore (v)})=F(κ,Δ,{overscore (v)})−F(κ,κ,{overscore (v)}) where κ is small to medium and Δ is medium to large, constructs a band-pass filter. Intervals C and E are segmented out by BP(1,2,100%).     3. LP(κ,{overscore (v)})=BP(κ,∞,{overscore (v)})=F(κ,∞,{overscore (v)})−F(κ,κ,{overscore (v)}) where κ is medium, constructs a low-pass filter. An interval F is segmented out by LP(4,100%).        
 
         [0046]     Also note that if the data in  FIG. 2  is sampled with a sampling interval  1 , then a union of intervals segmented by mutually exclusive sets HP(0,100%), BP(1,2,100%), and LP(4,100%) span all the data points for which the base signal is at 100% of its full operation range. For example, suppose the data are sampled at a sampling interval ΔT. The union of the mutually exclusive sets HP(i×ΔT,{overscore (v)}), BP((i+1)×(ΔT+1), j×ΔT, {overscore (v)}), and LP((i+j+2)×ΔT,{overscore (v)}), where i, j are integers, spans all the data points for which the base signal is larger than or equal to {overscore (v)}. Φ is used to generically refer to a union of all time intervals picked by the filter.  
         [0047]     To develop the fault detection principal, the time-axis is repeatedly decomposed into three disjoint sets of intervals characterized by high-pass, band-pass, and low-pass conditions using average speed of base compressors as the base signal and various filtering parameters. Each filtering results in a decomposition of time similar to the illustration in  FIG. 2 . For each repetition of the filtering process, the three disjointed data sets were analyzed based upon a set of statistical techniques consisting of sensor (variable) selection, calculation of distance between no-fault and full-fault data sets, and fault pattern discovery. The data filtering process resulted in selection of an opening associated with the PMV and a suction pressure as the variables carrying most information about system refrigerant charge leakage. An optimal filter was found to be BP(κ,Δ,{overscore (v)}) with κ≦2, 10≦Δ≦20, and {overscore (v)}≈50% of full compressor speed. Combining these results led to the low system refrigerant charge fault detection principle stated above.  
         [0048]     The fault detection principal for low system refrigerant charge is converted into an algorithm for low system refrigerant charge detection and implemented by calculating a low system refrigerant charge indicator over a batch of data of fixed length. The batch of data contains the most recent data points that have passed through the filter up to a point determined by its fixed length. As more data points become available, they replace the oldest data points, keeping the length of the data batch fixed.  
         [0049]     To illustrate, consider  FIG. 2  again. Suppose that the filter used in the calculations is BP(1,2,100%), segmenting intervals C and E. Also suppose that the sampling time is 1 and the fixed data length is 3. The batches of data in chronological order for this example will be {9,10,15}, {10,15,16}, and {15,16,17}. A time calculated as an average (or median) time of the batches member points is assigned to each data batch. In this example, the times associated with {9,10,15}, {10,15,16}, and {15,16,17} will be 11.3, 13.7, and 16 respectively using averaging over the time points in the batch. A low system charge indicator within each batch is calculated by finding the fraction of points within the batch for which the PMV opening is above a threshold or the suction pressure is below a threshold. The calculated low system charge indicator is assigned to the time associated with the batch.  
         [0050]      FIG. 3  is a flow chart detailing an example charge leakage calculation. In this flow chart: 
        Φ(t) is a binary indicator that assumes one or zero. Φ(t) is equal to one only if t belongs to an analysis period of interest, a set of all intervals picked by the designed filter, and under a cooling mode operating condition (excluding compressor protection conditions).     A batch of data, B, always keeps the most recent time points for which Φ(t)=1. A fixed number of points NB in the batch are set by a user. Typical values for the MMS application are 24×60 or 12×60.     C(B) denotes a number of elements of B.     δ represents a sampling time. δ should be larger than the data collection sampling time (in MMS case 1 minute), but may be selected larger to accelerate computations. 
 
 Development of an Air Filter Plugging Indicator. 
         
         [0055]     An air filter plugging algorithm uses simple physical relations to estimate a mass flow of air through each heat exchanger. Due to lack of sufficient airside measurements, these relations are based on several simplifying assumptions as described below: 
        Log-mean temperature and energy balance equations are applied to estimate the mass flow of air in the heat exchanger.     The mass flow of air is inversely proportional to air filter resistance.     A portion of the heat exchanger where refrigerant superheats has large UA.     Refrigerant flow is proportional to an opening associated with the PMV.     Air latent cooling is negligible.        
 
         [0061]     Based on these assumptions, the inverse of air mass flow in each heat exchanger is estimated via:  
       F   =       (       T   a     -     Tc   j       )       PMV   ⁡     (       Tc   1     -     Tc   j       )             
 
         [0062]     As illustrated in  FIG. 1 , T a  is an indoor air temperature, Tc j  is a middle of an indoor coil temperature (measuring a saturation temperature of the refrigerant in the heat exchanger), Tc 1  is a temperature of the superheated refrigerant flowing out of the indoor unit heat exchangers  26 , and PMV is an actuation signal of the pulse modulated expansion valves  34 A and  34 B.  
         [0063]     In  FIG. 1 , the measurement points T a ,Tc 1 ,Tc j  are respectively marked/numbered by TA-A, TC1-A, and 11 for indoor unit  14 A and TA-B, TC1-B, and 13 for indoor unit  14 B.  
         [0064]     Because these assumptions ignore several sources of variability such as latent cooling, an analysis is performed over longer periods to “average out” possible unknown effects. Similar to the development of the low system refrigerant charge indicator discussed above, batches of data are processed for calculation of the air filter plugging indicator. From this point of view, the algorithm is identical to the calculation of the low system refrigerant charge indicator for data batches of certain size.  
         [0065]     Because the air filter plugging algorithm estimates the inverse of air mass flow in each indoor heat exchanger  26  while the indoor fan  28  runs at various speeds, it is recommended that calculation of the air filter plugging indicator and its comparison to a threshold be performed for each fan speed separately. Although the fan speed data is not directly available, a set-point for the speed (FANTAP) is available and may be used as a proxy for the fan speed.  
         [0066]     The data filtering process for calculation of the air filter plugging indicator is much simpler than the data filtering process for low system refrigerant charge indicator calculation. The calculation of the air filter plugging indicator involves capturing sensor data that satisfies basic regularity principles needed for application of log mean temperature and energy balance in the cooling mode operations. The filter Φ is set to one at each time point if T a &gt;Tc j ,Tc 1 &gt;Tc j , PMV larger than a threshold (a value of about 100 is appropriate), and FANTAP is equal to a mode for which the calculation is performed.  
         [0067]      FIG. 4  is a flow chart detailing an example air filter plugging calculation. In this flow chart: 
        Φ(t) is a binary indicator that assumes one or zero. Φ(t) is equal to one only if t belongs to an analysis period of interest, a set of all intervals picked by the designed filter, and under a cooling mode operating condition (excluding compressor protection conditions).     A batch of data, B, always keeps the most recent time points for which Φ(t)=1. A fixed number of points NB in the batch are set by a user. Typical values for the MMS application are 24×60 or 12×60.     C(B) denotes a number of elements of B.     δ represents a sampling time. δ should be larger than the data collection sampling time (in MMS case 1 minute), but may be selected larger to accelerate computations.          
         [0072]      FIG. 5  and  FIG. 6  show an algorithm response to induced low system refrigerant charge conditions for a first example MMS and a second example MMS respectively. Both example MMS&#39;s include one outdoor unit and five indoor units which were designed to meet the cooling and heating demands of two office rooms and two conference rooms. In calculating the results, the thresholds that were optimized for the first system were directly applied to a data set of the second system without any further tuning. In other words, the first system data set was used for “algorithm training” while the second system data was used for validation.  
         [0073]     In regard to the low system refrigerant charge, the plots illustrated in  FIG. 5  and  FIG. 6  show the value of a low system refrigerant charge indicator as a function of the number of data batches processed by the algorithm. The data batches are ordered chronologically. The low system refrigerant charge indicator for each batch of data is calculated by finding a percentage of points within the batch for which an opening associated with the PMV is larger than a pre-determined threshold or a suction pressure is smaller than a pre-determined threshold. The variation in intensity of the plot line represents an actual charge averaged over all the points in the batch of data. For instance, if 40% of the points within the batch have 30% charge loss while the rest have full charge (0% charge loss), the average actual charge loss for that batch of data will be 12%, which is referred to as an average actual fault (AAF). From this plot, one may deduce that a value of indicator above 10-15% would flag low charge condition (&gt;25% AAF).  
         [0074]      FIG. 7  plots a receiver operating characteristic (ROC) of the low system refrigerant charge detection algorithm. ROC is a widely used tool for assessment of the performance of detection algorithms independent from a detection threshold. ROC plots detection rate (hit rate) of an algorithm as a function false alarms (false positives) generated by the algorithm. Detection rate or hit rate measures the probability of the low system refrigerant charge indicator raising an alarm given a failure event actually happens. False alarm or false positive measures the probability of indicating a failure when no failure is actually present. An ideal ROC curve will have 100% detection rate for any positive false alarm.  
         [0075]     Conceptually, ROC may be calculated by changing the detection threshold from its minimum to its maximum possible value, calculating a false alarm-detection rate pair for each chosen detection threshold, and plotting the calculated pairs. In other words, each point on the ROC curve is associated with one threshold value. When generating the ROC in  FIG. 7 , presence of a faulty condition for a batch of data is defined as AAF&gt;15% for that data batch. The absence of a fault for the data batch is defined as AAF&lt;5%.  
         [0076]      FIG. 8  and  FIG. 9  show an algorithm response to induced air filter plugging at a high wall unit and a 4-way unit respectively. The presented results are based on data collected from a first system at a high fan speed. The plots show a value of an air filter plugging indicator as a function of a number of data batches processed by the algorithm (ordered chronologically). A variation in intensity of the plot line represent the actual average fault (AAF), defined similarly above. For instance, if 40% of the points within the batch have 50% plugging while the rest have no plugging (0% plugging), the AAF for that batch will be 20%. From this plot, one may deduce that a value of the air filter plugging indicator above 15 would flag plugging conditions.  
         [0077]      FIG. 10  plots a ROC for the air filter plugging detection algorithm. When generating the ROC in  FIG. 10 , presence of a faulty condition for a batch of data is defined as AAF&gt;25%. The absence of fault for the data batch is defined as AAF&lt;5%.  
         [0078]     Although preferred embodiments of this invention have been disclosed, a worker of ordinary skill in this art would recognize that certain modifications would come within the scope of this invention. For that reason, the following claims should be studied to determine the true scope and content of this invention.