Abstract

Shale gas reservoirs have pores of various sizes, in which gas flows in different patterns. The coexistence of multiple gas flow patterns is common. In order to quantitatively characterize the flow pattern in the process of shale gas depletion development, a physical simulation experiment of shale gas depletion development was designed, and a high-pressure on-line NMR analysis method of gas flow pattern in this process was proposed. The signal amplitudes of methane in pores of various sizes at different pressure levels were calculated according to the conversion relationship between the NMR relaxation time and pore radius, and then, the flow patterns of methane in pores of various sizes under different pore pressure conditions were analyzed as per the flow pattern determination criteria. It is found that there are three flow patterns in the process of shale gas depletion development, i.e., continuous medium flow, slip flow, and transitional flow, which account for 73.5%, 25.8%, and 0.7% of total gas flow, respectively. When the pore pressure is high, the continuous medium flow is dominant. With the gas production in shale reservoir, the pore pressure decreases, the Knudsen number increases, and the pore size range of slip flow zone and transitional flow zone expands. When the reservoir pressure is higher than the critical desorption pressure, the adsorbed gas is not desorbed intensively, and the produced gas is mainly free gas. When the reservoir pressure is lower than the critical desorption pressure, the adsorbed gas is gradually desorbed, and the proportion of desorbed gas in the produced gas gradually increases.

1. Introduction

Deep shale gas in China is abundant and owning an immense potential for exploitation, which is becoming a key area for shale gas production in future [1, 2]. The shale gas reservoir is a kind of unconventional gas reservoir, with ultralow porosity and ultralow permeability, being complex in pore structure [3]. A large number of pores with diameter smaller than 50 nanometers are well developed in shale gas reservoir [46]. In addition, microfractures and micron-sized pores also exist, with a storage space system of multiscale. The gas in shale reservoirs is in a free or adsorbed state, with a complicated shale gas flow mechanism. According to the concentration of the gas, with Knudsen number (), the gas flow can be divided into four types, including continuous flow, slip flow, transitional flow, and free molecular flow [7, 8]. The is the ratio of the average molecular free path to the characteristic scale. The scale of pores in shale gas reservoirs varies significantly, resulting in different gas flowing patterns in different scale pores. However, there is a lack of the experimental method for quantitatively characterizing the pattern of shale gas flow in shale reservoirs.

NMR core analyzing method has been widely used in the field of oil and gas seepage mechanics, including the NMR porosity testing [9, 10], movable fluid saturation testing [11], water-driving-oil NMR analysis [1217], gas-driving-oil NMR analysis [18], and NMR pore structure characterization [1922]. When the study objects change from conventional reservoirs to unconventional reservoirs such as tight oil and shale oil reservoirs, the pore sizes are getting smaller and smaller, and the pore structures are becoming more and more complex. Particularly, large errors might occur when conducting the water-driving-oil and gas-driving-oil experiments on tight oil or shale oil cores with conventional sampling method due to the very small flow rate. In contrast, with the online NMR method under high-temperature and high-pressure, the variation of fluid relaxation in the core can be monitored in real-time, allowing a more accurate measurement of the variation of oil and water in the core. Therefore, the application of NMR core analysis technology will be further expanded when studying the seepage mechanics for unconventional oil and gas reservoirs [23].

During shale gas depletion development experiment, combined with high-temperature and high-pressure online NMR analysis technology, the occurrence and recovering regulation of adsorbed and free gas can be analyzed. Some key parameters can be calculated from the experiment, such as the critical desorption pressure and the recovering proportion of adsorbed gas [24]. The above experiments have all suggested that the gas production from shale gas reservoirs decreased rapidly at the initial stage and became slow at the middle and late stages, which is consistent with the gas production in the actual production process of shale gas wells [25]. The critical desorption pressure of the shale gas from the Longmaxi Formation in the Sichuan Basin is 12-15 MPa. When the reservoir pressure is lower than this critical desorption pressure, the recoverable degree of adsorbed gas increases significantly. When the produced shale gas is mainly free gas, the recoverable degree of adsorbed gas is low.

