Energy and Water Savings during Backwashing of Rapid Filter Plants

: This paper describes an analysis of the effects of adjusting the intensity of ﬁlter backwash to the water temperature. The consequences of the lack of such adjustment for the life of ﬁlter beds, the amount of water used for backwashing, the amount of energy used for backwashing and the quality of the ﬁrst ﬁltrate are presented. In order to determine the losses and proﬁts resulting from controlling the intensity of backwash water depending on its temperature, an analysis was carried out at a water treatment plant in southern Poland. Laboratory measurements were used to determine the granulation and speciﬁc gravity of sand grains ﬁlling the ﬁltration beds. On the basis of measurements on a semi-technical scale, the magnitudes of ﬁlter bed expansion were determined for average monthly wash water temperatures. They were ﬁrst calculated from the Richardson–Zaki equation, using different formulae for the value of the exponent of the power in this equation. Due to signiﬁcant differences in the density and shape of grains covered with a permanent deposit after several years of ﬁlter operation, a satisfactory match between the formulae known from the literature and the results of expansion measurements was not obtained. Therefore, an new formula for the bed expansion was developed based on the Richardson–Zaki equation. A good ﬁt of this formula to the experimental results was obtained. Monthly average values of water temperature were compiled, and on this basis the required amount of backwash water and energy was computed. The computations were made for 25% of ﬂuidized bed expansion. Possible energy and water savings were estimated, as well as further gains from keeping the required expansion of the porous bed constant regardless of the wash water temperature.


Introduction
Backwashing of rapid filters used in water treatment plants was traditionally carried out in the USA without the use of compressed air [1][2][3][4]. Since flushing with water alone is not very effective [5][6][7][8], it was necessary to use various procedures to support this effectiveness. However, for a long time in Europe, backwashing started with vigorous mixing of filter bed grains with compressed air, during which sudden changes in pressure gradient occurred and grains rubbed against each other [9][10][11]. First, however, the raw water inflow is closed and the water table in the filter chambers is allowed to drop to a height only a few centimeters above the top surface of the filter bed. As a result, the energy supplied to mix the grains of the filter bed is not dissipated in the water above this bed and the distances between grains remain small during "air washing", which favors frequent collisions leading to deposit abrasion. The most effective is the air-water flushing [12,13], during which the compressed air supply is closing and the backwash water supply is opening simultaneously. The compressed air supply has to be completely closed before the water table above the filter bed approaches the edge of the overflow. This is necessary to prevent the loss of part of the grain size of the filter bed during backwashing. Afterward, backwashing is carried out with water only. In this final stage of washing, the grains do not collide vigorously with each other and there is no effective removal of the sediment retained on the surface of the grains, but rather only the displacement of the water suspension that has formed during the two previous washing processes [14,15]. Backwashing is carried out with clean water, and the inlet to the filter is only closed when the backwash water flowing out is clean. In spite of the clarity of the water outflowing from a filter after backwashing, some of the removable sediment and durable deposits remain on the surface of the filtration bed grains, which have adhered so strongly to the grain surface that they are not only not removed during subsequent backwashing but are even difficult to remove partially during treatment of these grains with hydrochloric acid or a concentrated solution of hypochlorous acid [16]. Therefore, such treatments are not applied in practice. After closing the backwash water inflow, the fluidized grains inevitably fall to the bottom of the chamber and rub against each other [17][18][19]. This causes some of the removable sediment to be stripped off into the water filling the filter bed, as well as partially into the wash water remaining between the surface of the fixed filter bed and the edges of the overflow troughs [17,18]. When the filter is put back into operation, this water suspension causes poor quality in the first filtrate [15,17] and contributes up to 60% to 90% of all solid phase particles that enter the filtrate during the full filtration cycle. Several methods are used to counteract this unfavorable phenomenon [11,18,20]. Thus, it can be concluded that backwashing of rapid filters requires the consumption of treated water and energy and, moreover, it contributes to the deterioration of the quality of the first filtrate and, if carried out with too high intensity, may cause loss of grains from the filtration beds. In addition, it sometimes even contributes to a building catastrophe, resulting in the need for a general reparation of some filters. This last issue is not the subject of this article. Instead, attention will be focused here on all other problems mentioned.
The expansion of the filter bed at the same backwash intensity is higher in winter when the water viscosity is higher and lower in summer. These differences are greater the finer the grains in the bed are and the lower the density of the mineral of which they are composed. Thus, using the same backwash intensity throughout the year may lead to excessive energy expenditure and additional consumption of water and the danger of losing grains of the filtration bed in winter, or to keeping too small expansion of the bed in summer and thus to obtaining only partial washing out of contaminants contained in the bed. The selection of the right backwash intensity depending on the average temperature of the water in the individual months of the year is important both for water and energy demand and for the quality of drinking water. In addition, reducing the amount of heavily contaminated wash water reduces the costs associated with the treatment of these waters and the treatment of the resulting sludge [21][22][23][24]. Proper adjusting of the backwash intensity is important for preventing media grain loss [25,26].
The aim of the research was to develop and experimentally verify a computational algorithm allowing us to select the backwash intensity of rapid filters and granular activated carbon (GAC) filters, guaranteeing minimum water demand during fluidization and minimum energy consumption while maintaining the required level of bed expansion, guaranteeing an appropriately high degree of bed purification and at the same time not causing the loss of fine grains. The developed algorithm seeks to select parameters for the backwash intensity of filter in particular months of the year. It was assumed that the algorithm would be a tool to quickly estimate the annual gains for the filter station resulting from adjusting the backwash intensity. An analysis of one of the larger drinking water filter stations in southern Poland was performed to answer the question of how much reduction in water and energy demand could be expected by adjusting the backwash intensity of filter to the average water temperatures in the different months of the year in the climate prevailing in Poland.

