Search

The sealing capacities of shale, roof, and floor strata are essential for the preservation of shale gas. Sealing capacity development and its influence on shale gas composition and organic matter pore development are important but have not been addressed by existing studies. In this paper, the Upper Ordovician Wufeng-the first member of the Lower Silurian Longmaxi Fm. in the Yangtze region is selected as the study object. We dissected the J1 well in the Jiaoshiba block inside the Sichuan Basin and the Y6 well in the Youyang block outside the Sichuan Basin to analyze the gas components and carbon isotopes in recovered shale gas. In this study, two microscope techniques (i.e., FIB-SEM and FIB-HIM), porosity measurement, TOC content analysis, permeability tests, isothermal adsorption measurements, are utilized. The analytical results show that both the roof and floor strata have vertical sealing effects on the organic-rich shale due to their different physical properties compared to the shale. The shale overburden pressure increases with burial depth, and the adsorbed gas increases with increasing TOC content and thickness: all of these can lead to a decrease in permeability and an increase in self-sealing capacity. This paper determined the formation mechanism of the sealing capacity for shale, roof and floor strata. Furthermore, the control patterns of "sealing capacity strength" on the shale gas source and organic matter pore development are clarified. The main findings are as follows: when the shale, roof, and floor strata have a strong sealing capacity, only a small amount of liquid hydrocarbons can be discharged from the shale, but the natural gas formed due to liquid hydrocarbon cracking is abundant. Under such circumstances, shale gas has a mixed origin of liquid hydrocarbon cracked gas and kerogen cracked gas; the carbon isotopes are inverted; abundant pyrobitumen pores with high porosity and good connectivity are developed. When the sealing capacities are weak, the amount of liquid hydrocarbons discharged from the shale are high, and the amount of natural gas is low. The shale gas is mainly derived from kerogen cracking with weak or no carbon isotope inversion, and kerogen pores are small with poor connectivity. These findings will advance the study of shale gas formation mechanisms and provide a scientific basis for selecting shale gas exploration area. • Due to different physical properties, the roof and floor strata provide a seal for the organic-rich shale. • As the burial depth, TOC content, and the shale thickness increase, the shale has a stronger sealing capacity. • Sealing capacity plays an important role in the characteristics of shale gas and organic matter pore development. [ABSTRACT FROM AUTHOR]

• CH 4 and (sub and supercritical) CO 2 adsorption experiments on shales were performed. • The adsorption ratio of CO 2 over CH 4 for tested shale samples ranges from 1.66–8.32. • Langmuir, D-A, D-R and Ono-Kondo models were fitted to the adsorption data. • The favorable models used to fit CH 4 and CO 2 adsorption on shales are recommended. A thorough study of the adsorption behavior of shales to methane (CH 4) and carbon dioxide (CO 2) is critical for carbon dioxide sequestration in shale gas reservoirs and enhanced shale gas recovery. This paper discusses the results of an adsorption study of CH 4 and CO 2 , in single gas environment, on a set of shale samples taken from the Sichuan and Ordos Basins. The results indicated that shale exhibit higher affinity to CO 2 as compared to CH 4 under similar pressure and temperature conditions, and the preferential adsorption ratio of CO 2 over CH 4 varies between 1.66 and 8.32. Furthermore, the experimental data were modeled using Langmuir, Dubibin-Astakhov (D-A), Dubinin-Redushckevich (D-R) and Ono-Kondo models. The accuracy of the models in quantifying CH 4 and CO 2 adsorption on shale was compared using an error analysis technique. For CH 4 , both of Langmuir and D-A models performed satisfactorily with comparable accuracy. For subcritical CO 2 , the performance of the D-A model is the best. However, for supercritical CO 2 , Ono-Kondo model fitted the experimental data the best, which is the proposed model to use. The results of this study can provide the basis for the estimation of the CO 2 sequestration capacity in shale formations. [ABSTRACT FROM AUTHOR]

