In recent years, in reaction to the increasing usage of urban areas, the excavation of underground spaces has been developed. One of the most challenging issues encountered by engineers is the construction of subway stations as large underground spaces at shallow depth with soft surrounding soils. In this paper, Naghshe Jahan Square subway station located in Isfahan, Iran, has been simulated by geomechanical finite difference method (FDM). This station is located under important historical structures. Therefore, the ground displacement and surface settlement induced by the excavation of the subway station should be strictly controlled. Many of such problems are affected by selected excavation method. For these reasons, different underground excavation methods associated with construction have been studied. In this study, sequential excavation method and large-diameter curved pipe roofing method are used and the numerical results of the two methods are compared. The presence of groundwater table obliges us to choose special techniques for the stability of the ground around the subway station during construction;hence compressed air and ground freezing techniques are utilized in the simulations of the subway station. Finally, after choosing appropriate support systems, the large-diameter curved pipe roofing method with 1.5 m spacing between curved pipes is proposed.