The purpose of this study was to investigate the impact of different low-energy virtual monochromatic images (VMIs) in dual-energy CT on the performance of radiomics models for predicting muscle invasive status in bladder cancer (BCa). A total of 127 patients with pathologically proven muscle-invasive BCa (n = 49) and non-muscle-invasive BCa (n = 78) were randomly allocated into the training and test cohorts at a ratio of 7:3. Feature extraction was performed on the venous phase images reconstructed at 40, 50, 60 and 70-keV (single-energy analysis) or in combination (multi-energy analysis). Recursive feature elimination (RFE) and the least absolute shrinkage and selection operator(LASSO) were employed to select the most relevant features associated with BCa. Models were built using a support vector machine(SVM) classifier. Diagnostic performance was assessed through receiver operating characteristic curves, evaluating sensitivity, specificity, accuracy, precision, and the area-under-the curve (AUC) values. In the test cohort, the multi-energy model achieved the best diagnostic performance with AUC, sensitivity, specificity, accuracy, and precision of 0.917, 0.800, 0.833, 0.821, and 0.750, respectively. Conversely, the single-energy model exhibited lower AUC and sensitivity in predicting the muscle invasion status. By combining information from VMIs of various energies, the multi-energy model displays superior performance in preoperatively predicting the muscle invasion status of bladder cancer.