In shale plays, the in-situ effective stress is continuously being perturbed because of the production which impacts field operations. Therefore, mechanical properties of the formation and its response to stress field changes, with respect to mineral variations, is of great significance for the proper design of horizontal drilling and hydraulic fracturing. In this paper, clay-rich, calcareous and siliceous black shales are classified based on their XRD analysis, from the Wufeng and Longmaxi Formations in the Sichuan Basin of the Upper Yangtze Plate. Uniaxial quasi-static and cyclic loading tests, acoustic emissions, and CT scanning are combined to study mechanical behavior and fracturing mechanisms of the shale layer samples that own different mineral components. The results show that the strength and deformation response of the samples is a function of their mineral composition. The damage from cyclic loading is gradually accumulated in the samples which ultimately leads to fatigue and sample failure. The compressive strength of the samples under cyclic loading is lower than in the uniaxial quasi-static conditions. Furthermore, the fracturing of rock under uniaxial quasi-static loading mainly occurs after the peak strength is achieved while in cyclic loading, when a certain stress level is reached, fracturing begins, the damage accumulates up to the peak stress and the rock suddenly fails. Based on CT images from the samples after rupture, failure modes of the clay-rich, calcareous and siliceous black shales mainly constitute extensional and shear fractures, with some conjugate shearing. In comparison with uniaxial quasi-static loading, the fracture volume of the clay-rich, calcareous and siliceous black shales after cyclic loading increases by 12.2%, 25.9% and 32.9%, respectively. Overall, this study explains how cyclic loading that can be applied in the field via changing the injection pressure and other stimulation parameters, can effectively improve the complexity of the shale fracture network and lead to higher productivity of the formation. • Quasi static loading and cyclic loading are employed in the mechanical test. • The compressive strength of the cyclic loading is lower than that of uniaxial quasi-static loading. • The fracture volume of the cyclic loading is higher than that of uniaxial quasi-static loading. [ABSTRACT FROM AUTHOR]

The casing deformation problems during hydraulic fracturing process in shale gas wells seriously affect multi-stage hydraulic fracturing treatment and cause significant impact on efficient exploitation of shale gas in Sichuan, China. But by now, the mechanism of these casing deformation problems has not been completely understood. The tests of the casing deformation show that the shear failure is the main type of casing deformation. Combined with the detected seismic data, the fault slip caused by hydraulic fracturing is identified as the main reason for casing deformation in shale gas wells. A deep understanding of fault reactivation and slip is significant important. Hence, a semi-analytical model for calculating induced stress along the fault caused by hydraulic fracture is established in this paper. Based on the induced stress along the fault, semi-analytical models for calculating two types of fault slip are established, including the reactivated zone located at one fault tip and the reactivated zone located away from fault tips. Additionally, the effects of fluid pressure, dip angle of fault, scale of fluid stimulated area and the friction coefficient on fault slippage are discussed. At last, the method for calculating fault slippage proposed in this paper is applied in field and the calculated result is consistent with the field data. The semi-analytical model proposed in this paper provides a reasonable way for identifying casing deformation in shale gas wells or other wells which need hydraulic fracturing treatments. [ABSTRACT FROM AUTHOR]

• The desorption mechanism of oil from the nano-micro pores was studied. • The oil is desorbed when the particle self-rotation speed is faster than the critical speed. • When the inlet flow rate increased, the particle self-rotation speed increased, which could strengthen the desorption process. To realize the harmless treatment and resource recycling of oily shale cuttings, the most important step is to desorb the oil phase from the pores of the cuttings. The desorption treatment must overcome the capillary force and viscous resistance inside the pores. This paper discusses a cyclone oil desorption technology. The particles of the oily shale cuttings in the swirl field generated by the cyclone will have high-speed self-rotation, thus providing sufficient centrifugal forces that separate the phase oil from the pores. The critical speed of such self-rotation for the desorption of oil phase in pores of particles was analyzed using mechanical formula, and the distribution of particle self-rotation speed in the cyclone was calculated and simulated. We collected shale cutting samples from the Sichuan shale gas field in China. Experiments were conducted to observe the changes in the porosity and oil content of the cuttings before and after cyclone desorption treatment. The results show that the particle self-rotation in the cyclone can effectively achieve the oil phase desorption from the particle pores. For the shale cuttings in the experiments, when the particle self-rotation speed is faster than the critical self-rotation speed of 507 rad/s, the oil phase in the particle pores can be desorbed. When the inlet flow rate is between 55 m3/h − 95 m3/h, the particle self-rotation speed is faster than the required critical self-rotation speed. As the inlet flow rate increases, the particle self-rotation speed in the cyclone also increases, which generates better desorption effects. [ABSTRACT FROM AUTHOR]