Calculation Algorithm
Amirtharajah [7] and Kawamura [27,28] and later Sakkas and Lekkas [29] derived by means of analytical methods the equation for the backwash gradient, which is equivalent to the Camp-Stein equation for the mixing gradient. The gradients developed by Amirtharajah [7] and Kawamura [27,28] result in similar optimal porosity of fluidized filter bed [30], while using the approach of Sakkas and Lekkas [29] gives a slightly higher value [31]. However, many years of practice result in applying lower expansion of fluidized filter bed. Hence, in the algorithm adopted in the calculations, it was assumed that the expansion p defined by Equation (1) should be 25%. where: p-expansion of the bed during fluidization, V-volume of the filter bed during fluidization, V 0 -bed volume before fluidization.
For an expansion value p of 25%, a bed porosity ε was calculated from Equation (1) and the mass balance Equation (2): where: V-bed volume during fluidization, ε o -known bed porosity before backwash, ε-bed porosity during fluidization. It is known that the expansion of the filter bed depends on the temperature of the backwash water, especially for granular activated carbon beds and for anthracite and sand beds of fine grain size [32,33]. The finer the grain size and the lower the grain density, the lower the free sedimentation velocity v s of the filter grains and the more laminar character of their fall, and thus the greater the influence of the dynamic viscosity coefficient µ on this sedimentation. The µ value strongly depends on the water temperature T. For eighty years, the relationship between the free sedimentation velocity v s and the intensity of filter bed washing v and the porosity ε (v) has been traditionally described by the Richardson-Zaki Equation (3) [34]: where: v-apparent velocity of the fluid in the fluidization zone, v s -free sedimentation velocity of individual grains of the bed, n-power exponent.
In recent years, there have been attempts to describe the expansion of a fluidized filter bed [35] which are not based on Equation (3). However, the determination of the free sedimentation velocity of grains depending on the nature of the movement can be determined from the equations summarized in Table 1 without the knowledge of parameters that are difficult to determine, which makes Equation (3) of great practical importance. From Equation (3), it is possible to calculate the desired values of backwash intensity for the already known porosity ε of the bed during fluidization, providing that the values of the free sedimentation velocity v s and the exponent "n" are determined. Equation (3) is physically meaningful only for backwash intensities exceeding the minimum fluidization velocity and not exceeding the value of v s . Both the free sedimentation velocity v s and the exponent "n" in Equation (3) depend in the general case on the Reynolds number Re s . This relationship is more pronounced for fine grains and for grains of low-density material such as granular activated carbon. Free sedimentation velocities v s should be determined from the finest grain fractions up to d 90 . It is commonly assumed that preferably 90% and no less than 80% of grains by weight should be fluidized during backwashing by water alone. The expansion is calculated for each fraction separately and then a weighted average is determined taking into account the mass contribution of the individual fractions. This way of calculating the total expansion is a commonly used simplification, although the suspended grains mix with each other during backwash [36]. In monograph [37], several equations for the value of the exponent "n" in the Richardson-Zaki Equation (3) were compiled, in the general case depending on the material and grain size and the value of the Reynolds number Re s for free sedimentation of grains. The free sedimentation velocity v s of single grains depends on the fluid flow regime. It was calculated from the equations summarized in Table 1 [38]. In Table 1, the sedimentation velocities v s were calculated using the equations valid for the ranges of the Reynolds and the Archimedes numbers, defined by Equations (4) and (5). where: Re s -Reynolds number for free sedimentation of grains, Ar-Archimedes number, d-equivalent diameter of the grain, ρ g -the density of the grain, ρ w -density of water, ν-kinematic viscosity, g-acceleration of gravity.
For the backwash intensities v determined from Equation (3), the amount of backwashing water was calculated according to Equation (6) W i = t p ·v i ·A (6) where: W i -volume of water used for a single backwash of one filter in the 'i-th' month of the year, v i -intensity of bed washing in the "i-th" month of the year (depending on temperature), A-the area of a single filter bed, t p -backwash time with water only, increased by 50% of the water and air backwash time when simultaneously closing the air inlet and opening the water inlet.
The amount of water used to backwash the entire filter station throughout the year can be determined from Equation (7). where: W-volume of water used for backwash of all filters in the station during the year, FN-number of filters in the station, f -average backwashing frequency of one filter per day.
In compiling the energy gains, it was also necessary to first count the energy consumed in fluidizing the filter bed during the backwash with water alone. The hydraulic gradient H/L was determined from Ergun's Equation (8), in which the first expression describes the laminar linear component of the flow resistance through the porous bed, as in the Kozeny-Carman equation, while the second term describes the non-linear hydraulic gradient component as a function of backwash intensity v. where: h-pressure height loss, L-thickness of the filter bed, µ-dynamic viscosity, ψ-spherical coefficient. During backwashing of filter media by water alone, the resistance to the flow is almost linear in a function of backwash intensity, until the minimum fluidization velocity is reached. Thereafter, a disturbance occurs before the minimum fluidization velocity is reached and then the flow through the fluidized bed produces an almost constant hydraulic resistance. Increasing the backwash intensity increases the bed expansion, but hardly changes the magnitude of flow resistance, which can be explained by the balance of forces acting on the individual grains. For any backwash intensity, the gravity force minus the Archimedes buoyancy force is equal to the hydrodynamic force. The latter force does not change with backwash intensity, which is possible due to the increase in bed expansion. The amount of pressure loss h min caused by the flow through the fluidized bed can be calculated from Equation (9). It describes the equilibrium of forces at the moment of lifting all grains of the bed for a backwash intensity equal to the minimum fluidization velocity.
Excluding the small variations in pump-motor efficiency resulting from speed adjustments, the power consumed by the pump just to maintain fluidization during backwash with water alone is therefore proportional to the intensity of the backwash. The energy dissipation to overcome the flow through the pipelines and underdrain must also be added, which requires individual calculations for each filter station. Laminar flow in this case is beyond the scope of interest. The flow in the transient zone is more similar to turbulent flows for larger economic diameters of pipes, especially considering that these pipes are operating under higher economic flow velocities.
The annual amount of energy E used to backwash all filters in the station can be determined from Equation (10): where: E-energy, N p -power at pump shaft during filter flushing. where: Q p -pump capacity during backwash, H p -pump head during backwash, η p -pump efficiency. The pump capacity Q p during the backwashing by water alone is calculated from Equation (12) in which W i is defined by Equation (6).
The pump head during backwash can be calculated using Equation (13) where:

