Non aqueous phase liquid (NAPL) pollution exists in porous media such as soil. SVE technology can be used to remove this pollution in soil. However, few domestic and international studies have paid attention to the changes of soil temperature in the field, which we believe can be useful information to optimize the layout of heating wells. In this research we established partial differential equations of soil heat transfer using the COMSOL multi-field coupling tool to simulate the field distribution of the change in soil internal temperature in the process of SVE to obtain the change of effective heating area with time under certain initial heating conditions. At the same time, we used liquid ethylbenzene to represent NAPL pollutants, and designed a simulation of soil temperature field distribution under the movement of liquid ethylbenzene under external pressure. It was found that under the action of Darcy’s velocity field, the utilization efficiency of the SVE system for the heat source was significantly improved, that is, the temperature distribution of the soil was more uniform. However, the temperature of the heated area increased slowly because the extraction well took away the heat energy. The heat source power should be increased or prolonged to improve the effect of Darcy’s field. Through a coupled simulation, we obtained a variation relationship of the soil temperature field in 1800 min under the action of one extraction well and four heating wells. These data will provide the basis for our next step in designing an algorithm to optimize the distribution of heating wells.