Deep shale gas reservoirs, compared to shallow and intermediate shale gas reservoirs, usually exhibit lower fracture complexity, smaller stimulated volume, poorer conductivity, and faster production decline when the same injection scheme and construction parameters are used. As a result, specific injection scenarios for deep shale reservoirs must be investigated. In this paper, a series of laboratory fracturing tests were conducted to investigate the influence of injection scenarios on deep shale hydraulic fracture initiation and propagation. The black shale of the Lower Silurian Longmaxi formation in the Sichuan Basin was adopted, and four different injection schemes (constant fluid rate injection, cyclic injection, shut-in intermittent injection, and low-frequency pulse pressurization injection) were utilized in these experiments. The energy evolution during hydraulic fracturing with different injection schemes was analyzed by an array of acoustic emission (AE) sensors. Furthermore, the internal fracture geometry of the shale specimens after failure was determined using high-resolution CT scanning. Two parameters, fracture bulk density and fractal dimension were defined to quantitatively characterize the post-pressure fracture complexity. The results indicate that the rock breakdown pressure rises rapidly as the injection rate increases, for conventional monotonic rate injection. Although increasing the injection rate can serve the purpose of increasing the net pressure in the seam and enhancing the fracture complexity, it makes it easier to approach the equipment pressure limit. Unlike conventional hydraulic fracturing, which mainly induces tensile fractures near the wellbore, cyclic fracturing is more likely to cause fatigue damage to the rock, resulting in the formation of subcritical microcracks, which gradually grow into macroscale fractures. In addition, the cyclic injection method can reduce the rock breakdown pressure by approximately 24% compared with the conventional continuous injection method. The shut-in intermittent injection method can also reduce the breakdown pressure of shale by 22%, but this method is not conducive to enhancing the fracture complexity. Pulse fracturing creates the most complex fracture morphology among the four injection schemes tested. The rapid changes in wellbore pressure over a short time can easily result in the formation of multiple fractures around the borehole. Our experimental results and theoretical analysis can provide theoretical support for the optimization of injection schemes for deep shale hydraulic fracturing. • Laboratory tests were conducted to investigate the influence of injection scenarios on deep shale hydraulic fracturing. • All the experiments were monitored by an array of AE sensors. • The internal fractures geometry of the shale specimens after failure was determined by high-resolution CT scanning. • The cyclic fracturing can reduce the breakdown pressure by about 24% compared with the conventional continuous injection method. [ABSTRACT FROM AUTHOR]

