In this paper, we present a two-dimensional pore-scale numerical model to investigate the main mechanisms governing biofilm growth in porous media. The fluid flow and solute transport equations are coupled with a biofilm evolution model. Fluid flow is simulated with an immersed boundary-lattice Boltz-mann model while solute transport is described with a volume-of-fluid-type approach. A cellular automaton algorithm combined with immersed boundary methods was developed to describe the spreading and distribution of biomass. Bacterial attachment and detachment mechanisms are also taken into account. The capability of this model to describe correctly the couplings involved between fluid circulation, nutrient transport and bacterial growth is tested under different hydrostatic and hydrodynamic conditions (i) on a flat medium and (ii) for a complex porous medium. For the second case, different regimes of biofilm growth are identified and are found to be related to the dimensionless parameters of the model, Damköhler and Péclet numbers and dimensionless shear stress. Finally, the impact of biofilm growth on the macroscopic properties of the porous medium is investigated and we discuss the unicity of the relationships between hydraulic conductivity and biofilm volume fraction.