Sustainable Optimization for Wastewater Treatment System Using PSF-HS

Abstract: The sustainability in a river with respect to water quality is critical because it is highly related with environmental pollution, economic expenditure, and public health. This study proposes a sustainability problem of wastewater treatment system for river ecosystem conservation which helps the healthy survival of the aquatic biota and human beings. This study optimizes the design of a wastewater treatment system using the parameter-setting-free harmony search algorithm, which does not require the existing tedious value-setting process for algorithm parameters. The real-scale system has three different options of wastewater treatment, such as filtration, nitrification, and diverted irrigation (fertilization), as well as two existing treatment processes (settling and biological oxidation). The objective of this system design is to minimize life cycle costs, including initial construction costs of those treatment options, while satisfying minimal dissolved oxygen requirements in the river, maximal nitrate-nitrogen concentration in groundwater, and a minimal nitrogen requirement for crop farming. Results show that the proposed technique could successfully find solutions without requiring a tedious setting process.


Introduction
Sustainability is the capacity to maintain any system in terms of environment, economics, and society.In this regard, the sustainability of wastewater treatment systems is very important because it is strongly related with environmental pollution, economic expenditure, and public health.
This study is mainly interested in a wastewater treatment system in which effluent contamination dumped to a river stream should satisfy river quality standards; life cycle cost, instead of initial construction cost, should be minimized under the available budget, and groundwater contamination should satisfy drinking-water quality standards.
This optimal design of a wastewater treatment system can be complicated and venturous.The treatment system for the wastewater discharged to the river should consider various factors, such as waste sources and properties, treatment processes, transport of the wastes in the river, and waste effects on biota and human beings.While each component in the system can be analyzed separately, the problem can be also considered in its totality for managing it optimally [1].
In order to make the river biologically healthy, a certain level of dissolved oxygen (DO) is required and wastewater treatment, which is a procedure to make wastewater into an effluent that is turned back to the water cycle with less environmental issues [2], should be implemented.So far, there has been research on optimal DO control issues for various rivers, such as the Willamette River in the USA [3], the Nitra River in Slovakia [4], and the Yamuna River in India [5].Additionally, more broadly modeling and optimization techniques have been applied to wastewater treatment systems [6,7].However, Haith [1] proposed one of the most challenging and practical problems for the wastewater treatment optimization.
The real-scale problem by Haith was previously tackled by the harmony search (HS) algorithm [8].Results showed that while a calculus-based technique, generalized reduced gradient (GRG), found worse solutions or even diverged, HS robustly found better solutions.When compared with another meta-heuristic technique, genetic algorithm (GA), HS also found better solutions than those from GA in terms of best and average solutions.Nonetheless, the ranking of performance among metaheuristic algorithms is not fixed because of the no free lunch theorem [9][10][11].
Actually, the solution quality of the calculus-based GRG technique is very sensitive to the initial solution.Only with good starting solutions can good final solutions be obtained.On the contrary, soft-computing approaches, such as HS and GA, are free from the task of assigning a good initial solution because they handle multiple initial solutions.
However, those meta-heuristic algorithms also have their own inherent drawback: they have to carefully choose good values for algorithm parameters, such as crossover rate, mutation rate, harmony memory considering rate (HMCR), and pitch adjusting rate (PAR), in order to guarantee good final solutions [12].Thus, researchers have, so far, tried to find ways to avoid the tedious and boring task of setting proper algorithm parameter values by introducing dynamic or adaptive techniques (this technique was named as parameter-setting-free (PSF) in HS) for various engineering applications, such as chlorine injection [13], sewer network design [14], structural design [15,16], energy consumption of air-handling units [17], water network design [18], hydrologic parameter calibration [19], and economic dispatch [20].
The goal of this study is, for the first time, to introduce the parameter-setting-free harmony search (PSF-HS) technique to the sustainable wastewater treatment optimization for river conservation in order to eliminate the tedious parameter setting process, which helps any researcher or practical engineer to more easily optimize the design of the wastewater treatment system.

