Phthalic acid esters (PAEs) are emerging pollutants that need to be analyzed precisely. Chromatography-based determination of PAE content in soils are frequently affected by matrix effect, which may limit the quantification of different kinds of PAEs from different types of soil. Here we optimized a QuEChERS protocol combined with gas chromatography-mass spectrometry (GC-MS) for simultaneous determination of 16 PAEs in different soils. PAEs in different type of soils (fluvo-aquic soil, red soil, and black soil) were extracted with acetonitrile followed by GC-MS detection based on quantitative ion internal standard method. All 16 PAEs showed excellent linear relationships with mass peak areas (R2 > 0.99). The limits of detection (LOD) and limits of quantitation (LOQ) of all the samples were in the range of 0.91-66.97µg/kg and 2.7-200.9µg/kg, respectively. The accurate test at 0.5, 0.1, and 1.0mg/kg spiking level recorded recovery rate between 80.11% and 100.99% with relative standard deviations (RSDs) ranging from 0.37 to 8.50% in tested matrices. No significant matrix effect was observed for most tested PAEs. This is a simple method with high sensitivity and strong stability, which is suitable and reproducible for quantifying large number of PAEs in different types of soil.