Abstract

Pore structures with rich nanopores and permeability in tight gas reservoirs are poorly understood up to date. Advanced techniques are needed to be employed to accurately characterize pore structures, especially tiny pores which include micron and nanopores. In this study, various experimental techniques such as scanning electron microscopy (SEM), nuclear magnetic resonance (NMR) , nitrogen adsorption method, and NMR cryoporometry (NMRC) are combined to interrogate the complex pore systems of the tight gas reservoir in the Linxing formation, Ordos Basin, China. Results show that tight gas sandstones are primarily comprised of residual interparticle and clay-dominated pores. Clay and quartz are two dominate minerals while pyrite occupies a nontrivial amount as well. The permeability of tight gas sandstones is very low, exhibiting an extremely poor positive correlation with porosity. While pore types and relative pore contents are more influential factors on the permeability, accurate characterization of pore size distribution is critical for the permeability of tight gas sandstones. Therefore, complementary characterization methods are carried out, indicating that neither small pores with (around peak 1 in NMR distribution) nor large pores with (around peak 3 in NMR distribution) control the permeability by analyzing the connectivity of the pores in various size ranges, but rather pores averaging approximately (around peak 2 in NMR distribution) have sufficient connectivity to host and transmit hydrocarbons. The pore size of tight gas sandstones is dominated by the clay-rich mineral assemblage. The study shows that the NMRC technique can be a very promising method, especially when referred to as a promising “roadmap” on how to interrogate tight formations such as the tight gas sands or even shale especially for the nanopore characterization.

1. Introduction

With the gradual depletion of conventional oil and gas systems, unconventional reservoirs have become the focus of the energy market. Of these, a tight gas sandstone reservoir is an important member [1, 2]. The definition of a tight gas reservoir is “a reservoir that cannot be produced at economical flow rates nor recovery economic volumes of natural gas, unless the well is stimulated by a large hydraulic fracturing or produced by the use of a horizontal wellbore or multilateral wellbores” [3, 4]. In addition, low permeability reservoirs developed in the past are mostly sandstones.

Whether a gas reservoir qualifies as tight or not depends on both physical and economic factors, as with tight oil sandstones, tight gas is associated with low porosity (<10%) and commonly extremely low permeability (<0.1 mD) reservoirs produce mainly dry natural gas—i.e., methane [57]. Tight gas sandstones are characterized by a wide pore-size range, poor connectivity, and significant heterogeneity [5, 8]. Four types of pores observed by SEM include residual interparticle pores, grain dissolution pores, intercrystalline pores, and small microcracks [9]. Zhu et al. concluded that the large pores in the sandstone samples are more complex, which is also one of the primary causes of the sample heterogeneity [2]. Xiao et al. classified pore size of tight gas reservoirs into three scales: clay-associated pores and intraparticle dissolution pores are less than 0.5 μm, the narrow slits between grains and the quartz intercrystalline pores range 0.5-1.5 μm, and the interparticle-related pores are larger than 1.5 μm [10]. Lai et al. suggested that micropores and mesopores provide most of the pore volume, whose pore size is less than 50 nm [11]. According to Shi et al.’s results, almost half of the pore space is occupied by the clay and capillary bound water [12]. In short, pores with radii ranging around nanoscale are critical for the pore volume and permeability of tight gas sandstone.

