Transport of water and ions in partially water-saturated porous media. Part 1. Constitutive equations

I developed a model of cross-coupled flow in partially saturated porous media based on electrokinetic coupling including the effect of ion filtration (normal and reverse osmosis) and the multi-component nature of the pore water (wetting) phase. The model also handles diffusion and membrane polarization but is valid only for saturations above the irreducible water saturation. I start with the local Nernst-Planck and Stokes equations and I use a volume-averaging procedure to obtain the generalized Ohm, Fick, and Darcy equations with cross-coupling terms at the scale of a representative elementary volume of the porous rock. These coupling terms obey Onsager's reciprocity, which is a required condition, at the macroscale, to keep the total dissipation function of the system positive. Rather than writing the electrokinetic terms in terms of zeta potential (the double layer electrical potential on the slipping plane located in the pore water), I developed the model in terms of an effective charge density dragged by the flow of the pore water. This effective charge density is found to be strongly controlled by the permeability and the water saturation. I also developed an electrical conductivity equation including the effect of saturation on both bulk and surface conductivities, the surface conductivity being associated with electromigration in the electrical diffuse layer coating the grains. This surface conductivity depends on the CEC of the porous material.