The low magnetic field NMR spectrum distribution is directly related to the pore structure and can reflect the pore distribution of the sample to a certain extent [26, 27]. The key to study the pore size distribution with NMR test is to determine the surface relaxation rate and the shape of the pores. Generally, the surface relaxation rate is calculated by comparison with the capillary force curve obtained from the mercury intrusion capillary pressure experiment. By selecting an appropriate conversion coefficient (), the capillary pressure curve can be coincided with the pseudo capillary pressure curve. The derivation process of an equation for the conversion of the NMR transverse relaxation time value to the throat radius was proposed. On an assumption that a linear relationship exists between and the pore size distribution, the conversion coefficients of the NMR pore size distribution spectrum and the pore radius were calculated. Then, the NMR distribution was converted to the NMR capillary force curve, during which the value of the conversion coefficient () ranges from 33.3 to 250 [2830]. As for the double-peak spectrum, two capillary pressure curves were separately constructed on both micropore and macropore segments by using two different power functions, and the linear method and power function method were used to construct the capillary pressure curve separately [31]. Li et al. proposed that under the condition of pore diameters less than 5 μm, the -pore diameter relationship curves of the brine saturated sample and the oil-saturated sample are almost overlapped, indicating that the distribution is irrelevant with the fluid volume, but only associated with the pore size. In this case, no matter what fluid is saturated in the pores, the distribution is linearly related to the pore size distribution [32]. Hürlimann et al. determined the surface relaxation rate of sandstone and value through diffusion experiment [33].

This paper presents an on-line NMR analysis method of flow pattern in the experimental process of shale gas depletion development. Specifically, the relative amplitude of methane in pores of various sizes under different pressure conditions is calculated according to the conversion relationship between spectrum and pore radius, and then, the flow patterns of methane in pores of various sizes under different pore pressure conditions are analyzed as per the flow pattern determination criteria. Application of this method is exemplified in the first submember of the first member of the Silurian Longmaxi Formation (L11) in Well Ning 203.

2. Materials and Methods

2.1. Experimental Methods

The shale core sample was collected from the Long11 submember of the Silurian Longmaxi Formation in Well Ning 203 in the Changning area of the Sichuan Basin. The core sample, with a length of 73.62 millimeters and a diameter of 25.72 millimeters, was taken from depth interval ranging from 2384.69 m to 2384.86 m. The porosity of the core sample is 6.07%, and the Klinkenberg permeability is 1.02 mD. To truly reveal the seepage characteristics of natural gas in the reservoir, methane gas with a purity of 99.99% was used during the experimental process. A high-temperature and high-pressure online NMR core analyzer (ReCore-2515), which is independently developed by the Seepage Fluid Mechanics research laboratory of Research Institute of Petroleum Exploration and Development, PetroChina, was used in this experiment. The experimental setup is shown in Figure 1.

The experiment includes four main steps. Firstly, dry the core sample at 105°C for 48 hours, cool it to room temperature in a drying dish, and put it into the holder. Secondly, test the air tightness of the system, open the inlet valve of the holder and close the outlet valve, evacuate the core until it reaches the vacuum requirements, and then close the vacuum pump valve. Thirdly, the methane gas is injected into the core, pressurized to 20 MPa, and fully saturated for 24 days to simulate the original occurrence state of adsorbed methane gas. The change of methane gas signal in the rock core is detected in real time during the saturation process. And the final equilibrium pressure is set at 18.26 MPa. When the inlet and outlet pressures no longer change after 96 hours, the core is considered to return to its original state. Due to the tightness and adsorption characteristics of the shale reservoir, this process might last up to tens of days. Finally, open the outlet, start the depletion development experiment, and record the NMR spectrum, the pressure at each measuring point, and gas production data at the outlet valve in real-time.

The experiment includes four main steps. (1) The core was dried at 105°C for 48 h and then cooled to 25°C in a drying dish before being put into a holder. (2) Check the tightness of the system. Open the inlet valve of the holder and close the outlet valve to vacuum the core. Close the inlet valve after reaching the vacuum requirement. (3) Methane gas is injected into the core in a constant pressure mode to restore the original reservoir occurrence. Record the pressure at both ends of the core and close the gas source after saturation to 18.26 MPa at the inlet and outlet of the core. When the inlet and outlet pressures do not change during 96 h, the core is restored to the original occurrence state. Due to the density and adsorption properties of shale reservoirs, this process can last for tens of days. (4) Open the outlet value, start the depletion development experiment. The data of spectrum, pressure at each measuring point, and gas production at outlet were recorded during the experiment.

2.2. NMR Pore Size Distribution

Low-magnetic-field NMR test results can provide abundant information on the fluid in porous media. In general, the relaxation time of fluid in porous media includes bulk relaxation, surface relaxation, and diffusion relaxation. In this study, the Carr-Purcell-Meiboom-Gill test sequence has been used, and only the surface relaxation needs to be considered [34, 35]. The relationship between the information of fluid in pores and relaxation time of the NMR test is shown in Equation (1): where is the relaxation time, ms; is the surface relaxation rate, μm·ms-1; and is the ratio of specific surface area to pore volume, μm-1.

