Abstract
To mitigate the climate change issue, the Government of the Hong Kong Special Administrative Region (HKSAR) proposed the “Hong Kong Climate Change Plan 2030+” with an ambitious goal of reducing the carbon intensity to 65%–70% below 2005 levels by 2030. To achieve this target, the HKSAR Government has developed two District Cooling Systems (DCS) to promote the use of district energy system, enhancing the overall energy efficiency for cooling purpose. One of the government-owned DCSs is located at the Hong Kong-Zhuhai-Macao Bridge Hong Kong Port (HZMB-HKP), producing a total design cooling capacity of 24.16 MW to passenger clearance building and ancillary buildings there. This paper demonstrates how the DCS at HZMB-HKP could enhance overall energy efficiency for cooling purpose, through proper chiller sequencing control, and reduce CO2 emission significantly. An optimal seasonal chiller operation mode is also suggested to be adopted at the DCS at HZMB-HKP from cost-effectiveness perspective.
Introduction
Various DCSs and their installed cooling capacity.
Considering the rapid growth of energy consumption as well as the drastic increase of the carbon dioxide (CO2) emission in the territory, the Government of the Hong Kong Special Administrative Region (HKSAR) has proposed the “Hong Kong Climate Change Plan 2030+” to collaborate with other countries globally, in which, DCS development in Hong Kong is also one of the major measures for the reduction of carbon intensity to 65%–70% below 2005 levels by 2030. 4 The current installed cooling capacity of the DCSs, owned by the HKSAR Government, is about 308.16 MW including 284 MW for Kai Tak Development and 24.16 MW for Hong Kong-Zhuhai-Macao Bridge Hong Kong Port (HZMB-HKP). 5 To further reduce the CO2 emission, a comprehensive analysis of energy performance of the DCS at the HZMB-HKP is conducted based on real operational data and seasonal optimized chiller sequencing control pattern is also recommended for further enhancement of the DCS’ overall energy efficiency.
State-of-the-art DCS development in the HKSAR
There are currently two government-owned DCSs implemented in the HKSAR, namely the DCS at Kai Tak Development Area and the DCS at the HZMB-HKP.
District cooling system at Kai Tak development area
The one at Kai Tak Development Area is the first-ever DCS project in the HKSAR, which has been commissioned since January 2013. There are total of three plant houses including North Plant, South Plant, and seawater pump house. The current installed cooling capacity for this system is around 284 MW. Annual estimated electrical energy consumption of air-conditioning at non-domestic floor area in Kai Tak Development Area was around 243 million kWh with adoption of conventional air-cooled air-conditioning system (AACS) and 197 million kWh with adoption of water-cooled air-conditioning system (WACS), 1 while the energy consumption for DCS was only around 158 million kWh. Instead of adopting traditional AACS and WACS, the annual estimated amount of energy reduction for utilizing DCS was about 85 million kWh (35% reduction) and 39 million kWh (20% reduction), respectively, contributing to 59,500 tons and 27,300 tons of CO2 emission, respectively, using the carbon dioxide emission factor of 0.57kgCO2/kWh. 6 It is worth mentioning that an additional DCS is under construction phase, tailor-made to serve the increasing cooling demand of the extension of the Kai Tak Sport Park and the increase in the scale of the new Acute hospital. 7 Upon full utilization of the additional DCS, it is estimated that an annual amount of 37,000 tons of CO2 can be further reduced. 7
District cooling system at Hong Kong-Zhuhai-Macao Bridge Hong Kong Port
The DCS at the HZMB-HKP is the second DCS in the Territory, commissioned in 2018, designed to serve all facilities in the 130-hectare reclaimed island including Passenger Clearance Building (PCB) and ancillary buildings for various authorities. The simplified layout of the DCS at HZMB-HKP is shown in Figure 1. The design cooling capacity of this DCS is around 24.16 MW. The central chiller plant is located at the basement of the PCB, in which the chilled water is distributed through the underground piping network to the heat exchanger at discrete user buildings scattering at the island. In the meantime, seawater is utilized as a cooling medium for cooling the refrigerant, R134a, in the condensers of the chiller. The seawater is being pumped from the West side of the island to the condensers for refrigerant cooling and then being discharged at the North side of the island as indicated in Figure 1. The detailed system configuration, and specifications for the DCS at HZMB-HKP is highlighted in Configuration and specifications of DCS at HZMB-HKP of this paper. Simplified layout of district cooling system at HZMB-HKP.
Methodology
All data used for analysis and optimization of the chiller sequencing are the real operational data obtained from the Supervisory Control and Data Acquisition (SCADA) system of the HZMB-HKP for this study. The data are logged at 15-min intervals and 1 year of operational data in the period of 2020–2021 was collected for the analysis. Since the SCADA system can only record yearly operational data, including seawater temperature, only the data starting from February 2020 to February 2021 can be obtained for this study. In view of a lockdown due to outbreak of COVID-19 and only government authorities were operated at the HZMB-HKP. Hence, the cooling load calculated by using the operational data can be considered as a base load, without taken occupant load into account. According to the Hong Kong Observatory (HKO), Spring, Summer, Autumn, and Winter can be divided as the period of March–May, June–August, September–November, and December–February of the next year, respectively. 8 For the DCS energy analysis, it is generally accepted to consider the centralized water-side equipment including the condensing water pumps, chillers, and chilled water pumps.9, 10 The equipment considered for the chiller sequencing analysis in this paper includes the seawater pumps, seawater-cooled chillers, primary chilled water pumps, and the main secondary chilled water pumps.
Based on the analysis of the real operational data obtained from the SCADA system at the HZMB-HKP and seasonal weather in Hong Kong, optimal chiller sequencing pattern is then recommended.
Configuration and specifications of DCS at HZMB-HKP
Configuration of the DCS at HZMB-HKP
Configuration of the DCS at HZMB-HKP can mainly be divided into three loops, namely condensing water loop, primary chilled water loop, and secondary chilled water loop. For the condensing water loop, once-through seawater cooling is adopted, in which, the seawater is being pumped from the West side, entering the condensers, and being discharged at the North side of the island as indicated in Figure 1.
There is no heat exchanger provided in between the condensers and seawater inlets and certain filtration and screening processes are required before entering the condenser. Those processes include arresting small debris contained in seawater through the traveling band screens and filtration through automatic backwash strainers. After minimizing the debris contained in the seawater, the seawater enters the condensers for refrigerant cooling purpose.
For the chilled water loops, a constant-primary variable-secondary pumping system is adopted for chilled water distribution in which, the primary chilled pumps run at constant speed to maintain stable flow through the chillers and the secondary chilled water pumps, delivering the chilled water to the heat exchanger of PCB and ancillary buildings, run at various speed based on cooling load demand.
The schematic diagram of the condensing water loop, primary chilled water loop, and the main secondary chilled water loop is indicated in Figure 2. Schematic for condensing water loop, primary chilled water loop and the main secondary chilled water loop for DCS at HZMB-HKP.
Specifications of the equipment of DCS at HZMB-HKP
DCS configuration at HZMB-HKP.
Centrifugal pumps for DCS
There are three types of pumps adopted in the energy performance analysis and system coefficient of performance (SCOP) calculation, namely the seawater pumps (SWP), primary chilled water pumps (PCHWP), and secondary chilled water pumps (SCHWP). Since the SCADA system only takes log of supply water pressure, return water pressure and water flow rate, power consumption of the pumps 10 is evaluated by adopting equation (1), in which N is the electric power consumption in kW, H is the difference between the supply water and return water pressure in m, Q is the water flow rate in m3/s, ρ is the water density in kg/m3, g is the gravitational acceleration in m/s2, and η is the efficiency of the pump. The full-load operating efficiency of the seawater pumps, primary chilled water pumps, and secondary chilled water pumps ranges from 80.0% to 85.1%, which were evaluated during testing and commissioning stage of HZMB-HKP.
Seawater-cooled chillers
There are two types of seawater-cooled chillers (SWCCH) adopted for the DCS, that is, 435RT and another with 1000RT cooling capacity. The Coefficient of Performance (COP) characteristic curve for 435RT and 1,000RT SWCCH are shown in Figure 3. All the COP data are evaluated during testing and commissioning stage of HZMB-HKP. (a) Variation of COP with cooling load of 435RT SWCCH (b) Variation of COP with cooling load of 1,000RT SWCCH Variation of COP with cooling load at different seawater temperature for 435RT SWCCH and (a) Variation of COP with cooling load of 435RT SWCCH (b)Variation of COP with cooling load of 1,000RT SWCCH 1,000RT SWCCH.

