Sand is often produced with oil and gas even though sand control practices are used. To successfully transport solids in pipelines, the fluid velocity must exceed the critical velocity. Solids transport models are used to predict this velocity. However, for the same input field or laboratory condition, the models' critical velocity predictions may differ by orders of magnitude. Furthermore, none of the models provide information regarding the confidence in their velocity predictions. This paper introduces a systematic methodology, which selects the appropriate solids transport models for predicting the operational critical velocity envelope to ensure solids transport in the pipe to within a predetermined confidence level. Given a field or laboratory condition as input, the approach generates: (1) velocity predictions and uncertainty bounds of all models to within a predetermined confidence level, (2) recommendation of appropriate models for calculating the critical velocity at the input condition, and (3) impact of the uncertainty of input condition, experimental data, and model on the overall uncertainty bounds. The methodology uses data clustering, model evaluation, optimization, and uncertainty propagation approaches. Based on three Case Studies, the velocity envelopes suggested by the methodology, at the 90% confidence level, are consistent with experimental observations. For input laboratory conditions, which are used for validation purposes, the contributions of model, input condition, and experimental data uncertainties to the velocity prediction uncertainty were between 50 -90%, 3 -20%, and 6 -40%, respectively. If the parameters of the selected models were fine-tuned, for the same input conditions, the above uncertainties contribute about 40 -90%, 4 -15%, and 1 -50%, respectively, to the velocity prediction uncertainty. For an input field condition (the fourth Case Study), the above values become about 35 -70%, 7 -25%, and 25 -45%, respectively; and about 25 -30%, 30 -32%, and 40 -42%, respectively, if the parameters of the models are fine-tuned prior to uncertainty estimation. The results of these Case Studies suggest that for laboratory conditions, most of the uncertainty in a model's critical velocity prediction comes from the model, while for field conditions, the uncertainties of the model, input condition, and data may have nearly equal contributions.