According to Equation (1), small pores have larger values. Therefore, small pores have shorter relaxation times than large pores. The relationship between the spectrum and the pore radius can be expressed by: where is the shape factor of pores (3 for spherical pores and 2 for columnar pores), dimensionless, and is the pore diameter, μm.According to Equation (1) and Equation (2), the relationship between the pre diameter () and the relaxation time () can be obtained:

where is conversion coefficient set as 33 μm·ms-1 according to data range proposed in the literature [16].

2.3. Determination of the Flowing Pattern

Gas molecules frequently collide with each other during the movement process. The gas transportation process comes from the thermal motion of molecules. The collision causes the molecules to change continuously the direction and rate of movement, making the molecular movement path very tortuous. The collision also allows molecules to exchange continuously the energy and momentum, being the necessary prerequisites for the balance of the system. The path that the actual gas molecules travel between two collisions is called the free path. For gas molecules, the average distance between two adjacent collisions, known as the average free path of the molecule, is a physical quantity of only statistical significance, as shown in Equation (4) [36, 37]:

where is the effective diameter of the molecule, μm; is the gas pressure, MPa; is the Boltzmann’s constant, J/K; and is the temperature, K.

Knudsen number () represents the ratio of the average free path of molecules (λ) to the characteristic length () of the object in the flow field.

3. Results and Discussion

3.1. Experimental Results

Figure 2 illustrates the spectrum distribution at six stages during the experiment, including at the beginning and when the experiment was carried out for 7, 23, 43, 62, and 121 days, respectively. At the beginning of the experiment, there are two peaks in the spectrum, and the signal amplitude of the right peak is about 1 times higher than that of the left peak, indicating that there are relatively many macropores in the shale core sample. The first 23 days of the experiment is equivalent to the initial stage of the experiment. As the average pore pressure of the core decreased, the signal amplitude corresponding to the right peak also decreased, but the signal amplitude corresponding to the left peak was almost unchanged, indicating that the gas recovered at the initial stage was mainly free gas in macropores, and the free gas and adsorbed gas in the mesopores and macropores were less recovered. This gas production law is in good accordance with the actual gas well production performance. When the experiment entered the middle and late stages, that is, 23 days after the experiment, the signal amplitude of the right peak and the left peak began to decrease simultaneously, indicating that the gas in the pores of all scales has been recovered. At present, it is generally believed that the adsorbed gas mainly exists in pores of nanoscale (i.e., pores with a pore size of smaller than 100 nanometers) [36]. According to Figure 2, the signal of the left peak varies slightly as the average pore pressure is higher than 12.6 MPa, indicating that the adsorbed gas has not been widely desorbed. When the average pore pressure is less than 12.6 MPa, the signal amplitude of the left peak begins to gradually decrease with the decrease of the average pore pressure, indicating that the critical desorption pressure is around 12.6 MPa, and the large-scale desorption of adsorbed gas only occurs when the pore pressure is lower than the critical desorption pressure. When the average pore pressure is greater than 12.6 MPa, the proportion of adsorbed gas in the total produced gas is less than 10%. With the production of methane from shale, the proportion of adsorbed gas in the total produced gas gradually increases, and it reaches 27.2% when the average pressure is 4.48 MPa, as shown in Figure 3.

3.2. Discussion of the Flowing Pattern of Shale Gas during the Depletion Development

According to Equation (3), the relaxation time at the horizontal axis in Figure 4 was converted into pore diameter. The average pressure of pores in the shale gas depletion development experiment dropped from 18.26 MPa to 4.48 MPa, and the NMR spectrum corresponding to the two pressure states is shown in the following figure. According to Equation (3), the molecular free paths were, respectively, calculated  MPa and  MPa. The Knudsen number is obtained by the ratio of molecular free path to characteristic scale. The flow state can be divided by the interval of . That is, according to the degree of gas concentration, the flow can be divided into four regions based on . This comprises continuous flow (), slip flow (), transitional flow (), and free molecular flow () [8, 9].