Replacements of the Filter Media during Plant Operation
Over a longer period of time, abrasion of the grains occurs during air and water backwashing as well as rounding of the grain edges, and thus a change in the sphericity ψ. The sphericity ψ is defined as the ratio of the sphere surface area having the same volume as the grain to the real grain surface area, and ranges from 0.98 to 0.78 for natural mineral grains and may be even lower than 0.70 for artificially crushed grains [39]. For sand grains, the sphericity coefficient depends on the geological origin of the sand. Fluvial sands have more rounded shapes than those of glacial origin. According to a study in South Africa [40], only 0.3% of the weight of the activated carbon filter media was lost through abrasion, while 19% was lost during transport and almost 20% during thermal regeneration. Therefore, abrasion should not change the grain size classification of minerals from which filter beds are made. However, the formation of a biofilm by bacteria and the sticking together of grains, followed by the carrying away of low-density agglomerates formed in this way, is a common occurrence. The same applies to the loss of grains as a result of selecting too high backwash intensity for winter conditions. Moreover, with time, the grains of the filtration bed become partially covered with non-removable sediment, which changes the density of grains to such an extent that it may significantly influence the amount of bed expansion during washing [31]. Therefore, it is advisable to check periodically the filter grain size, or to measure on a technical scale the expansion for a known backwash water temperature.