Anti-corrosion measures taken for DCS at HZMB-HKP
The DCS at HZMB-HKP utilizes seawater as a cooling medium for refrigerant, in which corrosion is a critical problem that may depreciate the chiller operating performance. Several anti-corrosion measures are implemented on seawater contact parts of the DCS in this case study which are stated as follows: The condenser tubing of all chillers at HZMB-HKP is made of titanium which has excellent corrosion-resistant property Zinc metal is installed inside the surface of the water box of condensers of all chillers at HZMB-HKP to further enhance the anti-corrosion performance Fouling is a common problem in seawater condensing loop. Salt deposition and other impurities might easily build up inside the condenser tubes. From time to time, accumulating scaling and fouling will highly affect the performance of the condensers. An automatic tube cleaning system is integrated to the condensers of all chillers at HZMB-HKP to cater the problem of fouling. The basic principle is to circulate certain amount of soft sponge balls in the condensing water loop through the condenser tubes to minimize the formation of fouling and scaling, slowing down the depreciation of seawater condensers. Preventive maintenance is also carried out monthly, quarterly and yearly for the condensers and relevant parts to ensure a proper operation of condensing water system. The works include condenser tube and water box cleaning, coating repainting, replacing zinc anodes, and making good joint of condenser covers. The seawater contact parts of the seawater pumps are made of various types of stainless steel. The volute casing, impeller, shaft is made of duplex steel, and the shaft protecting sleeves are made of cold rolled (CR) steel. Both the duplex steel and CR steel are renowned for corrosion-resistant in the industry. There is anti-corrosion protection applied to the pipelines and fitting in the seawater piping network at HZMB-HKP. Cement mortar lining is also applied to all the straight pipes, with all the fittings incorporated with fusion epoxy coating, to minimum the possibility of corrosion.
Energy performance indicators for DCS
The following three indicators are generally accepted for assessing general energy performance of the DCS.
Coefficient of Performance of chiller plants
The Coefficient of Performance of chiller plants (COPPlant) has been widely used for assessing energy performance of the chiller systems and pumps, that is, SWCCH, PCHWP, and SWP. The governing equation is stated in equations (2) and (3), in which, CC is the instantaneous cooling demand in kW, WPlant, WSWCCH, WPCHWP, and WSWP are the power consumption of DCS plant, seawater-cooled chillers, primary chilled water pumps, and seawater pumps in kW, respectively. 9 The COP curve of different seawater temperature for 1,000RT SWCCH and 435RT SWCCH are indicated in Figure 3.
Water transfer factor of network distribution
The Water Transfer Factor of the cooling network distribution (WTFDistribution) is adopted for evaluating the energy performance of the cooling distribution network in the chilled water system, in which, CC and WDistribution represent instantaneous cooling demand in kW and power consumption used for distributing the chilled water (SCHWP) in kW, respectively. 9 The governing equation is stated in equation (4).
System coefficient of performance
The System Coefficient of Performance (SCOP) is the overall energy performance indicator for DCS, in which the CCEffective is the instantaneous cooling capacity accounting for the cooling loss of pipework, with governing equations stated in equation (5). The energy consumption of chilled water pumps will eventually be converted into heat gains in the water flow during distribution process and 80% of the pump energy is usually assumed.
9
By combining the equations (2) to (6), the SCOP for assessing the energy performance of DCS is indicated at equation (7)
Results and recommendations
Annual cooling load and distribution of power consumption
The peak cooling load of each month is indicated in Figure 4, calculated by using equation (8), in which, m is the mass flow rate of the chilled water through the SWCCH, cp is the specific heat capacity of water (4.198 kJ/kg°C), TCHWR and TCHWS are chilled water return temperature (°C) and chilled water supply temperature (°C), respectively. Since a control of chiller sequencing is only based on the peak cooling load of each season, only seasonal peak cooling load is indicated in Figure 4. The peak cooling load for Spring (March–May 2020), Summer (June 2020–August 2020), Autumn (September 2020–November 2020), and Winter (December 2020–February 2021) are 7,596 kW, 13,027 kW, 9,488 kW, and 3,056 kW, respectively. The distribution of the power consumption for the annual maximum cooling load extracted from the SCADA system is shown in Figure 5, in which the SWCCH, SWP, PCHWP, and main SCHWP accounts for 62% (2,378 kW), 9% (331 kW), 14% (559 kW), and 15% (594 kW) of total power consumption for the DCS at HZMB-HKP, respectively. Annual peak cooling load profile of the DCS at HZMB-HKP. Distribution of power consumption at annual peak cooling load of HZMB-HKP.

