We demonstrate how a network model can predict porosity and permeability changes in a porous medium as a result of biofilm buildup in the pore spaces. A biofilm consists of bacteria and extracellular polymeric substances (EPS) bonded together and attached to a surface. In this case, the surface consists of the walls of the porous medium, which we model as a random network of pipes.
Our model contains five species. Four of these are bacteria and EPS in both fluid and adsorbed phases. The fifth species is nutrient, which we assume to reside in the fluid phase only. Bacteria and EPS transfer between the adsorbed and fluid phases through adsorption and erosion or sloughing. The adsorbed species influence the effective radii of the pipes in the network, which affect the porosity and permeability.
We develop a technique for integrating the coupled system of ordinary and partial differential equations that govern transport of these species in the network. We examine ensemble averages of simulations using different arrays of pipe radii having identical statistics. These averages show how different rate parameters in the biofilm transport processes affect the concentration and permeability profiles.