Shale gas reservoir evaluation and production optimization both require geological models. However, currently, shale gas modeling remains relatively conventional and does not reflect the unique characteristics of shale gas reservoirs. Based on a case study of the Fuling shale gas reservoir in China, an integrated geological modeling workflow for shale gas reservoirs is proposed to facilitate its popularization and application and well improved quality and comparability. This workflow involves four types of models: a structure-stratigraphic model, reservoir (matrix) parameter model, natural fracture (NF) model, and hydraulic fracture (HF) model. The modeling strategies used for the four types of models vary due to the uniqueness of shale gas reservoirs. A horizontal-well lithofacies sublayer calibration-based method is employed to build the structure-stratigraphic model. The key to building the reservoir parameter model lies in the joint characterization of shale gas “sweet spots.” The NF models are built at various scales using various methods. Based on the NF models, the HF models are built by extended simulation and microseismic inversion. In the entire workflow, various types of models are built in a certain sequence and mutually constrain one another. In addition, the workflow contains and effectively integrates multisource data. Moreover, the workflow involves multiple model integration processes, which is the key to model quality. The selection and optimization of modeling methods, the innovation and development of modeling algorithms, and the evaluation techniques for model uncertainty are areas where breakthroughs may be possible in the geological modeling of shale gas reservoirs. The workflow allows the complex process of geological modeling of shale gas reservoirs to be more systematic. It is of great significance for a dynamic analysis of reservoir development, from individual wells to the entire gas field, and for optimizing both development schemes and production systems.