A general procedure which computes X-ray absorption near-edge structure (XANES) spectra starting from the individual configurations obtained from molecular dynamics (MD) simulations has been developed. This method allows one to perform a quantitative analysis of the XANES spectra of disordered systems and to estimate the effect of disorder on the low-energy region of the X-ray absorption spectra. As a test case, a Ni2+ aqueous solution has been investigated. The configurational averaged XANES theoretical spectrum obtained from the MD configurations compares well with the experimental data. This confirms simultaneously the reliability of the procedure and of the structural results obtained from the MD simulation.