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