BTEX compounds are among the light hydrocarbons of gasoline which have been noticed due to their toxicity and long term damaging effects in the environment. Understanding and prediction of how contaminations propagation are of great importance for designing an appropriate remediation method. In this paper, propagation pattern of BTEX compounds in saturated and unsaturated zones of soil is investigated by the finite element method. To investigate the accuracy and reliability of flow and transport modules of MOFAT software two different numerical analyses have conducted based on experimental and in-situ data and the obtained results have compared to results from previous studies. In the following, effect of physical factors of soil and contamination on the propagation pattern is studied based on flow module. The obtained results show that as soil hydraulic conductivity increases effect of NAPL’s physical factors in propagation rate are become more tangible. It is also observed that propagation depth of BTEX compounds is very sensitive to van Genuchten parameter α. The propagation rate of soluble components of BTEX compounds in aqueous, gas and NAPL phases have investigated by a larger model and by Simultaneous utilization of flow and transport modules.