Theoretical optimal chiller sequencing control at HZMB-HKP
Since the seawater-cooled chillers account for around 62% of the total power consumption of DCS, an optimal chiller staging strategy can effectively enhance the overall energy efficiency for the HZMB-HKP. Chiller sequencing control is then proposed.
Taking into account the average seawater temperature for Spring, Summer, Autumn, and Winter of around 24°C, 29°C, 27°C, and 19°C, respectively, the first option is sole operation of the chiller with 1000RT capacity (hereafter called “baseline”). The second option is to deal with the cooling demand with both 1,000RT and 435RT chillers (hereafter called “S+B”). If the cooling capacity of the 435RT chillers is insufficient, one more 1,000RT chiller will be operated and so on. The third option is almost the same as the second pattern, in which only difference is maximum cooling capacity of each chiller (hereafter called “S+B (MAX COP)”). The maximum cooling capacity of both 435RT and 1,000RT chillers which are set according to the maximum COP, indicated in the COP curve in Figure 3, (i.e., when the maximum cooling capacity with the highest COP of one chiller is reached, one more chiller will then be operated to cater the cooling demand). The 435RT chiller will also be put into operation before the start of the 1,000RT chillers. The minimum output (Turn-down ratio) of all chillers of DCS at HZMB-HKP is 25% of full cooling capacity as provided by the manufacturer. This study of “Baseline,” “S+B,” and “S+B (MAX COP)” chiller sequencing control have taken the turn-down ratio into account. Figure 6 and Figure 7 show the operating pattern for 435RT chiller and 1,000RT chiller, respectively. The starting point of each line represent the minimum cooling output of the operating chiller and the region between dark-blue dotted lines indicate the number of chillers that are operating and share equal amount of cooling load. All the above-mentioned three configurations in this study have taken the minimum turn-down ratio into account, in which, each chiller is operating above the 25% of the full cooling capacity. Minimum turn-down ratio and operation pattern of 435RT Chiller for different cooling load distribution of DCS at HZMB-HKP. Minimum turn-down ratio and operation pattern of 1000RT Chiller for different cooling load distribution of DCS at HZMB-HKP.

