## Abstract

Saturated sand box experiments were conducted to explore the effect of various physical factors on the transport and deposition of suspended particles in porous media. Red quartz powder and natural quartz sand were employed in the study and acted as suspended particles and porous media, respectively. Particles were injected into the sand box in two modes, i.e., pulse injection and continuous injection. Tests were performed at various particle concentrations, flow velocities, deposition rate coefficient and longitudinal dispersion coefficient by both injection modes. The breakthrough curves were described with the analytical solution of a convection–dispersion equation, in which first-order deposition kinetics were taken into account. Different behavior of suspended-particle transport and deposition in porous media was observed under different injection modes and experimental conditions. The results show that effluent concentration was approximately linear with the initial particle concentration. The deposition rate coefficient depends strongly on particle size and flow velocity, and the transport and deposition process was very sensitive to it. Furthermore, the longitudinal dispersion coefficient increases with increasing flow rate, and particles are easier to transport through pores as the longitudinal dispersion coefficient increases. This study shows the importance of particle concentration, flow velocity, deposition rate coefficient and longitudinal dispersion coefficient in the transport and deposition process of suspended particles.

## INTRODUCTION

Research on suspended particle transport and deposition in porous media is required to assess groundwater recharge efficiency and to protect drinking water supplies from contaminants (Sen 2011; Ward & Dillon 2012), to explore oil in low-permeability reservoirs and to microbially enhance oil recovery (Patel *et al.* 2015), to predict the stability of filter layers in the dam and to quantify colloid-facilitated transport of various contaminants (Rinck-Pfeiffer *et al.* 2000). Previous studies of particle transport, deposition and release mechanisms have mainly focused on colloids and microorganisms with sizes smaller than 1 μm (Bradford & Torkzaban 2013), and suspended particles with sizes larger than 5 μm have only attracted significant attention recently (Wang *et al.* 2000; Kim & Lawler 2012; Liu *et al.* 2014, 2015).

For particles with a size larger than 10 μm, hydrodynamic force, gravity and inertial effects are the dominant forces that control its fate (Kasel *et al.* 2013). For particles with smaller size (between 0.1 and 10 μm), factors such as electrostatic forces and Brownian motion can also contribute to particle deposition (Lanphere *et al.* 2014; Pugliese *et al.* 2015). The boundary between colloid and suspended particle is quite blurry, but it obeys the rule that as particle size decreases, particle interactions become increasingly significant relative to external forces. Diverse correlations have been developed between experimental parameters (i.e. attachment coefficients, straining rate coefficients and detachment coefficient) and variables of physical, chemical and biological particles (Alvarez *et al.* 2007; Bedrikovetsky *et al.* 2011; Yuan *et al.* 2013).

Experimental and theoretical work has demonstrated that physical factors, such as particle concentration, flow velocity, deposition rate and longitudinal dispersion coefficient, may affect the transport and fate of suspended particles in saturated porous media (Ahfir *et al.* 2007, 2009; Bedrikovetsky *et al.* 2011). Published literature on this field was mainly focused on colloid particles with a size smaller than 10 μm; the transport and deposition mechanisms for bigger particles have not been systematically studied. In other words, the effects of these physical factors on the transport and deposition of suspended particles in natural porous media have not been studied thoroughly, especially in experiments. Therefore, a systematic study of the interactions between these factors is essential for understanding the process of suspended-particle transport and deposition in a saturated porous medium.

The objective of this work is to quantify the effects of particle concentration, flow velocity, deposition rate coefficient and longitudinal dispersion coefficient on the transport and deposition of suspended particles in a saturated porous medium. A series of experiments was conducted with a self-developed sand layer transportation–deposition testing system, and two injection modes, i.e., pulse injection and continuous injection, were adopted. The theoretical result of the convection–dispersion model fits the experimental data fairly well. In addition, the transport and deposition mechanisms of particles under various factors were investigated.

## THEORY

### Mathematical model

*et al.*2007): where

*C*is the particle concentration (M/L

^{3}),

*t*is the time (T),

*x*is the travel distance of particles (L),

*D*

_{L}is the longitudinal dispersion coefficient (L

^{2}/T),

*u*is the average flow rate of water in the pores (L/T),

*K*

_{dep}is the deposition rate coefficient (T

^{−1}).

It should be noted that the process of particle detachment is neglected in this equation.

Two modes of injection were studied in this work, i.e. pulse injection and continuous injection. The former way can describe instantaneous pollution in the stratum and the latter one can represent typical dispersal of pollution, as well as groundwater recharge with constant particle concentration.

*et al.*2000)

*et al.*2013) with where erfc(

*x*) = is the complementary error function.

*S*) and particle concentration flow with water (

*C*) can be expressed as where is the porous medium bulk density, is porosity.

*v*is the interstitial particle velocity.