Economy of the Filter Backwash System
When comparing the backwash system of water rapid filters with and without temperature-dependent intensity control, a choice must first be made as to how to determine the backwash intensity. If you choose the intensity of backwash for summer conditions, you take the risk of losing some of the grains in winter. If, on the other hand, the backwash intensity is set for winter conditions, there is no risk of losing some of the grains in summer, but at high backwash water temperatures, the expansion is too small, so the backwash is carried out incorrectly, which is likely to affect the quality of the filtrate. Paradoxically, setting backwash for winter conditions can also indirectly lead to grain loss through insufficient bed expansion, which leads to an excessive biofilm sticking to the grains. There are many possibilities to set backwash parameters for conditions intermediate between winter and summer. For system comparison, it is assumed here that backwash is actually carried out year-round at the expansion set for summer conditions. A case study was then carried out in which the backwash intensity was adjusted to the average monthly temperatures.

Rapid Filter Station and Granular Activated Carbon Filters at the Rudawa River Water Treatment Plant
The analyses were carried out at one of four water treatment plants in the city of Krakow in southern Poland. The plant draws water from the Rudawa River. The maximum capacity of the station is 55,000 m 3 /d, and the current production varies between 22,000 m 3 /d and 28,000 m 3 /d. The raw water in the first stage is preliminarily oxidized with sodium permanganate, then it is coagulated in the rapid mixing chamber and flocculated in the slow mixing chambers. After flocculation, the water flows into horizontal settling tanks where floc sedimentation takes place. The water thus prepared goes first to rapid filters and then to filters with granular activated carbon. The filter station is currently equipped with 10 filter chambers filled with quartz sand and 8 chambers filled with granulated activated carbon. The height of the sand beds is 100 cm and the height of the granular activated carbon beds is 220 cm. The surface area of each of the sand filters as well as the filters filled with activated carbon is 30 m 2 . The capacity of one sand filter varies between 150 and 200 m 3 /h and the contact time in the GAC filter is 18 min. The sand filters are washed at an average frequency of once a day except on days with heavy rainfall. Backwash is carried out in air, water-air and water systems in sequence. First, the bed is agitated with air and then, while simultaneously closing the air inflow and opening the backwash water inflow, the backwash is conducted in the water-air system. The process ends with backwash with water only. GAC filters are only backwashed once a week not because of clogged pores, but because of the bacteriology of the treated water.
A specially designed optical device the Raven Environmental Products, Saint Louis, Missouri, USA SID-10200 Sludge Interface Detector was used to measure the expansion of the filter bed in the Rudawa River water treatment plant. The device allows the height of the bed to be determined during fluidization from a higher level.
Laboratory tests were carried out on a column ( Figure 1) with a square cross-section of 20.5 cm × 20.5 cm. It was assumed that the column cross-section is large enough that the wall effect accompanying fluidization should not significantly affect the experimental results. In addition, a uniform distribution of backwash flow during fluidization was maintained across the cross-section of the column by using perforated plates in the subfilter chamber and a stainless-steel wedge wire grating with triangular cross-section and narrow slots of 0.25 mm, spaced parallel every two millimeters.
The column was filled with sand taken from one of the rapid filters at the Rudawa River water treatment plant under analysis. Two sand samples were taken independently from the filter, dried in the laboratory, sieved, and weighed, and then two independent sieving curves were produced based on the measurements taken. As the curves are very close to each other, the sand samples taken were considered to be representative for this filter. The sand samples were mixed together and the sieving curve calculated after mixing at known weight ratios. The representative sieving curve is presented in Figure 2.  The height of the column was 2.0 m. The height of the unfluidized sand bed was equal to 42 cm. Backwashing of the sand bed took place from a much higher placed small tank with an overflow into which water was pumped from a large lower tank where tap water was heated to the required temperature. Water flowing through the overflow in the upper tank was returned to the lower tank. The flow was measured using two flow meters in series.