The SCOP of each chiller system combination is evaluated with polynomial regression model, which can be used to indicate the approximate SCOP for each cooling load. The chiller sequencing control for DCS is formulated based on the actual power consumption of each combination of chiller systems for all four seasons.
Chiller sequencing control for spring
The peak cooling demand in Spring was about 7,596 kW. The power consumption and the related SCOP of the chiller operation configurations for Spring are indicated in Figure 8 and Figure 9, respectively. In Figure 8, it shows that the configuration “S+B (Max COP)” performs the best theoretically as it had the least power consumption in Spring, and the highest SCOP of this configuration can be up to 4.2. Hence, the combination of operating both 435RT and 1,000RT chillers is recommended for Spring. Power consumption of different chiller operation pattern in Spring SCOP of different chiller operation patterns in Spring

Adopting this theoretical chiller sequencing control, the total energy consumption was about 1,505MWh, contributing to a reduction of 8.5% energy consumption when comparing with the “Baseline” configuration for Spring. The related CO2 reduction was around 80 tons, with the territory-based CO2 emission factor of 0.57 kgCO2/kWh. 11
Chiller sequencing control for Summer
The peak cooling demand in Summer was about 13,027 kW. The power consumption and the related SCOP of the chiller operation configurations are indicated in Figure 10 and Figure 11, respectively. In Figure 10, for the cooling load ranging from 0 to 5,000 kW, configuration “S+B (Max COP)” and “S+B” combination have a similar power consumption, which was lower than that of the “Baseline” operation mode. However, when the cooling load increased to 5,000 kW or above, the power consumption for the “S+B (MAX COP)” chiller system configuration kept the lowest and the related highest SCOP along that region was about 3.6. It is recommended to adopt the “S+B (MAX COP)” configuration in Summer theoretically. The total energy consumption of adopting this energy-saving sequencing control was about 7,380MWh, accounting for an energy reduction of 7% comparing with the “Baseline” operating mode. Power consumption of different chiller operation pattern in Summer SCOP of different chiller operation patterns in Summer