The Longmaxi Formation has abundant shale gas resources. High-yield wells have been obtained in medium-shallow and deep areas in the Sichuan Basin. Quantitative research on the brittle-ductile change in rocks is a basic geological component for caprock evaluation, and it can also be applied to shale gas exploration. This paper fully considers the difference between the lower shale (N1) and upper shale (N2) of the Longmaxi gas-bearing shale. Through total organic carbon (TOC), mineral components, field emission scanning electron microscopy (FE-SEM), and triaxial mechanics experiments, the two groups of shale samples, N1 and N2, were compared and analyzed, and the brittle-ductile transition was quantitatively studied. Based on previous work, a method for the fast calculation of the key over consolidation ratio (OCR) was proposed. The results show that N1 is mainly organic-rich siliceous shale, N2 is primarily organic-poor clay shale and mixed shale, and N1 is generally more brittle than N2. The key OCRs of N1 and N2 are 2.16 and 2.60, respectively. The key confining pressures at which N1 and N2 become fully ductile are 87.87 MPa and 70.28 MPa, respectively. Therefore, the key depths of the brittle zones of N1 and N2 (H b1 and H b2) and the key depths of the ductile zones of N1 and N2 (H d1 and H d2) are obtained. On this basis, a comprehensive model for the Longmaxi shale gas' brittle-ductile transition exploration zone in the Sichuan Basin and H b1 and H b2 maps of the Longmaxi gas-bearing shale have been established. Through comprehensive vertical and horizontal analysis, Longmaxi shale gas exploration risks in different blocks of the Sichuan Basin have been clarified. Therefore, this study of the variation in the brittle-ductile transition has an important guiding role in refining the exploration of Longmaxi shale gas in the Sichuan Basin. • There are obvious differences between the lower shale and the upper shale. • Difference of brittle-ductile transition makes shale gas exploration more complicated. • The new method can quickly calculate the threshold value of over consolidation ratio. • The new model provides more comprehensive guidance for shale gas exploration. [ABSTRACT FROM AUTHOR]

The upper Permian Longtan Formation (LTF) shale is a new shale gas exploration and development focus in China, which was formed in a transitional depositional environment. Because of the difference in the sedimentary environment between marine-continental transitional shale and marine shale, the enrichment characteristics and mechanism of organic matter (OM) are very different. Because of a weak research foundation, the enrichment mechanism of OM in the LTF shale in the Sichuan Basin is still unclear. In this paper, shale samples of the upper Permian LTF from the Well F5 were collected, and the abundance, source and sedimentary environment of organic matter were studied by using paleoclimate, detrital input, redox, paleoproductivity and sedimentation rate parameters. The results show that the TOC content of the LTF shale is high (average = 3.04%) and reaches the high to overmature stage (average = 2.46%), indicating a favorable source rock. The LTF shale was deposited in a humid and hot climate, which was conducive to the growth and prosperity of a large number of plants. The source of organic matter was mainly terrestrial plant fragments. The sedimentary environment of the LTF shale was oxic, it had a high sedimentation rate, which reduced the oxidation rate of OM, promoted the rapid deposition of OM, and ensured the preservation of organic matter. In summary, the enrichment of OM of LTF shale in the Sichuan Basin was affected by the sedimentary environment, which was complicated. The main controlling factors were detrital input, redox conditions and sedimentation rate. • The geochemical characteristics of Permian Longtan Formation (PLTF) shale in eastern Sichuan are analyzed. • The sedimentary paleoenvironment of PLTF shale in eastern Sichuan is reconstructed. • The main controlling factors of organic matter enrichment of PLTF in eastern Sichuan are studied for the first time. [ABSTRACT FROM AUTHOR]

