Abstract

Based on the life cycle cost (LCC) approach, this paper presents an integral mathematical model and particle swarm optimization (PSO) algorithm for the heating system planning (HSP) problem. The proposed mathematical model minimizes the cost of heating system as the objective for a given life cycle time. For the particularity of HSP problem, the general particle swarm optimization algorithm was improved. An actual case study was calculated to check its feasibility in practical use. The results show that the improved particle swarm optimization (IPSO) algorithm can more preferably solve the HSP problem than PSO algorithm. Moreover, the results also present the potential to provide useful information when making decisions in the practical planning process. Therefore, it is believed that if this approach is applied correctly and in combination with other elements, it can become a powerful and effective optimization tool for HSP problem.

1. Introduction

Humanity faces serious energy and environment problems at present. The environment is increasingly threatened. For instance, with the increase of greenhouse gas emissions in the atmosphere the environments have already reached concerning levels in terms of their potential to cause climate change. Air pollution, acid precipitation, and stratospheric ozone depletion are other serious environmental concerns. The severity of climate change impacts shows the increasing trend if significant action is not taken to reduce greenhouse gas emissions [1]. An important action to address energy and environmental challenges lies in the intelligent and efficient use of energy, including reducing energy waste and using low-carbon fuels.

In China, heating utilities have been developed rapidly, but the energy consumption of production and transport is still too much, which accounts for 21.5% of building energy consumption; building energy consumption accounts for 20.9% of social total energy consumption [2]. With the perfection of the systematic reform, the adjustment of energy structure in China and the requirement of environmental protection, heating energy structure had been changing, and it had been promoting the development of heating mode. It has very important significance to analyze, evaluate, and select heating mode correctly which suits its local characteristics. With the speeding up of urbanization, more and more heating systems will be built due to the importance of infrastructure in urban area. The research on optimal plan of heating system is very imperative for saving project investment, decreasing heating energy consumption, and improving enterprise benefit.

Sustainable development of heating system requires application of planning procedures, which includes optimization of both demand and supply sides of heating. Because the heat source site selection and heating pipe network optimizing plan have an important role in the HSP, there are many scholars concerning this subject and lots of optimization methods have been proposed. The methods of HSP can be classified into three separate categories [3]: planning by models, planning by analogy, and planning by inquiry. The planning by models can be based either on econometric or optimization models. Econometric models utilize mathematical or statistical methods and relay on statistical data. Optimization model allows for the identification of best possible solution—minimization or maximization of objective function, with the predefined set of constrains which describes the space of acceptable solutions. The planning by analogy utilizes the simulation of heating system. That kind of HSP is usually used for the verification of planning results which were achieved by other planning methods [4]. The planning by inquiry is used in the case when other aforementioned methods are not reliable. Good example of planning by inquiry is DELPHI method, which is based on the questioning of group of heating, ventilating, and air conditioning (HVAC) experts or municipal planners and statistical evaluation of their answers [5]. All the methods of HSP listed earlier have a limited transparency, especially for decision makers who do not have good mathematical background. Those methods do not give opportunity to create decision makers preference model or define that model a priori. Hence, many scholars have carried out extensive and deep research on optimization method of HSP. Shen and He [6] investigated optimal planning method of central heating system of water boiler and then put forward optimal planning model and solving method. According to the method, it can be determined the size, location, and so forth of regional heating plant and intermediate heat exchanger station, but there were no further discussions about how to design the pipe network. Wang et al. [79] investigated design method of central heating system with double “duct-station,” proposing two-step optimization method, but this method was only applicable to double “duct-station” system. For solving the problem, the study in [8] used fully stratified sequence method [1012] simultaneously taking into account only one objective function for each layer of heating source layout optimization and pipe network system. However, for HSP problem, there are kinds of complex logic even iterative relationships between/among objective functions of layers. This method usually can provide the optimal solution of each layer, but it cannot ensure that the solution of the objective function for last layer is the optimal solution of the whole system exactly. Shi and Li [13] first applied genetic algorithm (GA) for solving the heating source location problem in the study. This method described the cost of the heat source and heating substation as the function of heat load and described the cost of the heating network as the function of heat load and pipe length simply. So the calculated result by this method and actual situation often put in certain error. Shi et al. [14, 15] and Mu et al. [16] put forward the relatively consistent mathematical model for heating system optimization, based on the life cycle cost method, but formula or method for some of parameters was not given clearly and integrally in the model. It seems that limitation is inevitable in the process of the practical application of these methods. But we noticed the life cycle cost (LCC) and particle swarm optimization (PSO) algorithm in the more extensive research areas.