Chiller sequencing control for Autumn
The peak cooling load in Autumn was 9,488 kW. The power consumption and the related SCOP curve of the chiller operation configurations are indicated in Figure 12 and Figure 13, respectively. From Figure 12, it can be observed that the “Baseline” chiller configuration consumed the least power for most of the cooling load range in Autumn. Power consumption of different chiller operation pattern in Autumn SCOP of different chiller operation patterns in Autumn

Although the “S+B” sometimes overlap with the “Baseline” chiller operation mode, Figure 13 shows that the Baseline mode has a better performance regarding the overall system energy performance including the relevant SWP, PCHWP, and SCHWP, having largest SCOP of 4.0. Therefore, it is recommended to adopt the “Baseline” operation mode in Autumn. The total energy consumption for utilizing this strategy was around 4663 MWh with CO2 emission of about 3264 tons. Comparing with the “S+B” configuration, 7.7% (205 tons) of CO2 emission could be reduced.
Chiller sequencing control for Winter
The peak cooling load for Winter was 3,056 kW. The power consumption and the related SCOP of the chiller operation configuration are indicated in Figure 14 and Figure 15, respectively. In Figure 14, the “S+B (MAX COP)” performed the best, having the least power consumption in almost all cooling load range in Winter. In Figure 15, the highest SCOP of this control strategy was around 3. Hence, the “S+B (MAX COP)” is recommended to be adopted for Winter theoretically. The related energy consumption is about 650 MWh, contributing to a reduction of 49 tons of CO2 emission. It should be noted that the SCOP of the chiller operation configurations in Winter was the least among Spring, Summer, and Autumn in the year, which may be highly related to low part-load ratio of the chiller operation. Power consumption of different chiller operation pattern in Winter SCOP of different chiller operation patterns in Winter

