Influence of Pore Size Distribution on the Electrokinetic Coupling Coefficient in Two-Phase Flow Conditions

• Main Authors: Jan Vinogradov, Rhiannon Hill, Damien Jougnot
• Format: Article
• Language: English
• Published: MDPI AG 2021-08-01
• Series: Water
• Subjects: electrokinetic coupling coefficient; zeta potential; sandstones; partial water saturation; CO₂ geo-sequestration; bundle of capillary tubes model
• Online Access: https://www.mdpi.com/2073-4441/13/17/2316

Streaming potential is a promising method for a variety of hydrogeophysical applications, including the characterisation of the critical zone, contaminant transport or saline intrusion. A simple bundle of capillary tubes model that accounts for realistic pore and pore throat size distribution of porous rocks is presented in this paper to simulate the electrokinetic coupling coefficient and compared with previously published models. In contrast to previous studies, the non-monotonic pore size distribution function used in our model relies on experimental data for Berea sandstone samples. In our approach, we combined this explicit capillary size distribution with the alternating radius of each capillary tube to mimic pores and pore throats of real rocks. The simulation results obtained with our model predicts water saturation dependence of the relative electrokinetic coupling coefficient more accurately compared with previous studies. Compared with previous studies, our simulation results demonstrate that the relative coupling coefficient remains stable at higher water saturations but vanishes to zero more rapidly as water saturation approaches the irreducible value. This prediction is consistent with the published experimental data. Moreover, our model was more accurate compared with previously published studies in computing the true irreducible water saturation relative to the value reported in an experimental study on a Berea sandstone sample saturated with tap water and liquid CO₂. Further modifications, including explicit modelling of the capillary trapping of the non-wetting phase, are required to improve the accuracy of the model.