Abstract: Objective: Qingdao Metro Line 6 Huashan 1st Road Station is in complex and heterogeneous stratum formed by diabase intrusion into granite formation, different construction methods adopted for the station excavation will cause distinct mechanical responses, surrounding rock deformation and land subsidence, thus the optimal construction method needs to be selected. Method:Based on the engineering project of the above station (located in the complex and heterogeneous stratum of diabase intrusion into granite formation), a 3D numerical simulation analysis is carried out on four construction methods, including CD (center diaphragm) method, CRD (cross diaphragm) method, doubleside heading method and primary support archcover method. The land subsidence, tunnel hole surrounding rock deformation and station plastic zone distribution due to different construction methods are compared and the optimal construction method is selected. Results & Conclusion: After station excavation, the surrounding rock exhibits asymmetric deformation, with greater land subsidence and surrounding rock deformation on the diabase side compared to the granite side; the doubleside heading method and primary support archcover method demonstrate better control over land subsidence and surrounding rock deformation; in the aspect of deformation and destruction area development, the doubleside heading method and primary support archcover method outperform the CRD method, and the CRD method outperforms the CD method. Therefore, the doubleside heading method and primary support archcover method are effective in controlling deformation, and both methods are deemed feasible based on numerical analysis. The primary support archcover method is considered the optimal construction method. Monitoring results indicate that the primary support archcover method effectively controls surrounding rock deformation and land subsidence.