INTRODUCTION

Since the industrial revolution, humans have depended on energy by exploiting oil and coal for stable economic growth, future resource shortages will inevitably occur. Impacted by geographical factors, Taiwan does not have many natural resources, forcing the importation of resources from other countries. Because of complex accessibility to resources and the slow development of new resources, energy saving is important for Taiwan.

Electricity is the mother of industry and of the most important energy for modern life. With rapid economic development of industry and commerce, electricity consumption will also increase in the future. Industrial and residential electricity consumption will increase. Taiwan has a subtropical climate. Taiwan is characterized by hot and humid summers, causing demand for air conditioning to increase which has promoted the development of refrigerating air conditioning industry. Due to the popularity of air conditioning equipment, air temperature changes are closely related to peak loads for power systems. For electric power loads, power consumption of air conditioning equipment accounts for 30% of the total power consumption during the summer. For central air conditioning systems, chillers are the primary cooling source and comprise 60% of the total power consumed by the system. Reducing energy consumption for air conditioning equipment and increasing energy use efficiency is crucial. In air conditioning system planning, the maximum load for the system is often decided by the highest outside air temperature and maximum room load. The chiller rated capacity depends on the system load. Chillers seldom operate under the maximum load of air conditioners and often work under various loads. The optimization of chiller operations under various loads is one of the main energy-saving methods for central air conditioning systems (Shih, 2004).

PRINCIPLES AND METHOD

Biological neural networks exemplified by human brain operation and thinking are mysterious. Artificial neural networks were created by human beings to simulate biological neural network information processing systems. Neural network theory began in the 1950s. Prior to 1980, expert systems were the most popular artificial intelligence networks available. Artificial neural network theory was ignored and has not matured since. After 1980, artificially intelligent expert systems had bottleneck complications. Extensive attention has been paid to artificial neural network systems. By far, new architectures and theories for artificial neural networks have been proposed. With the increase of the computer speeds, artificial neural networks have more powerful functions and are in wide use. Their characteristics are as follows (Chiang, 2005):

Artificial neural network: Artificial neural networks are information processing systems that simulate biological neural networks. These are defined as:

... a mathematical model or computational model that is inspired by the structure and/or functional aspects of biological neural networks. A neural network consists of an interconnected group of artificial neurons and it processes information using a connectionist approach to computation. Artificial neurons are the simple simulation of the biological neurons and can obtain information from the external environment or other artificial neurons and transmit results to external environment or other artificial neurons.

Artificial neural networks are divided into two learning models (Chiang, 2005; Yeh, 2001):

The supervised learning model is represented by Back-propagation Neural Network (BPNN). The basic theory for this network uses the steepest gradient descent method to minimize error functions. Back-propagation neural networks have functions for learning and parallel information processing which include learning and recalling processes:

The modification process of weights and biases has forward and backward propagations. These are referred to as back propagation neural networks.

Back-propagation neural network model: Back-propagation neural network model is shown in Fig. 1, including (Jiang et al., 2004):

Fig. 1:	Back-propagation neural network model structure

BPNN input layer receives the signal input. The hidden layer and the output layer perform two operations: the linear sum and non-linear conversion of non-linear conversion function. The linear sum process is the product and sum of the first layer input (xi) and the connection weight (wji) plus biases (bi), the result is unit j plus sum uj. The non-linear conversion is used to obtain output value yj of unit j. The back-propagation neural network model flow chart is shown in Fig. 2.

The Equation is as follows (Huang and Huang, 1998):

Air conditioning energy-saving principle: The Mollier refrigerant system recycle diagram is shown in Fig. 3.

Qe	:	Refrigerating capacity
W	:	Compressor output power
C	:	(Performance coefficient):

If reducing condensation temperature and pressure, because:

and Q_e>Q_e>, W'<W, P_C>P_C' So C.O.P' >C.O.P

Fig. 2:	Back-propagation neural network model flowchart

Fig. 3:	Mollier recycle of reducing condensation temperature

System energy consumption decreases when condensation temperature and pressure reduce. Reducing cooling water input temperature can save chiller energy consumption. This study uses neural networks, i.e., back-propagation neural networks in a supervised learning model to train and simulate chiller parameters.

Experimental method: The experimental method employed historical data for the five parameters that affected cooling tower capacity. Cooling tower fan motor frequency, high pressure compressor, cooling water inlet temperature, cooling water pump frequency and temperature differences between cooling water inlet and outlet were used as input parameters (p) and the temperature differences between cooling water input and wet bulb were used as the output variable (t). The back-propagation neural network can be used to analyze this correlation. Further, the historical data from the last day were used to test the effectiveness of neural network modules. The network construction steps were as follows:

Step 1: We used MATLAB to convert the input parameters and output variables into vector RxQ matrix and the converted results were stored for test.

Step 2: Load data with MATLAB working space, normalize input, output objectives and set them to the zero mean and standard deviation to 1.

Step 3: Check converted data, we divided the data into training, verification and testing subsets; one quarter of the total data was used as the data verification set, the other quarter was used as the data testing set and last half was used as the data training set.

Step 4: We used Levenberg-Marquardt to set parameters, such as learning efficiency, learning objective and maximum number of training recycles.

Step 5: We used Levenberg-Marquardt for training and created a network for training. This study tested a two-layer network. This hidden layer within the network contained tangent sigmoidal transfer functions and the output layer had linear transfer function. According to this assumptions, the hidden layer would use five neurons, because there was only a single objective and, therefore, one output neuron was set.

Step 6: The plot error bar for training, verification and testing was used to check training progress.