The nuclear magnetic resonance spectrum of methane saturated shale has obvious bimodal characteristics. Methane is mainly adsorbed on the surface of nanopores in shale. The adsorption is mainly controlled by surface relaxation. The relaxation time is short between 0.1 and 1.0 ms, and the main peak is 0.4 ms, corresponding to the left peak on the spectrum of NMR. Free methane exists in larger shale pores, which is not bound by pore wall, and has a long relaxation time which corresponds to the right peak on spectrum. The relaxation time is between 1 and 100 ms, and the main peak is 10 ms. Thus, the signal quantity of adsorbed/free state of methane gas in shale can be determined by spectrum. However, it is generally believed that the adsorbed gas mainly exists in pores of nanoscale (i.e., pores with a pore size less than 100 nanometers) [34]. Therefore, by calculating the ratio of the amplitude of the left peak to the total amplitude in the spectrum at the pressure of 18.26 MPa, the proportion of the adsorbed gas in the total gas in the initial state can be determined as 33.3%.

The free path of the gas molecule was calculated according to equation (3), and then, the number was obtained by dividing it through the pore diameter. Figure 5 shows how the number varies with the pore size at average pore pressure conditions of 18.26 MPa and 4.48 MPa, respectively. According to the criterion for determining flowing pattern, there are three flowing patterns, which include continuous flow, slip flow, and transitional flow during the shale gas depletion development experiment. When the average pore pressure of cores is 18.26 MPa, the pore size is less than 1.39 nm for transitional flow, 1.39–12.47 nm for slip flow, and more than 12.47 nm for continuous medium flow. When the average pore pressure is 4.48 MPa, the pore size is less than 5.81 nm for transitional flow, 5.81–57.1 nm for slip flow, and more than 57.1 nm for continuous medium flow. In Figure 4, the ratio of the area enveloped by the signal amplitude line for each flow pattern under each average pore pressure and the -axis to the area enveloped by the whole signal amplitude line, and the -axis is the proportion of the flow pattern under the average pore pressure.

Between Figures 4 and 5, it can be inferred that (1) when the pressure is 18.26 MPa, the flow pattern is mainly slip flow and continuous medium flow, with the former observed in pores with diameter of 1.39–12.47 nm, accounting for 6.1% of the total flow, and the latter in pores with  nm, that is, the Darcy flow, accounting for 93.9% of the total flow; and (2) when the pressure is 4.48 MPa, transitional flow occurs in the pores with  nm, accounting for 0.3% of the total flow; slip flow occurs in pores with diameter of 5.81–57.1 nm, accounting for 41.5% of the total flow; and continuous medium flow occurs in pores with  nm, accounting for 58.2% of the total flow.

Continuous medium flow is dominant at high pore pressure. However, it is important to note that continuous production in the shale reservoir results in a decrease in pore pressure, increase in Knudsen number, and an expansion of the pore size range of the slip flow zone and transitional zone, with both zones also being characterized by an increase in the proportion of diffusion flux to total flux. In the slip flow zone, there is a thin gas layer near the wall, that is, the Knudsen diffusion layer, and the flow in this zone needs to be simulated by the slip boundary conditions near the wall. In the transitional flow zone, the Knudsen diffusion layer is wider to trap most of the fluids, and the flow in this zone needs to be described by rarefied gas dynamics. The gas flow patterns obtained by NMR analysis are consistent with the findings in previous studies [36, 38], proving that the propose method is scientific and practical.

At present, there is little research on shale gas flow pattern analysis. For shale, a porous medium with large pore space scale span, usually only static flow pattern analysis under a certain pressure, is carried out. However, this is limited to flow pattern analysis under several characteristic scales. This study provides a dynamic analysis method of flow pattern in the process of depletion development. It is expected to be applied to the flow pattern analysis of depletion development of deep shale gas reservoirs under high reservoir pressure.

4. Conclusions

The on-line NMR analysis method of gas flow pattern can quantitatively analyze the gas flow patterns occurred under different pressure and pore size conditions in the process of shale gas depletion development and clarify the contribution of each flow pattern to the total flux in such process.

The shale gas depletion development experiment reveals three flow patterns: continuous medium flow, slip flow, and transitional flow. When the pore pressure is high, the continuous medium flow is dominant. With the decrease of pore pressure, the Knudsen number increases, and the pore size range of slip flow zone and transitional flow zone expands.

In the process of shale gas depletion development, there is a critical desorption pressure. When the reservoir pressure is higher than the critical desorption pressure, the adsorbed gas is not desorbed intensively, and the produced gas is mainly free gas. When the reservoir pressure is lower than the critical desorption pressure, the adsorbed gas is gradually desorbed, and the proportion of desorbed gas in the produced gas gradually increases.

Data Availability

Data available on request.

Conflicts of Interest

The authors declare no competing interests.

Acknowledgments

This paper is supported by the National Major Research Program for Science and Technology of China (2017ZX05037-001) and the Research Fund for Basic Research and Strategic Reserve Technology of Institutes Directly under CNPC (2018D-500806).