Equation (8) is commonly referred to as the classical filtration model and has been widely used in describing the transportation of particles, colloids and microorganisms in saturated porous media (Lee *et al.* 2017).

### Sensitivity analysis of the model parameters

The transportation and deposition process of particles through porous media is affected by many parameters, thus it is of great importance to conduct a sensitivity analysis of the model parameters (Ouisse *et al.* 2012). In our study, parameters () under pulse injection and parameters () under continuous injection in the model were subjected to sensitivity analysis. The sensitivity analyses were conducted by fixing the remaining parameters and investigating the change trend of particle concentration in the pores () when a parameter changes. This procedure was conducted for and , separately, which may affect the process of particle transport and deposition. Referring to previous studies (Tufenkji & Elimelech 2004; Tufenkji 2007), the initial values of and were taken as 4.1 × 10^{−5} m^{2}/s, 2 × 10^{−2} 1/s, 3 × 10^{−3} m/s, 5 g and 1.8 × 10^{−}^{5} m^{3}/s, respectively. On this basis, particle concentration in the pores () was calculated at *x* = 0.6 m and *t* = 200 s for continuous injection by making changes of ±10%, ±20%, ±30%, ±40% and ±50% on initial values of the four parameters. The peak concentration () was calculated at *x* = 0.6 m for pulse injection in the same way.

The relationship between peak concentration () and the change of each parameter under pulse injection is shown in Figure 1. Under pulse injection mode, it can be seen that the change of peak concentration increases with the change of injected particle mass (*m*) and flow rate (*u*). The change of peak concentration ranges from −50% to 50% for the change of injected particle mass and −95.6% to 394.2% for the change of flow rate. Meanwhile, change of peak concentration decreases with the change of deposition rate (). The change of peak concentration ranges from 429.1% to −79.1% for the change of deposition rate. Unlike these parameters, the change of peak concentration almost remains constant with the change of longitudinal dispersion coefficient ().

The relationship between particle concentration in the pores () and the change of each parameter under continuous injection is shown in Figure 2. Under continuous injection, it can be seen that the change of particle concentration increases with the change of injected particle concentration (), longitudinal dispersion coefficient () and flow rate (*u*). The change of particle concentration varies from −50% to 50% for the change of injected particle mass, −20.6% to 19.0% for the change of longitudinal dispersion coefficient and −98.4% to 291.2% for the change of flow rate. Meanwhile, the change of particle concentration decreases with the change of deposition rate (). The change of particle concentration varies from 394.7% to −78.8% for the change of deposition rate. Similar to the pulse injection, the particle concentration decreases with the decrease of the flow rate until 0, beyond which it becomes almost constant, while the change of particle concentration increases rapidly as the flow rate increases.

These above behaviors indicate that the deposition rate is negatively correlated with particle concentration. However, the injected particle mass, the injected particle concentration and the flow rate are positively correlated with the particle concentration. The analysis also suggests that the particle concentration is highly sensitive to the flow rate and the deposition rate, while the longitudinal dispersion coefficient effect is negligible.

## MATERIALS AND METHODS

### Set up

In our study, a self-developed sand layer transportation–deposition testing system was used. A schematic diagram of the experimental apparatus is illustrated in Figure 3. All experiments described in this paper were run on this set, which mainly consists of a dynamic system, observation system and data collection system. The studied porous media are placed in a two-dimensional Plexiglas box and the particles are distributed in the water contained in the water tank. The Plexiglas box is sealed with silicon and strong enough to withstand tremendous water pressure. Six connectors are placed on the sides of the box and can be connected with the sensor or pipe. Interior dimensions of the sand box are approximately 1.2 cm (width) by 50 cm (height) by 60 cm (length parallel to flow). The water in the set was driven by a peristaltic pump and could maintain a constant flow rate. With excellent transparency, transportation and deposition phenomena can be viewed in real time in the Plexiglas box and the sample can be obtained from the outflow pipe.

### Characteristics of the medium and injected particulate matter

The medium used in the study was natural quartz sand and the SiO_{2} was not less than 99.6%. The grain size distribution was 15% at 0.18–0.25 mm and 85% at 0.25–0.5 mm.

Red quartz powder was used in the experiments and three size distributions were chosen. Red quartz powder is physically and chemically stable and easy to observe. In addition, its composition is similar to the particles that are clogging in the engineering. The median diameter of particles was chosen as 12.9 μm when we studied the effect of particle concentration. The detailed parameters of sand and quartz powder are shown in Table 1.

Name . | Diameter (μm) . | Median diameter (μm) . | Density (g/cm^{3})
. | Porosity (%) . | Non-uniform coefficient . | Curvature coefficient . |
---|---|---|---|---|---|---|

Quartz sand | 250–500 | 360 | 2.65 | 32 | 1.3 | 1.1 |

