The molecular models of nanopores for major rock constituents in deep shale were constructed. The microscopic adsorption behavior of methane was simulated by coupling the grand canonical Monte Carlo and Molecular Dynamics methods and the effect of rock constituents was discussed. Based on the illite and kerogen nanopore models, the discrepancies in microscopic water distribution characteristics were elucidated, the effects of water on methane adsorption and its underlying mechanisms were revealed, and the competitive adsorption characteristics between water and methane were elaborated. The results show a similar trend in the microscopic distribution of methane between different shale rock constituents. Illite and kerogen slit pores have no significant difference in methane adsorption capacity. The adsorption capacity per unit mass of kerogen is greater than that of illite due to the smaller molar mass of the kerogen skeleton and its large intermolecular porosity. Illite has a greater affinity for water than methane. With increasing water content, water molecules preferentially occupy the high-energy adsorption sites and then overspread the entire pore walls to form water adsorption layers. Methane molecules are adsorbed on the water layers, and methane adsorption has little effect on water adsorption. Kerogen is characterized as mix-wetting. Water molecules are preferentially adsorbed on polar functional groups and gather around to form water clusters. In kerogen with high water content, methane adsorption can facilitate water cluster fusion and suppress water spreading along pore walls. In addition to adsorption, some water molecules dissolve in the kerogen matrix.

Deep marine shale gas reservoirs are extremely rich in the Sichuan basin in China. However, due to the in situ conditions with high temperature and high pressure (HTHP), in particular reservoir pressure being usually much higher than the test pressure, it is difficult to accurately clarify the adsorption behavior, as seepage theory plays an important role in shale gas reserves evaluation. Therefore, three kinds of sorbent, including illite, quartz and kerogen, and two simulation methods, containing the grand canonical ensemble Monte Carlo method and molecular dynamics method, are synthetically used to determine the methane adsorption behavior under HTHP. The results show that both absolute adsorption and excess adsorption decrease with the increase of temperature. When the pressure increases, the absolute adsorption increases quickly and then slowly, and the excess adsorption first increases and then decreases. The superposition of wall potential energy is strongest in a circular hole, second in a square hole, and weakest in a narrow slit. The effect of pore size increases with the decrease of the pore diameter. Under HTHP, multi-layer adsorption can occur in shale, but the timing and number of layers are related to the sorbent type.

Although great success has been achieved in the shale gas industry, accurate production dynamic analyses is still a challenging task. Long horizontal wells coupling with mass hydraulic fracturing has become a necessary technique to extract shale gas efficiently. In this paper, a comprehensive mathematical model of a multiple fractured horizontal well (MFHW) in a rectangular drainage area with a rectangular stimulated reservoir volume (SRV) has been established, based on the conceptual model of “tri-pores” in shale gas reservoirs. Dimensionless treatment and Laplace transformation were employed in the modeling process, while the boundary element method was used to solve the mathematical model. The Stehfest numerical inversion method and computer programing techniques were employed to obtain dimensionless type curves, production rate, and cumulative production. Results suggest that 9 flow stages can be observed from the pseudopressure derivative type curve when the reservoir and the SRV are large enough. The number of fractures, SRV permeability, and reservoir permeability have no effect on the total production when the well is abandoned. As SRV and reservoir permeability increases, the production rate is much higher in the middle production stage. Although the SRV scale and its permeability are very important for early and intermediate production rates, the key factors restricting the shale gas production rate are the properties of the shale itself, such as adsorbed gas content, natural fractures, and organic content. The proposed model is useful for analyzing production dynamics with stimulated horizontal wells in shale gas reservoirs.

The microscopic pore throat structure of shale reservoir rocks directly affects the reservoir seepage capacity. The occurrence and flow channels of shale gas are mainly micron–nanometer pore throats. Therefore, to clarify the microstructural characteristics and influencing factors of the deep organic-rich shales, a study is conducted on the marine shale from the Upper Silurian to Lower Ordovician Wufeng–Longmaxi Formation in the southern Sichuan Basin. Petrographic lithofacies division is carried out in combination with petro-mineralogical characteristics, and a high-resolution scanning electron microscope, low-temperature nitrogen and low-temperature carbon dioxide adsorption, and micron-computed tomography are used to characterize the mineral composition and pore structure qualitatively and quantitatively, upon which the influencing factors of the microstructure are further analyzed. The results show that with the increase in burial depth, the total organic carbon content and siliceous mineral content decrease in the Wufeng formation to Long-1₁ subsection deep shale, while clay mineral content increases, which corresponds to the change in sedimentary environment from anoxic to oxidizing environment. Unexpectedly, the total pore volume of deep shale does not decrease with the increase in burial depth but increases first and then decreases. Using total organic carbon (TOC), siliceous mineral content showed a good correlation with total pore volume and specific surface area, with correlation coefficients greater than 0.7, confirming the predominant role of these two factors in controlling the pore structure of deep shales. This is mainly because the Longmaxi shale is already in the late diagenetic stage, and organic matter pores are generated in large quantities. Clay minerals have a negative correlation with the total pore volume of shale, and the correlation coefficient is 0.7591. It could be that clay minerals are much more flexible and are easily deformed to block the pores under compaction. In addition, the longitudinal heterogeneity of the deep shale reservoir structure in southern Sichuan is also controlled by the thermal effect of the Emei mantle plume on hydrocarbon generation of organic matter and the development of natural microfractures promoted by multistage tectonic movement. Overall, the complex microstructure in the deep shales of the Longmaxi Formation in the southern Sichuan Basin is jointly controlled by multiple effects, and the results of this research provide strong support for the benefit development of deep shale gas in southern Sichuan Basin.

The presence of a reasonable flowback system after fracturing is a necessary condition for the high production of shale gas wells. At present, the optimization of the flowback system lacks a relevant theoretical basis. Due to this lack, this study established a new method for evaluating the conductivity of artificial fractures in shale, which can quantitatively characterize the backflow, embedment, and fragmentation of proppant during the flowback process. Then, the mechanism of the stress sensitivity of artificial fractures on fracture conductivity during the flowback stage of the shale gas well was revealed by performing the artificial fracture conductivity evaluation experiment. The results show that a large amount of proppant migrates, and the fracture conductivity decreases rapidly in the early stage of flowback, and then the decline gradually slows down. When the effective stress is low, the proppant is mainly plastically deformed, and the degree of fragmentation and embedment is low. When the effective stress exceeds 15.0 MPa, the fragmentation and embedment of the proppant will increase, and the fracture conductivity will be greatly reduced. The broken proppant ratio and embedded proppant ratio are the same under the two choke-management strategies. In the mode of increasing choke size step by step, the backflow proppant ratio is lower, and the broken proppant is mainly retained in fractures, so the damage ratio of fracture conductivity is lower. In the mode of decreasing choke size step by step, most of the proppant flows back from fractures, so the damage to fracture conductivity is greater. The research results have important theoretical guiding significance for optimizing the flowback system of shale gas wells.