The lamellation fractures in marine shales are important to the accumulation and preservation of shale gas, whose characteristics and controlling factors are still unclear. This paper focus on the characteristics descriptions and controlling factor analyses based on core observations, thin sections, and scanning electron microscope (SEM) experiments. The marine shales of Paleozoic Wufeng Formation and Longmaxi Formation in the Fuling Area of Eastern Sichuan Basin are taken as an instance. The lamellation fractures are occurring along the direction of the lamellations, partially of them curved, bifurcated, and converged. The apertures of fractures typically range from 1 to 500 μm. Only a few (<1%) lamellation fractures are filled with quartz, calcite, pyrite, and bitumen. After a series of analyses, developments of lamellation fractures are controlled by total organic carbon (TOC), brittle minerals, laminae, and pyrite. (1) The organic matters control the formation and evolution of lamellation fractures by generating gas and resulting in overpressure. An increase in the TOC content of 1% can increase the lamellation fracture density by 0.26 cm-1. (2) Brittle minerals, especially authigenic quartz, is also an important controlling factor to the development of lamellation fractures. When BI increases by 10%, the lamellation fracture density increases by 19%. (3) Besides, siliceous and carbonate laminae in the shale are more suitable for lamellation fractures to develop. When the laminae density is smaller than 4 cm−1, it is easier to develop lamellation fractures with higher laminae density. However, too densely-distributed laminae will suppress the formation of lamellation fractures. (4) Higher pyrite content is good for the development of lamellation fractures. Lamellation fracture density will increase by 9%, with pyrite content increased by 1%. For the shale at the bottom of the Wufeng Formation - Longmaxi Formation, lamellation fractures develop extensively due to the high TOC and pyrite content, great brittleness, and moderate laminae density. The density of lamellation fractures decreases from the bottom to the top. • Quantitative characterize about development characteristics of the lamellation fractures. • The lamellation fractures is controlled by brittle mineral, TOC, laminae and pyrite. • TOC is the main one that affects the development of lamellation fractures. • Clarify vertical distribution characteristic of lamellation fractures in the Wufeng-longmaxi Fm. [ABSTRACT FROM AUTHOR]

The Lower Cambrian Qiongzhusi Formation source rock is one of the main sources of oil and gas in the Sichuan Basin. Owing to the controversial thermal history and the complex structural pattern in the Sichuan Basin, it is difficult to characterize the differences in maturity evolution and the timing of hydrocarbon generation periods of Qiongzhusi hydrocarbon source rocks in different tectonic zones. Based on the comprehensive analysis of the previous thermal history and denudation studies from the tectonic units of the basin, the thermal histories and denudation amount in main period of tectonic movement of a five-part tectonic development are obtained in this paper. Although the thermal history of different tectonic units varies, generally, the heat flow increased slowly in the early Paleozoic, and it experienced a peak heat flow at the end of the Middle Permian and decreased subsequently. Based on new thermal history and denudation data, the maturation of the Lower Cambrian Qiongzhusi Formation source rocks in the Sichuan Basin were modeled and their evolution can be classified into three patterns, which are: (1) continuous heating pattern in southern Sichuan, (2) Silurian-Permian heating cessation pattern in western Sichuan and northwest central Sichuan, and (3) Silurian-Permian and Triassic heating cessation pattern in the central uplift, northern structural zone and western part of the eastern structural zone. In addition, the key hydrocarbon generation periods of different tectonic units were determined in this study. Tectonic and sedimentary burial, together with high heat flow values during the Late Permian, were the main factors controlling levels of thermal maturity and thermal evolution models of the Qiongzhusi Formation. The present-day thermal maturity of the Qiongzhusi Formation depends on the thermal history and maximum burial depth before uplift in Late Cretaceous. However, the maturity evolution pattern affecting the key hydrocarbon generation stage was mainly related to the burial histories of source rocks in the Paleozoic. The thermal modelling of this study provides new evidence for the source rock evolution and new insight into Lower Paleozoic shale gas. • Thermal history, denudation amount and process of different tectonic units are summarized. • Maturity evolution and key HC generation periods of the Qiongzhusi Formation are obtained. • Three thermal evolution patterns of source rocks are classified. [ABSTRACT FROM AUTHOR]

The current fluid imbibition models used in shale usually assume a pore shape to be circular, whereas the pore shapes significantly deviate from that. One approach to correct this deviation is to introduce a pore shape factor into the models although it is challenging. This paper provides a way to calculate the specific shape factors of irregular inorganic and organic pores in shale based on SEM images analyses. Then, combining with the percentages of the organic and inorganic pores, a group of weighted-average pore shape factors of Longmaxi shale samples was calculated (0.587～0.667). Another group of pore shape factors (0.359～0.615) was also matched from the shale fluid imbibition curves, close to the previous values, which confirmed the validity of the proposed method. This work may provide an effective support in exploring the fluid imbibition or flow mechanisms in shale, especially considering the strong pore wall effects under nano-scale. • A method is provided to calculate the shape factors of irregular inorganic and organic pores. • The pore shape factor of organic pores is commonly higher than that of inorganic pores. • Ignoring the pore shape factor may cause an overestimation in fluid imbibition volume prediction. [ABSTRACT FROM AUTHOR]

