Approximation for Lambert W function near zero

I don't know how simple you need it, and since you never said anything on how accurate you want your approximant to be (i.e., to how many correct decimal places should the approximant match the Lambert function?),

$$W_0(z)\approx\ln(1+z)\frac{1+\frac{123}{40}z+\frac{21}{10}z^2}{1+\frac{143}{40}z+\frac{713}{240}z^2}$$

should be good enough, which has a maximum error of around $1.6\times 10^{-4}$ for $z\in[0,e]$.

The rational portion here is a Padé approximant; probably one might do better with a minimax rational approximation, but I don't have the patience and inclination to derive it since your question's rather vague to begin with.