Effects of Cemented Porous Media on Temporal Mixing Behavior of Conservative Solute Transport

The cementation of porous media leads to the variation of the pore space and heterogeneity of the porous media. In this study, four porous media (PM1, PM2, PM3, and PM4) with the different radii of solid grains were generated to represent the different cementation degrees of the porous media. The direct simulations of flow and conservative solute transport in PM1–4 were conducted to investigate the influence of the cemented porous media and Peclet number (Pe) on the temporal mixing behavior. Two metrics, scalar dissipation rates (SDR) and dilution index, were employed to quantify the temporal mixing behavior. It was found that the spatial velocity variability of the flow field was enhanced as cementation degree increased. The results of the coefficient of velocity variation ( C V U ) increased from 0.943 to 2.319 for PM1–4. A network consisted of several preferential flow paths was observed in PM1–4. The preferential flow enhanced the mixing of the conservative solute but had a negative influence on the mixing of the solute plume when the cemented solid grains formed several groups, and there were some stagnant regions where the flow was almost immobile. As the Pe increased, for PM1–3, the exponent of the best-fitting power law of the global SDR decreased. At the case of Pe = 400, the slope of the global SDR reduced to around −1.9. In PM4 where the preferential flow was enhanced by the cemented solid grains, the slope of the global SDR increased as the Pe increased. The global SDR results indicated that the temporal mixing behavior followed a Fickian scaling ( S D R ∝ p v − 1.5 ) in the early stage (Pv < 0.05), while the mixing behavior turned to be non-Fickian in the late stage. The transition time from the Fickian scaling to the non-Fickian scaling was found to be sensitive to the cementation degree of the porous media.

Multiscale Roughness Influence on Conservative Solute Transport in Self-affine Fractures

In this article, the influence of multiscale roughness on transport of a conservative solute through a self-affine fracture was investigated. The fracture roughness was decomposed into two different scales (i.e., a small-scale stationary secondary roughness superimposed on a large-scale non-stationary primary roughness) by a wavelet analysis technique. The fluid flow in the single fracture was characterized by Forchheimer's law and exhibited nonlinear flow features such as eddies and tortuous streamlines. The results indicated that the small-scale secondary roughness was primarily responsible for the nonlinear flow features. Numerical simulations of asymptotic conservative solute transport showed non-Fickian transport characteristics (i.e., early arrivals and long tails) in breakthrough curves (BTCs) and in residence time distributions (RTDs) with and without consideration of the secondary roughness. Analysis of multiscale BTCs and RTDs showed that the small-scale secondary roughness played a significant role in enhancing the non-Fickian transport characteristics. Removing small-scale secondary roughness delayed the arrival time and shortened the tail. The peak concentrations in BTCs decreased as the secondary roughness was removed, implying that the secondary roughness could also enhance the solute dilution. Fitting the one-dimensional (1D) Fickian advection-dispersion equation (ADE) to the numerical BTCs resulted in considerable errors that decreased with the small-scale secondary roughness being removed. The 1D mobile-immobile model (MIM) provided a better fit to the numerical BTCs and inclusion of the small-scale secondary roughness in numerical simulations resulted in a decreasing MIM mobile domain fraction and an increasing mass exchange rate between immobile and mobile domains.