Pore structures of tight sandstones control the flow capacity, water saturation, and the producible pore volume, which are key factors affecting the gas production. Lyu et al. and Huang et al. emphasized the importance of the pore throat size for the mobile fluid percentage of tight formations [13, 14]. Fluid bound to pores increases with the decrease in detrital quartz content and the increase in clay mineral content, whereas connected pores allowing mobile fluids increase with the increase in authigenic quartz content and the decrease in detrital quartz content [15]. The correlation of micropore (2.5-10 nm) abundance with porosity is altered by diagenesis and compaction, causing the weak correlations between porosity and permeability. Shao et al. concluded that large numbers of small pores contribute to high storage capacity, whereas small numbers of large pores improve the flow capability in tight formations [16]. However, Sakhaee-Pour and Bryant demonstrated that samples exhibit greater production with an increased fraction of the intragranular porosity [17]. This is particularly true when the permeability is less than 1 mD, and the permeability is dominated by nanopores and micropores, following the size classification scheme by Loucks et al. [18]. Hence, nanopores are increasingly important for reservoir storage capabilities with decreasing permeability [9]. Based on the empirical permeability models, Kolodzie found that the permeability is best correlated with the throat size corresponding to a mercury saturation of 35% (r35), and Zhao et al. indicated that a throat radius corresponding to 30% mercury saturation (r30) is the dominant throat radius [19, 20]. Further, some researchers reported that pore throat radii corresponding to a mercury saturation of 10% (r10) are the best predictor of permeability for tight gas sands [21]. Permeability estimation models can be obtained by pore distribution and porosity. A part of pore types strongly affects coal permeability, including movable meso-macropore porosity and movable cleat porosity [22]. Different degrees of diagenesis and compaction cause various types of sandstones and determine the dominant pore size, which also emphasize the necessity for characterizing the variable pore types and sizes influencing porosity and permeability of tight gas sandstones [23].

In summary, previous studies emphasized the stronger relationship between permeability and pore sizes of tight reservoirs than that between permeability and the total porosity [24, 25]. Therefore, attentions should be paid to a comprehensive description of pore sizes and pore types (pore structures). However, detailed characterization of pore features in tight gas sandstones is less well documented in the literature. Various current techniques have been applied to pore structure characterization of porous media, including fluid invasion methods (mercury injection porosimetry [13] and nitrogen adsorption method [26]), direct imaging methods (SEM, transmission electron microscopy) [27], and noninvasion methods (X-ray computed tomography, NMR) [2830]. Among the techniques above, mercury injection porosimetry and nitrogen adsorption, along with SEM, can detect tiny pores. However, pore structures can be distorted by mercury injection porosimetry with high pressure, while ink-bottle pore throat structures are probably not well constrained [20, 31].

The focus of this study is to explore a new technique to probe small pore features and their size, shape, and connectivity because nanopores (<50 nm) comprise a nontrivial fraction of porosity in tight sandstones. NMRC is a fairly new method that has been gradually applied to many porous media, such as silica gel, cements, bones, and fibers [3234]. For example, Hansen et al. reported highly consistent results between pore size distributions of industrial silica, determined by nitrogen adsorption isotherm and NMRC methods. Cement and cement paste were widely tested by the NMRC method. The experimental results of Korb et al. showed that pore sizes ranging from 2 nm to 500 nm can be measured by NMRC precisely and reproducibly [35]. Pores with a radius of 2 nm were identified in bones by Gun’ko et al. using NMRC, and also, a large fraction of pore volume with pore sizes in the mesopore range was indicated in bone by NMR [36, 37]. Recently, NMRC has been utilized by the petroleum industry to characterize pore structure [3840]. Z. P. Wang and T. Wang concluded that the pore information determined from NMRC demonstrated its advantage compared to other methods, with the ability to identify closed as well as open micropores accurately below 20 nm [41]. A number of studies on petroleum reservoir samples found that the NMRC method produced results consistent with the mercury intrusion method or the nitrogen adsorption curve, especially for shale and coal [38, 42, 43]. The NMRC results of shale powder reflect the information on not only open pores but also closed ones [42]. In addition, NMRC usually yields larger pore volume than the nitrogen adsorption method, since drying for BET analysis induces pore shrinkage/collapse [38]. However, as noted previously, few studies have been done on tight gas sandstones using NMRC.

The purpose of this study is to explore the application of NMRC on tight gas sandstone samples with different magnitudes of permeability collected from Linxing tight gas sandstone reservoirs, China. The focus is on the pore structures, especially nanopores and their influences on permeability. NMR was employed to obtain the full pore size distribution and also used to obtain the pore morphology when in concert with SEM. NMRC and nitrogen adsorption method were used to characterize pores with small size. The ultimate aim of the study is to explore the application of NMRC to nanopores of tight gas sandstones and obtain more detailed and accurate information of pore structures.