The growing research on the features of shale pore structures has deepened the understanding of shale oil and gas storage capacity. Nevertheless, few studies have focused on the differences of the specific parameters that are characteristic of the nanopore structure in shales with different lithofacies. Such parameters include the full-scale ranges of the pore specific surface area and pore volume distribution, porosity and fractal dimension. In this paper, a dozen Longmaxi Formation shale samples from the southern Sichuan Basin were studied to determine the differences in these parameters. A set of laboratory experiments were performed, such as geochemical experiments, X-ray diffraction, N 2 and CO 2 isothermal adsorption, mercury injection capillary pressure tests and scanning electron microscopy experiments. The results of this study indicate that shale samples can be classified into 5 types of lithofacies based on their mineral composition. The surface porosity of organic matter particles is 8–16 times greater than that of mineral particles. The major contributors to the pore volume and surface area of shale samples are mesopores and micropores, respectively. The full-scale pore size distribution is mainly from 0.4 to 0.6 nm, 2–7 nm, and 10–20 nm for the pore volume and from 0.3 to 1.0 nm and 1.5–7 nm for the pore surface area. The average values of the porosity, fractal dimension, pore volume and pore surface area of the M-2 shales are almost equal to those of the S group shales. This result indicates that the M-2 shale has a large gas storage capacity similar to that of shales with siliceous lithofacies. The development of organic matter pores that contribute to more than 50% of the porosity leads to the high porosity and fractal dimension of the M-2 shales and S group shales. In addition to organic matter, illite minerals also play an important role in the constitution of the majority of the storage space in the CM-1 shale. • Full-scale pore structure of shales with different lithofacies were contrasted. • The smaller the pore size, the stronger the heterogeneity of shales. • Porosity minerals and OM accounted for in different shales were quantified. • Except for OM, illite is another important factor to pore development. [ABSTRACT FROM AUTHOR]

The Lower Paleozoic Ordovician Wufeng Formation and the Silurian Longmaxi Formation shale in the southeastern part of Sichuan Basin in China are the key strata for shale gas exploration. In this paper, based on a large number of outcrops, cores and natural gamma spectroscopy logs combined with scanning electron microscopy, geochemical parameter testing and mineral component analysis, the geological conditions and exploration potential of shale gas accumulation were studied in this area, and the relationship between structure and shale gas preservation is discussed. The vertical reduction in the U/Th ratio measured in the outcrop sandstone lenses and core horizontal bedding and observed in the energy spectrum logs indicates that the sedimentary water body transformed from a low-energy anoxic reducing environment to a high-energy oxidizing environment. The TI index and the kerogen carbon isotope indicate that the organic matter type of the shale is mainly type I kerogen. The TOC content of high-quality shale is generally above 2.50%, and the R o value is between 2.00% and 3.08%. High organic matter abundance, moderate thermal evolution and a large thickness of high-quality shale ensure an adequate supply of hydrocarbon resources. Multiple types of inorganic pores, primary organic pores and pores formed during hydrocarbon generation provide abundant reservoir space for the enrichment of shale gas. The tectonics have an important influence on the preservation of shale gas. Moderate fracture development, a late tectonic uplift time, small-scale tectonic uplift, and high formation pressure coefficient are conducive to the preservation of shale gas, and the gas content and productivity will be higher. Overall, a structural stability zone with a high formation pressure coefficient is a favorable area for shale gas exploration. • Comprehensive analysis the sedimentary environment based on three scales: outcrop, core and energy spectrum logging. • Later and small-scale tectonic uplift, moderate fracture development, good roof preservation are conducive to gas enrichment. • A structural stability zone with a high formation pressure coefficient is a favorable area for shale gas exploration. [ABSTRACT FROM AUTHOR]

