Superconvergence analysis of finite element method for Poisson–Nernst–Planck equations