Life cycle cost (LCC) has been applied since the 1960s when the United States’ Department of Defense stimulated the development and application of LCC to enhance its cost effectiveness. Defense systems, such as an aircraft or a special land vehicle, are ideal for LCC analyses since the Department of Defense mainly controls the entire life cycle [17]. LCC may be defined as “the cost of acquisition, ownership, and disposal of a product over a defined period of its life cycle” [18, 19]. LCC is a standard engineering economic approach used for choosing among alternative products or designs that approximately provide the same service to the customer [20]. In many cases it may not be necessary to perform a complete LCC analysis, but rather to estimate the differences between the alternatives for the major cost elements [21]. The LCC process may also provide information, for example, in the assessment of the economic viability of products and projects, in the identification of the cost drivers and cost efficiency improvements, and in evaluations of different strategies for product operation, maintenance, and inspection, and so on [22].

There are two popular swarm inspired methods in computational intelligence areas: ant colony optimization (ACO) and particle swarm optimization (PSO). ACO was inspired by the behaviors of ants and has many successful applications in discrete optimization problems. The particle swarm concept originated as a simulation of simplified social system. The original intent was to graphically simulate the choreography of a bird block or fish school. However, it was found that particle swarm model can be used as an optimizer. A substantial review of the properties of the global optimization problems has been given by Parsopoulos and Vrahatis [23]. As one of the global optimization problems, PSO has been widely used in various kinds of planning problems, especially in the area of substation locating and sizing [2427]. But in area of heating supply, PSO is mainly applied in heating load forecasting [28, 29], but rarely used in HSP.

The main objective of this paper is to discuss the usefulness of the PSO algorithm for solving the HSP problem. Therefore, based on the LCC approach, an integral mathematical model is presented and PSO algorithm is introduced and improved for solving the problem. In the end, the results of the case study suggest the effectiveness of improved particle swarm optimization (IPSO) application to the optimal planning method for heating system.

2. Mathematical Formulation

2.1. Problem Definition and Assumptions

LCC is related to the systems engineering process, because economic considerations are very important in the process of creating systems. Life cycle economic analyses should be done early in the system or product life cycle, because the outcome of the systems engineering process cannot be influenced very much when the design is completed. Thus, LCC involves evaluation of all future costs related to all of the phases in the system life cycle including design, construction and/or production, distribution, operation, maintenance and support, retirement, and material disposal, and so on [30].

Cost models may range from simple to complex and are essentially predictive in nature. Parameters, such as the system’s physical environment, usage demand, reliability, maintainability, labor, energy, taxes, inflation, and the time value of money, may have a great influence on the life cycle costs [17].

The main objective of this paper is to discuss the usefulness of the PSO algorithm for owners in making sustainable heating system investment decisions and to improve their decision-bases for municipal administration. Therefore, we apply LCC approach to describe the HSP problem.

Moreover, HSP considered in this study works under the following definition and assumptions.(i)A heat consuming installation can connect with any heat source but cannot connect with two or more heat sources at the same time.(ii)The indirect connection between heat consuming installation and heat source is not allowed.(iii)A heat source must be connected with more than one heat consuming installation; otherwise, it will be closed.(iv)Any connection between any two heat sources is not allowed.(v)The location of heat consuming installation is fixed.(vi)A heat source can be sited in a given region.(vii)The elevation difference between heat consuming installation and heat source is ignored.(viii)Heating system planning and optimization can be achieved by changing the number and the heating capacity of heat source and the distance between the heat source and heat consuming installation.(ix)The measure between heat source and heat consuming installation is simplified to the Manhattan (or city block) distance.(x)There is no functional difference between any two heat sources and their products.

2.2. Notation

The notations used in the mathematical formulations are given as follows.

IndicesOptional heating sourceHeating equipmentHeat consuming installationHeat load distributing segment.