2. Experiment Section

2.1. Samples

The Linxing area is mainly located in the eastern edge of the Ordos Basin, China, and the structure is located around the Yishan slope, Ordos Basin, and Shanxi western flexure belt. Core samples (Figure 1) were collected between 1923 and 1929 meters depth of the LX-35 Well, Taiyuan Group, which belongs to tight gas sandstones. Table 1 indicates that the permeability of these low-porosity samples decreases with the increase in depth, from 0.88 mD to 0.05 mD. Contact angle results in Table 1 show a strong water wettability, which is related to the high clay content as shown in Table 2. X-ray results show that clay plus quartz exceeds 95% of the total mineral content, each exceeding 44%. Sample LX1 contains 3.2% pyrite, while minor feldspars also appear in samples, including K-feldspar and plagioclase.

In addition, the sample was divided into two parts: one part is for the NMR test and the other part is prepared for three separated assessments. Tight sandstone powder is for nitrogen adsorption, a tight sandstone fragment is for NMRC, and a tight sandstone slice is for SEM.

2.2. NMRC Setup

The sketch in Figure 2 is NMRC12-010V from Suzhou Niumag Analytical Instrument Corporation. The setup is mainly composed of two systems: a temperature controlling system and an NMR system. The instrument has a resonance frequency of 11.537 MHz, a magnet temperature controlled around 32.00 ± 0.02°C, and a probe coil diameter of 10 mm. The temperature range is -35°C to 40°C (-308~-233 K) with an accuracy of ±0.1°C.

2.3. Experimental Procedure

Nitrogen adsorption-desorption isotherm measurement has been frequently described in publications [4, 42]. Hence, the experiment procedure is not discussed in our manuscript. The experimental parameters and procedures of NMR and NMRC procedures are introduced as follows.

2.3.1. NMR Measurement

The procedures are described as follows. First, dry weight and core size (length, diameter) of the samples are measured. Second, samples are evacuated for 2 hours and then saturated with distilled water under 25 MPa for 48 hours. After the saturation, wet weights are measured. The water porosity of each core is calculated by dry weight, wet weight, and its volume. Third, raw NMR data of water-saturated cores are performed. The NMR CPMG parameters are listed as follows: echo spacing, 0.1 ms; waiting time, 300 ms; echo numbers, 6000; and numbers of scans, 32.

2.3.2. NMRC Measurement

As a comparison, the samples for NMRC are cut from NMR samples with the length of ~0.5 cm. Detailed experimental procedures are as follows: step (1): fragment samples were evacuated for 8 hours and then saturated with distilled water for 48 hours (20 MPa); step (2): NMR tubes filled with samples were inserted into the middle of the probe coil; and step (3): the temperature system lowered the sample temperature to – 30°C to freeze all fluids in pores and then increased it to 10°C at 0.1°C increments. At the same time, the fluid content was detected by NMR CPMG spin echo sequencing as a function of changing temperature. The NMR CPMG parameters are listed as follows: echo spacing, 0.1 ms; waiting time, 500 ms; echo numbers, 2500; and numbers of scans, 128.

2.4. NMRC Principles

During the NMRC measurement, porous media are saturated with a fluid, cooled to a base temperature, and then heated at increments of 0.1°C. The phase change of fluids in pores with temperature follows the Gibbs-Thomson equation (equation (1)) [44, 45]. where is the crystal diameter, since the pore is filled with the fluid or crystal; thus, also represents the diameter of the pore; is the melting point of crystal with diameter of ; is the melting point of the crystal with infinite diameter; is the surface energy at the crystal-liquid interface; is the solid density; is the contact angle between fluid and pore surface; and represents the bulk enthalpy of fusion (per gram of material). Equation (1) can be simplified to where is the geometric constant, , is a constant related to the solid-fluid system, and is a constant related to surface energy and reflects the contact angle between fluid and pore wall. is the melting point depression.

