Estimation of urban seismic damage using numerical simulation needs an automatic modeling method of surface layers and residential buildings. This study focuses on modeling of surface layers and shows a method of constructing models by interpolating boring data. An important property of the modeling method is robustness, that means that the method works for boring data with inconsistent soil layers. To satisfy this, we developed the method using artificial layers. We applied the method to a test site and checked its robustness. This test also showed that the method gave realistic models. Finally, we applied the method to the estimation of urban seismic damage and discussed the usefulness by comparing the result with one obtained by a conventional method.