In this paper, we present a general analytical model for light scattering by arbitrary Vogel spiral arrays of circular apertures illuminated at normal incidence. This model suffices to unveil the fundamental mathematical structure of their complex Fraunhofer diffraction patterns and enables the engineering of optical beams carrying multiple values of orbital angular momentum (OAM). By performing analytical Fourier-Hankel decomposition of spiral arrays and far field patterns, we rigorously demonstrate the ability to encode specific numerical sequences onto the OAM values of diffracted optical beams. In particular, we show that these OAM values are determined by the rational approximations (i.e., the convergents) of the continued fraction expansions of the irrational angles utilized to generate Vogel spirals. These findings open novel and exciting opportunities for the manipulation of complex OAM spectra using dielectric and plasmonic aperiodic spiral arrays for a number of emerging engineering applications in singular optics, secure communication, optical cryptography, and optical sensing.