Abstract

Natural gas hydrate is a kind of environmentally friendly unconventional energy with large reserves, which attract increasing attention in recent years. The microscopic pore structure and mineral composition of natural gas hydrate reservoir have a great influence on its seepage characteristics during exploitation. In this study, computed tomography and scanning electron microscope are used to obtain digital images of three clayey-silt natural gas hydrate reservoir samples in the Shenhu area of South China Sea, and then, the pore structure and the mineral composition of the samples are obtained after image processing. The result indicates that the clayey-silt samples show strong hydrophilic characteristics, small particles, good sorting properties, variable pore distribution, small average pore and throat radius, large porosity, and a large content of submicron pores. The mineral composition of the reservoir is mainly quartz and contains a certain amount of clay minerals, mainly illite, while the mineral composition of reservoir rich in microbial fossils is mainly carbonate. The results indicate that the pore structure and mineral composition of natural gas hydrate reservoirs in Shenhu area lay a foundation for the efficient development of the clayey-silt reservoir.

1. Introduction

Natural gas hydrate (NGH) is a kind of clean unconventional energy formed of water and natural gas under low temperature and high pressure [1, 2]. It is generally distributed in submarine sedimentary strata, land permafrost, and some deep-water lake sediments, with large reserves and shallow burial depth [36]. Many field production tests have been tried to exploit such resources from onshore and offshore natural gas hydrate deposits worldwide, including West Siberia [7], Alaska North Slope [8], Krishna-Godavari Basin [9], and Nankai Through [10]. In the Shenhu area of South China Sea, two production tests were carried out for gas hydrate reservoirs in 2017 and 2020 [1113], respectively. During the development process, the gas production efficiency of the South China Sea hydrate reservoir fluctuates greatly, and the main controlling factors for the decline of the gas production efficiency in the later stage of development are still unclear. The reservoir samples and data obtained during the exploitation provided the basis for subsequent research on detailed seepage and development issues [11, 1316]. Clarifying the complex mineral composition and pore structure of reservoir sediments is one of the key scientific issues to improve the recovery efficiency of hydrate reservoirs [1, 9, 17].

The lithology of natural gas hydrate reservoirs in the Shenhu area is clayey silt with complex mineral composition. The overall mineral sources of the reservoir vary greatly and contain a certain amount of clay minerals [18, 19]. Meanwhile, the reservoir pore structure is complex, the sources of pores are diverse, and the size of pores also spans multiple scales. Effective characterization of pore structure is helpful to clarify the reservoir productivity [20]. Existing reservoir mineral and pore information is only obtained indirectly through geophysical logging, and there is a lack of direct characterization of gas hydrate reservoir sediments in the South China Sea using laboratory experiments. Digital rock technology, as a method of characterizing reservoir samples graphically, has been widely used in unconventional reservoirs such as shale gas and tight oil [2124]. The scanning electron microscope (SEM) and microcomputed tomography (micro-CT) techniques are usually used to obtain high-resolution images for digital rock construction [2527]. The mineral and pore structure features of in different rock samples can be visualized and quantified [2831]. In terms of natural hydrate reservoirs, some works have also applied digital rock technology to study the hydrate synthesis-decomposition process [3234], hydrate occurrence [35, 36], and seepage characteristics [3740]. However, there is a lack of precise characterization of the mineral composition and micro- and nanoscale complex pore structure of gas hydrate reservoirs in the South China Sea by directly using high-precision imaging instruments, such as CT scanning and SEM imaging.

In this study, the micro-CT is used to obtain the three-dimensional grayscale images of the hydrate reservoir samples, and then, the pore and throat distribution are analyzed based on the processed binary image. Subsequently, two-dimensional nanoscale pore structure and mineral composition images are obtained by SEM, and the submicron pore characteristics and mineral composition of the sample are quantitatively analyzed.

2. Samples and Methods

The three NGH reservoir samples studied come from the Shenhu area. The location map and the strata of the reservoir are detailed in the study by [11, 39]. The lithology of the samples is mainly clayey silt with a median particle size of 12 μm [11]. The wettability angles are determined by sessile drop method which measures the angle between the bubble surface and the core surface using the KRÜSS DSA100E apparatus. As shown in Figure 1, the average wettability angles of three samples on both sides show strong hydrophilic characteristics (53.9°, 46.5°, and 48.3°).

2.1. Micro-CT