Recommendations for an optimal chiller sequencing control based on different design criteria
Cases considered for energy efficiency and cost-effectiveness analysis.
Optimal chiller sequencing control base on energy efficiency
Figure 16 showed the total energy consumption of the studied chiller sequencing control for each season. As the most energy-efficient operational mode, it is proposed to adopt configuration number 4, that is, “S+B (MAX COP)” in Spring, Summer, and Winter and “Baseline” in Autumn, leading to an annual potential energy saving 785,316kWh. Seasonal energy consumption for each chiller sequencing control.
Optimal chiller sequencing control base on cost-effectiveness
In this section, a simple payback calculation for different chiller combinations is performed to figure out the most cost-effective chiller sequencing control. The installation cost of a SCADA system for the DCS at HZMB-HKP is around HKD15 million, while the potential cost for reprogramming and interfacing SCADA system to chillers would be around HKD3 million on top of the original installation cost of SCADA system. The average electricity tariff rate in the studied period was around HKD1.22/kWh.
Net energy saving of various configurations comparing with “Baseline” control of chiller in kWh
Positive (+): The configuration consumed less energy than that of “Baseline” Negative (−): The configuration consumed more energy than that of “Baseline”
Summary of payback period of four studied configurations
Summary of recommended chiller sequencing configurations
Percentage energy reduction and recommendations of chiller sequencing control for each season regarding the energy efficiency and cost-effectiveness
Remark: Negative (−) = consume less energy comparing with baseline mode Positive (+) = consume more energy comparing with baseline mode
Summary of seasonal optimal energy efficient chiller sequencing control for DCS at HZMB-HKP
Remark: S = 435RT chiller, B = 1,000RT chiller
Summary of seasonal optimal cost-effective chiller sequencing control for DCS at HZMB-HKP.
Remark: S = 435RT chiller, B = 1000RT chiller
In fact, there are certain constraints for this study. The regular scheduled preventive maintenance for the SCADA system may lead to data loss for some period of time. Also, as the SCADA system does not log the instantaneous power consumption data directly, corresponding power consumption for all HVAC pumps and chillers, it can only be calculated based on theoretical equations and assumptions from academic thesis. Hence, it would be better to conduct retro-commissioning individually for each premise to tailor-made the most technically feasible and cost-effective energy saving solution.
Proposed configurations in different dimensions
Optimal chiller sequencing control is one of the most essential components in both design and operation and maintenance stage to minimize the CO2 emission from reducing energy consumption. It is recommended to apply chiller configuration 4, adopting “S+B (MAX COP)” for Spring, Summer, and Winter and “Baseline” for Autumn to the DCS at HZMB-HKP if the aim of the design is to achieve higher system energy efficiency. With consideration shortest payback period, it is suggested to adopt configuration 3 with payback period of 18 years, utilizing “S+B (MAX COP)” for Spring, Summer, and Winter and “Baseline” for Autumn for the DCS at HZMB-HKP. The payback calculation was performed based on data collected during the COVID-19 lockdown and it is anticipated that the annual energy cost saving will increase in the post-COVID-19 period as the higher PLR of the chillers contribute to a higher COP. The whole study demonstrates a simple method that assist the system designers with the decision making on various chiller sequencing controls with different design criteria.
Practical application on energy assessment of various chiller sequencing design for new DCS development
District energy systems are a popular means of enhancing regional energy efficiency. This study demonstrates a simple approach for evaluating energy performance of various chiller sequencing control options which can assist operators to determine optimal operation patterns of chillers base on different operational criteria. The following items should be considered in this study: This case study has not considered the factor of depreciation of chillers in the energy performance evaluation and SCOP calculation. It is recommended to use the instantaneous data of power consumption with the polynomial regression model and equation to predict the operating SCOP as the factor of depreciation can be reflected by SCOP from time to time. As the instantaneous data of power consumption have not been logged by the SCADA system, it is not technically feasible to obtain a practical SCOP curve with different sea water temperatures. It is recommended to group the data points with same condensing temperature and set up a COP curve for that particular temperature in the future DCS development. It is also feasible to formulate an accurate and practical chiller sequencing strategy for different condensing temperatures through a regression curve of the operating condensing temperatures. The DCS pumping systems at HZMP-HKP consists of a constant-primary pumping system and a variable-secondary pumping system. Primary chilled water pumps are required to maintain a stable chilled water circulation through evaporator to avoid freezing problem, it is considered as the WPlant in the SCOP calculation. For DCS utilizing only one type of pumping system, such chilled water pumps should only be considered as WDistribution for the SCOP calculation as its aim is to deliver chilled water to the terminal units. The seawater cooling system at the condensing water loop is a once-through system. In case energy evaluation is to be conducted on a DCS adopting an indirect seawater cooling system for condensing water loop, the power consumption of the circulating pumps between condensers and seawater heat exchangers should be considered in the SCOP calculation. Thermally driven cooling system (absorption cooling), utilizing waste heat or heat from renewable energy sources, is also commonly adopted for the DCS. The grades of energy indicate the order state of molecules stored in these forms of energy. The energy sources for an absorption cooling system are considered as low-grade energy, while energy sources for an electrically driven cooling system is a kind of high-grade energy. Because of different energy grades and sources for absorption and electric cooling system, the energy input cannot be directly adopted for SCOP calculation and determining the chiller sequencing control. Hence, the energy evaluation method in this case study cannot be adopted for the DCS utilizing absorption cooling system including combined cooling, heat and power application
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