The Exponential "n" in Richardson-Zaki Equation
The analysis was conducted for the filter station at the Rudawa River water treatment plant. The tests were carried out directly at the station and at the stand described in the previous paragraph and built in the laboratory of the Cracow University of Technology.
The average monthly raw water temperatures at the Rudawa water treatment plant are shown in Figure 3. The plant is supplied with river water through retention tanks with a retention time of 7 days, which results in an equalization of the daily raw water temperature entering the station. Measurements were carried out to determine the grain density of the filter bed in the Rudawa water treatment plant (WTP). The results show that the density of sand grains in the studied bed together with a layer of sediment permanently bound on their surface was 2087 kg/m 3 and was significantly lower than the density of different quartz varieties, which usually ranges between 2648 and 2651 kg/m 3 . Such low density values of the grains meant that they were covered by a large amount of permanently bound sediment. The biofilm consists of approximately 90% extracellular polymeric substances and 10% bacteria. Its density is close to that of water, such this adhesion reduces the density of the sand covered by biofilm. Most of the equations proposed in the literature for calculating the "n" value of the power exponent appearing in the Richardson-Zaki Equation (3) are applied to grains of silica sand or other materials, but do not cover a significant amount of permanently bound sediment. It was therefore necessary to develop an equation that would describe the expansion of grains with such low density and very strongly curved edges. To calculate the exponent "n" in Equation (3), Function (14) was developed by fitting the coefficients in this equation using the method of least squares applied to the results of expansion measurements. Figure 4 shows measured values and the results of the bed expansion calculations. Calculations were performed using the equations in Table 1 and Equations (1)-(5) and Equation (14) fitted to the experimental measurements for temperatures of 13, 19, and 25 • C. The coefficients of determination R 2 for the filter bed expansion determined, considering Equation (14)   To smooth the results of the calculations shown in Figure 4, trend lines were drawn for them together with the equations summarized in Table 2. High values of R 2 with a relatively high number of measurements, and therefore a high number of degrees of freedom, indicate a good fit of the trend lines. Based on the proposed algorithm, considering the equations in Table 1, Equations (1)-(5) and the newly developed Equation (14) fitted to the experimentally measured points in Figure 4, expansion curves for temperatures of 3, 5, 10, 13, 19, 25, and 28 • C were determined and are shown in Figure 5.  (3) with the exponent of the power "n" described by Equation (14) and the free sedimentation velocity according to the nature of motion calculated according to the equations presented in Table 1.