CT scans are performed on three NGH samples using a Zeiss MicroSCT-200micron CT scanner by iRock Technologies. The voltage and power used for scanning are 150 kV and 10 W, respectively. The pixel number and resolution are and ~0.6 μm, respectively.

2.2. Preparation for Electron Microscope

Considering the lithology of the hydrate reservoir, dispersion happens when it contacts with water; meanwhile, the clay cementation and shrinkage will produce microcracks after drying. Therefore, targeted measures should be taken during the preparation of electron microscope samples to keep the drying temperature at 60°C to prevent the clay from losing water. During the glue injection process, using double injection is applied on both sides so that the epoxy resin can better fill the pores. The dry grinding method is used to polish samples to prevent particles from falling off after being exposed to water. The specific electron microscope sample preparation process is shown in the flow chart of Figure 2. The diameter of the accomplished round thin sample is 20 mm.

2.3. MAPS

Due to the limitation of the instrument itself, the resolution of the CT scan is up to the submicron level, resulting in the low accuracy in imaging nanolevel pores in the clay. Certain limitations in mineral identification also arise, viz., accurately identifying mineral types remain unattainable. Therefore, it is necessary to further supplement the two SEM methods, including modular automated processing system (MAPS) as well as the quantitative evaluation of minerals by scanning electron microscopy (QEMSCAN), to obtain the micro-nano pores and accurate mineral composition of hydrate reservoirs.

MAPS is aimed at samples requiring large-area observation (containing multiscale pore structures) [41]. A series of continuous and overlapping edges in a large number of high-resolution samples in a selected area is scanned. These small images will be stitched together to obtain a super high-resolution and large-area two-dimensional backscattered electron image. The MAPS image can be zoomed in or out arbitrarily like Google MAPS to observe multiscale rock characteristics. The test instrument used in this study is Helios NanoLab 650 of iRock Technologies.

2.4. QEMSCAN

QEMSCAN is a comprehensive automatic mineral detection method, which scans the surface of a sample with a high-energy electron beam accelerated along with a preset raster scan mode and produces a color map of the mineral assemblage’s distribution characteristics [42]. The instrument emits an X-ray energy spectrum and provides information on the element content at each measurement point. Through the combination of backscattered electron grayscale image and X-ray intensity, the content of the element can be obtained and converted into a mineral phase. The test instrument used in this study is QEMSCAN 650F of iRock Technologies.

3. Results and Discussion

3.1. Micropore Structure Distribution

In this section, the microscale pore structures of samples are obtained by CT scan. Figure 3(a) shows the original CT scan grayscale image of sample 1, and small particles are visible with good sorting properties. The presence of large particles is rare, while it is rich in clay and debris. Most of the pores are developed in clay and microorganisms such as diatoms. Figures 3(b), 3(c), and 3(d) show the pore structure images obtained from the original gray image after post-processing (intercepting, filtering, and binarization) [43]. The corresponding CT porosity can also be extracted from the binarization pores (Table 1). It is more distinct that only a few connected macropores appeare in sample 1 from the binarization 3D image, yet most pores are isolated micropores. The CT porosity is 17.3%, and the connected porosity is 13.7%.

Based on pore structure space shown in Figure 3, Figure 4 shows pore and throat distribution curves obtained from nuclear magnetic resonance (NMR) and mercury injection simulation [4447], respectively. The pore radius of sample 1 ranges from 0 μm to 12 μm, presenting a single-peak distribution with a maximum peak at 1 μm and an average pore radius of 1.29 μm (Table 1). The throat radius is in the range of 0~2.9 μm, and the distribution characteristic of the throat radius is multipeak. The maximum peak is located at 1.2 μm, and the average throat radius is 0.94 μm (Table 1).

Figure 3(e) shows the original CT grayscale image of sample 2. Compared with sample 1, sample 2 has the same small particles and good storability. Large particles are rarely distributed, while more amounts of clay and debris are presented. However, the pores are mostly elongated pores with good pore connectivity and supplemented by macropores of microorganisms’ origin. Figures 3(f), 3(g), and 3(h) show the pore structure image obtained from the original gray image after post-processing (intercepting, filtering, and binarization) [43] of sample 2. The binarization 3D image shows it is obvious that more large pores in sample 2 with better connectivity, and fewer isolated pores are observed. The CT porosity is 19.0% and the connected porosity is 18.2%.

The pore radius of sample 2 ranged from 0.4 μm to 2.5 μm, presenting a single-peak distribution with a maximum peak at 1.2 μm (Figure 4) and an average pore radius of 1.15 μm (Table 1). The throat radius is in the range of 0.3~2.5 μm, and the distribution characteristic of the throat radius is multipeak. The maximum peak is located at 1.2 μm, and the average throat radius is 1.04 μm (Table 1).