The product of , , and can be represented by ,which is the thermodynamic constant. The value of is 58 (with the unit of ) for distilled water when used to convert the melting curves.

From equation (3), the pore size has the one-to-one relationship with the melting point depression [46]. The process of increasing the total fluid content with increasing temperature can be regarded as a progressive accumulation process of pore volume from small pores to large ones, in which the fluid volume can be measured by the probe (Figure 3). The relationship between fluid content and temperature in the process can be represented by either solidification or melting curves. Instead of tracking the solidification curve, the melting curve is selected as it yields greater reproducibility, which can be transferred to pore size distribution by the Gibbs-Thomson equation (equation (1)). Pores with radii ranging from 1 nm to 500 nm can be accurately measured by the NMRC method within the temperature range of -30°C to 0°C (±0.1°C). In addition, water wettability is a favorable situation for the accuracy of NMRC results when using distilled water as the filling fluid.

3. Results

3.1. Pore Types and Minerals

A TESCAN-VEGA SEM was employed to document the microscopic morphology, and the results are shown in Figure 4. Figures 4(a) and 4(d) show a full microscopic view of sandstone, which reveals most interparticle pores ranging between 1 and 5 microns.

Pores are filled with various minerals, including irregular platey kaolinite (Figure 4(f)), thick platey dickite (Figure 4(i)), and filamentous or flaky illite (Figures 4(h) and 4(i)). Clay-dominated pores are distributed in filamentous and platey illite and range between 1 and 4 microns (Figures 4(c) and 4(h)). Although pyrite was only detected in sample LX1 by X-ray, it can be observed by SEM in all three samples (Figures 4(b), 4(e), and 4(g)). As can be seen from Figure 4, residual interparticle pores and clay-dominated pores are the main pore types observed in tight gas formation. SEM observations indicate that microfractures are in low abundance, which differs from what others have seen in tight oil sandstones [13].

3.2. NMR Distributions at Full Saturation

As can be seen from Figure 5, NMR distributions of three tight gas sandstone samples all show three peaks, marked as P1, P2, and P3. Three peaks, respectively, range around 0.3 ms, 10 ms, 300 ms, respectively, except peak P3 of core LX3 distributed around 100 ms. Areas under P1 or P2 are both larger than that of P3 for all three samples, respectively, corresponding to three scales of pore size. Therefore, small pores are enriched developed in tight gas sandstones. The increase in peak intensities for P2 and P3 correlates with larger pores and increasing permeability.

The values of NMR spectra have one-to-one relationship with pore radii. Three samples share a similar size range of small pores represented by P1 and P2, whereas pores represented by P3 for LX3 exhibit a smaller size compared to samples LX1 and LX2. In addition, three peak amplitudes of LX3 are greater than those of cores LX1 and LX2, indicating greater porosity and hence greater capacity to host fluids. Peaks P2 and P3 of LX3 are well connected, while two peaks of other two samples are disconnected. Because the permeability and porosity of LX3 are both larger than those of LX1 and LX2, as listed in Table 1, we suggest that the contradiction may be explained by the differences in connectivity between peak P2 and peak P3. Peak P2 and peak P3 of LX3 are connected, indicating that pores are well connected.

3.3. Nitrogen Adsorption-Desorption Isotherms

Nitrogen adsorption-desorption isotherms of three samples belong to type IV of IUPAC definition [47]. As shown in Figure 6, the presence of mesopores (2-50 nm) is mainly reflected. Three adsorption isotherms start at the relative pressure value of 0.02, reflecting the presence of micropores (<2 nm). Because of capillary condensation or micropore filling during the adsorption, two isotherms do not coincide with the desorption curve located above the adsorption curve, which always forms a hysteresis or lag loop at high relative pressure, referred to as adsorption equilibrium isotherm of type IV. The characterization of lag loops reflects the distinct pore structure of tight gas sandstones [47, 48].