Problem Description
A city with the population of 100,000 people disposes of wastewater into a river.The wastewater of 40,000 m 3 /day is first treated with two operations of settling and biological oxidation to eliminate organic material.However, these treatments are not enough to satisfy the water quality standard of 5 ppm (or 5 mg/L) dissolved oxygen (DO) during summertime, which was set by regulation agencies for maintaining aquatic life.Thus, the city is forced to construct additional treatment systems to enhance the effluent quality.
The expanded wastewater treatment system consists of two additional treatment systems (filtration and nitrification) and one diverted irrigation system as shown in Figure 1.In the system, we have to decide five decision variables, including (1) wastewater amount discharged directly into the river (Q 1 10 3 m 3 /day); (2) wastewater amount treated by filtration operations (Q 2 10 3 m 3 /day); (3) wastewater amount treated by nitrification operations (Q 3 10 3 m 3 /day); (4) wastewater amount treated by both nitrification and filtration operations (Q 4 10 3 m 3 /day); and (5) wastewater amount diverted for irrigation and fertilization (Q 5 10 3 m 3 /day).Table 1 shows the effluent quality after each treatment.
Secondary+Nitrification+Filtration ( 4 Additionally, there exists one more decision variable r which stands for irrigation rate (cm/week).For the diverted wastewater amount 5 Q , the following mass balance equation can be considered: where T is irrigation period (13 weeks in this study), P is the precipitation amount (cm) during the irrigation period, and ET is the evapotranspiration amount (cm) during the irrigation period ( ET is assumed to be equal to P in this study).If the nitrogen concentration of the diverted amount 5 Q is n (20 ppm in this study), the total nitrogen amount becomes rTn 1 .0 (kg/ha).If crop nitrogen uptake amount is NC (170 kg/ha in this study), the surplus nitrogen amount of 5 , and the nitrate-nitrogen concentration n c (ppm) in the percolation amount becomes as follows: Here, the city-regulated n c is less than or equal to 10 ppm for drinking water usage.
Additionally, the nitrogen amount in 5 Q should be enough to cover the crop's nitrogen requirement as follows:  Additionally, there exists one more decision variable r which stands for irrigation rate (cm/week).For the diverted wastewater amount Q 5 , the following mass balance equation can be considered: where T is irrigation period (13 weeks in this study), P is the precipitation amount (cm) during the irrigation period, and ET is the evapotranspiration amount (cm) during the irrigation period (ET is assumed to be equal to P in this study).If the nitrogen concentration of the diverted amount Q 5 is n (20 ppm in this study), the total nitrogen amount becomes 0.1rTn (kg/ha).If crop nitrogen uptake amount is NC (170 kg/ha in this study), the surplus nitrogen amount of Q 5 becomes 0.1rTn ´NC, and the nitrate-nitrogen concentration c n (ppm) in the percolation amount becomes as follows: Here, the city-regulated c n is less than or equal to 10 ppm for drinking water usage.Additionally, the nitrogen amount in Q 5 should be enough to cover the crop's nitrogen requirement as follows: For the discharged amount (Q 1 `Q2 `Q3 ) into the river, dissolved oxygen Cpxq at a distance (x km) from the point of wastewater discharge can be expressed as the following differential equation: where u is river flow velocity (7.9 km/day in this study); k 2 is the re-aeration rate (0.5/day in this study); C s is saturation DO (8.0 ppm in this study); B and N are the remaining carbonaceous biochemical oxygen demand (CBOD) and nitrogenous biochemical oxygen demand (NBOD) (ppm) at distance x; and k 1 and k n are rate constants (0.35/day and 0.2/day, respectively).The first term on the right-hand side of Equation ( 4) denotes an oxygen increase due to re-aeration and the second and third terms denote an oxygen decrease due to oxidation of carbonaceous and nitrogenous material, respectively.From the following relationship: Equation ( 4) becomes: where B 0 and N 0 are, respectively, river CBOD and NBOD right after discharge.Additionally, Equation ( 7) can be rewritten in the format of a first-order non-homogeneous ordinary differential equation as follows: We can set an integral factor µ as follows: µ " e Then, if we multiply µ and both sides of Equation ( 8), it becomes as follows because Or, it becomes: If we integrate Equation ( 11), it becomes: If we divide both sides of Equation ( 12) by µ, it becomes: Sustainability 2016, 8, 321 Here, e k 2 u x qpxqdx becomes: Thus, Equation ( 13) becomes: If x " 0 and D " r r C `r C, Equation (15) becomes: where C 0 is the river DO right after discharge.Finally, we can obtain the analytic solution as follows.
Actually, the name of the equation is Streeter-Phelps [21], which is also known as the DO sag equation, and first derived in 1925 based on data from the Ohio River.The model conveniently shows how DO decreases along the river reach by degradation of the biological oxygen demand (BOD).
If river water and wastewater are completely mixed at the discharge point, the initial C 0 , B 0 , and N 0 can be calculated using a weighted average.Since the river flow is 110,000 m 3 /day, river DO is 8.0 ppm, river CBOD is 2.0 ppm, and river NBOD is 5.0 ppm, C 0 , B 0 , and N 0 become as follows: The above-mentioned problem originally came from Haith [1].However, because it does not reflect full aspects of the optimization procedure, this study modifies the problem by adjusting cost coefficients to obtain more reasonable optimization results.

Optimization Formulation
The objective function of this wastewater treatment problem is the total design cost of the system, which should be minimized while satisfying all the design constraints, as follows [1,8]: The objective function consists of three costs, such as filtration cost C f t , nitrification cost C nt , and irrigation cost C ir .While many design optimization problems just consider initial capital (or construction) cost, this study considers life cycle costs, including operation and maintenance costs, crop sales benefit, and land rental fee.
The filtration cost C f t ($10 3 /year), which consists of capital costs (first term in right hand side) and operation and maintenance (O and M) costs (second term), is expressed as follow: The nitrification cost C nt ($10 3 /year), which consists of capital costs and O and M costs, is expressed as follow: C nt " 6.9Q 3 0.68 `5.3Q 3 0.42 The irrigation cost C ir ($10 3 /year), which consists of capital cost of transmission pipe, capital cost of the storage system, O and M costs of the storage system, capital cost of irrigation system, O and M costs of the irrigation system, and net benefit of crop farming (crop sales minus land rent), is expressed as follow: C ir " 21.9Q where A is the irrigated area (ha), which can be calculated as 70Q 5 {r.Constraints for this problem are: 8p1 ´e´0.063xq `C0 e ´0.063x ´2.33B 0 pe ´0.044x ´e´0.063xq ´0.67N 0 pe ´0.025x ´e´0.063xq ě 5 (28) 6.54 ď r ď 13.07 The constraint in Equation (25) represents that total wastewater amount is equal to 40,000 m 3 /day; the constraint in Equation (26) represents that the amount for nitrification should be greater than, or equal to, that for both nitrification and filtration; the constraint in Equation (27) represents that any individual amount cannot exceed the total amount as a non-negative value; the constraint in Equation (28) represents that DO concentration at any point (x km downstream from the discharged point) should be greater than or equal to 5 ppm as regulated.In reality, we cannot test all the points because the number of them is astronomical.Thus, we instead test eight points with an even interval (5 km or 10 km) as previous studies did, which contain the DO sag curve; and the constraint in Equation (29) represents that the nitrate-nitrogen concentration c n should be less than, or equal to, 10 ppm as in Equation ( 2) and the total nitrogen amount should be greater than, or equal to, the crop nitrogen requirement as in Equation (3).

Parameter-Setting-Free Harmony Search
The basic harmony search algorithm starts with a solution pool named harmony memory (HM), where randomly generated and feasible solutions are initially filled to the size of the harmony memory (HMS; 30 in this study) as follows [8]: Next, a new harmony (solution vector) x New is improvised based on HM as follows: x where ∆ is the pitch adjusting amount, which is calculated as px U pper i ´xLower i q{1000 ˆup´1, 1q in this study (up´1, 1q is uniform random number between ´1 and 1), HMCR denotes a harmony memory considering rate (0 ď HMCR ď 1) and PAR denotes pitch adjusting rate (0 ď PAR ď 1).In memory considering operation, any value in HM can be selected using uniform random index intpHMS ˆup0, 1qq `1 (intp¨q is integer function which returns the integer portion of a number).
If the newly generated harmony x New satisfies all the constraints and is better than the worst harmony x Worst stored in HM, the former replaces the latter as follows: x New P HM ^xWorst R HM (32) The basic HS algorithm repeats Equations ( 31) and (32) until a termination criterion is satisfied.In this study, the PSF-HS algorithm has one more matrix, named operation type matrix (OTM), as follows: OTM retains the operation (memory consideration, pitch adjustment, or random selection) by which each variable value is obtained.At the early stage of computation (for example, number of iterations is less than 4 ˆHMS), HMCR and PAR are fixed into 0.5 and OTM are being updated.Then, after the early stage, based on OTM, HMCR, and PAR are iteratively calculated, instead of using fixed values, as follows: where ctp¨q denotes the function which counts the number of designated operation.Since each decision variable has its own HMCR and PAR, Equation (31) should be slightly modified as follows:

Computation Results
In previous research [8], the basic HS algorithm, which was applied to the optimal design of the wastewater treatment system for river conservation, found good solutions while GRG got stuck in local optima, or even diverged.
To make matters worse, the original approach had a tendency to choose only Q 1 (no treatment) and Q 5 (irrigation) without choosing Q 2 (filtration) or Q 3 (nitrification).In other words, although we elaborate the filtration and nitrification treatment formulation, it is useless if it is not chosen in optimization process.Thus, a couple of cost functions (Equations ( 22) and ( 23)) were modified in order for every decision variable to have certain value.
When the PSF-HS approach (iterations = 5000 and 10 runs) in this study was applied to the modified real-scale wastewater treatment problem, it was able to find good solutions with the best solution of 175.18 ($10 3 /year) and the average solution of 207.87 ($10 3 /year).PSF-HS successfully produced the solutions without demanding the tedious algorithm parameter setting process.Table 2 shows the computational results of 10 runs from PSF-HS.
As seen in Figure 3, all HMCR's have converged into higher values (more than 0.95) after 2000 iterations.Meanwhile, as seen in Figure 4, all PARs, except that from irrigation rate r , came down to lower values between 0.0 and 0.05 after 2000 iterations.The PAR for irrigation rate r soared to around 0.5 after fluctuating between 0.05 and 0.4.

Conclusions
For river sustainability, a wastewater treatment system, which has three different options of treatments, such as filtration, nitrification, and diverted irrigation, was optimally designed using the PSF-HS technique which does not require the tedious task of algorithm parameter value setting.
When the PSF-HS technique was applied to the system, it could find the lowest cost design while As seen in Figure 3, all HMCR's have converged into higher values (more than 0.95) after 2000 iterations.Meanwhile, as seen in Figure 4, all PARs, except that from irrigation rate r, came down to lower values between 0.0 and 0.05 after 2000 iterations.The PAR for irrigation rate to around 0.5 after fluctuating between 0.05 and 0.4.

Conclusions
For river sustainability, a wastewater treatment system, which has three different options of treatments, such as filtration, nitrification, and diverted irrigation, was optimally designed using the PSF-HS technique which does not require the tedious task of algorithm parameter value setting.
When the PSF-HS technique was applied to the system, it could find the lowest cost design while satisfying the minimal DO requirement in the river, maximal nitrate-nitrogen concentration in groundwater, and the minimal nitrogen requirement for crop farming by automatically adapting its parameter values.This user-friendliness can help any designer/engineer to more easily use this model in real-world design process.
For future study, more realistic and up-to-date wastewater treatment systems can be considered and more robust PSF-HS techniques are expected to be developed.Additionally, because the optimization model in this study searches for the minimal budget solution in terms of life cycle costs (construction, operation, and management costs), while satisfying environmental criteria, the general format of this model can possibly be applied to other sustainability problems, such as a waste-disposal plant design.n nitrogen concentration of the diverted amount Q 5 (20 ppm

Figures 3 and 4
Figures 3 and 4 show the trends of HMCR and PAR, which have been automatically adapted, as the iteration goes.As seen in Figure 3, all HMCR's have converged into higher values (more than 0.95) after 2000 iterations.Meanwhile, as seen in Figure 4, all PARs, except that from irrigation rate r , came down to lower values between 0.0 and 0.05 after 2000 iterations.The PAR for irrigation rate r soared to around 0.5 after fluctuating between 0.05 and 0.4.

Figures 3 and 4
Figures3 and 4show the trends of HMCR and PAR, which have been automatically adapted, as the iteration goes.

Table 1 .
Effluent quality after each treatment.

Table 2 .
Optimization results by PSF-HS.