Pore network extraction, which is defined as the transformation from irregular pore space to a simplified network in the form of pores connected by throats, is significant to microstructure analysis and network modeling. A physically realistic pore network is not only a representation of the pore space in the sense of topology and morphology, but also a good tool for predicting transport properties accurately. We present a method to extract pore network by employing the centrally located medial axis to guide the construction of maximal‐balls‐like skeleton where the pores and throats are defined and parameterized. To validate our method, various rock samples including sand pack, sandstones, and carbonates were used to extract pore networks. The pore structures were compared quantitatively with the structures extracted by medial axis method or maximal ball method. The predicted absolute permeability and formation factor were verified against the theoretical solutions obtained by lattice Boltzmann method and finite volume method, respectively. The two‐phase flow was simulated through the networks extracted from homogeneous sandstones, and the generated relative permeability curves were compared with the data obtained from experimental method and other numerical models. The results show that the accuracy of our network is higher than that of other networks for predicting transport properties, so the presented method is more reliable for extracting physically realistic pore network.

Characterization of shale cores with low and anisotropic permeability is complicated, due to the presence of multiscale pore structure and thin layers, and defies conventional methods. To accurately reproduce the morphology of multiscale pore structure of the shale core, a novel core-scale reconstructing method is proposed to reconstruct 3D digital-experimental models by means of the combination of SEM, EDS images, nitrogen adsorption and pressure pulse decay experiment result. In this method, the multiscale and multicomponent reconstructing algorithm is introduced to build the representative multiscale model for each layer, which can describe the complex 3D structures of nano organic pores, micro-nano inorganic pores, micro slits and several typical minerals. Especially, to reproduce the realistic morphology for shale, the optimization algorithm based on simulated annealing algorithm uses the experimental data as constrain conditions to adjust and optimize the model for each layer. To describe the bedding characteristics of the shale core, bedding fractures are constructed by analysis of the mineral distribution in the interface of two layers, and then the representative models for different layers are integrated together to obtain the final core-scale digital-experimental model. Finally, the model is validated by computing its morphological and flow properties and comparing them with those of the actual 3D shale sample. This method provide a way for systematically and continuously describe the multiscale and anisotropic pore structure (from nm-cm) of the shale core, and will be helpful for understanding the quality of the shale reservoir.

Organic-rich shales have abundant micronanopores and slits. Most micronanopores exist in organic matter and pyrites. Pore size distribution significantly influences shale gas flow. In micronanopores of organic matter and pyrites, pore size distribution has a self-imitating or fractal property. This property is helpful in upscaling calculation. In this work, intermingled fractal model for organic-rich shale is investigated. First, on the basis of scanning electron microscopy (SEM) images of shale samples in Longmaxi shale formation, fractal units are built-in representative elementary surfaces of organic matter and pyrites. Second, intermingled fractal models (IFMs) are then built considering the organic matter fractal pore units, the pyrite fractal pore units, and slits. Influences of fractal parameters on apparent gas permeability of IFMs are investigated. Third, an upscaled shale model is successfully built using the mirror image method. This upscaling method is verified by comparing the apparent gas permeability between upscaled models and experimental results. The upscaling method based on fractal theory can significantly shorten the calculation time and provide a promising way for apparent gas permeability calculation in the upscaled model.

Food security is a goal and means of global sustainable development, and an important component of China’s national security. Based on grain production data from 2000 to 2020, 31 provinces (cities, autonomous regions) in China were used as research units to analyze the spatiotemporal differences and driving forces of grain production in China using a combination of local correlation index, center of gravity transfer model, and geographic detector. The results as follows, ① During the research period, China’s total grain production showed a decrease followed by an increase, and the yield per unit area of grain showed an increasing trend. Corn has become the "largest staple food" in China; ② During the research period, the focus of China’s grain production continued to shift northward, with Heilongjiang, Henan and Shandong provinces becoming the main grain production areas, with Henan being the province where China’s grain production center was located. Among the factors affecting grain yield, the effective irrigation area (0.971) has the strongest explanatory power. Finally, countermeasures and suggestions were proposed from five aspects, stabilizing grain production, reducing grain inventory pressure, implementing regional grain security responsibilities, improving grain circulation efficiency, promoting high-quality grain engineering construction, adjusting grain production structure, strictly implementing farmland protection responsibilities, scientifically applying pesticides, fertilizers and other production materials, emphasizing agricultural infrastructure construction and stabilizing the number of agricultural labor force.

In this study, shale cores from 20 wells in the S1l¹_1-1 sublayer of Longmaxi Formation buried in shallow shale (3500 m) in the southern Sichuan Basin, China were collected to compare their pore structures and gas-bearing properties using multiple experiments. Results showed that the deep layer has relatively lower brittle mineral content, which is disadvantageous in terms of the higher requirements it imposes on hydraulic fracturing. Results also showed that the most important factor controlling the differential enrichment of S1l¹_1-1 shale gas in southern Sichuan Basin is porosity. Moreover, the porosity composition of shallow shale and deep shale has significant differences: the porosity of shallow shale is dominated by organic pores, while for deep shale, both organic and inorganic pores are important. The inorganic pores provide significant storage space for free gas in deep shale; their contribution warrants more attention. We also found that the difference in organic porosity of the shallow and deep shale samples resulted from large differences in pore development ability, while the highest inorganic porosity was concentrated near the optimal mineral composition when the content of quartz plus feldspar plus pyrite was about 70%. This study revealed the primary factor controlling the difference in gas content between shallow and deep shale and detailed the characteristics of microscopic pore structure, providing a basis for the exploration and development of deep shale gas in the Wufeng-Longmaxi Formation in the southern Sichuan Basin.

Low-temperature nitrogen adsorption is a widely used method for the research and evaluation of gas shale’s pore structure. The existing interpretation method, utilizing gas adsorption isotherms to obtain pore size distribution (PSD), is always based on the one-dimensional geometry model, while the void space of gas shale has strong multi-dimensional characteristics. It is necessary to investigate the nitrogen condensation and evaporation behavior in multidimensional structures. In this study, a series of two-dimensional and three-dimensional models based on ink-bottle pores were constructed. A hybrid molecular simulation approach combining grand canonical Monte Carlo (GCMC) and molecular dynamics (MD) is proposed to simulate the low-temperature nitrogen adsorption isotherms. Three aspects have been analyzed in detail. Compared with the conventional understanding that the threshold of cavitation in the ink-bottle pore only relates to throat diameter, this study discloses a wider and more comprehensive range of conditions of cavitation that considers both the throat length and diameter. As pore spaces of shale samples consist of many complex interconnected pores, the multi-stage ink-bottle pore model is more suitable than the single ink-bottle pore model to similarly reproduce the wider cavitation pressure range. A more universal parameter is proposed that quantitatively unifies the influence of cavity diameter and length on condensation pressure and has good applicability in cavities with different shapes. This work quantitatively studies the nitrogen adsorption isotherms of three-dimensional complex nanopore structures using molecular simulation and provides a reasonable explanation for the low-temperature nitrogen adsorption isotherms of gas shale.