Parameters Number of heat source Number of heating equipment which could be installed at the heating source ; Number of heat consuming installation Number of heat load distributing segments Life cycle fixed cost of the heat source Life cycle fixed cost of the heating equipment which is in the heat source Variable production and transport discounted costs within life cycle of heating equipment to satisfy the heat load distributing segment of heat consuming installation , which is in the heat source ; , where is the variable production discounted cost within life cycle of specific heat load; is the transport discounted cost within life cycle per specific heat load Continuous variable, the load of the heat load distributing segment of heat consuming installation , which is supplied by heating equipment of the heat source Load of the heat load distributing segment of heat consuming installation Maximum supply capacity of heating equipment of the heat source Maximum supply capacity of heat source Major repair depreciation discounted costs within life cycle of heat source Labor discounted cost within life cycle of heat source Coefficient of sum; , where is the standard discount rate, and is the life cycle Price of fuel Calorific value of fuel Thermal efficiency of heat source Water and electricity consumption costs of specific heat load Duration of heat load distributing segment Sulfur content in fuel Standard emission charge for Pipe network discounted cost per specific heat load, which is supplied by heat source to heat consuming installation The discounted cost of pipe segment The major repair depreciation discounted costs within life cycle of pipe segment The heat loss discounted costs within life cycle of pipe segment The heat load-bearing of pipe segment The power consumption discounted cost within life cycle of pipe segment per specific heat load The transport discounted cost within life cycle per specific heat load, which is supplied to heat consuming installation The major repair depreciation discounted costs within life cycle of heat consuming installation The heat load of pipe network for heat consuming installation The power consumption discounted cost within life cycle of heat consuming installation per specific heat load The labor discounted cost within life cycle of heat consuming installation Investment of meters diameter double-pipe per meter length, where and are coefficients of pipe laying Rate of major repair depreciation Rate of gross fixed capital formation Conversion coefficient of the units Specific frictional resistance Equivalent length of local resistance for pipe segment Heating period Electricity price for industrial uses Efficiency of circulating water pump Supply/return water temperature of pipe segment Conversion coefficient of the units Heat transfer coefficient Local heat loss coefficient of pipe fittings Annual costs of heat loss Annual mean supply water temperature of pipe segment Annual mean return water temperature of pipe segment Annual mean temperature Comprehensive coefficient of investment Heat load of pipe network Correction factor Conversion coefficient of the units Pressure difference between supply and return water of pipe network for heat consuming installation Average annual wages of operating personnel and manager Number of operating personnel and manager per 1 MW heat load Conversion coefficient of the units.

Decision Variables 1, if the equipment is installed or set up in the heat source ; 0, if the equipment is not installed or set up in the heat source 1, if the heat source is set up; 0, if the heat source is not set up.

2.3. Mathematical Model of HSP

In this study, the problem is summarized into a multisource, multifacility, single-commodity, multiraw material plant location problem, and a mixed 0-1 integer planning model has been formulated. The cost model of the heat source and the heat-transmission network concerned in the optimization model are considered in this study. The objective function of heating system planning problem is to minimize the total heat production cost. The proposed mathematical model formulation for HSP problem can be found as follows.

Minimize subject to Objective function (1) minimizes the discounted costs within life cycle of heating system as the general objective; it is an index of dynamic economy evaluation, where , , and are composed of respective discounted costs together. Constraint (2) is each heat consuming installation’s heat load, which is heat consumption for each user and the requirements of the heating quantity and quality. Constraint (3) means that each of the heating equipment in the heating system bear heat load should not exceed the maximum heating capacity. Constraint (4) means the maximum heating capacity of heating source, which is allowed under the restrictions of objective conditions. Constraint (5) means that only open heating source first can install equipment in it. In the model, there are two decision variables, in which is related to heating source, and is related to heating equipment.

Because the piecewise function of heat load duration curve is introduced in the process of solving the model, this model can be applied to any form of heating system.

2.4. Formulation of Heating System Cost Model
2.4.1. The Heating Source Cost Model

