Journal of Chemistry

Volume 2015 (2015), Article ID 313105, 10 pages

http://dx.doi.org/10.1155/2015/313105

## Parameter Identification of the 2-Chlorophenol Oxidation Model Using Improved Differential Search Algorithm

^{1}Department of Environmental Engineering, Anhui Jianzhu University, Hefei 230022, China^{2}Key Laboratory of Water Pollution Control and Wastewater Resource of Anhui Province, Hefei 230601, China^{3}School of Resources and Environmental Engineering, Hefei University of Technology, Hefei 230009, China

Received 10 December 2014; Accepted 23 January 2015

Academic Editor: Jian Lu

Copyright © 2015 Guang-zhou Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Parameter identification plays a crucial role for simulating and using model. This paper firstly carried out the sensitivity analysis of the 2-chlorophenol oxidation model in supercritical water using the Monte Carlo method. Then, to address the nonlinearity of the model, two improved differential search (DS) algorithms were proposed to carry out the parameter identification of the model. One strategy is to adopt the Latin hypercube sampling method to replace the uniform distribution of initial population; the other is to combine DS with simplex method. The results of sensitivity analysis reveal the sensitivity and the degree of difficulty identified for every model parameter. Furthermore, the posteriori probability distribution of parameters and the collaborative relationship between any two parameters can be obtained. To verify the effectiveness of the improved algorithms, the optimization performance of improved DS in kinetic parameter estimation is studied and compared with that of the basic DS algorithm, differential evolution, artificial bee colony optimization, and quantum-behaved particle swarm optimization. And the experimental results demonstrate that the DS with the Latin hypercube sampling method does not present better performance, while the hybrid methods have the advantages of strong global search ability and local search ability and are more effective than the other algorithms.

#### 1. Introduction

In environmental technology, supercritical water oxidation (SCWO) is an innovative technology. Supercritical water exhibits some properties very different from those of water at standard conditions. One of its main applications is the oxidative treatment of aqueous wastes containing organic compounds in the so-called “supercritical water oxidation” (SCWO) process [1]. Its potential advantage lies in the fact that it can rapidly and efficiently decontaminate diverse organic wastes [2]. In wastewater treatment, supercritical water oxidation technology is considered to be an effective clean technology and is applied to the organic wastewater treatment [3–5].

The key technical issue in the design and application of SCWO process systems is to understand the principle of kinetic reaction and identify the key reaction mechanisms and parameters in the kinetic model [2]. In modeling, the mathematical model is designed to simulate the system behavior for a given set of model parameters representing the actual response properties of the system. The more accurately the model matches the experiment observations, the more representative the model parameters are assumed to be. Therefore, the above parameter identification or parameter estimation is essentially an inverse problem [6].

In many applications, inverse problems can be nonlinear and ill posed which makes them difficult to solve numerically. The basic idea is to carry out the optimization by constructing the objective function. The best set of parameters is selected from within reasonable ranges by adjusting the values until the discrepancies between the model generated values and observations are reduced to a minimum in the weighted least squares sense [7]. In addition, in the process of model identification, along with calibration to obtain a set of optimal parameter sets giving the best performance of a simulation model, there are other important issues such as uncertainty and identify ability of the parameters which should be taken into account for a proper model evaluation.

Generally speaking, the optimization methods for solving the inverse problems are classified into two groups. One sort is traditional methods: (1) gradient-based optimal methods such as Levenberg-Marquardt algorithm and the conjugate gradient descent method; (2) no gradient-based algorithms such as simplex. And the other is search-based heuristic optimization methods such as genetic algorithm (GA), simulated annealing (SA), and tabu search algorithm (TA), which are used in dealing with complex nonlinear inverse problems. They have been proved to have the ability to identify global or near-optimal solutions in the global scope by a large number of engineering optimization cases.

In recent years, many swarm intelligence optimization algorithms, such as ant colony, particle swarm optimization, artificial fish swarm, and artificial bee colony, have been put forward unceasingly. Many swarm intelligence algorithms have been applied in the parameter identification from diverse engineering domains. For the chemical kinetic models, chaos genetic algorithm, clonal selection algorithm, hierarchical differential evolution algorithm, adaptive differential evolution algorithm, and ABC-simplex and DNA genetic algorithm have been applied to estimate the kinetic parameters [8–13].