Step 7: We introduced the data sets into the network and conducted linear regression for network outputs and the corresponding objectives. In this experiment, there was one input and one linear regression required. The historical data of p2 and t2 from the last days were used to normalize the input and the output objectives. By setting them as the zero mean and the standard deviation was 1.

Step 8: We conducted the network simulation analysis for the total data sets and did not normalize the network output.

Step 9: We analyzed and simulated the error between the output value and the actual value and performed a linear regression.

Step 10: We used constructed modules and simulated parameters.

EXPERIMENTAL RESULTS AND ANALYSIS

This study used the actual measured data from air conditioning systems from the research institute and back-propagation network analysis to find out the chiller energy consumption weights and construct the network modules.

Fig. 4:	Network training progress of simulation process

The data from the second day were simulated to test modules and the simulated cooling water pump frequency was used to compare the simulated cooling water inlet temperature from the second day with actual results.

Result of training network modules: Figure 4 demonstrates the rationality of simulated results. The errors of the data testing set and data verification data have similarity. Further, no significant over-fitting occurred. The training process will stop when error verification was performed at 22 times epochs. The scatter plot for the output results is shown in Fig. 5. The five input and output parameters had a close relationship. The R-value represented the correlation coefficient of output and objective and was the measure of verification. If the value was 1, the objective and input had a perfect correlation. In this study, the R-value was 0.986, so the network consisted of the five input and output parameters was feasible. The network modules were used to simulate the data for the next day, sharp climate changes were excluded, because the data of two consecutive days had a correlation. If the R-value was higher, then the modules used the data from the first day to effectively simulate the data from the second day. The results are shown in Fig. 6. If the R-value was 0.853, then there was a slight difference. Nevertheless, the simulated data from the second day was effective.

Simulated temperature difference result using simulated fan motor frequency: The neural network module was stored as pro-model 1. The fan motor frequency and the temperature difference was represented by t and p. After relevant simulation training, the results are shown in Fig. 7 and 8. The training results indicated that the correlation R-value was in the accepted range.

Fig. 5:	Scatter plot of simulated output results

Fig. 6:	Scatter plot of output results using the simulated data of the last day in the training network

The module is stored as pro-model 2. The temperature difference was gradually increased based on the historical data for temperature differences of 1.7°C. The pro-model 2 simulation frequency results are shown in Table 1. The simulated frequency was 58.394Hz when the temperature difference was set to 2.2°C. If increasing the temperature difference, then the simulated value of the water tower fan motor frequency had no significant change. The simulated frequency of pro-model 2 can propagate back to the simulated temperature difference of pro-model 1, as shown in Table 2. The simulated temperature difference was 1.6979°C when the frequency was 56.684 Hz.

Fig. 7:	Scatter plot of simulated output results

Fig. 8:	Network training progress in simulation process

Table 1:	Analysis list of simulated frequency and temperature difference

Table 2:	Mean error

Fig. 9:	Cooling water inlet temperature comparison

If increasing the frequency, then the temperature difference had no significant change. Thus, the fan motor frequency was adjusted to 56.684 Hz when temperature difference was higher than 1.6979°C and increased of the fan motor frequency had no effect. As a result, the cooling water tower motor frequency would not work under full load of 60 Hz; reducing operation frequency can effectively reduce cooling water tower fan motor power consumption.

Adjusting cooling water pump frequency to simulate cooling water inlet temperature: The constructed network modules were used to simulate water pump frequency for the following day. The value was used to estimate cooling water inlet temperature and compared the estimated value to the actual value. Figure 9 shows the Bode diagram for the actual cooling water inlet temperature and the simulated cooling water inlet temperature. From this figure, two values were closer and the error difference was less than 0.5°C. The module had a high feasibility of simulating cooling water inlet temperature. Average error of estimated values of 12 days are shown in Table 2. The results indicated most errors were lower than 0.5°C. Because the climate changes frequently in Taiwan, the effectiveness of the network-based modules based on historical data was decided by climate stability. From Table 2, the fourth day and the tenth day had the greatest average errors. If two days had significant climate changes, then the effectiveness of the simulation prediction was lower. Therefore, network modules for different climates should be constructed.

CONCLUSIONS

In this analysis of air conditioning system power consumption, chillers are of the primary cooling source. Industrial semiconductor and wafer manufacturers require ten thousand ton refrigeration loads with chillers running year round. Air conditioning power consumption is considerable. The cooling water tower efficiency is closely correlated to wet bulbs. As the temperature of wet bulb changes, the fan motor frequency is adjusted to optimize the cooling effect to cool the water tower to maximize efficiency.

This study used back-propagation neural networks to analyze and classify chiller historical data. Present findings showed that the correlation and the feasibility of the simulated results. The cooling tower fan motor frequency can be adjusted to simulate temperature differences between cooling inlet water and wet bulbs. The known impacts of the cooling water tower fan motor frequency and cooling water pump frequency on cooling water inlet temperature. Further control cooling water inlet temperature can optimize the chiller energy efficiency. Future studies can use more sampling points to find greater regularity. In combination with PLC, the control equation can be used in field tests to achieve energy saving efficiency for air conditioning systems.

ACKNOWLEDGMENT

This study was financially supported by the National Science Council under the “Chiller Energy Saving Simulation and Measurement Optimization Research” of the University-Industry Cooperation Plan (Project No: NSC 99-2622-E-167- 005-CC3.). The authors would also like to thank Skytech Science and Technology Engineering Co., Ltd. for its additional funding and technical support.