A three-dimensional time-dependent quantum mechanical approach is used to calculate the reaction probability (PI) and the integral reaction cross section (sigma(R)) for both channels of the reaction He + HD+(nu = 0, 1, 2, 3; j = 0) -> HeH(D)(+) + D(H), over a range of translational energy (E-trans) on two different ab initio potential energy surfaces (McLaughlin -Thompson -Joseph - Sathyamurthy and Palmieri et al.). The reaction probability plots as a function of translational energy exhibit several oscillations, which are characteristic of the system. The vibrational enhancement of the reaction probability and the integral reaction cross section values are reproduced qualitatively by our calculations, in accordance with the experimental results. The isotopic branching ratio for the reaction decreases in going from nu = 0 to nu = 1 and then becomes nearly nu-independent in going from nu = 1 to nu = 3 on both the surfaces.