Figure 3(i) is the original CT grayscale image of sample 3. Compared with samples 1 and 2, sample 3 also shows better sorting with smaller particles while the rare occurrence of large particles. However, the pores mainly developed into fractures and round pores with good connectivity, and there are some microorganism-derived pores. Figures 3(j), 3(k), and 3(l) show the pore structure image obtained from the original gray image after post-processing (intercepting, filtering, and binarization) [43] of sample 3. It is even much clearer that the fractured macropores and circular macropores with good connectivity are in the majority in sample 3, from the binarization 3D image, while isolated pores still account for a certain part. The CT porosity is 18.7%, and the connected porosity is 16.5%.

The pore radius of sample 3 ranges from 0 μm to 12.5 μm, presenting a bimodal distribution corresponding to the bimodal pore radius values of 1 and 6.5 μm (Figure 4), respectively, and an average pore radius of 3.04 μm (Table 1). The throat radius is in the range of 0.3~19 μm, and the distribution characteristic of the throat radius is multipeak. The maximum peak is located at 0.8 μm, and the average throat radius is 1.19 μm (Table 1).

In summary, under the CT scan images of micron level, the pore and throat distributions of three samples are basic distribution within 19 μm, the average radius of pore and throat is between 1 and 3 μm, and porosities are in the range of 17~19%.Namely, although the overall pore size of the hydrate sample is small, a large number of pores still constitute a considerable pore space.

3.2. Micro- and Nanopore Structure Distribution

In this section, the micro-nano scale pore structures of samples are obtained by MAPS. Figure 5(a) is the electron microscope grayscale image of sample 1 obtained by MAPS after sample preparation. A total of 1,230 single scanning images are combined into an overall image, with a pixel number of 8 and a field of vision of  mm. The pixel number of the overall image is , and the field of view is  mm. The pixel resolution of the image is 195 nm.

To eliminate the large pores generated in the sample preparation process and improve the calculation efficiency, the original grayscale image of sample 1 is first intercepted, and then, the binary segmentation is employed to divide the pores and the matrix, obtaining a segmented image of the binarized two-dimensional pore distribution of sample 1 (Figure 5(b)). The corresponding porosity can also be obtained from the binarized pore image, which is 16.78%. Figure 6 is the two-dimensional pore size distribution of sample 1 calculated from Figure 5(b). The simulation calculation results show that the sample pores have an average radius of 1.543 μm and a median radius of 1.341 μm, and the pores are distributed in a normal distribution.

Figure 5(d) is the electron microscope grayscale image of sample 2 obtained by MAPS after sample preparation. A total of 1,152 single scanning images are combined into an overall image, with a pixel number of and a field of vision of  mm. The pixel number of the overall image is , and the field of view is  mm. The pixel resolution of the image is 195 nm.

To eliminate the large pores generated in the sample preparation process and improve the calculation efficiency, the original grayscale image of sample 2 is first intercepted, and then, the binary segmentation is employed to divide the pores and the matrix, obtaining a segmented image of the binarized two-dimensional pore distribution of sample 1 (Figure 5(e)). The corresponding porosity can also be obtained from the binarized pore image, which is 24.14%. Figure 6 is the two-dimensional pore size distribution of sample 2 calculated from Figure 5(e). The simulation calculation results show that the sample pores have an average radius of 4.1 μm and a median radius of 1.23 μm, and the pores are distributed in a normal distribution.

Figure 5(g) is the electron microscope grayscale image of sample 3 obtained by MAPS after sample preparation. A total of 1,023 single scanning images are combined into an overall image, with a pixel number of and a field of vision of  mm. The pixel number of the overall image is and the field of view is  mm. The pixel resolution of the image is 195 nm.

To eliminate the large pores generated in the sample preparation process and improve the calculation efficiency, the original grayscale image of sample 3 is first intercepted, and then, the binary segmentation is employed to divide the pores and the matrix, obtaining a segmented image of the binarized two-dimensional pore distribution of sample 1 (Figure 5(h)). The corresponding porosity can also be obtained from the binarized pore image, which is 20.77%. Figure 6 is the two-dimensional pore size distribution of sample 3 calculated from Figure 5(h). The simulation calculation results show that the sample pores have an average radius of 2.86 μm and a median radius of 0.5 μm, and the pores are distributed in a normal distribution.