The heating source cost model is aimed to resolve the calculation problem of , in the objective function (1), and , which is a part of in the objective function (1). It consists of fixed costs and variable costs, the former refers to all necessary costs of heating source, so long as open a heating source or install a piece of heating equipment, the latter only associated with the size of the heat load and running status of equipment. Consider Equation (7) is the formulation of heating source fixed costs, and (8) is the formulation of heating source variable costs.

2.4.2. The Heating Network Cost Model

The heating network cost model is aimed to resolve the calculation problem of , which is a part of in the objective function (1), and also to optimize the direction of heating network and the pipe diameter. Heating network (heat consuming installation included) cost consists of the heating network operation cost and heat consuming installation costs. Heating network operation cost consists of major repair depreciation discounted cost, power consumption discounted cost, pipe network heat loss discounted cost, and labor discounted cost. By dividing the discounted cost within life cycle of pipe segment allocation to the total heat load bearded by itself directly and evenly, the transport discounted cost within life cycle per specific heat load can be obtained, which is supplied by heat source to heat consuming installation . Consider Equation (9) is the heating network transportation cost model. The pipe segment cost model is composed of (10)–(14), where (10) is the pipe network discounted cost per specific heat load supplied by heat source to heat consuming installation ; (11) is the investment cost of pipe segment ; (12) is the major repair depreciation discounted cost of pipe segment ; (13) is the power consumption discounted cost of pipe segment per specific heat load; and (14) is the heat loss discounted cost of pipe segment . Consider The heat consuming installation cost model is composed between (15) and (19), where (15) is the transport discounted cost within life cycle per specific heat load supplied by heat consuming installation ; (16) is the investment cost of heat consuming installation ; (17) is the major repair depreciation discounted cost of heat consuming installation ; (18) is the power consumption discounted cost of heat consuming installation per specific heat load; and (19) is the labor discounted cost of heat consuming installation .

3. PSO and Its Improvement

3.1. PSO Algorithm

The PSO is proposed by Kennedy and Eberhart [31, 32] in 1995, and the motivation for the development of this algorithm was studied based on the simulation of simplified animal social behaviors, such as fish schooling and bird flocking. Similar to other population-based optimization methods such as genetic algorithms, the particle swarm algorithm starts with the random initialization of a population of particles in the search space [33]. However, unlike in other evolutionary optimization methods, in PSO there is no direct recombination of genetic material between individuals during the search. The PSO algorithm works on the social behavior of particles in the swarm. Therefore, it provides the global best solution by simply adjusting the trajectory of each individual toward its own best location and toward the best particle of the entire swarm at each time step (generation) [31, 34, 35]. The PSO method is becoming very popular due to its simplicity of implementation and ability to quickly converge to a reasonably good solution.

3.2. Formulation of General PSO

Specifically, PSO algorithm maintains a population of particles, each of which represents a potential solution to an optimization problem. The position of the particle denotes a feasible, if not the best, solution to the problem. The optimum progress is required to move the particle position in order to improve the value of objective function. The convergence condition always requires setting up the move iteration number of particle.

The position of particle move rule is shown as follows: where represents the velocity vector of particle in time; represents the position vector of particle in time; is the personal best position of particle ; is the best position of the particle found at present; represents inertia weight; , are two acceleration constants, called cognitive and social parameters, respectively; and and are two random functions in the range .

The flow chart of general PSO is shown in Figure 1.

3.3. Improvement of Particle Swarm Optimization (IPSO) for HSP Problem

For HSP problem and its model in this paper, the value of LCC depends mostly on the distance between heating source and heat consuming installation, and the number of heating source . It is necessary to make corresponding improvements on PSO, in order to solve this problem more accurately and effectively.

The evolution of the solution set begins with an initial solution set in the PSO; initial solution set is composed of initial particles. Each solution location is represented by an -dimensional vector; represents the number of variables of each solution, and it represents the number of heating sources in particularly in this paper.

The position coordinate of heating source has two components, which is represented by two -dimensional vectors, where direction coordinates are represented by vector , and direction coordinates are represented by vector . Therefore, direction component for the position vector of particle in time can be represented by , and the rest can be done in the same manner.

In the same way, the velocity for location change of heating source has two components, which is represented by two -dimensional vectors, where direction component for the velocity vector is represented by vector , and direction component for the velocity vector is represented by vector . Therefore, direction component for the velocity vector of particle in time can be represented by , and the rest can be done in the same manner.