Viscosity Coefficient
The dynamic viscosity of water varies with temperature according to Equation (18) [41], in which µ(T) and µ(T 0 ) are the dynamic viscosities corresponding to the instantaneous temperature T and the reference temperature T 0 , respectively. where: T-temperature in • C, T 0 -reference temperature in • C.
For backwash below the minimum fluidization velocity, the hydraulic resistance of the flow depends almost linearly on the backwash intensity. However, once the sand is fluidized, the flow is no longer laminar linear, as the current lines are not fixed in time and the local flow velocities at the grain-water interface are different from zero. For granular activated carbons in work [42], expansion was described for different temperatures of washing water using the results of measurements carried out for a reference temperature. It was assumed that the expansion is identical for the same product of the dynamic viscosity raised to the power of "z" and washing intensity v, thus for the same product v · µ z . Equation (18) shows that µ(T) = µ(T 0 )(1 + 0.0337 · T 0 + 0.00022 · T 0 2 )/(1 + 0.0337 · T + 0.00022 · T 2 ). For the reference temperature T 0 = 13 • C, the value µ(t 0 ) = 1.2028 mPa· s. Taking the expansion curve p(T 0 = 13 • C) as a starting point, the trend line v = v(p) for backwash intensity was first calculated. Measurements and the calculations are technically irrelevant for p expansions below 6%. To obtain a better fit of the trend line, these data were not considered. After omitting the measurements for p < 6%, a very good fit was obtained with a high R 2 value for the approximation v(p) with a power function. This function was used to calculate v(p) for T 0 = 13 • C, for p values the same as those recorded in measurements taken for T = 19 • C and 25 • C. The equation v(T) = v(T 0 ) · [µ(T 0 )/µ(T)] z was then used for the same expansions to predict calculated v(p) values for 19 and 25 • C. The results of calculations and measurements are compared in Figure 6 for z = 1. The calculated v(p) values are close to the measurements even for z = 1. However, it was decided that measurements of p made for only three close temperatures (one of which is the reference temperature) were not sufficient to evaluate the applicability of the equation v(T) = v(T 0 ) · [µ(T 0 )/µ(T)] z . Another widely accepted prediction method p(v,T) was used to calculate the potential for water and energy savings during the backwashing of filters at the Rudawa water station (south of Poland).

Calculations for an Expansion of 25%
Since the 25% expansion lies in the middle of the recommended expansion range of the filter bed, calculations of the expected water and energy savings were carried out for this value of expansion. For this purpose, from the equations in Table 2 and from Equation (14), the values of the backwash intensity v(T) at which 25% expansion will occur were calculated for temperatures of 3, 5, 10, and 28 • C. The results are quite accurately approximated by the linear Equation (19), which is presented in Figure 7. In Equation (19), the temperature T is expressed in • C and backwash intensity v in m/h. v = 0.5928·T + 32.095 (19) Figure 7. Intensity of rapid filter backwashing for 25% expansion, as a function of temperature.

Analysis of Annual Water Savings in a Rapid Filter Plant
Based on the developed mathematical algorithm, the annual quantity of water used for the backwashing of rapid filters was determined for the analyzed filter station in the water treatment plant, once with the assumption that the intensity of backwash was maintained at the same level throughout the year, and the second time that it was adjusted to the average monthly temperature. In the case of the adjustment of the backwash intensity, it was assumed that it was adjusted to the average monthly temperature so as to maintain a constant expansion of 25% throughout the year. In the case of the unadjusted backwash, it was assumed that the backwash intensity was the same throughout the year and that it guaranteed an expansion of 25% in the warmest month. Table 3 shows the results of the water consumption calculations for the backwashing with intensity adjustment. Table 4 shows the results of the calculations of water consumption for backwashing without monthly adjustment of intensity. Finally, Table 5 compares the monthly and annual results and the water savings resulting from adjusting the backwash intensity. These gains amount to about 854,100 m 3 .  From Tables 3-8, it can be concluded that the greatest water and energy savings occur in winter when water temperatures are at their lowest. Comparisons were made for the case of adjusting the backwash intensity once to the summer water temperature (July) and then to the average water temperature in each month. Hence, in July, the backwash intensity is identical in both cases, so there are no water or energy savings in this month.