The Upper Ordovician Wufeng Formation and Lower Silurian Longmaxi Formation marine shale are the major shale gas exploration and development target in the Middle-Upper Yangtze region. Although three national shale gas demonstration areas have been built and several high-yield shale gas wells have been discovered in the Wufeng Formation and Longmaxi Formation, shale gas in these areas has not been fully commercialized because the enrichment degree of shale gas varies greatly in different tectonic settings. By analyzing the formation processes and distribution characteristics of the organic-rich shale, reconstructing the burial history and hydrocarbon generation history, and comparing the shale gas enrichment and preservation conditions in different tectonic settings in South China, this paper points out that (1) The craton depression and foreland basin on the cratonic margin are favorable zones for the development of organic-rich shale. Because the Yangtze plate is small, the foreland basin is greatly affected by the input of terrigenous clastic and the thickness and area of organic-rich shale in the craton basin is larger than that in the foreland basin. The organic-rich shale is mainly distributed in south Sichuan deep-water shelf and west Hubei-east Chongqing deep-water shelf, and the enrichment interval is the Dillcellograptus complexus-Cystograptus Vesicolus graptolite zone from the Wufeng Formation to the bottom of Longmaxi Formation are favorable intervals for shale gas enrichment. (2) Since the Late Jurassic, significant differential intracontinental tectonic deformation, uplift and denudation occurred in the Middle-Upper Yangtze region. Moderate compression and uplift were beneficial to the development of micro-fractures and bedding fractures in shale, resulting in good preservation conditions of shale gas reservoirs with a high-pressure coefficient. In contrast, strongly deformed areas, where faults and high-angle large fractures were developed, are not conducive to the preservation of shale gas reservoirs. (3) The craton depression of the Middle-Upper Yangtze region, where organic-rich shale was deposited, is favorable to shale gas enrichment because it provides a good material basis, and also a relatively stable and compression-resistant tectonic environment for later tectonic deformation and uplift-derived destruction. Based on the analysis and evaluation of tectonic construction–transformation conditions under the macro-control of plate tectonics, it can well explain the existing commercial shale gas discoveries and failure cases and can guide the optimization of future exploration targets. • ►Deep water in the craton depression controlled the distribution of organic-rich shale. • Stable tectonic settings are favorable for shale gas enrichments. • Two depressions in Yangtze Craton are the most favorable areas for shale gas accumulation. [ABSTRACT FROM AUTHOR]

Numerous studies have been conducted on the isotopic geochemical characteristics of shale gas. However, a systematic study of the relationship between the isotopic geochemical characteristics and gas content has not been addressed to date. In this study, we investigated the geochemical and isotopic compositions of shale gas samples, and the gas content of the Lower Silurian Longmaxi marine shale in the Huangjinba area, southern Sichuan Basin, China. In addition, we analyzed the quantitative relationship between the geochemical parameters and the gas content. The Longmaxi shale gas is a dry gas with a low wetness (0.12%–0.67%), dominated by methane (89.49%–99.00%) with a small amount of ethane (0.11%–0.61%). The average values of δ13C methane , δ13C ethane and δ13C propane are −29.8‰, −35.0‰, and −33.2‰, characterized by a partial isotopic reversal (δ13C methane > δ13C propane > δ13C ethane) with a maturity level (Ro ≈ 2.30%–2.51%). The total gas content of the Longmaxi shale varies widely from 0.42 m3/t to 3.88 m3/t. Correlation analyses revealed that the ethane isotope (δ13C ethane) and the magnitude of isotopic reversal (δ13C methane - δ13C ethane) exhibit strong correlations with the total gas content and the free gas content. This study demonstrates that carbon isotopic reversal of shale gas largely depends on the secondary cracking, and thermal maturity and residual liquid hydrocarbons content are important factors affecting the magnitude of carbon isotopic reversal. Based on the quantitative study conducted in this paper, we provide a new proxy for predicting the gas content of high maturity shale using an isotopic method. • Carbon isotopic reversal of shale gas largely depends on the secondary cracking. • Sufficient residual liquid hydrocarbon is necessary for the carbon isotopic reversal. • Wetness of 1.4% in high maturity shale gas indicates the peak of secondary cracking. • δ13C methane –δ13C ethane can be used to predict the gas content of high maturity shale. [ABSTRACT FROM AUTHOR]