According to the IUPAC method, lag loops of three samples belong to the H3 type, in which the adsorption and desorption isotherms slowly rise; then, the amount of adsorption increases rapidly at higher pressure and the desorption loops are initially narrower. These characteristics of adsorption and desorption isotherms reflect mesopores and may be caused by stacking platey particles. These pores are mainly slit pores. Based on the SEM results of Section 3.1, these platey minerals are probably illite, kaolinite, and dickite (Figures 4(c), 4(f), 4(h), and 4(i)).

Pore size distributions were also obtained using the BJH (Barrett-Joyner-Halenda) method, especially for micropores (2-50 nm). Figure 7 shows pore volumes of LX2 and LX3 evenly distributed between 2 and 100 nm, while that of LX1 ranges between 20 and 110 nm. The shape of curves of pore volume versus pore sizes for three samples is very similar for pore sizes greater than 70 nm. As noted above, the biggest difference in pore development falls in the 2-30 nm range wherein LX2 and LX3 differ significantly from LX1.

3.4. NMRC Results

NMRC can directly obtain the pore volume versus pore size curves according to the linear relationship between pore volume and signal intensity. Figure 8 shows the pore size (radius) distributions of three samples by NMRC, in the form of multiple peaks across a pore diameter range of 2 nm to 500 nm. The main pore volumes for three samples display somewhat different patterns in Figure 8; LX1 exhibits peaks distributed between 1 and 100 nm, and LX2 and LX3 have peaks ranging between 2 and 500 nm. The first five peaks of LX1 and LX2 are all located around 2 nm, 4 nm, 8 nm, 16 nm, and 65 nm, respectively. Pores of LX1 and LX2 distributed at 2 nm show the largest pore volume compared to pores observed at the other peaks, whereas seven peaks of LX3 are evenly spread from 2 nm to 500 nm. Similar to the results of nitrogen adsorption, pore volume versus pore size curves of LX2 and LX3 exhibit somewhat greater pore volume compared to LX1, but all three samples display similar peak positions especially for LX1 and LX2 apart from pores ranging around 1-2 nm. Pore volumes of nanopores (pore with ) of LX1 tend to be dominant, while pore volumes of micropores (pore with radii ranging 2-50 nm) and mesopores (pore with ) of LX2 and LX3 are more prominent compared to LX1 [49]. In summary, the pore volume of LX1 is mainly contributed by nanopores, whereas pore volumes of LX2 and LX3 are mainly contributed by micropores and mesopores.

4. Discussion

4.1. Application of Three Methods to Nanoscale Pore Characterization of Tight Gas Sandstones

NMRC, nitrogen adsorption, and NMR were employed to characterize the pore structure of tight gas sandstones. NMR distribution at full saturation can reflect the full pore size distribution of porous media. However, the conversion factor between NMR and pore radii is difficult to determine. Radii of pores described by nitrogen adsorption ranged from one nanometer to tens of nanometers, while NMRC can describe pores with radii ranging from 2 nm to 500 nm. A comparison of results from nitrogen adsorption and NMRC is shown in Figures 7 and 8, respectively, indicating that pore size distributions determined by NMRC exhibit higher pore volumes relative to those from nitrogen adsorption. Hence, NMRC has a practical application to tight gas sandstones, especially for the nanopore characterization.

4.2. Permeability, Porosity, and Pore Size Distribution

