  We study the quasi-normal modes of the charged scalar perturbations in the
background of the Einstein-Maxwell-aether black hole through three methods (WKB
method, continued fraction method, generalized eigenvalue method). Then we
propose the specific treatment for the generalized effective potential with
$\omega$-dependence and the complete procedure of transforming calculation
continued fraction method into finding the zero point of the corresponding
complex function numerically. These methods are valid because the results from
different methods are consistent. We also investigate the allowed region of the
second kind aether black hole among the system parameters ($c_{13},c_{14},Q$).
Finally we show the existence of quasi-resonances of massive perturbation for
Einstein-Maxwell-aether black hole even with large aether parameter.
