Gas reservoir numerical simulation is an important method to optimize the development strategy of shale gas reservoirs which has been influenced by the multi-stage fracture. The regular fracture network model was used to build a conventional numerical simulation, in which it was difficult to show the true situation of fracture propagation. However, the physical parameters not only affect the production, but also influence the stimulation effect; moreover, the quality of the fracturing effect also affects the production which causes the input and out parameters to be inaccurate. To solve this problem, the process simulation must be completed from geology to engineering to gas reservoir. The main controlling factors of production are identified with geological and engineering factors such as horizontal stage length, the volume of fracturing fluid, well spacing, production allocation, and proppant mass. Therefore, on the basis of the integrated simulation method of a hydraulic fracturing network simulation and an unstructured grid high-precision numerical simulation, this paper builds an integrated numerical simulation of a shale gas reservoir coupled with geology and engineering to optimize the development strategy with production as the target. Taking four wells of a platform as an example, the EUR (estimated ultimate recovery) has increased by 25% after the optimization of the development strategy.