Despite the growing interest in the potential electrochemical applications of both aluminium and ionic liquids in batteries, the microstructure of mixtures of trivalent salts and these dense ionic environments is completely unknown. In this work, the solvation of Al3+ cations in highly dense ionic solvents is investigated. For this purpose, molecular dynamics simulations of mixtures of a protic ionic liquid, ethylammonium nitrate (EAN), with aluminium nitrate (Al(NO3)3), both in bulk and confined between graphene walls, are performed. Several structural quantities of the system are calculated for different salt concentrations, such as densities, radial distribution functions, structure factors, coordination numbers and hydrogen bonds for the bulk mixture and ionic density profiles for the confined ones. Moreover, vibrational density of states is calculated for the salt cations, both in bulk and when close to the walls. The results obtained are analyzed and compared to those for mixtures of EAN with monovalent and divalent salts, in order to probe the influence of the salt cation charge on the system's properties. Finally, ab initio density functional theory calculations were performed in order to analyze the structure of the Al3+-ligand complexes, and their predictions for the Raman spectrum are compared both to the corresponding experimental one and the one coming from molecular dynamics simulations. According to our calculations, [Al(NO3)6]3- octahedral complexes do not significantly change the microstructure of the mixtures relative to those of Mg2+-based ones.