Quantitative lithofacies modeling is important to understand the depositional and diagenetic history, and hydrocarbon potential of unconventional resources at a regional scale. The complex heterogeneous nature and large data dimensionality of unconventional mudstone reservoirs increase the challenge of lithofacies interpretation by conventional qualitative methods. Quantitative shale lithofacies, which are meaningful, mappable, and predictable at core, well log, and regional scales, can be defined based on mineralogy and Total Organic Carbon (TOC) derived from core analysis and advanced geochemical spectroscopy logs (e.g. Pulsed Neutron Spectroscopy, PNS). However, access to numerous and widespread core samples and geochemical log responses is typically limited by cost and time. We apply different mathematical techniques to ubiquitous conventional well log suites calibrated to rock types, defined by the limited number of wells with high-quality core and geochemical logs. The documented interrelationships between lithofacies and conventional logs are propagated with a quantified degree of accuracy in wells without advanced log or core data. Our study addresses issues of different approaches of quantitative lithofacies classification and