Numerical modeling of discrete spatial heterogeneity in seismic risk analysis: application to treated ground soil foundation