This study is to simulate the transport phenomena of LNAPLs and DNAPLs spilled in saturated and unsaturated zones. The model system consists of four heterogeneous layers with different permeabilities and goundwater flow at the bottom layer. In the study, LNAPL and DNAPL are assumed to be n-hexane and 1,1-dichloroacetone, respectively. Results show that the transport phenomena and morphology of the NAPLs change with the extent of water saturation, groundwater flow and buoyancy. In the unsaturated zone, the DNAPL migrates vertically two times faster than the LNAPL due to the gravity. Also, the lateral deformation of the DNAPL becomes smaller than that of the LNAPL in the interface between the layers. When the spilled NAPLs reach at the top of the water-table, the LNAPL forms the continuous film on the water table and spreads horizontally with water flow. In contrast, the DNAPL builds up the discontinuous plumes on the water-table and creates the several nuclears under the water table. Also, the profiles of volume fraction of DNAPL and LNAPL in vertical direction are estimated. Results would provide a tool to understand the transport phenomena of NAPLs spilled in the underground and provide a basis for the selection of the remediation technologies of both soils and ground-water polluted with NAPLs.