1. Introduction
When carrying out physical modeling to study the characteristics of gas and water, gas–water driving (GWD) models are generally used, with semi-permeable baffle (SPB) models being used less frequently to study the hydraulic properties of rock. Furthermore, nuclear magnetic resonance (NMR) and magnetic resonance imaging (MRI) techniques have recently been used to study both free and combined water in this context.
Many articles in the literature have studied gas–water driving characteristics using physical models [
1,
2,
3,
4,
5,
6,
7], considering factors such as the design of the driving fluid and displacement conditions, measurement of the driving pressure, the volume of discharge water, and the calculation of residual water. These articles have generally focused on tight sandstone or heavy oil. If necessary, numerical simulations may also be carried out according to the experimental conditions.
Some scholars have also conducted studies involving the visualization of driving experiments [
1,
8]. Lysova [
8] analyzed the influence of the exothermic reaction by determining the residual water distribution using magnetic resonance imaging (MRI), while Masri [
1] placed CT imaging devices on the core displacement chamber to analyze the spatial changes in water saturation after displacement. These works allowed the residual water to be easily observed in a single picture, but they lacked an analysis of the change in residual water over different periods of the experiment.
There have been many articles focused on determining the distribution of core residual water according to the results of Centrifugal NMR experiments [
9,
10,
11,
12,
13,
14,
15], through which the fluidity of the reservoir fluid can be analyzed.
Semi-permeable baffle experiments are mostly used in the field of material research [
16,
17]. The wettability of porous materials has been previously studied [
18], which has reference significance regarding the wettability of rock samples, as semi-permeable experiments considering porous materials are similar to those using geological samples.
Few previous studies have focused on visualization of the dynamic displacement process. Therefore, in this paper, we attempt to visualize such a process through a gas-driven water displacement–magnetic resonance imaging (GWD-MRI) experiment. Furthermore, considering the lack of gas–water characteristic analysis in different physical models, we analyzed the residual water results obtained with different models.
The three types of natural gas charging simulation experiments considered in this paper have varying characteristics, as described below.
The centrifugal nuclear magnetic resonance (Cen-NMR) experiment uses different speeds to centrifuge the core, such that the water in the pores is driven out of the core under the action of the centrifugal force. Generally, the speed ranges from 1000 rpm to 10,000 rpm, such that the centrifugal force produced is in the range of 0.1–5 MPa and the time required for centrifugation is short. After centrifugation, the samples are tested by NMR spectroscopy and weighed at the same time, allowing for an analysis of the water content and water-bearing pore distribution characteristics under different centrifugal force conditions [
10,
12].
The semi-permeable baffle (SPB) experiment involves the use of a ceramic baffle, which only allows the liquid phase to pass through at the end of the core; this can be used to measure the capillary data. The core chamber is filled with natural gas at a certain pressure, which will drive out the water in the core under the action of pressure. The driven liquid phase enters the measuring cup through the semi-permeable ceramic plate, while the gas phase remains in the core chamber [
19]. Regarding the experimental conditions, the gas chamber pressure is generally 0–1 MPa and the experiment equilibrium time is very long (generally >6 h). Therefore, this type of experiment mainly simulates low-pressure trap charging.
The Gas-driven water displacement–magnetic resonance imaging (GWD-MRI) experiment uses a displacement device combined with a magnetic resonance imaging (MRI) device. According to the design, the core driving experiment can simulate different driving pressures, different driving times, and different types of fluid. Under the gas-driven water model, the driving pressure depends on the core petrophysical properties to select the appropriate displacement pressure (in the range of 0–5 MPa), and a driving time between 1–48 h must be appropriately selected. The displacement is very close to migration layer charging. In order to effectively visualize the drive process, this paper uses a technical method involving the combination of a drive device and a magnetic resonance imaging (MRI) device.
According to the above, the Cen-NMR experiment takes a short time, has a large driving pressure range, and can analyze the residual water distribution, making it suitable for analyzing the characteristics of short-term high-intensity charging. In contrast, the semi-permeable baffle experiment takes a long time and has a low driving pressure, which makes it better suited for simulation of the charging process in the underground trap, with the final water saturation generally well-matching the current water saturation of a typical gas field. Finally, GWD–MRI can be used to observe the spatial distribution characteristics of the residual water. Gas generally drives water under low-pressure conditions. In the driving process, the end of the core will discharge water and gas, such that the core driving experiment can better reflect the charging of the migration layer.
According to the experimental time and whether the core is damaged, we conducted GWD-MRI, Cen-NMR, and SPB experiments in sequence. Combining these three types of experiments allowed us to analyze the gas–water characteristics under migration layer charging, short-term high-intensity trap charging, and slow trap charging, respectively.
3. Result and Discussion
The displacement-MRI experiment was mainly conducted to analyze the gas infiltration characteristics in the migration layer, the semi-permeable baffle experiment was mainly conducted to analyze the trap charging characteristics, and centrifugal-NMR was mainly conducted to analyze the short-term charging characteristics.
3.1. Migration Layer Charging Results and Discussion
We analyzed the nuclear magnetic resonance imaging (MRI) results while carrying out gas-driven water displacement (GWD), which was completed for three samples.
3.1.1. Sample A10 Experimental Results and Discussion of Start-Up Pressure
Sample A10 (4320.9 m from Well A-2) had a permeability of 0.34 mD and a porosity of 8%. After a drive pressure of 0.3 MPa for 5 min and a drive pressure of 0.6 MPa for 20 min (for a cumulative driving time of 25 min), the nuclear magnetic imaging (MRI) results showed that the residual water was evenly distributed, indicating that natural gas cannot charge low-permeability cores under the condition of 0.6 MPa. After 0.3 MPa and 0.6 MPa driving, the water saturation of the samples decreased, which is considered to be the result of water evaporation caused by heating the core throughout the 40 min required for imaging.
Then, the drive pressure was increased to 1 MPa for 10 min (for a cumulative driving time of 35 min). The water content in the front of the core decreased, while that at the center and the end of the core did not change, and the water saturation was 93%, indicating that 1 MPa pressure allowed for charging into the core. When we continued to use 1 MPa driving pressure for 15 min (for a cumulative driving time of 50 min), the water content in the front of the core continued to decrease, and the water content in the center and end of the core began to decrease, with a water saturation of 87%. With a cumulative driving time of 60, 70, and 80 min, the water content decreased at all different core positions. Finally, when the cumulative driving time was 95 min, the water content at the front of the core was very low and the water content at the end of the core had also decreased significantly, with a water saturation of 75%.
An important feature of the migration layer charging infiltration is that there is a start-up pressure gradient. When the driving pressure gradient is lower than the starting pressure gradient, the natural gas cannot drive the formation water; only when the driving pressure gradient is greater than the start-up pressure gradient can the natural gas drive the formation water. The water saturation of sample A10 was basically unchanged at driving pressures of 0.3 MPa and 0.6 MPa. However, when the driving pressure was 1 MPa, the water saturation began to decrease. After the experiment had completed, the final water saturation under 1 MPa was 75% (
Figure 4i). Samples A25 and A45 had a higher permeability, and no start-up pressure gradient was observed during the physical modeling experiments. For a migration layer of about 0.3 mD, the starting pressure was approximately 0.6 MPa (
Figure 4c), the core length was 3.47 cm, and the starting pressure gradient was 0.15 MPa/cm (
Figure 5). Some scholars have also observed a start-up pressure gradient in oil-driven water displacement experiments [
4], where the start-up pressure gradient for oil-driven water is relatively high: generally, for a 0.1 mD core, the start-up pressure was about 0.6–1 MPa/cm [
4], which was affected by the oil viscosity. From this physical modeling experiment, it can be considered that, when the permeability was greater than 1 mD, the start-up pressure gradient could not be observed.
3.1.2. Sample A25 Experimental Results and Discussion of Driving Frontier
Sample A25 (4324.9 m from Well D-1), with a permeability of 1.39 mD and porosity of 8.2%, was tested at only one driving pressure: 0.3 MPa. After 5 min of driving, the water content in the front of the core decreased significantly, and the water content in the center and the end of the core also began to decrease, with a water saturation of 93%. This result indicates that 0.3 MPa was sufficient for charging into a medium-permeability core. After driving for a further 5 min (for a cumulative driving time of 10 min), the water content in the center and the end of the core decreased slowly, while the water content of the front end of the core decreased significantly, with a water saturation of 87%. The water content at the end of the core was similar to that in the original sample, but the water content in the center and the front of the core decreased significantly, reflecting that the natural gas had driven the water at the front of the core to the end of the core. With a cumulative driving time of 25 min, the water content throughout the core was equal, indicating that the natural gas had completely infiltrated the core, displacing most of the movable fluid; the water saturation was 74% at this time. Continuing to use 0.3 MPa for driving, the water content throughout the whole core decreased, indicating that natural gas could still carry fluid away from the core. After a cumulative driving time of 50 min, the water saturation was 71%.
The simulation experiments demonstrated that the charging front of the migration layer was basically flat (i.e., planer; see
Figure 4d–f and
Figure 6b–d). Some scholars believe that charging occurs in a branch-like manner [
23,
24,
25]; however, the core displacement imaging experiment cannot prove whether the charging occurs in a branch-like way or not. The results support the idea that the driving front is flat, with gradual changes in water saturation. When sample A10 was driven for 60 min at a driving pressure of 1 MPa, the core was still divided into three sections: a complete displacement section (the left third of
Figure 4f), the driving front (the middle third of
Figure 4f), and an unchanged section (the right third of
Figure 4f). Here, the core length was 3.36 cm and the driving front was about 1–1.5 cm. Similarly, when sample A25 was driven for 5 min under a pressure of 0.3 MPa, according to the MRI results, the core could also be divided into three sections: A complete displacement section (
Figure 6b, left third), the driving front (
Figure 6b, middle third), and an unchanged section (
Figure 6 b, right third). Here, the core length was 3.47 cm and the driving front was, again, about 1–1.5 cm. The physical modeling experiments demonstrated that the core driving front was about 1–1.5 cm in both cases. It can be considered that the gas-driven water displacement is approximately planar (
Figure 7). The completely displaced area is mainly dominated by the residual water of the water film, the driving front is dominated by a gradual decrease in movable water, and the unchanged section is basically full of water. In the physical simulation, due to the use in small core sections, no obvious rock heterogeneity affected the charging of the migration layer.
3.1.3. Sample A45 Experimental Results and Discussion of Carrying Effect
Sample A45 (4106.86 m from Well F-1) was displaced with a driving pressure of up to 3 MPa. After using 0.3 MPa driving for 1 min, the residual water content of the core was evenly distributed and the water saturation decreased to 86%, reflecting that natural gas quickly infiltrated the core, thus displacing the movable water. Continuing to use 0.3 MPa pressure for another 1, 2, and 11 min (for a cumulative driving time of 15 min), the core water saturation was 76% and the water saturation was equal and decreased slowly, reflecting the carrying effect as the main displacement method. Subsequently, the driving pressure was gradually increased to 3 MPa. For simplicity of calculation, 0.3 MPa driving was carried out for 10 min, 0.5 MPa driving for 5 min, 0.8 MPa driving for 10 min, 1 MPa driving for 5 min, 2 MPa driving for 10 min, and 3 MPa driving for 5 min. The water saturation decreased by 4%, 1%, 4%, 2%, 7%, and 2% in these stages, respectively, indicating that the decrease in water saturation caused by the carrying effect was related to the gas flow rate. At the same time, high driving pressures can displace water from smaller pores, but the driving effect was not directly observable.
After the natural gas completely infiltrated the entire core, while the displacement pressure remained unchanged, the water saturation continued to decrease (
Figure 6e–g and
Figure 8b–d). The experimental results prove that the gas charging of the migration layer was influenced by the carrying effect by the flooding gas, leading to the water content decreasing throughout the entire core. Representative images are those for Sample A25 in
Figure 6e–g. Here, with a core driving pressure of 0.3 MPa, the water saturation of the core dropped from 76% to 70%. The MRI results show that the water content dropped uniformly in the core. Similarly, the images for Sample A45 in
Figure 8b–d were obtained with a driving pressure of 0.3 MPa, showing a water saturation drop from 86% to 76%; again, the MRI results indicate that the water content in the core decreased uniformly.
After completing MRI testing of the three samples at different driving times, NMR T2 spectrum tests were also carried out. It can be seen from the T2 spectra that the T2 relaxation time of the movable fluid in the core was mainly 10–500 ms, while that of the bound water was mainly <10 ms. The movable fluid gradually decreased during the driving process, but was not completely displaced. There was no change in bound water. Under low-pressure driving for 60–90 min, the water content of low-permeability cores generally decreased by less than 25%, the water content of medium-permeability cores decreased by about 30%, and the water content of the high-permeability cores decreased by around 40%. The NMR T2 spectra also revealed that a small part of the movable fluid tended to remain in the core. According to the equal distribution of water content observed in the MRI analysis, the remaining movable fluid decreased through the carrying effect of flooding gas.
According to the analysis of the modeling experiments, the carrying effect should not be ignored, which can reduce the water saturation by about 10% (
Figure 9d). According to the water saturation change plot in
Figure 9d, when the driving time was 20 min for Sample A25, the natural gas had completely infiltrated the core and the water saturation was 76% at this time. As the driving continued to 50 min, the core water saturation slowly dropped to 70%. According to the trend in water saturation, there was still room for it to further decline. Under the condition of the driving pressure of 0.3 MPa, the final water saturation can be expected to be around 60–65%; that is, the carrying effect can reduce the water saturation by about 10–15%.
After the natural gas had infiltrated the entire core, the driving pressure continued to increase and the water saturation decreased (
Figure 8f–k). When the core driving pressure increased from 0.5 MPa to 2 MPa, the water saturation dropped from 76% to 60%, with an associated reduction of 16%, and the water content of the core decreased uniformly throughout the whole core. Our analysis suggests that, in the displacement process, as the driving pressure increases, there is one effect through which the water in smaller pores is driven by high-pressure gas, and another effect through which flowing gas carries the water out from core. The carrying effect generally reduced the water saturation by about 10% in the analysis of Sample A25 (
Figure 9d) and, so, an increase in the driving pressure could produce a water saturation drop of about 5%. It can be seen that, when the natural gas completely infiltrated the core, the water saturation reduction produced by the carrying effect was about two times greater than that due to the pressure increase.
3.2. Slow-Charging Results and Discussion
3.2.1. SPB Experiment Results
The main feature of trap charging is that the natural gas will not go away and only water can discharge, with the semi-permeable baffle (SPB) model being the best experimental fit for this type of charging. We completed an SPB charging experiment using six core samples, in which the charging pressure was 0–1.4 MPa and the balance time at each pressure was 4–6 h.
After analyzing the six samples, we found that the water saturation change was almost unaffected by the physical properties of the rock (
Figure 10), as the water saturation of different cores decreased in basically the same way with an increase in charging pressure; however, the physical properties of the rock did affect the final water saturation. Sample A1 (3449 m, Well A-1) had a permeability of 0.1 mD, a porosity of 7.6%, and a water saturation of 55% at a charging pressure of 1 MPa. When the charging pressure increased again, the water saturation did not decrease and remained stable at 50%. Samples A20 (4008 m, Well B-2) and A26 (5106.9 m, Well D-1) had permeability of 0.29 mD and 0.37 mD, and porosity of 7.2% and 8.5%, respectively. At a charging pressure higher than 1.2 MPa, the water saturation was basically unchanged and remained at 45%. Samples A4 (3823.1 m, Well A-1), A7 (3614.9 m, Well A-2), and A17 (3752.7 m, Well B-2) had aa permeability of 1.9 mD, 8.1 mD, and 2.6 mD, respectively, and the water saturation at above 1.1 MPa was basically unchanged, being stable at about 40%. Overall, the final water saturation of high porosity and permeability cores tended to remain at 40%, while that for low permeability cores basically remained at around 45–50%.
3.2.2. Time Factor
The permeability of the reservoirs in the study area approximately ranges between 0.1–10 mD, and there is obvious gas−water contact. The charging pressure of a gas trap can be calculated according to the gas–water density difference:
where P is the trap charging power (MPa),
is the density of the formation water (g/cm
3),
is the density of natural gas under formation conditions (g/cm
3), g denotes gravitational acceleration (N/m), and Δ
h is the distance between the gas–water contact and the measuring point (km).
According to the formation conditions of the study area, the natural gas density is generally around 0.2 g/cm3. The typical height of the gas column of the main gas reservoirs is 150–200 m. The maximum charging pressure at the reservoir top is about 1.3–1.6 MPa, while the max SPB charging pressure is 1.4 MPa, which is similar to the maximum charging power of the reservoir in the study area. Even though the equilibrium time for the semi-permeable charging model was very long, there is still a huge gap compared with the historical period of the charging geology.
The gas–water contact drilled in Well A-4 is at about 3550 m and is obvious; furthermore, the water saturation at about 3510 m is 40–55% according to the measurements of the sealed cores (
Table 2). Considering the charging pressure cause by the difference in density between gas and water, the charging pressure in a small area above the gas–water contact is very small (only about 0.2 MPa). The water saturation in the SPB experimental model was 80–90%, which is very different to that observed in the sealed core test. This indicates that the time factor is the most important factor to control the water saturation of the reservoir. Conversely, the pressure factor in the modeling experiment had an almost negligible influence on the water saturation.
Considering the time factor, even if the charging pressure is low, the reservoir can still reach a high gas saturation. Therefore, for traps with relatively good physical properties, the water saturation of the reservoir above the gas–water contact can reach a minimal value. The conclusion that the charging pressure affects the water saturation in the experimental model is not suitable when considering the slow charging of traps in the process of geological accumulation. From the experimental results and analysis of sealed cores, slow trap accumulation is only affected by the time factor.
3.2.3. Final Water Saturation
Even though the SPB charging experiment model cannot truly describe the trap accumulation process, the minimum water saturation obtained in the SPB experiment matched the sealed core data well.
According to the SPB charging experiment, when the charging pressure is about 1.2 MPa, the water saturation of the reservoir is basically stable, an increase in the charging pressure will no longer cause the water saturation to decrease (
Figure 8a–c), and the water saturation will basically be maintained at about 40–50%. This indicates that a water saturation of 40–50% is the lowest water saturation that the reservoir can reach; that is, the maximum gas saturation in the study area is 50–60%. The water saturation for the sealed cores from well A-4 was in the range of 40–55% (
Table 2), confirming the final water saturation of 40–50% obtained through SPB experimental modeling (
Figure 11).
The water saturation data of the sealed cores was negatively correlated with the rock physical properties, but the correlation was poor; this also indicates that the physical properties of the rock have little influence on the water saturation of the reservoir, being more greatly affected by the charging time.
3.3. Rapid-Charging Trap Result and Discussion
3.3.1. Cen-NMR Experiment Result
We selected seven samples to carry out centrifugal experiments at 1500, 3000, 4000, 5000, 6000, 6500, and 7000 rpm. After each centrifugation step was completed, the NMR T2 spectrum was recorded (
Figure 12a–c).
In order to establish the relationship between centrifugal force and water saturation (
Figure 12d), it was necessary to convert the rotational speed into the pressure at the bottom of the cores. According to the speed, the length of the centrifugal arm, and the length of the core, the centrifugal speed can be converted into the pressure at the bottom of the core as follows:
where
Pc is the centrifugal charging pressure (MPa), Δρ is the difference in density between the two fluids (g/cm
3), L is the core length (cm),
Re is the maximum radius of the core (cm), and
n is the centrifugal rotation speed (rpm).
Converting the centrifugal speed to the pressure at the bottom of the core, it was found that A7, A25, and A45 have permeabilities greater than 2.5 mD. For these samples, when the centrifugal pressure was less than 0.6 MPa, water saturation rapidly decreased; meanwhile, when the pressure was above 0.6 MPa, the water saturation decreased slowly, and the water saturation was considered to be stable at 15–30%. Samples A5, A10, and A26 had permeabilities in the range of 0.3–1 mD. For these samples, when the pressure was higher than 1 MPa, the water saturation decreased slowly. According to the recorded curve, we considered the water saturation as stable at about 30%. The permeability of Sample A1 was 0.1 mD. Under the experimental conditions, there was no inflection point for the decrease in water saturation, and the residual water saturation at the same pressure was higher than that in the other samples. This reflects the fact that the inflection point of the decrease in water is influenced by the core permeability; namely, when the physical properties are better, the inflection point is lower.
3.3.2. Rapid Trap Charging Characteristics
Centrifugal modeling is most suitable for simulating the short-term gas charging process at the initial stage of episodic charging, while mid–late episodic charging is more suitably simulated using the semi-permeable baffle charging model. The NMR curves obtained after different centrifugal forces (
Figure 9) indicate that cores with permeability less than 0.1 mD had almost no large pores, cores with permeability around 0.5 mD had a relatively similar water content in large pores and small pores, and cores with permeability greater than 1 mD had a very high water content in large pores. After centrifugation, the water saturation and permeability presented a good relationship: when the permeability was higher, the water saturation was lower after centrifugation (
Figure 13). The results of the Cen-NMR experiments demonstrate that short-term rapid charging is completely affected by the physical properties of the rock (
Figure 13).
The water saturation decreased quickly when the centrifugal force was lower than 0.8 MPa, and decreased slowly when it was above 0.8 MPa (
Figure 12d). When the centrifugal force was above 2 MPa it became difficult to judge the change trend in water saturation due to experimental limitations; however, according to the curve extension for Samples A7 and A45, it can be seen that the water saturation of samples with a higher permeability will not change as the pressure increases. For samples with lower permeability, the water saturation will still decrease as the centrifugal force increases. This indicates that, when considering the rapid charging of a trap, it can be completely filled at about 2 MPa if the reservoir has good physical properties.
3.4. Final Water Saturation with Different Physical Models
The final water saturation of cores with different physical properties differed with the models used in the experiments. Only Sample A1 had permeability less than 0.1 mD. The semi-permeable baffle (SPB) and the Cen-NMR experiments were completed with this sample, but it was not submitted to MRI analysis. Under the same driving pressure for this ultra-low-permeability sample, the water saturation in SPB modeling was low, while that obtained in the Cen-NMR experiment was higher (
Figure 14).
Three sample had a permeability in the range of 0.1–0.5 mD; namely, A10, A20, and A26. Among them, Sample A10 completed the GWD-MRI and Cen-NMR experiments, while Sample A20, with similar properties, completed the SPB experiment. The final water saturation in the SPB experiment for the low-permeability samples was the lowest, while that for Cen-NMR was medium, and GWD-MRI gave the highest values.
The cores with permeability of 0.5–2 mD were represented by Samples A4 and A25. Sample A25 completed GWD-NMR and Cen-NMR experiments, while Sample A4, with similar permeability, completed the SPB experiment. The water saturation values obtained in the SPB, Cen-NMR, and GWD-MRI experiments were relatively similar for these medium-permeability samples.
The samples with permeability greater 2 mD were Samples A7, A17, and A45. Among them, Sample A45 completed the GWD-MRI and Cen-NMR experiment, while Sample A7, with similar permeability, completed the SPB experiment. For these samples, GWD-MRI yielded the highest water saturation, while the final water saturation in the SPB and Cen-NMR results was relatively close.
The water saturation values obtained in the displacement imaging experiments were generally 10% higher than those of the semi-permeable baffle and centrifugal-NMR experiments under the same driving pressure (
Figure 10b,d), indicating that the gas charging of the migration layer was weaker than that of the reservoir within the trap. For samples with permeability greater than 0.1 mD, the semi-permeable baffle and centrifugal nuclear magnetic water saturation values were relatively similar (
Figure 10b–d), indicating that trap charging leads to similar water saturation when the drive pressure is low (when not taking the charging time into consideration).