The vibrationally resolved absorption spectra of zinc phthalocyanine (ZnPc) aggregates (up to 70 monomers) are explored using the non-Markovian stochastic Schrödinger equation. Various types of local excitations, charge-transfer (CT) excitations, and exciton-phonon couplings are explicitly included in a comprehensive model Hamiltonian, which is parameterized by first-principles calculations. The absorption spectral simulations clarify that the two absorption bands in the Q-band region observed in experiments can be assigned to the contribution from the CT-mediated interactions, rather than the mixtures of different-type aggregates, as prevailingly assumed. Furthermore, the relative intensities of the two bands are found to be closely related to the intermolecular distance and molecular number in a ZnPc aggregate. From the investigation of the decoherence process after optical excitation, it is found that CT states can induce coherence regeneration as the time scale of charge separation is much faster than that of the vibration-induced decoherence. However, they would instead boost the decoherence process as the two time scales become comparable. The two different effects of CT states may suggest a novel way to regulate the decoherence process in excitation energy relaxation.