Based on the results of three samples from the NMR method, T2 values of main peaks are lower than 100 ms, and half of the total T2 signal is composed of the signal from water saturating small pores. The overall NMR T2 peak distributions reflect an increased relative pore volume with increasing permeability, as well as the total pore volume. However, P3 of LX1 and LX2 is located between 100-500 ms, while P3 of LX3 is located approximately between 50-200 ms. This indicates that LX1 and LX2 have wider pore size distributions than LX3. Pores associated with P3 of LX3 are smaller size than those of LX1 and LX2, whereas the pore connectivity of LX3 is greater than that of the other two samples. On one hand, the differences in connectivity among three samples may explain the variability in the permeability and pore size. Conversely, pores reflected by P1 are abundant but small, versus pores represented by P3 that are larger but rare. Therefore, we surmise that pores associated with the peak P2 (~7 ms) likely control the permeability. Based on our previous work and related studies [13], the conversion between NMR T2 and pore size for tight sandstone formations is about 50 nm/ms, where pores associated with P2 are about 350 nm. The combination of Nitrogen adsorption and NMRC provide a solid basis for characterizing these pores accurately across roughly two orders of magnitude in dimension.

The pattern of pore size distribution curves determined by NMRC shown in Figure 8 does not coincide with the sequence of permeability values. Based on the weighted mean calculation of NMRC pore sizes, the results of the weighted mean pore diameters are 25.8 nm, 41.31 nm, and 49.99 nm, respectively, for LX1, LX2, and LX3. These values of the weighted pore diameters are consistent with the sequence of permeability values. Because pores with radii ranging 1-10 nm exhibited by LX1 and LX2 are abundant, it results in a low relative proportion of larger pores as well as the poorer connectivity. Three tight gas samples share almost the same porosity and mineral composition but very different permeabilities. Figure 9 indicates that the correlation between permeability and porosity is weak. The significant difference in permeability probably reflects differences in the relative ratios of pore size to pore throat diameters. The pore size can be divided into three intervals, which is marked in Figure 10. The sequence of pores falls in NMR interval II of the pore size distribution in Figure 10 of three samples and follows the sequence of permeability. Pores with size belonging to interval II are probably potential ones determining permeability. Therefore, it is not the absolute amount of potential pores but the relative proportion of potential pores determining the permeability (Figure 10), in which potential pores means pores determining connectivity. Combined with the results of SEM and NMR , pores fall in interval II in Figure 10 with radii ranging between nanometer to several microns, which are mostly clay-dominated pores identified by SEM. As noted above, we can conclude that these potential pores determining the permeability are clay-dominated pores. Therefore, neither small pores with (around peak 1 in NMR T2 distribution) nor large pores with (around peak 3 in NMR distribution) control the permeability by analyzing the connectivity of the pores in various size ranges, but rather pores averaging approximately (around peak 2 in NMR distribution) have sufficient connectivity to host and transmit hydrocarbons.

5. Conclusions

Three tight gas sandstone samples from the Linxing area, Ordos Basin, were collected to study the pore structure using various experimental techniques, including SEM, nitrogen adsorption, NMR , and NMRC. The following conclusions are obtained. (1)Residual interparticle pores and clay-dominated pores are the two primary pore types for our target tight gas sandstone samples. Clay and quartz are the main minerals, accounting for almost excess 95% of the mineralogy content. The combined results of SEM and X-ray test show that pyrite is significantly distributed(2)The NMRC technique has a wide range of application to the pore structure characterization of tight gas sandstones, especially for pores with nanoscale radii. It can be a very promising method especially when referred to as a promising “roadmap” on how to interrogate tight formations, such as the tight gas sands or even shale(3)Porosity of tight gas sandstones is very low, exhibiting extremely poor correlations with the permeability. It is neither tiny pores nor large pores controlling the permeability of tight gas sandstone but relatively abundant pores, ranging around peak P2 of NMR distribution (~350 nm). Pore types and their relative amounts of pores are critical factors determining the permeability. These potential pores determining the permeability of tight gas sandstones are associated with clay packets and folia (clay-dominated pores), which contribute a nontrivial volume to tight gas sandstones

Data Availability

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

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The research is also supported by the Science Foundation of China University of Petroleum, Beijing (Grant No. 2462021XKBH007). In addition, DRC was supported by the US Department of Energy, Office of Science, Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences, Geosciences Program, under grant DE-SC0006878 to Ohio State University.