Abstract: Based on Python and Paraview software, the 3D suitability evaluation of 0~-100m underground space of Zhongshan city was realized by bilinear interpolation method and analytic hierarchy process. The data of 8192 bore holes was collected and analyzed, and five 3D data of single factor (spt blow count, water content, internal friction angle, compression modulus and compressive strength) was obtained. In addition, the geotechnical conditions, hydrogeological conditions, adverse geological conditions, resource sensitivity and development value were selected as theme layer, them contain 21 signal factor 2D data. The coupling calculation of 2D and 3D data is realized based on the Variable weight method. According to the sensitivity of different engineering types to single factor, the architectural engineering and tunnel engineering are evaluated respectively. This evaluation was wide area, large depth and many evaluation factors, which can provide reference for the suitability evaluation of underground space development and utilization in other cities.