Analysis of Annual Energy Savings in a Rapid Filter Plant
The energy savings resulting from adjusting the backwash intensity of the filter plant on a monthly basis were also analyzed. The differences in energy consumed were mainly due to the different amounts of water required for backwash, and to a lesser extent the different bed resistances during fluidization. Bed resistance to flow during fluidization was determined from Equation (9), and power and energy from Equations (10)- (11). The results of calculations for filters washed with the same intensity are summarized in Table 6, with the intensity adjusted monthly to temperature in Table 7, and a comparison of both calculations in Table 8.

Granulated Activated Carbon
In the water treatment plant for which an analysis of the benefits of adjusting the backwash intensity to the water temperature was carried out, all of the water after the sand filters is diverted, without ozone treatment, to eight filters with GAC which are in urgent need of replacement. A decision on the type of granulated activated carbon to be used has not yet been taken. Four types of new GAC produced by Chemviron Carbon were selected, for which equations previously developed by one of the co-authors of this article were used [42]. Equations (20)-(23) for backwash intensity as a function of water temperature and GAC bed expansion are presented in Table 9. Based on these equations and monthly average temperatures, the backwash intensity required to achieve bed expansion of 25% for the four selected GAC types during the year was calculated and is shown in Figure 8. Based on the backwash intensities thus determined, the volume of backwash water expected to be used each month of the year for all the GAC filters was calculated for the conditions at the Rudawa River water treatment plant and is shown in Figure 9. As in the case of rapid sand filters, the monthly and annual water consumption for backwashing GAC filters without and with backwash intensity control was calculated. Calculations of water savings achieved by using backwash intensity control of the filters were carried out for the four selected GAC types and are presented in Tables 10 and 11. Tables 10 and 11 compare the backwash water savings from adjusting the backwash intensity of the GAC filters first to the summer water temperature in July and then to the average water temperatures in each month. Hence, the backwash intensity in July was identical in both cases, so there were no water savings in that month. As expected, the highest savings were achieved in winter, when the required expansion could be achieved at a much lower intensity than the adjusted intensity for the summer season (July). Table 9. Approximation equations for backwash intensity as a function of water temperature t and bed expansion p taken from [42].

