Modeling added spatial variability due to soil improvement: Coupling FEM with binary random fields for seismic risk analysis