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 application of horizontal wells with multistage hydraulic fracturing technologies has made the development of shale gas reservoirs become a worldwide economical hotspot in recent years. The gas transport mechanisms in shale gas reservoirs are complicated, due to the multiple types of pores with complex pore structure and special process of gas accumulation and transport. Although there have been many attempts to come up with a suitable and practical mathematical model to characterize the shale gas flow process, no unified model has yet been accepted by academia. In this paper, a comprehensive literature review on the mathematical models developed in recent years for describing gas flow in shale gas reservoirs is summarized. Five models incorporating different transport mechanisms are reviewed, including gas viscous flow in natural fractures or macropores, gas ad-desorption on shale organic, gas slippage, diffusion (Knudsen diffusion, Fick diffusion, and surface diffusion), stress dependence, real gas effect, and adsorption layer effect in the nanoshale matrix system, which is quite different from conventional gas reservoir. This review is very helpful to understand the complex gas flow behaviors in shale gas reservoirs and guide the efficient development of shale gas. In addition to the model description, we depicted the type curves of fractured horizontal well with different seepage models. From the review, it can be found that there is some misunderstanding about the essence of Knudsen/Fick diffusion and slippage, which makes different scholars adopt different weighting methods to consider them. Besides, the contribution of each mechanism on the transport mechanisms is still controversial, which needs further in-depth study in the future.

Shale gas reservoirs are unconventional resources with great potential to help meet energy demands. Horizontal drilling and hydraulic fracturing have been extensively used for the exploitation of these unconventional resources. According to engineering practice, some shale gas wells with low flowback rate of fracturing fluids may obtain high yield which is different from the case of conventional sandstone reservoirs, and fracturing fluid absorbed into formation by spontaneous imbibition is an important mechanism of gas production. This paper integrates NMR into imbibition experiment to examine the effects of fractures, fluid salinity, and surfactant concentration on imbibition recovery and performance of shale core samples with different pore-throat sizes acquired from the Longmaxi Formation in Luzhou area, the Sichuan Basin. The research shows that the right peak of T₂ spectrum increases rapidly during the process of shale imbibition, the left peak increases rapidly at the initial stage and changes gently at the later stage, with the peak of the left peak shifting to the right. The result indicates that water first enters the fracture system quickly, then enters the small pores near the fracture wall due to the effect of the capillary force, and later gradually sucks into the deep and large pores. Both imbibition rate and capacity increase with increased fracture density, decreased solution salinity, and decreased surfactant concentration. After imbibition flowback, shale permeability generally increases by 8.70–17.88 times with the average of 13.83 times. There are also many microcracks occurring on the end face and surface of the core sample after water absorption, which may function as new flowing channels to further improve reservoir properties. This research demonstrates the imbibition characteristics of shale and several relevant affecting factors, providing crucial theory foundations for the development of shale gas reservoirs.

In order to account for big uncertainties such as well interferences, hydraulic and natural fractures’ properties and matrix properties in shale gas reservoirs, it is paramount to develop a robust and efficient approach for well spacing optimization. In this study, a novel well spacing optimization workflow is proposed and applied to a real shale gas reservoir with two-phase flow, incorporating the systematic analysis of uncertainty reservoir and fracture parameters. One hundred combinations of these uncertainties, considering their interactions, were gathered from assisted history matching solutions, which were calibrated by the actual field production history from the well in the Sichuan Basin. These combinations were used as direct input to the well spacing optimization workflow, and five “wells per section” spacing scenarios were considered, with spacing ranging from 157 m (517 ft) to 472 m (1550 ft). An embedded discrete fracture model was used to efficiently model both hydraulic fractures and complex natural fractures non-intrusively, along with a commercial compositional reservoir simulator. Economic analysis after production simulation was then carried out, by collecting cumulative gas and water production after 20 years. The net present value (NPV) distributions of the different well spacing scenarios were calculated and presented as box-plots with a NPV ranging from 15 to 35 million dollars. It was found that the well spacing that maximizes the project NPV for this study is 236 m (775 ft), with the project NPV ranging from 15 to 35 million dollars and a 50^th percentile (P50) value of 25.9 million dollars. In addition, spacings of 189 m (620 ft) and 315 m (1033 ft) can also produce substantial project profits, but are relatively less satisfactory than the 236 m (775 ft) case when comparing the P25, P50 and P75 values. The results obtained from this study provide key insights into the field pilot design of well spacing in shale gas reservoirs with complex natural fractures.