Parameter-robust discretization and preconditioning of multiple-network poroelasticity equations