Effect and minimization of errors in in-situ ground impedance measurements Roland Kruse, Volker Mellert Oldenburg University, Inst. of Physics, 6 Oldenburg, Germany Abstract The transfer function method is a procedure to measure the surface impedance of grounds in-situ. In this article, the influence of measurement errors on the predicted surface impedance is investigated numerically. Even small errors in the range of accuracy of common measurement equipment can lead to significant errors in the impedance. This is especially true for errors in the transfer function at frequencies below about 5 Hz and for highly reflecting grounds. To a lesser degree, errors in the measurement geometry contribute to the uncertainty of the estimated impedance. To minimize these effects, an improved geometry is suggested for the frequency range 4 Hz significantly reducing the average error. However, even with this optimized geometry the average error for high impedance grounds, like compacted silt, in this frequency range will be around 5 %. Therefore, the use of the transfer function method can t be recommended in this case unless the requirements in accuracy are very low for a specific application or particular favourable measurement conditions are given.
Keywords: Ground impedance; In-situ impedance measurement PACS 43.58.Bh Introduction The surface impedance of outdoor grounds is an important parameter for the prediction of the sound propagation, especially for small distances when the effect of meteorological factors is often low compared to the ground effect. Generally, the measurement of the ground impedance can only be done using in-situ methods because small samples for use in the impedance tube [] do not correctly represent the properties of the extended and sometimes layered / inhomogeneous outdoor ground. Commonly used for in-situ impedance measurements both in room acoustics and for outdoor use - are impulse-echo ( Adrienne ) methods [,3] and the two-microphone (transfer function) technique [4, Iterative Procedure ]. While the impulse-echo methods, because of their simple theoretical background, are easy the implement, they rely on a plane wave like reflection and are therefore not adequate for low frequencies and small or even grazing incidence angles. The two-microphone method, on the other hand, is commonly used in combination with a sound field model taking the spherical nature of the waves at low frequencies into account and is therefore, theoretically, well suited for that purpose. A variant of this method is used in the American standard ANSI S.8 [5]. In this standard, the magnitude of the transfer function between two microphones situated above the ground is measured using one of three predefined geometries. The user will then (visually) compare this level difference spectrum with spectra provided in the standard and choose the one that fits best over the whole frequency range. These spectra ( templates ) were derived using a sound field model and a surface impedance from either a one or a two- parameter absorber model. Thereby, the user will obtain absorber parameters by which the surface impedance can easily be calculated.
Currently, a new ANSI standard is proposed (working group S/WG ) to allow a direct deduction of the ground impedance without the use of templates. This procedure, while requiring a sophisticated data evaluation, would have the advantage of being independent of any absorber model. The use of the transfer function method in such a standard raises concerns about the robustness of the method, because the method may become more widespread and may be used by persons less experienced in this area of acoustics. The influence of selected measurement errors on the predicted surface impedance is investigated in this paper. Small errors in the measurement geometry can hardly be avoided. Distance measurements, e.g. with a tape measure, have limited accuracy and the ground surface is not well defined. Another source of errors is the fact that the model requires the knowledge of the position of the acoustical centre of the loudspeaker which can be different from its geometrical centre [6]. Errors may not only occur in the geometry but also in the measured transfer function T. Phase errors arise from the fact that standard microphones are not phase-matched. Errors in the magnitude of T can arise from slight differences in calibration and frequency response of the microphones. Both errors may also follow from unwanted reflections, e.g. from the measurement equipment or other reflecting objects nearby as well as from meteorological influences (wind, temperature) causing spatial or temporal changes in the speed of sound. The necessary measurement precision is derived from the error predictions. Subsequently, an optimization is done to select the best measurement geometry for a given frequency range and surface impedance. All impedances presented are normalized to the impedance of air and assume an exp(iωt) time dependence. Experimental procedure Measurement set-up Measurements of the two-microphone transfer function were done on a (grass covered) soccer field. The temperature was 7 C and the wind speed < m/s. The background noise
was below 45 db(a). The measurement set-up is shown in fig.. The sound source was a cm loudspeaker in a (3 cm)³ closed cabinet. The emitted signal was pink noise with a level of not less than 8 db(a) at the microphone positions to ensure the signal was well above the background noise. B&K 489 microphones (/, IEC 65 Type ) with windscreens were used. The transfer function T as defined in eq. was determined using B&K PULSE, p representing the sound pressure at the two microphone positions. It is equal to the ratio of the velocity potentials φ. T p( upper) φ( upper) = p( lower) φ( lower) = () Five measurements at different locations were done. Geometry B from ANSI S.8 was used as it showed the best overall performance on this type of ground: Source height (h s ) = cm, upper microphone height (h ru ) = cm, lower microphone height (h rl ) = 5 cm, source- receiver distance (R) = m. Sound field model and impedance deduction A widespread model for the reflection of spherical waves on a locally reacting impedance plane [7] was used for the calculation of the surface impedance from the two-microphone transfer function. This model assumes an exp(-iωt) time dependence. The velocity potential φ is a function of the lengths of the direct path r and reflected path r, the wave number k, the angle of incidence θ and the normalized surface admittance β. It is defined by eq. - 5. with e = r ikr + [ R p + ( R p e ) F] r ikr φ () F = + i π λe λ ² erfc( iλ) (3) R p is the reflection coefficient for plane waves (eq.4).
p sin( θ ) β = sin( θ ) + β R (4) The numerical distance λ is defined in eq.5 representing a common simplification of the definition by [8]. λ = ikr (sin( θ ) + β ) (5) The surface impedance Z was calculated from the measured transfer function T (eq.) and geometry by the Newton-Raphson algorithm which finds the zero of the function (predicted T observed T) and therefore gains the surface impedance Z which minimizes the difference between the observed transfer function and the transfer function predicted by the model. This very efficient method has been described in [9]. The calculation was started at the lowest frequency point ( Hz) with a seed value for the iteration of β =.5. The obtained solution is used as seed value for the next higher frequency point. While the iteration may, in general, sometimes converge on a wrong solution (local minima), such a behaviour was not observed even if the initial seed value was varied. For all calculations, Matlab R6b was used. The Matlab code can be found on the authors webpage. Effect of errors To determine the influence of errors in the geometry or transfer function on the estimated surface impedance, a two-parameter model (eq.4,5) [5] was (least square) fitted to the impedance using data from 4 4 Hz, averaged over /3 rd octave bands, to obtain an estimate Z of the impedance at low frequencies. Im( Z σ e Re( Z / ρc ) = (4) πγρ f / ρ c σ f c α e e ) = + (5) πγρ 8πγf http://www.physik.uni-oldenburg.de/aku/
ρ and c c are the density and speed of sound in air, γ the ratio of specific heats, σ e the effective flow resistivity and α e a parameter representing an effective rate of change of porosity with depth. Now, the transfer function T was calculated from this impedance using geometry B (or geometry B with slight errors). T was then slightly changed in amplitude (±. db) or phase (±.5 ) and the surface impedance calculated. The predicted impedance with measurement errors is then compared to the nominal impedance Z. These error estimated were chosen to account mainly for the differences in the microphone frequency responses which, in a free sound field, were determined to be. db resp..5, on average. Additional. db were added to allow for other (external) influences. Geometry optimization To find an optimized geometry for the lower frequency range, which is least sensitive to measurement errors, the following procedure was used: ) Bounds for the geometry were chosen to allow for an easy set-up and keeping the distances small to minimize the influence of meteorological factors. ) The achievable precision of the measurement (resp. the measurement errors) was defined. 3) The magnitude of the estimated surface impedance Z was calculated for all geometries on a 5 cm grid. 4) Because of the possible nonlinear relation between the errors and their effect on Z, it is not sufficient to calculate the error in Z only for the maximum error estimate. The calculation was done for four values within the error range, e.g. for the error in T, estimated to be ±. db, the values -., -.,.,. db were chosen. The resulting error E x was defined as the standard deviation of Z for the four error values, calculated at each frequency point with the predictand Z. 5) To allow for an efficient calculation, the effects of the errors were considered to be independent except for the magnitude and phase of the transfer function. This reduces
the number of calculations for each frequency point and geometry to 56 instead of 496 needed if all error combinations were considered. 6) Step 3 was repeated for three surface impedances Z (low, medium, high) from twoparameter models [5]: Pine forest floor (σ e = 7.5 kpas/m²,α e = 6 m - ), soccer field (σ e = kpas/m²,α e = 4 m - ) and compacted silt (σ e = 4 kpas/m²,α e = -5 m - ). The calculation was done for the frequency range 4 Hz divided into /3 rd octave steps, the number of frequency points N thereby being seven. 7) In analogy to the Gaussian error propagation law, the (relative) average total error ATE for each geometry was defined as the geometrical sum of the single errors E x divided by the true impedance Z - averaged over the frequency range f (eq.6). ATE = N f E T ² + E hru ² + E Z hrl ² + E hs ² + E R ² (6) The optimal geometry - observing the boundary conditions in table shall be the geometry with the lowest ATE. Results The estimated surface impedance of the soccer field (average of five measurements) is shown in fig.. While the course of the impedance is rather smooth and does agree well with impedance models for grass-covered ground [] for frequencies above about 4 Hz, the decrease in impedance at lower frequencies does not at all agree with models for porous absorbers. Effect of errors The fitted two-parameter model on which the error analysis is based is shown in fig.. The obtained absorber parameters σ e = kpas/m² and α e = 4 m - are comparable to the values σ e = 83 kpas/m² and α e = 4 m - which are stated for lawn in the ANSI standard.
In fig.3, the effect of an error of. db resp..5 in the transfer function is shown. Even small errors in T cause high errors in the estimated surface impedance at low frequencies. It should also be noted that the magnitude of the error does not necessarily increase continuously with decreasing frequency. In fig.4, the effect of an error in the microphone positions is presented. The effect is much smaller than the effect of an error in T, but errors in the upper microphone position can still cause significant deviations. Not shown but also observed is the fact that errors in the source height and source- receiver distance have only a minor effect on the surface impedance estimate (for geometry B and the assumed impedance of the soccer field). The large effect of transfer function errors will now be further analysed. In fig. 5 and 6, the relative error in the real and imaginary part of the impedance is shown in relation to magnitude and phase errors of T. The effect is not only large but also asymmetric: For most of the error combinations the real part of the impedance will be too low and the imaginary part too high. Furthermore, there s a strong interaction between the effect of magnitude and phase errors. From fig.5 and 6 as well as likewise diagrams for the combinations h ru, h rl and h s, R, the following estimates for the necessary measurement precision (geometry B, 5 Hz) are derived, if the error is to be smaller than % in both real and imaginary part and no interactions between the three error combinations are considered: T ±.3 db, phase(t) ±., h ru ± cm, h rl ± cm, h s ±.5 cm, R ± 5 cm A likewise analysis revealed that for high surface impedances or lower frequencies the necessary precision becomes even higher. In case of high surface impedances this is not only true for the transfer function but also for the geometry. Matlab code for generating these diagrams can be found on the authors website.
Geometry optimization The results of the optimization of the geometry for the frequency range 4 Hz for three different surface impedances can be found in tab.. For all impedances, the optimum values R = 3 m, h s =.5 m and h rl = 5 cm are the same, while the upper microphone height varies between m and.8 m. The effect of this optimization on the average total error ATE is shown in fig.7. In comparison with the predefined geometries from the ANSI standard, the average error can be reduced significantly for all three surface impedances to about half of the value of the best of these geometries. However, for the high surface impedance the ATE is still 53%. In fig.8, the effect of an error of. db resp..5 in the transfer function is shown, comparable to fig.3 for geometry B. The error is significantly reduced, esp. in the imaginary part. Discussion Taken into account the results of the error analysis, it seems that the often reported decrease of the predicted surface impedance at low frequencies is the result of measurement uncertainties. The investigated soccer field was large, flat and homogeneous, it therefore resembles one of the best measurement location one can hope to encounter during ground impedance measurements. On the other hand, small errors in the transfer function result in deviations from the expected surface impedance as observed in the measurement, with the real part being too low and the imaginary part too high. To a lesser degree, this is also true for errors in the geometry, esp. the upper microphone height. To obtain acceptable results at low frequencies with geometry B, which is well suited for frequencies above 4 Hz and ground with medium impedance, the necessary measurement precision for an error less than % is already high at 5 Hz and can only be achieved with very great care and high quality equipment including phase-matched microphones. Fortunately, the effect of errors in the geometry on the impedance is less asymmetric compared to the effect of errors in the transfer function and can therefore be reduced by
averaging multiple measurements on the same ground including repositioning of source and receivers. The optimized geometry for the frequency range 4 Hz does feature a significant advantage over the predefined geometries with respect to its error sensitivity, but the average error for high impedance surfaces is still high. The very high error for geometry C in this case does agree with the standards recommendation to use this geometry only for very soft grounds. It should be kept in mind that the necessary measurement precision itself depends on how accurate the surface impedance needs to be measured and therefore on the specific application. In general, the effect of a given percental error in the impedance is highest near the first ground dip, the frequency range in which the presence of the ground leads to a high attenuation compared to the free field situation. The position of the ground dip depends on the impedance and geometry. As an example, the ground dip above the soccer field with source and receiver at m height and the distance between them being m is located around 75 Hz. In this case, an error of -5 % in the surface impedance leads to an underestimation of the sound pressure at 4 Hz by 7 db, while at khz the sound pressure is overestimated by only 3 db. Below 4 Hz, the error decreases and reaches.3 db at Hz. For higher propagation distances, the error does decrease slower with decreasing frequency. Thus, the frequency range below khz should receive particular attention regarding errors in the predicted surface impedance. Conclusion The two-microphone method is an established procedure for the in-situ measurement of the surface impedance. A disadvantage is that it is sensitive to errors both in the measured transfer function and, to a minor degree, to errors in the geometry. The effect of these errors depends on the measurement geometry and the surface impedance. For low impedance surfaces, it is possible to measure with acceptable accuracy at frequencies > Hz even with the predefined geometries from the ANSI standard. Unfortunately, such highly absorbing
materials rarely occur in nature with the important exception of snow. For materials with higher flow resistivity, the use of the predefined geometries can t be recommended for frequencies below about 5 Hz. The optimization procedure, on the other hand, leads to a geometry which provides a significant higher error tolerance and should enable the researcher to obtain reasonable results down to Hz. However, for materials with a very high surface impedance like compacted silt or asphalt, even the optimized geometry cannot guarantee low errors. Taken in mind that the proposed source-receiver distance is 3 m and therefore the maximum allowed distance for the optimization it is possible that better geometries with larger distances exist. Larger distances, however, would be more affected by meteorological factors and require large surfaces. Therefore, the use of the two-microphone method is not recommended in such situations unless the requirements in accuracy are low for a specific application or very favourable measurement conditions (large, flat surfaces, no wind) and high quality equipment including phase- matched microphones are available. Acknowledgement The authors would like to thank their student Ping Rong for his helpful comments on the data processing and interpretation and his valuable Matlab programming. References. ISO 534-. Determination of sound absorption coefficient and impedance in impedance tubes - Part : Transfer-function method. International Organization for Standardization. 998. Mommertz E. Angle-Dependent In-situ Measurement of Reflection Coefficients Using a Subtraction Technique. Applied Acoustics (995); 46: 5-63
3. CEN/TS 793-5. Road traffic noise reducing devices - Test method for determining the acoustic performance - Intrinsic characteristics - In situ values of sound reflection and airborne sound insulation. European Committee for Standardization. 3 4. Allard JF, Champoux Y. In Situ Two-Microphone Technique for the Measurement of the Acoustic Surface Impedance of Materials. Noise Control Engineering Journal (989); 3(): 5-3 5. ANSI S.8-999 (R4). Template Method for Ground Impedance. American National Standards Institute. (4) 6. Fuhs S, Höldrich R, Tomberger G. Validierung des Entfernungsgesetzes und Korrektur der Gruppenlaufzeit und des akustischen Zentrums des Lautsprechers im Adrienne- Verfahren. DAGA 6. Braunschweig. 7. Nobile MA, Hayek SI. Acoustic propagation over an impedance plane. Journal of the Acoustical Society of America (985); 78(4): 35-336 8. Ingard U. On the Reflection of a Spherical Wave from an Impedance Plane. Journal of the Acoustical Society of America (95); 3(3): 39-335 9. Taherzadeh S, Attenborough K. Deduction of ground impedance from measurements of excess attenuation spectra. Journal of the Acoustical Society of America (999); 5(3): 39-4. Donato RJ. Impedance models for grass-covered ground. Journal of the Acoustical Society of America (977); 6(3): 449-45
Figures R = m Upper mic. Loudspeaker hs = cm θ Ground hru = cm hrl = 5 cm Lower mic. Fig. : Measurement set-up for the transfer function method. Geometry B from ANSI S.8.
5 5 Impedance -5 - -5 - -5 3 5 4 Frequency [Hz] Fig. : Predicted surface impedance of a grass covered soccer field. Real part ( ), imaginary part ( ), fitted two-parameter model ( ).
3 Impedance - - -3-4 3 5 4 Frequency [Hz] Fig. 3: Effect of errors in the transfer function on the estimated surface impedance for geometry B. Correct impedance ( ), transfer function increased by. db ( ), phase error of +.5 ( ).
3 Impedance - - -3 3 5 4 Frequency [Hz] Fig. 4: Effect of errors in the microphone position on the estimated surface impedance for geometry B. Correct impedance ( ), upper microphone cm too high ( ), lower microphone cm too high ( ).
Error in magnitude of the transfer function [db].5..5..5 -.5 -. -.5 -. -.5. -.. -.5.5.5. -. -.5.5.5. -. -.5 -.5. -. -.5 - - -.5.5 Error in phase of the transfer function [deg] Fig. 5: Effect of errors in the magnitude and phase of the transfer function on the real part of the predicted surface impedance for geometry B at 5 Hz. Relative error ( = %). Error in magnitude of the transfer function [db].5..5..5 -.5 -. -.5 -. -.5. -.. -.5.5.5. -. -.5.5.5. -. -.5 -.5. -. -.5 - - -.5.5 Error in phase of the transfer function [deg]
Error in magnitude of the transfer function [db].5..5..5 -.5 -. -.5 -. -.5.5.. -...5 -.5. -. -. - -.5..5 - -.5. -. -.5 - - - -.5.5 Error in phase of the transfer function [deg] Fig. 6: Effect of errors in the magnitude and phase of the transfer function on the imaginary part of the predicted surface impedance for geometry B at 5 Hz. Relative error ( = %). Error in magnitude of the transfer function [db].5..5..5 -.5 -. -.5 -. -.5.5.. -...5 -.5. -. -. - -.5..5 - -.5. -. -.5 - - - -.5.5 Error in phase of the transfer function [deg]
.5 Average total error.5 Z_low Z_medium Z_high Fig. 7: Average total error (relative) in the predicted surface impedance for three surface impedances and geometries A C [5] as well as the optimized geometry. Geometry A ( ), B ( ), C ( ), optimized geometry ( )..5 Average total error.5 Z_low Z_medium Z_high
4 3 Impedance - - -3-4 3 5 4 Frequency [Hz] Fig. 8: Effect of errors in the transfer function on the estimated surface impedance for the optimized geometry. Correct impedance ( ), transfer function increased by. db ( ), phase error of +.5 ( ).
Tables Range [m] Error estimate Transfer function - ±.db, ±.5 Lower microphone height.5 - ± cm Distance upper lower mic.. - ± cm Source height..5 ± cm Source-receiver distance - 3 ± 5 cm Tab. : Boundary conditions for the geometry optimization. Ranges of the geometry and estimates of the achievable measurement precision. Low impedance Medium High impedance impedance Lower microphone height 5 cm 5 cm 5 cm Upper microphone height m.95 m.8 m Source height.5 m.5 m.5 m Source-receiver distance 3 m 3 m 3 m Tab. : Results of the geometry optimization for the frequency range 4 Hz for three different surface impedances.