Harmony Search and NatureInspired Algorithms for Engineering Optimization
View this Special IssueResearch Article  Open Access
Parameter Estimation for Traffic Noise Models Using a Harmony Search Algorithm
Abstract
A technique has been developed for predicting road traffic noise for environmental assessment, taking into account traffic volume as well as road surface conditions. The ASJ model (ASJ Prediction Model for Road Traffic Noise, 1999), which is based on the sound power level of the noise emitted by the interaction between the road surface and tires, employs regression models for two road surface types: densegraded asphalt (DGA) and permeable asphalt (PA). However, these models are not applicable to other types of road surfaces. Accordingly, this paper introduces a parameter estimation procedure for ASJbased noise prediction models, utilizing a harmony search (HS) algorithm. Traffic noise measurement data for four different vehicle types were used in the algorithm to determine the regression parameters for several road surface types. The parameters of the traffic noise prediction models were evaluated using another measurement set, and good agreement was observed between the predicted and measured sound power levels.
1. Introduction
It is important to evaluate the primary impact of traffic noise for several different road surfaces in order to predict the interactive noise between road surfaces and vehicle tires. A number of noise prediction models [1–6] have been developed for environmental estimation of traffic noise levels in terms of vehicle and pavement types. For example, the vehicle types in the ASJ model [1, 2] are large vehicles, medium vehicles, light trucks, and cars, and the pavement types are densegraded asphalt (DGA) and permeable asphalt (PA). Thus, the ASJ traffic noise prediction model is restricted to four vehicle types, varying vehicle speeds, and two pavement types.
In order to expand the applicability of the ASJ model, the present paper introduces a parameter estimation procedure based on a harmony search (HS) algorithm as follows: (a) the traffic noise for a targeted road surface is measured, (b) the parameters of the noise prediction model are estimated via the HS algorithm, and (c) the resulting coefficients are evaluated using another set of measurements (consisting of vehicle speeds and vehicle types) for a different traffic volume on the targeted road surface. To validate the proposed traffic noise prediction technique, traffic noise measurement sets from three different surface types were used in this study: stone mastic asphalt surfaces (SMA), 30 mm transversely tined Portland cement concrete surfaces (30 mm trans.), and 18 mm longitudinally tined Portland cement concrete surfaces (18 mm long.), as shown in Figure 1. The measurement site is a 7.7 km long. 2lane highway along the side of the south bound of Jungbu inland highway in South Korea. This measurement section includes both asphalt and Portland cement concrete pavements.
(a)
(b)
(c)
This paper is organized as follows. Section 2 describes the ASJ model and vehicle characterization. Section 3 explains the application of the HS algorithm to estimate the parameters of ASJbased noise prediction models, and Section 4 presents the conclusions of this research.
2. ASJ Model
The Acoustic Society of Japan (ASJ) published the ASJ model [1, 2] for calculating road traffic noise. The procedure involves calculating the noise level generated by traffic as well as the attenuation during noise propagation. The octaveband power spectrum for nominal midband frequencies (63 Hz to 8 kHz) can be generated according to the ISO 96132 standard [7]. The ASJ model classifies vehicles into four types listed in Table 1.

The Aweighted overall sound power levels () of the noise emitted by interactions between vehicles and pavement are listed in Table 2 for a densegraded asphalt (DGA) surface. In terms of nominal midband frequencies, the individual Aweighted sound power level for each octave band () is calculated as follows: where is the overall sound power level (dB) and is the relative level (dB) at the th nominal midband frequency given by is a factor of the standard Aweighting frequency given by the following equation: where is the nominal midband frequency. Finally, is the correction factor defined by
 
