The lymphovascular invasion (LVI) status facilitates the determination of the optimal therapeutic strategy for superficial esophageal squamous cell carcinoma (SESCC), but in clinical practice, LVI must be confirmed by postoperative pathology. However, studies of the risk factors for LVI in SESCC are limited. Consequently, this study aimed to identify the risk factors for LVI and use these factors to establish a prediction model. The data of 516 patients who underwent radical esophagectomy between January 2007 and September 2019 were retrospectively collected (training set, n=361, January 2007 to May 2015; validation set, n=155, June 2015 to September 2019). In the training set, least absolute shrinkage and selection operator (LASSO) regression and multivariate analyses were utilized to identify predictive factors for LVI in patients with SESCC. A nomogram was then developed using these predictors. The area under the curve (AUC), calibration curve, and decision curve were used to evaluate the efficiency, accuracy, and clinical utility of the model. LASSO regression indicated that the tumor size, depth of invasion, tumor differentiation, lymph node metastasis (LNM), sex, circumferential extension, the presence of multiple lesions, and the resection margin were correlated with LVI. However, multivariate analysis revealed that only the tumor size, depth of invasion, tumor differentiation, and LNM were independent risk factors for LVI. Incorporating these four variables, model 1 achieved an AUC of 0.817 in predicting LVI. Adding circumferential extension to model 1 did not appreciably change the AUC and integrated discrimination improvement, but led to a significant increase in the net reclassification improvement (p=0.011). A final nomogram was constructed by incorporating tumor size, depth of invasion, tumor differentiation, LNM, and circumferential extension and showed good discrimination (training set, AUC=0.833; validation set, AUC=0.819) and good calibration in the training and validation sets. Decision curve analysis demonstrated that the nomogram was clinically useful in both sets. Thus, it is possible to predict the status of LVI using this nomogram scoring system, which can aid the selection of an appropriate treatment plan.