Discussion
It is well known that porous media bed expansion depends on water temperature, at least for fine grains and for GAC. However, most filter plants do not have a backwash intensity control system adjusting the backwash to the current temperature. There are several reasons for this situation. One of them is that it has not been reported how much water and energy can be saved. The study presented here gives an estimation of the saving opportunities. These possibilities are underestimated because they do not take into account the savings in the backwash water treatment or the pressure losses in the backwash water supply lines. Surprisingly significant changes over time can occur to filter media, including not only the changes in shape and size, but also the overall density of the grains, taking into account impurities permanently adhered to the grains surface. A new formula (14) for the exponential "n" of the Richardson-Zaki equation was developed here. This equation is valid for the sand from the filter plant studied here, which has been in operation for many years. It would be interesting to see if this formula is also correct for other sand media used for filtration purposes in surface water treatment plants. In the future, more possibilities for water and energy saving should be analyzed and the results obtained should be applied.
Analyses carried out in this article for a medium-sized water treatment plant located in the northern part of Europe have shown almost 10% water savings per year for backwashing of rapid filters and also GAC filters through the use of backwash intensity control. The amount of water saved in this way was 836,310 m 3 per year using rapid filters and, depending on the type of GAC, from 4570 to 7511 m 3 .
Although the analysis of energy savings due to backwash intensity regulation showed a gain of more than 10% in the year, relatively small benefits are expected, since a relatively low amount of energy is consumed for filter backwash.
In addition to adjusting the backwash intensity to the water temperature, other actions can be taken to save water and energy. When using anthracite-sand beds, and multi-layer beds in general, the aim is usually to lay down a coarse layer of material with a lower density on the bed surface. The pores in this layer are large, leading to a significantly longer filter cycles. A granular activated carbon (GAC) bed can be used simultaneously as a mechanical filter and sorbent, which also leads to a significant reduction in the amount of backwash water. A thin layer of GAC is also used as the top layer of multi-layer water rapid filters. GAC filters upstream of sand filters are increasingly used in Japan. The reason for this is the large number of microorganisms inhabiting the carbon dust from filter beds filled with granular active carbon. These microorganisms are extremely resistant to disinfection because they are partially hidden in the largest pores of the dust. Moreover, the most popular disinfectants based on chlorine turn out to be ineffective in these cases. Although water savings for backwashing the filtration-adsorption system in the order GAC-sand filters are obviously not the main reason for using sand filters after GAC, these savings are evident. In the traditional sand filter-ozonation-GAC systems, granular activated carbon is usually washed weekly due to the bacteriology of the treated water, and the fine-grained sand filter bed becomes clogged quickly, filtering the water suspension of high concentration flowing from the settling tanks. On the other hand, in the GAC-sand system, granulated activated carbon does not clog so quickly due to the large size of the pores and the sand filter can be backwashed rarely due to the low concentration of the inflowing suspension. The choice of the filter operation system also has a great impact on saving water and energy. Many years of research, summarized in [43], have shown that it is difficult to decide on the basis of numerous experiments which filtrate is better, that coming from the operation of identical filters with a constant or with a declining variable filtration rate [44]. The performed numerical calculations [43,45,46] based on a simple capillary model (unit bed element) concluded with the suggestion that the filtrate quality from an identical rapid filter operated in a variable declining rate (VDR) system was slightly better, when the filter cycles ended after obtaining the same head pressure losses of flow through the bed. However, according to the calculations, a slightly better filtrate was obtained in the system with a constant filtration rate, when the filter cycles were of the same length [43]. First, theoretically [47], and then experimentally [48], it has been shown that in the case of filter plants controlled by variable declining rate operation systems, in comparison with constant rate systems, longer filter cycles, and thus saving water for backwashing, are achieved in the VDRF operation systems if the ratio of the hydraulic load of the most recently backwashed filter to the average hydraulic load of the all filters in the plant is higher than 1.2. This VDR system can be used for the time-varying turbidity of raw water [49,50], but it requires adjustment when changing the water temperature and it is troublesome to use this system in large stations due to small fluctuations in the water table above the filters, which makes it difficult to determine the appropriate time of subsequent backwashing in stations [51]. Thus, financial benefits do not always determine the adopted technical solutions. Other energy and water saving possibilities such as appropriate internal pipe lining [52,53] or leakage protection [54] inside the filter station are of marginal importance.

•
An algorithm for selecting backwash intensity for rapid filters and GAC filters in a water treatment plant was successfully developed and applied to an example. • An algorithm has been developed for determining water and energy savings due to the use of backwash intensity control for rapid and GAC filters in the water treatment plant in individual months of the year. • A new formula was developed to determine the exponent "n" in the Richardson-Zaki equation for selecting appropriate backwash intensity. • It was found that the temperature dependence of backwash intensity at rational 25% expansion is almost linear. • An analysis of a municipal water treatment plant located in northern Europe with a current average capacity of about 25,000 m 3 /d showed a potential annual saving of about 71,045 m 3 of water used in backwashing the sand filter beds, representing about 8.5% of the total volume of water used for this purpose and about 0.78% of the total production of the water treatment plant. • Potential water saving during backwashing of GAC filters with monthly intensity adjustment ranged from 4570 to 7511 m 3 , depending on bed type. This represents from 7.57 to 13.7 per cent of the water used for this purpose.

•
Although the analysis of energy saving due to backwash intensity regulation showed a gain of 13.6% of total annual energy used for this purpose, relatively small energy benefits are expected as relatively little energy is used to backwash rapid and GAC filters.