is a velocity. 
2.1. Equivalent Sound Power Level for a Road
The Aweighted sound power level emitted by a specific type of vehicle moving along a road over a specified time period, , can be calculated via the following equation: where is the basic sound power level (dB) emitted by the vehicle (listed in Table 2), is the length of the road segment (in meters), is the mean speed for the vehicle type (in km/h), and is the hourly traffic flow for the vehicle type (in vehicles/h). The equivalent sound power level, , emitted by all vehicles moving along the road can then be calculated as follows: where is the Aweighted sound power level for each vehicle type. Here, , , , and are the Aweighted sound power levels for a large vehicle, a medium vehicle, a light truck, and a car, respectively, as listed in Table 2. Furthermore, the attenuation (e.g., geometrical divergence, atmospheric absorption, ground effect, and screening structures) was calculated based on ISO 96132 [7].
3. Application of the Harmony Search Algorithm
3.1. Harmony Search Algorithm
This section describes the procedure for estimating the parameters of noise prediction models, using a harmony search (HS) algorithm that employs a heuristic algorithm based on an analogy with natural phenomena [8–12]. The detailed procedure for applying a harmony search consists of four steps as follows.(1) The algorithm parameters are specified. These include the harmony memory size (HMS), initialized as the number of solution vectors in the harmony memory (HM), the harmony memory consideration rate (HMCR, between 0 and 1), the pitch adjustment rate (PAR, between 0 and 1), and the maximum number of improvisations (or stopping criterion), which terminates the HS program. The optimization problem is specified as follows: where is the objective function, is the set of decision variables , is the number of decision variables, and is the possible range of values for the th decision variable; that is, , where and are the respective lower and upper bounds for the th decision variable. To estimate the parameters of a noise prediction model, the following minimization function can be used: where is the predicted equivalent sound power level of (5), which can be calculated as follows: where , , , and are the Aweighted sound power levels for a large vehicle, a medium vehicle, a light truck, and a car, respectively. is the measured equivalent sound power level obtained from previous research [13, 14]. In this optimization problem, the Aweighted sound power levels can be defined as given in Table 3. Thus, the coefficients of , , , and must be determined via the HS algorithm. The slope is fixed in the ASJ models for both surface types (DGA and PA); therefore, the slope given in Table 3 is fixed at 30.(2)The HM matrix is initially filled with randomly generated solution vectors up to the HMS, together with the corresponding function values of the random vectors : (3) A column vector of the newly generated harmony, , is improvised utilizing the following three mechanisms: (a) random selection, (b) memory consideration, and (c) pitch adjustment. In the random selection, the value of each decision variable, , in the column vector is randomly chosen within the range of values with a probability of (1HMCR). HMCR (which is between 0 and 1) is the rate at which a single value is chosen from the historical values stored in the HM. The value of each decision variable selected by memory consideration is examined in terms of pitch adjustment. This operation uses the PAR parameter (which is the rate of the necessary pitch adjustment according to the neighboring pitches) with a probability of HMCR × PAR. Pitch adjustment is applied to each variable as follows: If the newly generated column vector is better than the worst harmony in the HM, based on evaluation of the objective function, the newly generated column vector is included in the HM, and the existing worst harmony is excluded from the HM.(4)If the stopping criterion (or maximum number of improvisations) is satisfied, the computation is terminated. Otherwise, Steps 3 and 4 are repeated.
 
is a velocity. 
3.2. Application to Parameter Estimation for Noise Prediction Models
In order to estimate the parameters of noise prediction models based on the objective function of (8), noise measurements for three different road surfaces were obtained from previous research, which was conducted on a test track [13, 14]. The vehicle velocities and hourly traffic flows are listed in Table 4.

To apply the HS algorithm to parameter estimation for the noise prediction models, the four coefficients were determined for each road surface type (stone mastic asphalt (SMA) surface, 30 mm transversely tined Portland cement concrete surface (30 mm trans.), and 18mm longitudinally tined Portland cement concrete surface (18 mm long.)), as shown in Figure 2 and Table 5, based on the training data from Table 4. In this way, the coefficients of the noise prediction models, which are dependent on the road surface type, can be updated via the HS algorithm. As a result, noise prediction models can be provided for various surface types by using the HS algorithm to update the ASJ model equations.
Another set of testing data (given in Table 6) was used to evaluate whether or not the updated noise prediction models (with the four coefficients estimated by the HS algorithm) provided results consistent with measured noise levels. The predictions and measurements are compared in Table 7; good agreement was noted for all three surface types. Using the original ASJ model, the prediction value is resulted in a same noise level regardless of surface types and in bad agreement when compared with the measured noise levels.