Thus, the update rule of velocity for each particle is indicated by (22)-(23), and the update rule of position for each particle is indicated by (24)-(25). Consider

The meanings of parameters are consistent with previous description.

3.4. Calculated Flow of  IPSO

The calculated flow of proposed IPSO is described as follows.

3.4.1. Initial Solution

The initial solution for HSP problem is obtained by random initial position of each heat source; a matrix is employed in recording the coordinates and the heat load-bearing information of heat source, and the calculated flow of initial solution is as follows. (1)Set up the number of heat source , and generate an empty matrix for the initial position of heat source.(2)Based on randomly and evenly distributed manner, generate the position coordinates of heat sources, into the matrix.(3)Call the decoding function; calculate the heat load-bearing and the cost for each heat source, into the matrix.(4)Calculate the LCC, the fitness value of the initial particle.

3.4.2. Decoding Function

In this paper, decoding function will call the matrix for current position and heat load of heat consuming installation, and then according to the matrix for the position of heat source, which is represented by current particle, divide the heating range of each heat source, and calculate the LCC.

Information matrix of heat consuming installation (heat_point) is a -line four-column matrix; the first column represents the serial number of heat consuming installation, the second column represents the coordinate of heat consuming installation, the third column represents the coordinate of heat consuming installation, and the heat load of heat consuming installation is represented by the fourth column. The calculated flow of initial solution is as follows.(1)Read matrix heat_point, and let .(2)Calculate the distance to all heat source from the heat consuming installation , into the vector .(3)By substituting into (10)–(14), calculate the cost of the heat consuming installation connected with each heat source.(4)Find out the minimum cost, and the heat consuming installation connected with the corresponding heat source.(5)If is the last heat consuming installation then stop; otherwise, go to Step 1.

3.4.3. The Evolution of Particle Swarm

After one generation of particles, a new generation is evolved as follows.(1)Call the decoding function; calculate the fitness value of the particle swarm.(2)Update the individual optimal solution and the global optimal solution .(3)Update the speed vector, by using (22)-(23).(4)Update the speed vector, by using (24)-(25).

3.4.4. Improvement Approach

The PSO’s convergence is fast, so it is liable to fall into local optimal solution. In order to improve the optimizing capability, we add modular arithmetic of velocity vector into each iterative operation. If the norm of velocity vector is less than the predetermined minimum value , then generate a random velocity, let the current particle swarm out of local convergence region, and search other solution spaces. However, after it falls into local optimal solution, the norm of velocity tends to 0 in the general PSO algorithm, the solution stabilized near the local optimal solution, and it cannot explore search space furthermore.

The flow chart of IPSO is shown as Figure 2.

4. Case Study

4.1. Basic Information of Case

This is a heating plan for a new area in China covering the area of 3.346 million square meters, and heat load is 167.3 MW in total. Based on the road network, the new area is divided into 29 heating districts (Figure 3), and the heating load of each district (Figure 4) is supplied by their small gas-fired boiler.

4.2. The Parameters of Algorithms

The role of the inertia weight , in (20), (22), and (23), is considered critical for the PSO’s convergence behaviour. The inertia weight is employed to control the impact of the previous history of velocities on the current one. Accordingly, the parameter regulates the trade-off between the global and local exploration abilities of the swarm. A large inertia weight facilitates global exploration, while a small one tends to facilitate local exploration. A suitable value for the inertia weight usually provides balance between global and local exploration abilities resulting in a reduction of the number of iterations required to locate the optimum solution. Initially, the inertia weight was constant. However, experimental results indicated that it is better to initially set the inertia to a large value, in order to promote global exploration of the search space, and gradually decrease it to get more refined solutions [32, 36]. Thus, an initial value around 1.2 and a gradual decline towards 0 can be considered as a good choice for . The parameters and , in (20), (22), and (23), are not critical for PSO’s convergence. However, proper fine-tuning may result in faster convergence and alleviation of local minima. A further study of the acceleration parameter in the first version of PSO is given in [37]. As default values, were proposed, but experimental results indicate that might provide even better results. Some work reports that it might be even better to choose a larger cognitive parameter, , than a social parameter, , and [38, 39], but was suggested by [35]; the strategy of acceleration parameter linear changing with iterations was proposed by Ratnaweera et al. [40], but acceleration parameter is the nonlinear function of the ratio -to-, which was proposed by Arumugam et al. [41]; Jie et al. [42] suggested to adjust the acceleration coefficient by measuring diversity.

