Abstarct: Today, numerical methods are used as an effective tool for predicting behavior of ground. Tendency to use the numerical methods have increased due to the ability to simulate different material, high speed and precision, as well as the ability to change the properties of the ground and structures materials to investigate the effect of each of these factors, in solving geomechanical problems. In recent years, the need for underground spaces has increased due to the expansion of urban area. It should be note that most cities are located in alluvial areas. Therefore, ground surface settlement and its impact on surrounding structures must be considered in excavation of large underground spaces such as Tunnels and Statins of subways. So it is important to analysis of the stability of different methods of subway station. Typically, subway stations are drilled by Cut and Field method or underground method. The selection of subway station drilling method is depended on various conditions such as station depth, traffic Considerations, economic aspects, and so on. Each method of drilling the metro station has advantages and disadvantages, but it is a suitable method that, in addition to having less disadvantages, provides underground space safety and long-term sustainability, and is also the best option in terms of cost. In this research, in order to select the methods for numerical investigation, at the first various methods of drilling the subway stations have been studied on variety of Sequential methods, Mechanical methods and variety of Pre-Supporting methods. The Sequential methods including Side Drift, Centural Diaphragm. The Mechanical methods including Mechanical Pre-Tunnelin methods and the use of Multi-face Shield Construction. Also the Pre-Supporting , including Carved Pipe Roofing Method and Celular Arc Method . Finally, the Concrete Arc Pre-Supporting System method and Side Drift method have been selected for comparison. Then, the FLAC3D software has been used to investigate the stability of subway station drilling methods. The data of the Marzdaran station of subway-Line 6 of Tehran have been used to make a numerical model. The Marzdaran station of subway-Line 6 of Tehran is built by Concrete Arc Pre-Supporting System method. The effect of some soil parameters and maintenance elements is also analyzed using numerical modeling. The results of the numerical modeling are compared with the results instrumentation of The Marzdaran station of subway-Line 6 of Tehran. Also, the approximate cost of drilling and maintenance operations in both methods is estimated and compared baxsed on the road, runway and airport prices. Finally, the most suitable method for drilling The Marzdaran station of subway-Line 6 of Tehran is presented in terms of sustainability and cost.