The production and transportation of high viscosity liquid/gas two-phase along petroleum production system is a challenging operation due to the lack of understanding the flow behavior and characteristics. In particular, accurate prediction of two-phase slug length in pipes is crucial to efficiently operate and safely design oil well and separation facilities. The objective of this study is to develop a mechanistic model to predict high viscosity liquid slug length in pipelines and to optimize the proper set of closure relationships required to ensure high accuracy prediction. A large high viscosity liquid slug length database is collected and presented in this study, against which the proposed model is validated and compared with other models. A mechanistic slug length model is derived based on the first principles of mass and momentum balances over a two-phase slug unit, which requires a set of closure relationships of other slug characteristics. To select the proper set of closure relationships, a numerical optimization is carried out using a large slug length dataset to minimize the prediction error. Thousands of combinations of various slug flow closure relationships were evaluated to identify the most appropriate relationships for the proposed slug length model under high viscosity slug length condition. Results show that the proposed slug length mechanistic model is applicable for a wide range of liquid viscosities and is sensitive to the selected closure relationships. Results revealed that the optimum closure relationships combination is Archibong-Eso et al. (2018) for slug frequency, Malnes (1983) for slug liquid holdup, Jeyachandra et al. (2012) for drift velocity, and Nicklin et al. (1962) for the distribution coefficient. Using the above set of closure relationships, model validation yields 37.8% absolute average percent error, outperforming all existing slug length models.