Though the above intelligence optimization algorithms can usually obtain better optimization results, a more effective and convenient optimization algorithm is still worth seeking. Differential search (DS) algorithm, proposed by Civicioglu in 2012 [14], is a new optimization algorithm with a short code, short run time, and simple principle. Subsequently, some applications show that it is an effective optimization algorithm and has better optimization performance [14, 15].

Aimed at the case where the sensitive and uncertainty analysis is easy to be ignored in parameter identification process, the main objectives of this paper are to use Monte Carlo method to carry out sensitivity analysis based on Monte Carlo analysis toolbox (MCAT); to develop two improved algorithms to enhance optimization performance; to analyze their optimization performance by a case and compare it with that of some other swarm intelligence optimization algorithms.

This paper is organized as follows. In Section 2, we introduce the sensitivity analysis and MCAT. Then, the improved DS algorithms are presented in Section 3. In Section 4, we apply MCAT to make sensitivity analysis and improved DS algorithms to test a case of parameter identification from a chemical kinetics model and analyze its performance. Finally, we give some conclusions in Section 5.

#### 2. Sensitive Analysis and Introduction of MCAT

##### 2.1. Sensitive Analysis

To provide guidance for a better understanding of the modeling and reduce the response uncertainties most effectively, sensitivity analysis (SA) of the input importance on the response variability can be useful [16]. Sensitivity analysis is valuable tool for identifying important model parameters, testing the model conceptualization, and improving the model structure. They help to apply the model efficiently and to enable a focused planning of future research and field measurement [17].

##### 2.2. Monte Carlo Analysis Toolbox

The Monte Carlo analysis toolbox (MCAT), developed by Wagener et al., is a collection of analysis and visualization functions integrated through a graphical user interface based on MATLAB language [18]. The toolbox can be used to analyze the results from Monte Carlo (MC) parameter sampling experiments or from model optimization methods that are based on population evolution techniques. A number of powerful functions are included in the toolbox to investigate the structure, sensitivity, and parameter and output uncertainty of models. In view of the problem that system identification techniques are often limited to linear systems and are subject to many statistical assumptions, MCAT uses the Monte Carlo simulation approach to solve the “equifinality” [19] instead of classical system identification techniques [20]. Monte Carlo parameter sampling provides an alternative measure to classical statistical modeling and uncertainty estimation methods. It is designed to help interpret the results of such modeling experiments.

The core of the toolbox is based on the concept of Regional Sensitivity Analysis [21–23], and its extension to the generalized likelihood uncertainty estimation (GLUE) technique was developed at Lancaster University by Beven and coworkers [19, 23]. For detailed description, see guide of MCAT.

#### 3. Improved Differential Search Algorithm

##### 3.1. Main Idea of the Basic DS Algorithm

DS is an algorithm developed for solution of numeric optimization problems [14]. It was inspired by migration of living beings which constitute superorganisms utilizing the concept of Brownian-like motion. Suffering from the impact of periodical change of regional resources, many species of the living beings show seasonal migration behavior during the year. In the migration movement, the migrating species of living beings constitute a superorganism containing large number of individuals.

Migration behavior allows them to move from one habitat to more efficient habitat. Of course, they change their position by moving toward more fruitful areas. The movement of superorganisms can be described by a Brownian-like rand-walk model [24]. In algorithm, evolving process of populations is corresponding to artificial superorganisms migration to global optimum solution of the problem. During the migration, the superorganisms judge whether some randomly selected positions are suitable for temporary stopovers. If the position is suitable to stop over for a temporary time during the migration, the members of superorganisms that made the discovery settle at this position; subsequently, they continue their migration from this position.

##### 3.2. Main Steps of the Basic DS

Step 1. Set up the size of population, dimension of the problem, and numbers of iterations. Step 2. Evaluate the initial population. Step 3. Execution cycle is as follows:(1)use the random shuffling method to produce the donor population;(2)compute the scale factor;(3)obtain the stopover site position population;(4)use the random process to decide the individuals to participate in the search process of stopover site;(5)check whether the elements of stopover site are beyond the limits of the search space; if beyond, a position in the search space will be randomly produced;(6)evaluate the individuals of the stopover site position population; if it is better, replace the original one; or it remains unchanged;(7)update the entire group. Step 4. If the numbers of iterations are satisfied, then stop; otherwise, continue step 3.