The initial water saturation of shale reservoirs is in general lower than the irreducible water saturation. However, the hydraulic fracture network region is highly saturated by fracturing fluid, which results in a two-phase flow inside fractured region. A long-time two-phase flow period is always observed in the fields even after a high flowback rate, which cannot be explained by existing models. In this paper, we develop a 2D gas-water flow model coupling multi-scale discrete fractures and continuum framework to describe the flowback process of a multi-stage fractured horizontal well in shale gas reservoirs. The model considers the non-linear flow mechanism of free gas, desorption of adsorbed gas and the co-existence of single gas flow and gas-water flow. The capillary pressure model and relative permeability model are both employed with different parameter values for matrix and fractures, respectively. The permeability of naturally fractured shale is described using a fractal model. To overcome non-convergence problem of finite element method for solving two-phase flow in cross fractures, the model is solved by finite element method and finite volume method. The model is verified through a simple two-phase flow simulation and a single gas phase simulation with some simplifications. Then the model is used for history matching as well as production prediction of a well in Changning block, Sichuan, China. The simulation results show that the high water saturation induced "skin factor" in the vicinity of hydraulic fractures significantly delays the gas production for a long time because some of the fracturing fluids become irreducible water and lower the gas phase permeability. Due to the great influence of fracture interference on water flow, the fracturing fluids concentrate on regions near hydraulic fractures and their junctions. The long-term gas production of the complex fracture network is slightly higher than that of the simple fracture network, which is different with the simulation results of single gas flow model. • A physical model where the gas-water coexists with single gas flow is built. • The viscous flow, Knudsen diffusion and adsorption / desorption phenomenon are considered of the gas flow in shale matrix. • The finite element method (FEM) and the finite volume method (FVM) are employed to solve the mathematical model. • The flowback and production process can be simulated continuously. [ABSTRACT FROM AUTHOR]

• Shale permeability upscaling is proposed considering multiple transport mechanisms. • A REV size is determined as 420 μm with experimental data of Sichuan basin. • Bulk gas in inorganic matter contributes roughly 70% to core-scale permeability. • Water film effect can be neglected at mean pore size larger than 40 nm. This paper presents a new upscaling method to derive the core-scale apparent gas permeability from an improved pore-scale permeability model and experimental data, with more rigorous incorporation of varying gas storage/transport mechanisms in nano/micro pores. First, in use of SEM images of a gas-rich shale field example in Sichuan Basin from our lab, pore network models of inorganic-matter (IOM) and organic-matter (OM) are characterized by using a digital-core technique. Next, an improved pore-scale real gas apparent permeability is modeled rigorously for both IOM/OM, respectively, with 1) bulk gas transport, gas adsorption, surface diffusion, pore-size confined phase behavior, and stress-dependent rock properties and 2) an additional reduction in inorganic pore sizes by water film adhered on pore surfaces. Core-scale permeability is then derived by assembling the permeabilities of stochastically distributed IOM/OM patches with different pore network models properties using the Monte Carlo sampling method. The new core-scale permeability model is validated by pulse-decay permeability experiment. Moreover, the representative elementary volume (REV) size is determined by analyzing the relative standard deviation of apparent gas permeability in cases with different sample sizes. The contributions of different gas transport mechanisms are discussed, and the impacts of stress-dependence for several field examples (i.e., Sichuan, Pierre and Barnett Basins) and water film with varying relative humidity (RH) on core-scale apparent permeability are analyzed. This work provides an effective approach to determine the core-scale shale permeability by directly using pore-scale experimental data, which is a common challenge in the unconventional resources. [ABSTRACT FROM AUTHOR]