But so far, the research on the most appropriate values for , , and has no unified conclusion. And how the variable values impact the solution to HSP problem is unknown. For HSP problem on kinds of values is unknown. So we set the , , and to common values in this study.

4.3. Analysis of Results

By applying PSO and IPSO algorithm, respectively, we solved the HSP problem in this paper. The parameters of PSO and IPSO are summarized in Table 1.

In this study, 29 kinds of schemes of heating (from one heat source to twenty-nine heat sources) were calculated for 10 times through reading initial conditions from the excel file successively, which contains the coordinates and heat load of heat consuming installation, preset maximum number of heat source. The results of LCC and the -value for the optimum and the average between PSO and IPSO at the same number of heat source are shown as Figures 5 and 6 and Table 2.

From analyzing the results, we can draw the following conclusion about the HSP problem.(1)The original plan (the heating load of each district is supplied by its small gas-fired boiler) is not an economic and reasonable plan for the case, and the LCC is the second highest in 29 schemes, which is only better than the scheme which plans to set up one heat source only.(2)From one heat source to twenty-nine heat sources, LCC is monotone decreasing until a minimum value first, then monotone increasing.(3)Only one minimum value of LCC that appeared throughout the change process, which is 1.4828 billion Yuan, the scheme of which plans to set up 8 heat sources, is the best choice for the case. (The detailed calculation results of this scheme are shown in Table 3.)

By observing the algorithms, the following is also concluded.(1)The optimal solution of IPSO is better than PSO. The optimum LCC which calculated by IPSO is not larger than PSO for all 29 schemes. The maximum -value is 2.4 million Yuan in the scheme which plans to set up 29 heat sources.(2)The real minimum LCC was not calculated by PSO. The minimum LCC calculated by PSO is 1.9 million Yuan larger than the minimum LCC calculated by IPSO.

Figure 7 compares the LCC convergence curves of two algorithms in three kinds of schemes, respectively. When the population size and the iteration number of PSO are same as those of IPSO during the HSP optimization process, although the PSO algorithm is faster for giving the optimization results, but the optimal results by IPSO are better than the searcher values by PSO. The main reason for current performance is that IPSO can avoid local optimal solution and then further expand the search space so as to find a better solution.

Hence, it can be concluded that the improvement approach is effective, and the proposed method IPSO has better significance in solving the HSP problem and competitive to PSO algorithm.

5. Discussion

Section 4.2 referred to the values of , , and which may influence the computational results. In Figures 8 and 9, the results obtained by IPSO were also proved. Algorithm calculation comparison with different parameters is shown in Figures 8 and 9, which is a visual display of the coordinates of heating sources with different parameters.

The parameters of each case in the figure are summarized in Table 4.

The influence aspect of the algorithm is worth further study, but because of the major goal of the present study, the more details were not presented here and will be discussed in a separate paper.

6. Conclusions and Prospects

(1)In this paper, we presented an integral mathematical model for solving the heating system planning (HSP) problem taking into account minimizing the cost of heating system for a given life cycle time.(2)According to the particularity of HSP problem, the particle swarm optimization (PSO) algorithm was introduced and improved, the new definition and update rule of velocity and position vector were proposed, and the improvement approach about generating a random velocity was adopted to avoid particle swarm into local optimal solution. Then an actual case study was calculated to check its feasibility in practical use. The results show that the IPSO algorithm can more preferably solve the HSP problem than PSO algorithm.(3)Although there is no more discussion about the influence of computational results by changing the values of algorithm parameters (, , and ), but the results of the case study still show the potential to provide useful information when making decisions in the practical planning process. Thus, it is believed that if this approach is applied correctly and in combination with other elements, such as the accurate prediction of heating load, the running efficiency of equipment, and the real operation situation, it can become a powerful and effective optimization tool for HSP problem.