We construct bivariate polynomial approximations of the Lerch function, which for certain specialisations of the variables and parameters turn out to be HermitePadé approximants either of the polylogarithms or of Hurwitz zeta functions. In the former case, we recover known results, while in the latter the results are new and generalise some recent works of Beukers and Prévost. Finally, we make ...