High-order BDF convolution quadrature for subdiffusion models with a singular source term