Particle 1 | 5–6.5 | 5.5 | 2.65 | – | – | – |

Particle 2 | 10–15 | 12.9 | 2.65 | – | – | – |

Particle 3 | 20–25 | 22.7 | 2.65 | – | – | – |

Name . | Diameter (μm) . | Median diameter (μm) . | Density (g/cm^{3})
. | Porosity (%) . | Non-uniform coefficient . | Curvature coefficient . |
---|---|---|---|---|---|---|

Quartz sand | 250–500 | 360 | 2.65 | 32 | 1.3 | 1.1 |

Particle 1 | 5–6.5 | 5.5 | 2.65 | – | – | – |

Particle 2 | 10–15 | 12.9 | 2.65 | – | – | – |

Particle 3 | 20–25 | 22.7 | 2.65 | – | – | – |

### Operating procedure

Generally, the injection mode of particles during column-type tracing experiments involves continuous injection and pulse injection (Tufenkji & Elimelech 2004). In groundwater recharge systems, water source heat pump engineering and pollutant transportation, the concentration of particles is substantially stable. Meanwhile, pulse injection may occur in pollution emergencies, thus continuous injection and pulse injection were both adopted in the study.

To prepare the tests, the quartz sand was washed thoroughly and then packed into the box after being dried for 24 h at 105 °C. The water tank in the upper stream was filled with particle suspension and the inlet valve was then opened, with flow rate recorded by an electromagnetic flowmeter. The effluent concentrations of suspended particles were measured every 10 seconds in tests, and the transportation and deposition conditions were observed simultaneously.

## RESULTS AND DISCUSSION

### Initial particle concentration

Figure 4 presents the experimental breakthrough curves and corresponding simulated curves for different particle concentrations under two typical ways of injection.

The experimental results show that the sand box was clean in the beginning and the particle concentration started to increase and then reached the peak as particles were transported with the water. The changing trend of particle concentration was different after that, i.e., particle concentration continued to decline under pulse injection and remained stable under continuous injection. All results above-mentioned suggest that particle concentration can influence the transportation and deposition process, and the deposit formations may also vary as the concentration changes.

### Deposition rate coefficient

The deposition rate coefficient of suspended particles is determined by fitting the experimental breakthrough curves with the theoretical solution of the convection–dispersion equation.

*et al.*2006), the deposition rate coefficient can be expressed as a function with the single-collector contact efficiency (

*η*) and the sticking efficiency (

*α*), as follows: where is porosity,

*u*is the average flow rate of water in the pores, is the average grain diameter.

The deposition of particles in porous media results from multi-factors which are dependent upon physical and geometrical parameters. In this study, the chemical mechanisms were neglected from the assessment of the deposition process. Thus particle straining and sedimentation are the principal mechanisms that lead to the deposition of particles, which are closely related to the density and size of suspended particles.

*et al.*2000; Ahfir

*et al.*2009; Kim & Lawler 2012). The relationship between flow rate and deposition rate coefficient is often described with a power law: where and

*n*are constants that change with particle size.

The various tests carried out under different particle distributions and flow rates indicate that the deposition rate coefficient increases with flow velocity (see Figure 5). These results are in accordance with particle behaviors observed by other researchers (Kim & Lawler 2012).

Figure 5 presents the effects of the flow rate and particle size on deposition rate coefficient. We can observe that increases with the size of injected particles at a given velocity. Factors such as constrictions, crevices, and cavities in the sand box are reasons that gave rise to the deposition of particles in porous media. Kretzschmar *et al.* (1997) found a power value of 0.31 for carboxyl latex colloids and 0.18 for humic-coated hematite colloids. For silt particles with large size distribution, the value varies between 0.56 and 0.70 (Wang *et al.* 2000; Ahfir *et al.* 2007, 2009; Kim & Lawler 2012). In the present study, the power value varies from 0.239 to 0.306.

Figure 6 presents the experimental breakthrough curves and corresponding simulated curves for different deposition rate coefficients under two typical ways of injection. As changes from 28(1/h) to 122(1/h) under continuous injection, the final value of effluent particle concentration changes from 2 g/m^{3} to 181 g/m^{3}. This indicates that the particle concentration of the effluent is very sensitive to the deposition rate coefficient under both injection modes. For pulse injection, the times needed to reach the peaks of particle concentration of different deposition rate coefficients close to each other. This indicates that the deposition rate coefficient can influence the deposition mass of particles but will not affect their transport velocity markedly.

Different mechanisms of deposition may be the reason for the difference observed in tests. Large particles are dominated by gravity under low flow velocity, thus they can be captured in constricted and cavity sites. As the flow velocity increases, the dominant factor turns into hydrodynamic forces gradually. According to these results, one can deduce that the particle size and flow velocity play important roles in the particle deposition rate coefficient, as well as the process of particle transport and deposition.

