Individual survival prediction and risk stratification are of vital importance to optimize the individualized treatment of metastatic leiomyosarcoma (LMS) patients. This study aimed to identify the prognostic factors for metastatic LMS patients and establish prognostic models for overall survival (OS) and cancer-specific survival (CSS). The data of LMS patients with metastasis between 2010 and 2015 were extracted from the Surveillance, Epidemiology, and End Results (SEER) database. The entire cohort was randomly divided into a training cohort and a validation cohort. The influences of primary tumor site, localized and distant metastases, and sites and number of metastases on the prognosis of metastatic LMS patients were firstly explored by Kaplan–Meier curves and log-rank tests. Furthermore, the effective therapeutic regimens and prognosticators for metastatic LMS patients were also analyzed by Cox analysis. In addition, two prognostic nomograms for OS and CSS were established, and their predictive performances were evaluated by the methods of receiver operating characteristic (ROC) curves, time-dependent ROC curves, calibration curves, and decision curve analysis (DCA). A total of 498 patients were finally collected from the SEER database and were randomly assigned to the training set (N = 332) and validation set (N = 166). No significant differences in OS were observed in patients with distant organ metastasis and localized metastasis. For patients who have already developed distant organ metastasis, the sites and number of metastases seemed to be not closely associated with survival. Patients who received chemotherapy got significantly longer survival than that of their counterparts. In univariate and multivariate Cox analyses, variables of surgery, chemotherapy, age, and tumor size were identified as independent predictors for OS and CSS, and distant metastasis was also independently associated with CSS. The areas under the curve (AUCs) of ROC curves of the nomogram for predicting 1-, 3-, and 5-year OS were 0.770, 0.800, and 0.843, respectively, and those for CSS were 0.777, 0.758, and 0.761, respectively. The AUCs of time-dependent AUCs were all over 0.750. The calibration curves and DCA curves also showed excellent performance of the prognostic nomograms. Metastasis is associated with reduced survival, while the sites and the number of metastases are not significantly associated with survival. The established nomogram is a useful tool that can help to perform survival stratification and to optimize prognosis-based decision-making in clinical practice.