Femtosecond laser-induced volume nanograting formation is numerically investigated. The developed model solves nonlinear Maxwell's equations coupled with multiple rate free carrier density equations in the presence of randomly distributed inhomogeneities in fused silica. As a result of the performed calculations, conduction band electron density is shown to form nanoplanes elongated perpendicular to the laser polarization. Two types of nanoplanes are identified. The structures of the first type have a characteristic period of the laser wavelength in glass and are attributed to the interference of the incident and the inhomogeneity-scattered light waves. Field components induced by coherent multiple scattering in directions perpendicular to the laser polarization are shown to be responsible for the formation of the second type of structures with a subwavelength periodicity. In this case, the influence of the inhomogeneity concentration on the period of nanoplanes is shown. The calculation results not only help to identify the physical origin of the self-organized nanogratings, but also explain their period and orientation.