##### 3.3. Improved DS Algorithms

###### 3.3.1. Construction of the Initial Population Using Latin Hypercube Sampling Method

Latin hypercube sampling (LHS) is a statistical method for generating a sample of plausible collections of parameter values from a multidimensional distribution. The sampling method is often used to construct the design of experiments. Latin hypercube sampling can be viewed as a compromise procedure that incorporates many of the desirable features of random sampling and stratified sampling and also produces more stable analysis outcomes than random sampling [25]. This sampling approach ensures that each of the input variables has all portions of its range represented; LHS is computationally cheap to generate and can cope with many input variables [26]. Consequently, LHS was used to produce the initial population instead of uniform distribution method so as to enhance the performance of algorithm. Its realization was simply described [26]: generate a sample size from the variables . The range of each variable is divided into nonoverlapping intervals on the basis of equal probability size . One value from each interval is selected at random with respect to the probability density in the interval. The values thus obtained for are paired in a random manner with the values of . These pairs are combined randomly with the values of to form triplets and so on, until a set of -tuples is formed. This set of -tuples is the Latin hypercube sample. In this paper, we use “lhsdesign” function of MATLAB to realize it.

###### 3.3.2. Combination of DS and Simplex

Simplex was proposed by Spendley et al. in 1962. Subsequently, it was improved by Nelder and Mead in 1965 [27]. It is a direct search optimization method, and has the fast searching ability to the local optimization without the need of gradient information. Its biggest weakness is sensitive to the choice of the initial points. For the DS algorithm, though it has the powerful global search ability, the convergence speed of the basic DS algorithm is slow and needs a large number of iterations.

Therefore, to own both the merits of the better local searching ability from the simplex method and the global searching ability from the DS algorithm, the hybrid algorithms of DS and simplex are developed.

The paper designed two methods: (1) firstly completely run the DS algorithm to obtain the optimization solution; then simplex method is employed to process local search based on the above solution which is considered as the initial solution; (2) at the end of each run of the DS algorithm, simplex method is employed to execute local search based on the optimal solution of the current generation. If the better solution can be obtained, the optimal solution of the current generation will be replaced. Therefore, it effectively absorbs the advantages of DS and simplex method. It not only can overcome the deficiency that DS has weak local search ability but also can effectively solve the sensitivity problem of simplex to the initial search points.

#### 4. Parameter Estimation for Kinetic Model of the 2-Chlorophenol

##### 4.1. The Kinetic Model

Supercritical water oxidation (SCWO) is an effective treatment technology for organic waste in environmental field. The removing rate expression for 2-chlorophenol (2-CP) is as follows [28]: where , , and represent 2-CP, , and concentration, respectively. is the removing rate of 2-CP, is the preexponential factor, denotes the activation energy, is the molar gas constant, and , , and are the reaction orders of 2-CP, , and , respectively. After the analysis and treatment, the following equations can be obtained [28]: where represents the conversion rate of 2-CP. The objective is to identify the five parameters (). This optimization problem is a complex one with high nonlinearity.

##### 4.2. The Objective Function

This parameter estimation problem is based on minimizing the objective function given by the sum of squared deviation. The objective function is defined as where is the sample number and is the real conversion rate of 2-CP in the th sample data while is the conversion rate calculated by the kinetic model.

##### 4.3. Sensitivity Analysis

We used MCAT to carry out the sensitivity analysis and uncertainty analysis. After the 50000 runs, the threshold of sum of squared errors (SSE) for the objective function was taken to 2, and then the following plots were obtained.

###### 4.3.1. Dotty Plots of Objective Function

The dotty plots shown in Figure 1 is the result of rejection threshold of sum of squared errors (SSE) set to 2. The “surfaces” of the dotty plot for , , and parameters have a clearly defined minimum, so they can be considered to be well identified. Followed by , however, the parameter is difficult to be defined.