Finally, the Aweighted sound power levels for the individual octave bands were estimated for the SMA, 30 mm trans., and 18 mm long. surface types, utilizing the parameters estimated via the HS algorithm and (1). For example, Figure 3 shows the Aweighted sound power levels for the octave bands for the three different surface types in the case of the large vehicle speeding at 80 km/h.
4. Conclusions
In this study, it was shown that the optimization problem related to updating the noise prediction models for several surface types could be solved using an HS algorithm process. The process involves (a) obtaining measurements for different road surfaces, (b) estimating the coefficients of the noise prediction models using this measurement set as training data, and (c) evaluating the estimated coefficients using another measurement set as testing data. When this procedure was utilized, an evaluation of the parameters of the traffic noise prediction model yielded good agreement between predicted and measured sound power levels.
Conflict of Interests
The authors of the paper do not have a direct financial relationship with the commercial identities mentioned in this paper.
Acknowledgment
This research was supported by a grant from a Strategic Research Project “Development of Traffic Noise Modeling and Noise Reduction Technology as Road Traffic Noise by Location” funded by the Korea Institute of Construction Technology.
References
 Acoustic Society of Japan, “ASJ prediction model for road traffic noise,” Bull. Acoustic Society of Japan, vol. 55, no. 1, pp. 281–321, 1999. View at: Google Scholar
 K. Yamamoto, “Road traffic noise prediction model “ASJ RTNModel 2008”: report of the research committee on road traffic noise,” Acoustical Science and Technology, vol. 31, no. 1, pp. 2–55, 2010. View at: Publisher Site  Google Scholar
 U. Sandberg, “Road traffic noise—the influence of the road surface and its characterization,” Applied Acoustics, vol. 21, no. 2, pp. 97–118, 1987. View at: Google Scholar
 C. Steele, “Critical review of some traffic noise prediction models,” Applied Acoustics, vol. 62, no. 3, pp. 271–287, 2001. View at: Publisher Site  Google Scholar
 D. S. Cho, J. H. Kim, T. M. Choi, B. H. Kim, and D. Manvell, “Highway traffic noise prediction using method fully compliant with ISO 9613: comparison with measurements,” Applied Acoustics, vol. 65, no. 9, pp. 883–892, 2004. View at: Publisher Site  Google Scholar
 T. Bennert, D. Hanson, A. Maher, and N. Vitillo, “Influence of pavement surface type on tire/pavement generated noise,” Journal of Testing and Evaluation, vol. 33, no. 2, pp. 94–100, 2005. View at: Google Scholar
 International Organization for Standardization, ISO 96132:1996. Acoustics—attenuation of sound during propagation outdoors—part 2: general method of calculation, 1996.
 S. Mun and Z. W. Geem, “Determination of viscoelastic and damage properties of hot mix asphalt concrete using a harmony search algorithm,” Mechanics of Materials, vol. 41, no. 3, pp. 339–353, 2009. View at: Publisher Site  Google Scholar
 S. Mun and Z. W. Geem, “Determination of individual sound power levels of noise sources using a harmony search algorithm,” International Journal of Industrial Ergonomics, vol. 39, no. 2, pp. 366–370, 2009. View at: Publisher Site  Google Scholar
 Z. W. Geem, “Particleswarm harmony search for water network design,” Engineering Optimization, vol. 41, no. 4, pp. 297–311, 2009. View at: Publisher Site  Google Scholar
 J. Fourie, R. Green, and Z. W. Geem, “Generalised adaptive harmony search: a comparative analysis of modern harmony search,” Journal of Applied Mathematics, vol. 2013, Article ID 380985, 13 pages, 2013. View at: Google Scholar  MathSciNet
 S. Lee, S. Mun, and H. Moon, “Determination of pavement rehabilitation activities through a permutation algorithm,” Journal of Applied Mathematics, vol. 2013, Article ID 252808, 5 pages, 2013. View at: Publisher Site  Google Scholar
 S. Mun, D.S. Cho, and T.M. Choi, “Influence of pavement surface noise: the Korea Highway Corporation test road,” Canadian Journal of Civil Engineering, vol. 34, no. 7, pp. 809–816, 2007. View at: Publisher Site  Google Scholar
 D.S. Cho and S. Mun, “Determination of the sound power levels emitted by various vehicles using a novel testing method,” Applied Acoustics, vol. 69, no. 3, pp. 185–195, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 DeokSoon An 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.