### Flow rate and longitudinal dispersion coefficient

The longitudinal dispersion coefficient is a simplification of hydrodynamic dispersion, in which mechanical dispersion and molecular diffusion are both included. Mechanical dispersion and molecular diffusion are caused by variation in fluid velocity in the pores and the random motion of molecules, respectively.

*et al.*2009): where is the molecular diffusion coefficient, is the tortuosity of the porous media, is the longitudinal dispersivity and

*n*is an empirical coefficient.

where is the median particle size, is a coefficient (approximately 1) and is a constant.

For materials with uniformity coefficient smaller than 2, the value of the dispersivity is approximately equal to the median grain size (Bernabé *et al.* 2016).

Figure 7 presents the longitudinal dispersion coefficient of suspended particles as a function of the pore water velocity for different particle sizes (i.e., 5.5, 12.9 and 22.7 μm). The experimental results indicate that the longitudinal dispersion coefficient increases with increasing flow rate. The fitted results show that experimental data satisfy a power law quite well. The results show that the longitudinal dispersion coefficient has an exponent relation to the flow velocity, and the indexes are close to 1. Based on this and considering the difficulty in controlling in tests, we chose flow velocity to be a substitution parameter for .

The longitudinal dispersivity may be affected by the materials of the porous media. Ahfir *et al.* (2009) found the difference of longitudinal dispersivity between glass beads and gravel media through laboratory study, and the results are shown in Figure 8. Due to the broad pore-size distribution, the pore velocity functions are more important. Meanwhile, a uniform pore distribution of glass beads with a single size is also an important factor that led to the difference.

Figure 9 shows the evolution of particle dispersivity in the porous media according to the median particle size (5.5, 12.9 and 22.7 μm). The result indicates that particle dispersivity increases with the suspended particle size and its value varies from 0.98 to 3.03.

Figure 10 presents the observed and theoretical effluent concentration curves with different flow rates under pulse injection and continuous injection, respectively. It can be observed that the theoretical result fits the experimental data fairly well. As the flow velocity changes from 0.2 cm/s to 0.5 cm/s under continuous injection, the final value of the effluent particle concentration changes from 10 g/m^{3} to 70 g/m^{3}, while for pulse injection, the peak value of effluent particle concentration changes from 53 g/m^{3} to 287 g/m^{3}. This indicates that the particle concentration of the effluent is highly sensitive to the flow rate for both injection modes. We can also get the conclusion that as the flow rate decreases, more particles are deposited in the pores. For both injection modes, there is a big difference between the times required to reach the peak of particle concentration for different flow rates. This indicates that the flow rate can both influence the deposition mass of particles and their transport velocity markedly.

## CONCLUSION

This research was undertaken to study the characteristics of particle transport in saturated porous media. An experimental and analytical investigation is presented, which was aimed at studying the influence of various physical factors on the transport and deposition processes of suspended particles in saturated porous media. The response to pulse and continuous injection provides an effective way of analyzing the process of suspended-particle transport and deposition in porous media.

Results from these experiments show that there exist remarkable differences in transport and deposition characteristics between the two injection modes. For pulse injection, effluent particle concentration has a peak point after injection, and declines rapidly soon after that, while for continuous injection, effluent particle concentration can keep stable after a rising process.

The deposition rate coefficient and longitudinal dispersion coefficient were estimated from the experimental results with the analytical solution of the convection–dispersion model. Sensitivity analysis under the two injection modes indicates that peak concentration and particle concentration are highly sensitive to flow velocity and deposition rate coefficient.

Particle concentration may affect the transport and deposition process in a simple way, i.e., effluent particle concentration is almost in linear correlation with initial concentration. In contrast, the deposition rate coefficient may affect the transport and deposition process of particles more significantly. The relationship between flow velocity and deposition rate coefficient can be described by a power law, and particle size can also affect the deposition rate markedly. It is worth noting that particle peak concentration and effluent particle concentration change significantly even though there is little change in deposition rate coefficient. The deposition rate coefficient can influence the deposition mass of particles but will not affect their transport velocity markedly. The longitudinal dispersion coefficient has an exponent relation to the flow velocity, and the indexes are close to 1. As flow velocity increases, the particle peak concentration and effluent particle concentration also increase. Meanwhile the time needed to reach this concentration is reduced for both injection modes. Flow rate can influence both the deposition mass of particles and their transport velocity markedly.

From this experimental contribution, it arises that the particle concentration, flow velocity, deposition rate coefficient and longitudinal dispersion coefficient are important factors affecting the transport and deposition process of suspended particles in saturated porous media.

## ACKNOWLEDGEMENTS

This research was supported by the High-level Scientific Research Foundation for the Introduction of Talent, the National Natural Science Foundation of China (Grant No. 41272272) and the National Natural Science Foundation of China (Grant No. 51478368).