1. Introduction
Wind turbines are one of the most important sources of renewable energy [
1]. They are widely promoted and used around the world to efficiently provide environmentally friendly energy [
2]. One of the main disadvantages of wind turbines is the generation of infrasound and low-frequency noise (ILFN), which can cause annoyance and nuisance to nearby residents as well as people at a considerable distance from the turbine [
3,
4,
5]. Therefore, the propagation of ILFN should be modelled at the stage of predicting the acoustic impact of wind turbines on the environment. Infrasound is defined as noise in the range of 1 Hz to 20 Hz [
6,
7], while low-frequency noise has no standardised definition. Different frequency range limits for this noise can be found in the literature. In the DIN 45680 standard [
8], the low-frequency noise band is specified from 8 Hz to 100 Hz. According to Leventhall [
9], low-frequency noise ranges from 10 Hz to 200 Hz, while according to ACGIH limits [
10], infrasound and low-frequency noise is noise in the range of 1 Hz–80 Hz and according to other researchers up to 100 or 250 Hz or even 500 Hz [
11,
12,
13,
14]. In Poland, the Leventhall approach has been applied to low-frequency noise. In contrast, this study considers ILFN in octave bands with centre frequencies from 4 Hz to 250 Hz.
Many methods for modelling wind turbine noise can be found in the literature. The most common are complex methods that use different models. These include algorithms based on the “ray tracing” method [
15,
16,
17,
18], where the ray paths from each turbine source to each receiver are calculated based on favourable weather conditions or atmospheric inversion conditions. Current research is aimed at validating the models [
16] and extending the models to include multiple reflections from the ground and the influence of atmospheric turbulence [
15].
ILFN can be modelled using a sequential approach that involves a Computational Fluid Dynamics (CFD) model, an aeroelastic HAWC model, and an acoustic model. This approach enables the analysis of how different factors and turbine design parameters affect the noise generated by turbines [
19,
20].
A theoretical method using Parabolic Equations (PE) [
21,
22,
23,
24,
25,
26], as well as the Fast Field Program (FFP) method [
24], is also used to model ILFN generated by wind turbines.
Keith et al. [
24] conducted a comparative analysis of sound pressure level (SPL) modelling results for infrasound and low-frequency sound. The results of calculations using the PE and FFP methods were compared with the results of calculations using the ISO 9613-2 method for the extended frequency range (below 63 Hz). These results were also compared with the results of long-term measurements carried out in Canada as part of the “Health Canada’s Community Noise and Health Study” project.
Keith et al. [
24], based on their study, clearly stated that sound speed profiles obtained from actual meteorological data can represent conditions not included in the Harmonoise weather classes. In such cases, SPL calculations using the Harmonoise class and actual meteorological conditions were similar at a distance of a few hundred metres from the turbine, while they differed by more than 20 dB at a distance of 10 km [
24]. The authors also concluded that calculations using Harmonoise weather classes at long distances can lead to serious errors and recommend the use of actual meteorological profiles.
In conclusion, Keith et al. [
24] state that the ISO 9613-2 extended the frequency calculation method can be used to calculate the annual average SPL of ILFN when the area to be assessed is within a few kilometres of the nearest wind turbines. When calculations need to be made for larger distances from the turbine or specific meteorological classes, the authors recommend using the FFP method using real atmospheric properties.
Bertagnolio et al. [
27] investigated the possibility of using the coupling of the aeroelastic model HAWC2 and the so-called Formula 1A developed by Farassat to model the ILFN generated by a wind turbine operating upwind. The calculation results obtained with the above model were compared with other models used for modelling wind turbine noise, namely those of Viterna [
28,
29] and Amiet [
30].
Bertagnolio et al. [
27] analysed the influence of different factors and parameters on the level of the calculated ILFN at a point downstream of the turbine. The following were analysed: the influence of the time window of the Fourier analysis, the influence of the time step, the influence of the blade loading model, the influence of the tower and wind uplift, and the influence of the inflow turbulence and disturbance caused by the operation of the preceding turbine. The authors have shown that the proposed calculation method gives results that are consistent with other existing calculation methods. The main conclusion from these analyses is that the field disturbance caused by the preceding turbine has a significant effect on the estimate of the ILFN level generated by the turbine under study when the intensity of atmospheric turbulence is low. When the turbulence intensity is high, the influence of the field disturbance caused by the operation of the preceding turbine on the calculated ILFN level is low.
The computational model proposed in [
27] has been verified experimentally in papers [
31,
32]. The authors compared the results of the numerical calculations with the measured results for two test wind turbines located at the DTU-Riso station. All meteorological input parameters to the model were determined from measured data recorded at the site. The authors found good agreement between the measured and calculated results when background noise is neglected. The results of the proposed model correctly reproduce the quantitative increase in the ILFN as a function of wind speed; although, in some cases, differences with the measured data are apparent in some frequency ranges [
32].
However, research is still underway to improve methods for predicting wind turbine noise, as exemplified by the PIBE project [
33], which was run in France from 2019 to 2023.
In papers [
34,
35], the authors have pointed out that the modelling of wind turbine noise requires the consideration of the wind turbine as an extended noise source together with aeroacoustic phenomena [
36,
37,
38] as the main mechanisms for wind turbine noise generation. In contrast, the PIBE project showed that for distances up to 800 m in all propagation directions and for distances up to 1000 m for propagation directions from −120° to 120°, the effect of atmospheric turbulence on A-weighted sound pressure level predictions is negligible [
39].
In a publication by Mascarenhas et al. [
40], produced as part of the PIBE project, the authors presented the results of comparisons of wind turbine noise predictions (using a model combining Amiet’s theory with a wide-angle parabolic equation) with field measurements. The results of the model and measurements were first compared at points close to the wind turbine and then at distances from 350 m to 1300 m from the source. Based on the results presented in this publication, the authors concluded that the model combining Amiet’s theory with the wide-angle parabolic equation can be used for calculations in the frequency range from 100 Hz to 4 kHz. The authors also found that the turbulence dissipation rate parameter must be correctly estimated in order to accurately estimate the noise of a wind turbine at the receiving point for the model used. On the other hand, in publication [
41], the authors presented a comparison of calculation results using a model combining Amiet’s theory with a wide-angle parabolic equation with measured results in the 50 Hz to 10 kHz band for measurement points located upwind and downwind. The authors found that the model underestimated the measured results by about 2 dB for a measurement point at a distance of 1318 m from the nearest turbine in the frequency range below 200 Hz and above 400 Hz. This paper also shows that accurate prediction results can be obtained using this model for 1/3 octave band spectra averaged over 10 min for measurement points located up to 1300 m from a wind turbine.
All the above models are very complex and require specialised knowledge and the collection of input data that are very difficult or even impossible to obtain at the forecasting stage. These are mainly meteorological data that affect the generation of ILFN by the wind turbine (turbulence intensity, turbulence dissipation rate parameter) and its propagation (temperature gradient, wind speed gradient, wind direction in different parts of the atmosphere). For the above reasons, the engineering calculations commonly used in the prediction phase can be subject to very large uncertainties.
For the reasons outlined above, it is therefore necessary to use and improve existing computational methods (ISO 9613-2 [
42], CNOSSOS-EU [
43], and Nord2000 [
44,
45]) in engineering applications for modelling the ILFN generated by wind turbines.
Therefore, this paper identifies values for the differences between the calculated value (using three commonly used calculation methods) and the measured value of the sound pressure level generated by wind turbines in the 1/1 octave bands with centre frequencies from 4 Hz to 250 Hz. The sources of the differences are also identified to further work on these models to minimise these differences to a level acceptable for engineering applications.
4. Analysis of Results
The absolute values of the differences between the calculated and measured value
using the ISO 9613-2 method of the sound pressure level generated by the wind turbine for an average wind speed of 3.3 m/s (
Table 6) range from 0.0 dB (31.5 Hz, 500 m) to 13.5 dB (31.5 Hz, 1500 m). For the average wind speed of 4.2 m/s, they range from 0.1 dB to 13.2 dB for the 4 Hz band at 500 m and 1500 m, respectively. In contrast, for the highest average wind speed analysed, 4.6 m/s, the absolute values of
range from 0.2 dB (8 Hz, 1500 m) to 10.9 dB (63 Hz, 1000 m).
At the measurement points 250 m, 500 m and 1000 m from the turbine for octave bands with centre frequencies of 4 Hz, 8 Hz, 16 Hz (
Table 3,
Table A1 and
Table A3), and for the band 4 Hz at a measurement point 1000 m from the turbine (
Table 5,
Table A9 and
Table A11), the difference between the calculated and measured value of the sound pressure level was not determined. This is because the ILFN emission level generated by the wind turbine was not determined at these measurement points for the aforementioned frequency bands because the recorded value of the background sound pressure level was higher than the sound pressure level recorded during turbine operation.
The measurement and calculation results obtained for each calculation method at all measurement points for an average wind speed of 3.3 m/s are shown in
Figure A1, while the values of the differences
between the calculated value
and the measured value
are shown in
Figure 5. In the ideal case of agreement between the model results and the measured values, each curve representing a given model should be flat and its value should be 0 dB. It can be seen that the
values for each method at a given measurement point form parallel lines. Only for the 250 m point does the line of difference obtained for the ISO 9613-2 method cross the line for the CNOSSOS-EU method, and for the 63 Hz–250 Hz frequency bands, the
values for the ISO 9613-2 method are smaller than for the CNOSSOS-EU method. For a measurement point 500 m away, the
values for the ISO 9613-2 and Nord2000 methods from the 31.5 Hz frequency band approach each other and for the 250 Hz band are very similar. For measurement points 1000 m and 1500 m away, the lines move closer together as the frequency of the band analysed increases, and in the 250 Hz band, the
values for each method are very similar.
The
values determined for the prediction results using the CNOSSOS-EU algorithm for favourable propagation conditions at an average wind speed of 3.3 m/s range from 1.9 dB (250 Hz, 500 m) to 11.9 dB (31.5 Hz, 1500 m). For an average wind speed of 4.2 m/s (
Table 7), they range from 0.1 dB to 11.5 dB for the 31.5 Hz band at 500 m and the 4 Hz band at 1500 m, respectively. In contrast, for an average wind speed of 4.6 m/s, the absolute values of
range from 0.2 dB (16 Hz, 500 m) to 9.6 dB (63 Hz, 1000 m).
The
difference values obtained at a given measurement point for each calculation method at an average wind speed of 4.2 m/s are shown in
Figure 6, while the results of the measurement and calculation are shown in
Figure A2. As in the case of an average wind speed of 3.3 m/s, the curves
are parallel and, for a measurement point 250 m away, the line of difference values obtained for the ISO 9613-2 method intersects the line for the CNOSSOS-EU method. As a result, the
values for the ISO 9613-2 method are smaller than those for the CNOSSOS-EU method in the 63 Hz, 125 Hz and 250 Hz frequency bands. For a measurement point at a distance of 500 m, the
values for the ISO 9613-2 and Nord2000 methods from the 31.5 Hz frequency band approach each other and for the 250 Hz band are almost identical. On the other hand, at 1000 m and 1500 m, the lines move closer together as the frequency of the band analysed increases, and at 250 Hz, the
values for each method are very similar.
The absolute values of the differences between the calculated and measured values
using the Nord2000 method of the sound pressure level generated by an operating wind turbine for an average wind speed of 3.3 m/s range from 0.0 dB (4 Hz, 1500 m) to 6.9 dB (31.5 Hz, 1500 m). At an average wind speed of 4.2 m/s, they range from 0.3 dB to 7.6 dB at the measurement point at a distance of 1000 m, for the 16 Hz and 31.5 Hz frequency bands, respectively. On the other hand, for an average wind speed of 4.6 m/s, the absolute values of
(
Table 8) range from 0.1 dB (125 Hz, 250 m) to 11.3 dB (4 Hz, 500 m).
The values of the differences
between the calculated value of
and the measured value of
obtained at measurement points at distances of 250 m, 500 m, 1000 m and 1500 m for each calculation method at an average wind speed of 4.6 m/s are shown in
Figure 7, while
Figure A3 shows the measured and calculated results for the same points. As with the two previous average wind speeds analysed, the curves
shown in
Figure 7 are parallel up to a frequency band of 63 Hz and then begin to converge as the frequency increases. Only at a point 250 m from the turbine does the
curve obtained for the ISO 9613-2 method cross the curve for the CNOSSOS-EU method. The
values obtained for the CNOSSOS-EU method in the frequency bands 63 Hz–250 Hz are higher than those obtained for the ISO 9613-2 method. For the 500 m measurement point, the
values for the ISO 9613-2 and Nord2000 methods in the 250 Hz band are almost identical, while for the 1000 m and 1500 m points, the
values for the 250 Hz frequency band are very similar for all three calculation methods analysed.
Non-parametric statistical tests were performed to check for statistically significant differences between the measured results and the predicted results obtained for the calculation methods used. These tests were carried out on 12 sets of data: all the measurement points (4 points) and all the average wind speeds analysed (3 average wind speeds) for the whole frequency range analysed.
First, the Kruskal–Wallis test (
Section 2.3.1) was performed at a significance level of
. The results of this test clearly showed that there were no statistically significant differences between the measurement results and the prediction results at the significance level chosen. The obtained test probability values ranged from 0.08 to 0.60. However, a non-parametric Tukey–Kramer multiple comparison test (
Section 2.3.2) was also performed to determine which of the calculation methods used produced results most similar to the measurement results. Analysis of the results of the Tukey–Kramer test showed that in eight cases the results obtained using the Nord2000 method were closest to the measurement results. The results obtained using the ISO 9613-2 method were closest to the measured results in three cases (250 m, 4.6 m/s; 500 m, 3.3 m/s and 4.2 m/s) and for the CNOSSOS-EU method in one case (500 m, 4.6 m/s).
The same test procedure was performed for the value of the differences
between the calculated and measured results. In this case, it was found that statistically significant differences for
occurred in 8 of the 12 data sets analysed (marked in bold in
Table 9). In the remaining four cases, there were no statistically significant differences. The test probability values obtained from the Kruskal–Wallis test performed on the
differences are shown in
Table 9.
A non-parametric test for multiple comparisons was also performed at a significance level of to determine for which method the differences in were closest to zero (the ideal case of agreement between the model results and the measured values). The results of this test are very similar to those of the previous test procedure performed for predicted and measured results. The differences closest to zero were obtained using the Nord2000 calculation method in the nine data sets analysed. The ISO 9613-2 method had difference results closest to zero in two data sets ( 500 m, 3.3 m/s; 500 m, 4.2 m/s), while CNOSSOS-EU had only one (500 m, 4.6 m/s).
Analysing the results of the tests carried out, it can be concluded that the ISO 9613-2 method performs best for distances up to 500 m from the turbine for all the average wind speeds analysed, while in other cases, the Nord2000 method gives much more accurate results.
Analysing the measured
and calculated
sound pressure levels shown in
Figure A1,
Figure A2 and
Figure A3, it can be seen that the shape of the octave spectrum is well reproduced by all three calculation methods analysed for all cases.
In order to confirm this objectively, a correlation analysis was carried out between the measured and calculated ILFN values generated by the wind turbine for all cases analysed. For this purpose, the Spearman correlation coefficient
was determined, as the acoustic data do not have a normal distribution [
74].
The values of the correlation coefficient
are shown in
Table 10 and they range from 0.79 to 1.00. This means that there is a strong and very strong positive correlation between the measurement and calculation results. Based on the values of the Spearman correlation coefficient, it is clear that all the computational models analysed reproduce the shape of the ILFN octave spectrum very well in all the cases analysed.
The values of the correlation coefficients determined between measured and calculated values have the same value for each calculation method in almost all cases analysed. Only at a distance of 250 m from the turbine, for an average wind speed of 4.6 m/s, the value of the correlation coefficient between measured and calculated values for the ISO 9613-2 method is slightly lower than for the other two methods.
The suspicions arising from the analysis of
Figure A1,
Figure A2 and
Figure A3 were objectively confirmed by the correlation analysis carried out between the measured and calculated results.
5. Summary and Conclusions
In this paper, a measurement verification of the prediction results of the ILFN generated by a wind turbine obtained by the ISO 9613-2, CNOSSOS-EU and Nord2000 methods has been carried out. The aim was to test the suitability of these models for the determination of sound pressure levels in octave frequency bands in the range from 4 Hz to 250 Hz.
For this purpose, a geometric-acoustic model of a wind farm operating in central Poland was created in the SoundPlan software. The SoundPlan software was also used to calculate the sound pressure levels using the above calculation methods.
The prediction results were compared with the actual measurements taken at the wind farm in question. Calculations and measurements were carried out for points located at distances of 250 m, 500 m, 1000 m and 1500 m behind the turbine at 0 m above ground level. The calculation points and the measurement microphones were placed on a measuring plate according to the guidelines of IEC 61400-11:2012/AMD1:2018 [
46] with a single windscreen. The analysis was carried out for octave bands with centre frequencies from 4 Hz to 250 Hz for three average wind speeds (3.3 m/s, 4.2 m/s and 4.6 m/s) recorded at 10 m above ground level.
The absolute values of the differences between calculated and measured results for the ISO 9613-2 method range from 0.0 dB (3.3 m/s, 500 m, 31.5 Hz) to 13.5 dB (3.3 m/s, 1500 m, 31. 5 Hz), for the CNOSSOS-EU method from 0.1 dB (4.2 m/s, 500 m, 31.5 Hz) to 11.9 dB (3.3 m/s, 1500 m, 31.5 Hz) and the Nord2000 method from 0.0 dB (3.3 m/s, 1500 m, 4 Hz) to 11.3 dB (4.6 m/s, 500 m, 4 Hz).
The differences
shown in
Figure 5,
Figure 6 and
Figure 7 obtained for each method form parallel lines over the whole frequency band analysed, regardless of the average wind speed used in the calculations. Overall, it can be concluded that the Nord2000 model overestimates the results by an average of 0.8 dB, while the ISO 9613-2 and CNOSSOS-EU models underestimate the results by an average of 3.4 dB and 3.8 dB, respectively.
It was also noted that the prediction results obtained using the ISO 9613-2 method are 3 dB higher than those obtained using the CNOSSOS-EU method for each octave frequency band analysed between 4 Hz and 250 Hz, but only at the calculation point 500 m downstream of the turbine regardless of wind speed.
The differences between calculated and measured values depend on several factors. The most important are the following:
A different height of the source location above ground level (105 m) than assumed in the calculation methods (according to the ISO 9613-2 and CNOSSOS-EU calculation methods, the maximum height of the source location should not exceed a value of 30 m);
A different value of the ground coefficient G used in the calculations than in reality, due to the angle of incidence of the sound wave on the ground;
The difference between the actual values of the sound absorption coefficients of the ground and the atmosphere and the values assumed for the calculations;
In the ISO 9613-2 and CNOSSOS-EU models, it is not possible to take into account the direction and speed of the wind;
It is not possible to take into account the wind speed and direction for the different atmospheric layers present in the propagation path of an acoustic wave.
Non-parametric statistical tests were performed at a significance level of to determine whether there were statistically significant differences between the predicted and measured results of the ILFN generated by the wind turbines. A total of 12 data sets were analysed. The calculation results obtained for the Nord2000 method were closest to the measurements in eight cases, for the ISO 9613-2 method in three cases and for the CNOSSOS-EU method in only one case.
The same test procedure was also applied to the value of the differences between the calculated and measured results. Once again, 12 sets of data were analysed. In nine cases the differences were closest to zero for the Nord2000 method, in two cases for the ISO 9613-2 method and only one case for the CNOSSOS-EU method.
As a result of the tests, it was found that the most accurate results were obtained using the Nord2000 calculation method, and the least accurate results were obtained using the CNOSSOS-EU model.
A correlation analysis was also carried out between the results of the calculations for each method and the results of the measurements . The Spearman’s rank correlation coefficients were determined for this purpose. All models showed a strong or very strong correlation with the measurement results, which means that they represent the shape of the spectrum in the analysed frequency range very well.
It is not easy to verify the calculation results obtained using the ISO 9613-2 and CNOSSOS-EU models in specific wind conditions by measurement. This is because these models do not take into account the wind speed and direction recorded during the measurements. The ISO 9613-2 method assumes that there are favourable propagation conditions in all directions from the source to the receiver. The CNOSSOS-EU method, on the other hand, allows calculations to be made for homogeneous or favourable sound propagation conditions, assuming that the same propagation conditions prevail in all directions from the source to the receiver. Favourable propagation conditions in all directions do not exist in reality. Homogeneous conditions are sometimes encountered.
Such limitations are not present in the assumptions of the Nord2000 method. The model allows the calculation of sound pressure levels under specific meteorological conditions. To perform such calculations, knowledge of wind direction and speed, air pressure, relative humidity and air temperature is required, as well as knowledge of parameters such as roughness length, temperature gradient, standard deviation of temperature gradient, structure parameter of turbulent temperature fluctuations, standard deviation of wind speed and structure parameter of turbulent wind speed fluctuations.
Taking into account all the analyses carried out, as well as the amount and availability of data input to the model for predicting ILFN generated by wind turbines, it was concluded that the Nord2000 model can be used for calculations with increased accuracy (on average, the results are overestimated by 0.8 dB), but with high labour intensity. This model requires input data that are difficult to obtain. These are mainly the meteorological data mentioned in the previous paragraph.
On the other hand, in the case of the ISO 9613-2 model, the collection of the required input data is not problematic, but the accuracy of the prediction results is lower than in the case of the Nord2000 model. The results obtained using the ISO 9613-2 model are underestimated by an average of 3.4 dB. This model can therefore be used for simplified calculations with lower accuracy.
In order to obtain accurate ILFN modelling results, it is necessary to build an accurate digital twin of the wind farm to be analysed, taking into account all important elements that affect the attenuation of sound during propagation outdoors (landform, land cover, land use). All parameters required by a given calculation method (turbine sound power level, turbine geometric dimensions, turbine operating schedule, meteorological conditions prevailing in a given area) should be determined with the highest precision available. A very important element influencing the accuracy of the calculations is the determination of the sound attenuation coefficient by the atmosphere for each frequency band analysed. For this purpose, its value should be determined based on Equations (3)–(5) given in the ISO 9613-1.
The research results presented in this article relate only to the measurement verification of ILFN calculations generated by a working wind turbine under specific meteorological conditions for three average wind speeds. The research results and analyses presented indicate that these computational models can be successfully used to predict ILFN from wind turbines. Therefore, it seems reasonable to carry out a measurement verification of the ILFN prediction results obtained with the ISO 9613-2, CNOSSOS-EU and Nord2000 methods for a longer period (e.g., one year), taking into account the different sound propagation conditions. This research may help to identify a model that should be used to determine long-term noise hazard indicators
, as well as indicators describing the harmful effects of noise, which include Ischaemic Heart Disease (IHD), High Annoyance (HA) and High Sleep Disturbance (HSD), in accordance with Commission Directive (EU) 2020/367 [
75].