In　this　paper,　ETD3-Pade　and　ETD4-Pade　Galerkin　finite　element　methods　are　proposed　and　analyzed　for　nonlinear　delayed　convection-diffusion-reaction　equations　with　Dirichlet　boundary　conditions.　An　ETD-based　RK　is　used　for　time　integration　of　the　corresponding　equation.　To　overcome　a　well-known　difficulty　of　numerical　instability　associated　with　the　computation　of　the　exponential　operator,　the　Pade　approach　is　used　for　such　an　exponential　operator　approximation,　which　in　turn　leads　to　the　corresponding　ETD-Pade　schemes.　An　unconditional　L-2　numerical　stability　is　proved　for　the　proposed　numerical　schemes,　under　a　global　Lipshitz　continuity　assumption.　In　addition,　optimal　rate　error　estimates　are　provided,　which　gives　the　convergence　order　of　O(k(3)　+　h(r))　(ETD3-Pade)　or　O(k(4)　+　h(r))　(ETD4-Pade)　in　the　L-2　norm,　respectively.　Numerical　experiments　are　presented　to　demonstrate　the　robustness　of　the　proposed　numerical　schemes.