3.3. Mineral Composition

In this section, the mineral compositions of samples are obtained by QEMSCAN. Figure 5(c) shows the QEMSCAN result of sample 1, with a scan accuracy of 5 μm and a scan area of  mm. The mineral of sample 1 is mainly quartz, whose proportion is less than 50%; the proportion of feldspar is 12.16%; the proportion of carbonate is 15.37%; and the content of clay is close to 20%. Clay is mainly composed of illite (72.65%) and contains a certain amount of mica (6.98%).

Figure 5(f) shows the QEMSCAN result of sample 2, with a scan accuracy of 5 μm and a scan area of  mm. The mineral of sample 2 is mainly quartz, whose proportion is less than 35%, the proportion of feldspar is 9.89%, and the proportion of carbonate is 13.74%, and high in clay, which is close to 30%. Clay is mainly composed of illite (72.65%) and contains a certain amount of mica (6.98%).

Figure 5(i) shows the QEMSCAN result of sample 3, with a scan accuracy of 5 μm and a scan area of  mm. The mineral of sample 3 is mainly calcite, whose proportion is 69%; the proportion of quartz is less than 16%; and the proportion of feldspar is 5%, and low in clay, which is close to 5.1%. Clay is mainly composed of illite (76.02%).

3.4. Typical Minerals and Pores

Observing the MAPS images at different scales and comparing the mineral analysis results of the three samples, the particles of quartz and feldspar in the reservoir are relatively large and poor in roundness. The foraminifer fossils are mainly calcite. Clay particles are small and fill between quartz and feldspar, mainly illite (Figures 7 and 8).

Observing the MAPS images of three samples at different scales, the pore types of clayey-silt reservoirs can be divided into microbial fossil pore, intercrystalline pore, intergranular pore, and dissolution pore (Figure 9). Among them, clay pores and microbial fossil pores take up a relatively large proportion of the total pores. There are many microbial fossils in the hydrate reservoir, most of which are foraminifera, which can be used as a typical feature of the reservoir (Figure 10).

3.5. Comparison of Minerals and Pores

A comprehensive comparison of the mineral of the three samples shows that the mineral content and distribution of samples are quite different (Figure 11). Sample 1 has similar mineral properties to sample 2, with relatively high quartz content and a certain amount of clay. The minerals of sample 3 are mainly carbonate, which should be related to the fossil minerals.

From the comparison of the two-dimensional pore distributions of the three samples, the pore distributions of different samples are also quite different (Figure 6). Under the nanolevel SEM image, all three samples have a large number of submicron pores. Submicron pores of samples 1 and 2 account for 38.9% and 45%, respectively. As the reservoir minerals of sample 3 are mainly carbonate, there are many submicron pores, accounting for 75% (Table 2).

4. Conclusions

In this study, micro-CT and nano-SEM are applied to obtain digital images of clayey-silt hydrate reservoir samples in the South China Sea, and the pore structure characteristics and mineral compositions of the reservoir samples analyzed based on the images. The following conclusions can be derived: (1)Hydrate samples show strong hydrophilic characteristics. Under the resolution of micro-CT, the clayey-silt reservoir samples have small particles and good sorting properties and are rich in clay and debris contents. The pore morphology characteristics are variable, and the pore types are mainly clay and microorganisms’ pores. The pore and throat radius are distributed below 19 μm, but the porosity is large, indicating that a large number of small pores constitute a relatively considerable pore space(2)Under the resolution of SEM, the pore types of clayey-silt reservoirs can be divided into microbial fossil pore, intercrystalline pore, intergranular pore, and dissolution pore. Among them, clay pores and microbial fossil pores account for a relatively large proportion of the total pores. The pore distributions of different samples vary greatly, but they all have submicron pores with larger content(3)There are obvious differences in the mineral compositions of the three clayey-silt samples. Samples 1 and 2 are mainly composed of quartz and contain a certain amount of clay, which is mainly illite. The mineral composition of sample 3 is mainly carbonate, which is speculated to be caused by fossil minerals

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The author(s) declare(s) that they have no conflicts of interest.

Acknowledgments

The authors are grateful to the National Natural Science Foundation of China (No. 51991365); the China Geological Survey Project (No. DD20211350); the Key Program of Marine Economy Development (Six Marine Industries) Special Foundation of Department of Natural Resources of Guangdong Province [2021]56; and the Guangdong Major Project of Basic and Applied Basic Research (No. 2020B0301030003).