It is strange to me that for a symmetry which involves $\dot{x}$, there seems to always appear a term with $\dddot{x}$ in the variation of the equations of motion, which doesn't makes much sense. I think that probably the procedure I am following is wrong.
I will show you an example: Consider the simple case of a free particle in one dimension, it's Lagrangian is:
$$L=\frac{1}{2}\dot{x}^2$$
It is obvious that the system conserves energy, so the symmetry that must be valid is $(\delta x=0 ,\delta t=\epsilon)$. I may rewrite these as symmetries that doesn't involve time variations (as the paper by E. L. Hill do):
$$(\delta_{*}x=\epsilon\dot{x},\delta_{*}t=0)$$
Now, I calculate the variation of the equations of motion, with the hope of finding that $\delta(\text{e.o.m})=\text{e.o.m}.$ Such result would mean that the equations of motion are invariant under the symmetry in consideration. So:
$$\delta{(\text{e.o.m})}=\delta(\ddot{x})=\ddot{\eta}$$
In this case $\eta = \dot{x}$ (remember that a variation is of the form $\delta{x}=\epsilon\eta$). So $\ddot{\eta}=\dddot{x}$. Hence:
$$\delta(\text{e.o.m})=\dddot{x}\neq\text{e.o.m}$$
This doesn't make sense because, for the system in consideration, the time translation is a Noetherian symmetry that gives conservation of Energy.
My question is: What is failing in this procedure? Is there a general way of showing that some Symmetry is indeed Noetherian?