We investigated models of the internal structure of initially homogeneous Moon differentiated as a result of partial melting, using data on seismic velocities according to the seismic models assume the zonal structure of the lunar mantle is a model of the Moon which was obtained with using the array processing methods of high velocities in the lower mantle. As a result of inversion of gravity (mass, moment of inertia), seismic (P- and S-waves velocities) and petrological (balance ratios) data, the Monte Carlo method was used to reconstruct the chemical composition and internal structure of the Moon. The phase composition and physical properties of the mantle were obtained with Gibbs free energy minimization method and equations of state in the five-component system CaO-FeO-MgO-Al2O3-SiO2. For all models, possible values of seismic velocities and concentrations of the main oxides in three zones of the mantle were obtained, satisfying the geochemical and geophysical constraints and the possible sizes of the Fe-10%S core were determined. It was found that the lunar mantle chemical composition (concentration of FeO, Al2O3 and CaO) differs depending on the mantle zone. Constraints on the values of seismic velocities in the lower mantle and the most probable size of the lunar core were determined: VP ≤ 8.45 km/s; Fe-10%S core radius is ∼360 km.