Multigrid method for nonlinear poroelasticity equations