Based on core observations, thin sections, X-ray diffraction (XRD), and seismic data, the lithofacies types in the organic-rich Longmaxi shale (Lower Silurian) in the Changning area of the southern Sichuan Basin were identified. The factors controlling the spatial variations in the shale lithofacies and the influences of the shale lithofacies on shale gas development were also analyzed. Results indicate that there are seven main types of shale lithofacies in the Long11 sub-member of the Longmaxi Formation, including siliceous shale (S-1), mixed siliceous shale (S-2), carbonate-rich siliceous shale (S-3), clay-rich siliceous shale (S-4), carbonate/siliceous shale (M-1), mixed shale (M-2), and argillaceous/siliceous shale (M-4). A vertical transition from the carbonate shale association + mixed shale association at the bottom of the sub-member to a siliceous shale association and mixed shale association + siliceous shale at the top generally appears in the Long11 sub-member. The shale lithofacies of the Long11 sub-member also laterally change from the central depression (low-lying area) to the geomorphic highland in the east and west parts of the Changning area. The spatial variations in shale lithofacies in the Long11 sub-member of the Changning area were mainly controlled by palaeogeomorphology and relative sea level. The geomorphic highland area is dominated by carbonate-rich siliceous shale and mixed siliceous shale, but the depression (low-lying area) is mainly dominated by mixed siliceous shale and argillaceous/carbonate shale.

There are a large number of natural fractures in shale reservoirs, which create great challenges to hydraulic fracturing. Activating the natural fractures in reservoirs can form a complex fracture network, enhance fracturing effects, and increase shale gas production. Reservoir geological conditions (low in situ stress, natural fracture distribution, and cement strength) and operation parameters (fracturing fluid viscosity and injection rate) have an important influence on fracture network propagation. In this article, a two-dimensional hydraulic fracturing fluid-mechanic coupling numerical model for shale reservoirs with natural fractures was established. Based on the global cohesive zone model, the influence of geological conditions and operation parameters on the propagation of the hydraulic fracture network and fracturing process is investigated. The numerical simulation results show that when the horizontal in situ stress difference, approach angle, and cement strength are low, it is easier to form a complex fracture network. Research on the construction parameters indicated that when the viscosity of the fracturing fluid is low, it is easier to form a complex network of fractures, but the length of the fractures is shorter; in contrast, the fractures are straight and long. In addition, increasing the injection rate is beneficial for increasing the complexity of the fracture network while increasing the initiation pressure and width of the principal fracture reduces the risk of sand plugging. This article also proposes an optimization solution for hydraulic fracturing operations based on numerical simulation results.