ORIGINAL RESEARCH article

Front. Earth Sci., 04 March 2026

Sec. Georeservoirs

Volume 14 - 2026 | https://doi.org/10.3389/feart.2026.1669602

The influence of gas injection and production rates on pore-scale gas-water movement in low-permeability heterogeneous gas storage based on microfluidic technique

  • 1. Exploration and Development Research Institute, PetroChina Changqing Oilfield Company, Xi’an, Shaanxi, China

  • 2. National Engineering Laboratory of Low-permeability Oil & Gas Exploration and Development, Xi’an, Shaanxi, China

  • 3. School of Human Settlement and Civil Engineering, Xi’an Jiaotong University, Xi’an, China

Abstract

Low permeability heterogeneous gas storage is one of the important types of gas storage currently, and clarifying its microscopic gas water movement laws is crucial for improving the operational efficiency of gas storage. This study investigated the impact mechanism of gas injection and production rates on the gas-water movement behaviors at the pore-scale in low-permeability heterogeneous gas storages. To achieve this, a series of microfluidics experiments was conducted using large-scale chips under different injection and production rate conditions. The dynamic evolution of gas-water distributions, migration pathways, and gas-bearing pore spaces was recorded and quantitatively analyzed. The results indicate that slow injection combined with fast production significantly improved the utilization efficiency of gas-bearing pore spaces. Slow injection promoted more uniform gas displacement and reduced fingering, while the fast production suppressed water occupation in pore throats and mitigated water-locking effects. Conversely, fast injection and slow production led to rapid early gas invasion but increased the formation of water-locked capillary valves, resulted in a higher proportion of bound pore-throats and reduced effective gas storage spaces during subsequent rounds. Microscopic observations showed that repeated injection–production processes progressively increase water-locked throats, limiting further gas migration. Parameter analysis based on capillary number demonstrated that the slow injection/fast production strategies favored the expansion of dynamic gas-bearing boundaries and lowered bound-throat ratios. Furthermore, pore-scale transport mechanisms are discussed using single pore-throat models, highlighted the roles of wettability, corner flow, and capillary resistance in controlling gas advance and retreat behaviors. The findings providing valuable insights into optimizing injection-production strategies for low-permeability heterogeneous gas storages, and offering theoretical support for improving gas storage efficiency and reducing water-locking effects.

1 Introduction

With the continuous growth of global energy demands and the requirements for clean energy transformation, natural gas as a relatively clean transitional energy source, has been increasingly valued for its development and utilization (; ; ). However, the relationship between supply and demand of natural gas is easily affected by seasons, that is, supply exceeds demand in summer, while opposite in winter (). Gas storage is an important way to regulate this relationship (). As a typical gas storage, low-permeability heterogeneous gas storages have complex pore structures, and the gas-water migration characteristics are affected by many factors (; ). Therefore, their gas-water movement laws at the pore-scale are complex. Although there are many factors that can affect the injection and production efficiency of gas storage, for the gas storage site, the parameter that can be directly and easily changed is the injection and production speed, which can be achieved by controlling the pipeline pressure of the injection station. Hence, study the impact of gas injection and production rates on the gas-water movement laws at the pore-scale in low-permeability heterogeneous gas storage reservoirs is of great significance for optimizing the operating parameters and improving the efficiency (; ). Due to the low physical properties and strong heterogeneity of the gas storages, the migration law of gas is complex (; Xu et al., 2025; ). During the gas injection processes, as the pressure increased, the porosity increased due to pressure sensitivity effects, and the gas migration mode gradually shift from Knudsen diffusion to viscous flow (; Zhao et al., 2024; ). During the gas production processes, the decrease in gas pressure will cause the diffusion and release of dissolved gas in the pores, further affected the gas-water movement law (; Zhao et al., 2021; ). To sum up, it is urgent to carry out the impact of injection and production rates on efficiency and gas-bearing spaces, in order to provide a theoretical basis for the efficient development of gas storage.

Low permeability heterogeneous gas storage is transformed from low permeability reservoir, so the research results of low permeability reservoir are also applicable to this type of gas storage, which also opens a new idea for literature research (Zhang et al., 2025; ). Currently, scholars have conducted extensive research on the gas-water movement law in low-permeability heterogeneous reservoirs and have made significant progress (; ; ; ). For example, in terms of revealing fluid distributions at the pore-scale, previous studies have used microscopic testing techniques, such as scanning electron microscopy and nuclear magnetic resonance, to analyze the microscopic pore structures and clarified the initial distribution states of fluid (; ; ). In order to explore the impact of gas injection and production rates on the fluid movement behaviors, scholars have designed experiments under different gas injection and production rate conditions to monitor the dynamic flow processes of fluid in the pores, and analyzed the impact of rate changes on fluid distribution, migration path, and spatial distributions (Zhang et al., 2022; ; ). Regarding the fluid-solid interactions in the process of gas-water flow, some scholars have also utilized comprehensive experimental and numerical simulation methods to clarify the dynamic issues involved (; ; ). In terms of establishing in-situ models to study fluid movement laws, previous studies have constructed models represented the low-permeability heterogeneous gas storages based on the experimental data and theoretical analysis, providing theoretical support for optimizing the operation of the reservoir (; ; ). Some scholars have also paid attention to the fracture propagation prediction under strong injection and production conditions (; Xia et al., 2024; ; ). These studies provide references for the mechanism of gas-water flow in low permeability heterogeneous gas storages.

In recent years, with the continuous popularization of new energy business, researchers have paid attention to the research of gas storage for other gases, especially hydrogen (; ; ; Wang H. et al., 2025). At the same time, in order to cope with global climate change, the research on carbon dioxide storage has also attracted much attention (; ; Zhou et al., 2017; ). These microscopic results can also provide a basis for the study of gas water relationship in gas storage. For example, some researchers used microfluidic methods to study the gas-water distribution mechanisms of hydrogen in the process of repeated injection production cycles from both qualitative and quantitative aspects, and analyzed the influence of pore-scale hydrogen storage capacity, gas water distribution and velocity on the utilization rate of gas-bearing spaces (; ; ; ). Other researchers have studied pore-scale flow behavior and hysteresis in underground hydrogen storage in sandstones, and studied the impact of both injection and production processes on the gas-water relationship, especially the gas storage capacity under a variety of gas water conditions, such as hydrogen/brine system and carbon dioxide/brine system (; ; Wang S. et al., 2025). Other researches explained the microscopic mechanism of fluid migration by studied the dynamic changes of contact angle at different time points during injection and production (; ; ). The research on the microscopic flow of these other gases provided methods and theoretical references for the fluid flow mechanism of low-permeability heterogeneous gas storages.

However, although there are many researches on the microscopic movement of gas and water at present, most studies focused on the impact of permeability changes, gas migration pathways, and reservoir heterogeneity on production performance characteristics at the macroscale (). The studies of the fluid distributions at the pore-scale, and the impact of gas injection and production rates are not enough yet. In addition, existing gas injection and production experiments of gas storages were mostly focused on salt cavern gas storage facilities, with relatively less researches on low-permeability heterogeneous gas storages. By the way, most existing experiments are conducted at low pressures, which greatly differ from real formation pressure conditions, making it difficult to reflect the characteristics of gas storages under real conditions (Zhang et al., 2016). Moreover, most of the microfluidic chips’ researches on the laws of gas and water movement is realized by using small chips at present, and there are relatively few large chips with more than 10 cm in size (). Based on this, this study used 20 cm × 20 cm large-scale chips to study the gas and water movement laws of gas storage. Firstly, experimental schemes under different gas injection and production rates are designed, and microfluidic gas reservoir chip-related devices are used to simulate the gas injection and production process under real conditions. Then, the dynamic flow process of gas and water in the pores were recorded, and the distribution, migration path, and migration rates of fluid were analyzed. Finally, combined with experimental data, using fluid mechanics, percolation mechanics and other theories, the impact mechanism of gas injection and production rates on the fluid movement laws were analyzed, and the pore-scale gas-water movement law models were constructed. This research hopes to provide new theoretical insights and technical support for the efficient and stable operation of low-permeability heterogeneous gas storages.

2 Methodology

This experiment used the microfluidic chip made of glass material. Compared with rock chips, microfluidic chip manufacturing processes have stronger controllability, better sealing, and relatively simpler processes. More importantly, the use of microfluidic chips enables repetitive experiments, providing a fundamental guarantee for injection and production experiments under different conditions. At the same time, the main mineral of sandstone in gas reservoirs is quartz, which has the same composition as glass, making microfluidic chips more representative for conducting experiments. In summary, this experiment conducted research using microfluidic chip, and the results obtained have representativeness.

2.1 Chip manufacturing processes

  • Before making chip, it is necessary to obtain the reservoirs’ pore structures. The pore structure of the chip is taken from real sandstone cores by CT methods, that is, obtain pore channels of different slices through X-ray method, and the pore distributions were generated through continuous flattening and iteration. The results of CT scan indicated that the radius of the channel ranges from 50 μm to 1,000 μm, with a throat radius of 50 μm and a pore radius of 100–1,000 μm. The channel of the chip is strongly hydrophilic, which conformed to the wettability of most sandstones.

  • Repeatedly concatenated the pore channels of different slices to obtain a pore size of 15 cm × 15 cm, and fabricated microfluidic chip on a 20 cm × 20 cm glass. Compared with traditional<10 cm microfluidic chips, large-sized microfluidic chips have stronger representativeness. The production of microfluidic chip for experiments is the foundation of all the experiments. This study used laser etching, ultrasonic cleaning, oven drying, and high-temperature bonding fabricated microfluidic chip.

  • After the chip production was completed, the pore channels were rinsed with alcohol to remove residues. After that, the pore channels were rinsed with deionized water for 10 min. Then, all fluids were sucked out, and the chip was dried using a constant pressure pump. After drying at 100 °C in a drying oven, the chip has been produced and completed.

2.2 Experimental process of microfluidic chip

  • The initial microfluidic chip consists of dry pore channels, and the first step requires saturated water. The chip was saturated with water, and the specific steps were consistent with the injection production experiment. Finally, we connected the gas injection port to a constant flow pump to ensure that the gas was injected at a constant rate. At the same time, the discharge outlet was sealed with water to simulate a closed geological environment, prevented gas leakage during the experimental processes, and ensure the accuracy of experimental results and analysis.

  • In the first stage of the experiment, we conducted a round of slow injection and fast production experiments. The gas injection rate is 10 μL/min, the gas production rate is 20 μL/min, and the total gas injection amount was 350 μL. This step simulated the process of natural gas accumulation and production. By injected gas slowly and collected it quickly, we created an initial water and gas-containing state, which was the state before the depleted gas reservoirs were transformed into gas storages.

  • Subsequently, a new round of slow injection and fast production experiments were conducted, with a gas injection rate of 10 μL/min and a gas production rate of 20 μL/min. The total injection volume was 350 μL, with a total of four rounds of injection. This parameter simulates the scenario of real gas storage injection and production, that is, long-term gas injection in summer and short-term gas extraction in winter.

  • After the experiment was completed, repeat step 1) to clean the chip, and follow step 3) to create an initial environment. Then, conduct a fast injection and slow production experiment with an injection rate of 20 μL/min and a production rate of 10 μL/min. The total injection volume was 350 μL, with a total of four rounds of injection. This experiment was a parameter reversal experiment of step 4), aimed at studying the influence of different injection and production rates on the gas-containing spaces of gas storages. The use of this method for comparison was also to compare with the current conditions of the gas storages, providing theoretical support for the reasonable setting of injection and production rates.

  • To further investigate the influence of injection and production rate on gas-water distributions, we recorded in detail the state of gas-water distributions in each round of experiments, including the diffusion range of gas and the displacement of water. Subsequently, we used image processing techniques to process the recorded data in order to observe the changes in gas-water distribution more intuitively. Finally, we conducted an in-depth analysis of the experimental data and compared the dynamic boundaries between slow injection and fast production, as well as fast injection and slow production, in order to reveal the specific impact of injection and production rate on the distribution of fluids.

3 Results and discussion

3.1 Determination of injection and production rounds

Clarifying the impact of injection and production rounds on the gas-containing spaces of gas storage can improve the experimental testing efficiency of this study, that is, to determine the stable operation of the gas storage after several rounds of injection and production, in order to achieve the best experimental results with the least number of experiments. Therefore, priority tests were given with lower injection and production volumes (100 μL). The results of the gas reservoir chip experiment indicate that with the gradual increase of injection and production rounds, the gas storage gas-containing spaces showed gradually expanding trends (Figure 1). This phenomenon indicated that during continuous gas injection and production processes, the permeation and diffusion range of gas in the gas reservoir chip expanded, gradually occupied more pore space. After four production rounds, the gas-containing pore space reached stable state, which means that in the subsequent injection production processes, the changes in the gas-containing pore spaces remained nearly the same, indicated that after a certain number of operations, the spatial distribution of the gas storage gradually reached dynamic equilibrium.

FIGURE 1

In addition to qualitatively identified changes in gas-bearing spaces, we also quantitatively calculated the gas bearing space change rate under different injection production rounds (Figure 2). The statistical data showed that after four rounds of injection and production, the increase in gas-containing pore spaces were less than 2%, and the subsequent rounds below to 1.1%, which met the principle of stable operation of the gas storage if it is below 3% on site. Therefore, in the subsequent experiments, we conducted four rounds of injection and production.

FIGURE 2

3.2 Chip characteristics at different injection and production rates

The gas reservoir chip images showed the brine-gas distributions under the conditions of slow injection and fast production from the injection beginning and production in the end at different rounds (Figure 3). The images illustrated that with the increased of injection and production rounds, the gas-containing space was gradually expanded. The expansion of the brine-gas boundary was obvious in the first two rounds, and then the expansion rate decreased. The gas-containing space is basically fixed, which is consistent with the experimental results that Section 3.1 have mentioned. This phenomenon clearly indicated that during the injection production processes, the gas-containing spaces were constantly increased, and the gas can effectively penetrated and occupied wider areas, demonstrated good fluidity and diffusion of the fluids. These operations are conducive to the uniform distribution and effective displacement of gas in the pores, thus improved the utilization efficiency of the pore spaces.

FIGURE 3

The fast injection and slow production experiments exhibited different characteristics. The experimental results showed that under the condition of fast injection and slow production, with the increase of injection and production rounds, the gas-containing spaces were relatively high at the beginning, but the growth rates were low (Figure 4). It was illustrated from the microscopic images that under the condition of fast injection, the gas moved rapidly, resulted in the rapid expansion of gas bearing spaces. However, under the condition of slow recovery, the wetting phase quickly occupied the tiny throats, resulted in the increase of capillary valves and serious water-lock, led to the limited gas injection efficiency subsequently. Therefore, fast injection and slow production are not conducive to the expansion of gas bearing space qualitatively, although gas could occupy the more pore spaces initially, the pore spaces could not be utilized efficiently during slow production.

FIGURE 4

Comparing the gas storage boundaries can clarify the qualitative change of the gas storage gas-containing spaces under different injection and production conditions. The figures revealed that the expansion of the boundaries is faster in the early stage of injection and production (Figure 5A). These phenomena resulted from the gas-water boundaries rapid expansions due to fast gas injection. At the second to fourth rounds of injection and production, the expansion rate of the gas-containing spaces of slow injection and fast production was significantly higher (Figure 5B). These phenomena indicated that during the rapid injection of gas, although the storage space increased, the diffusion and permeation rate of gas relatively slow, resulted in limited expansion of the storage gas-containing spaces. In contrast, the slow injection and fast production strategy were more conducive to the widespread distribution of gas and the effective expansion of gas storage dynamic boundaries.

FIGURE 5

3.3 Analysis of microscopic images and flow parameters

Microscopic images indicate that as the round of injection and production increased, the capillary valve caused by water-lock became more difficult to broken through. The gas storage space after the second round of injection and production became bound pore after the third round of injection and production (Figures 6A,B). This is because the water in the throat formed capillary valve, based on the capillary pressure formula, the breakthrough pressure required for the capillary valve is relatively high, resulted to the difficult for gas to enter the pore through the throat. Due to the existence of water-locked capillary valve, the proportion of invalid pores gradually increased. These invalid pores are pores that are occupied by water, making it difficult for gas to enter. After continuous gas injection, the actual effective gas-bearing space is not fully utilized due to the increasing proportion of invalid pores. With the increased of injection production rounds, the reason for the increase of invalid pores is that within the scope of effective gas storage space, more and more throats were filled with water, resulted in the increase of capillary valves (Figures 6C,D). Therefore, the increase in gas-bearing space after gas injection mainly occurs in areas that have not yet been affected by water-locking effects.

FIGURE 6

Capillary number (NCa) is an important parameter to characterize two-phase flow, which is the ratio of viscous force to capillary force. NCa is usually defined as NCa = Uμ/σ = Qμ/σLdφ, where U is the velocity, Q is the injection/production rate, μ is the injected fluid viscosity, μN2 = 3.3 × 10−5 Pa·s, μH2O = 1.0 × 10−3 Pa·s (; ). σ is the interfacial tension, and the value is 0.022 N/m in this research, because we define the system as nitrogen-water system here (). L, d, and φ are microchips’ length, depth, and porosity, respectively. Table 1 showed the relevant parameters and the calculation results. In the process of gas injection, the lower capillary number led to uniform displacement, uniform expansion of gas-bearing area and less fingering of gas. In the process of gas production, the higher capillary number resulted to the rapid recovery of non-wetting phase, which made a higher final recovery under the condition of slow injection and rapid production.

TABLE 1

No.Qi (μL/min)Qw (μL/min)NCaCycle numberInitial gas saturationsResidual gas saturations
InjectionProduction
15103.8 × 10−82.3 × 10−610.410.05
20.510.07
30.570.11
40.620.16
21057.4 × 10−81.1 × 10−610.450.14
20.500.18
30.520.20
40.530.25

Flow parameter calculation results.

Qi-Injection rate, Qp-production rate, NCa-Capillary number. No.1 represent slow injection and fast production; No.2 represent fast injection and slow production.

3.4 Parameter characteristics at different injection and production rates

The comparison of data from gas-containing spaces experiments under different injection and production rates showed that the expansion of gas storage dynamic boundaries was pronounced under slow injection and fast production conditions (Figure 7). Specifically, compared with fast injection and slow production, the increase in gas-bearing spaces under the slow injection and fast production strategy is larger, indicated that gas could penetrate and occupied the reservoir effectively and then formed wider gas-bearing areas. Meanwhile, in terms of total volume, the gas content under slow injection and fast production conditions was relatively high, which further confirmed the advantages of this strategy in expanding the gas-containing spaces of gas storage and improving gas storages efficiency. This discovery is of great significance for guiding the selection of injection and production strategies in the operation of actual gas storages. It suggested that while pursuing injection and production efficiency, we should consider the influence of injection and production rates on the gas-containing spaces and the gas storages. At the same time, it is worth noting that as the gas-containing spaces expanded, the proportion of confined water pores-throats in the gas storage area (i.e., the coordination pores and throats of single or continuous water-containing pores and throats are all gas) gradually increased (Figure 8).

FIGURE 7

FIGURE 8

Statistical data also showed that as the number of injection and production rounds increased, the proportion of bound water pores and throats growth, but the proportion of bound water pores and throats formed under slow injection and fast production conditions was relatively low (Figure 9). This difference was caused by capillary action, according to a variant of the capillary force formula:where r is radius, σ is the interfacial tension, θ is the contact angle, and Pc is the capillary pressure. During the gas injection processes, the water phase acted as a capillary resistance, while played a capillary dynamic role at the production rounds. Therefore, under slow injection conditions, the non-wetting phase of the gas was subjected to greater capillary resistance, resulted in uniform propulsion and less water-lock sites. Under fast injection conditions, non-wetting gas can easily break through the capillary blockage, formed discontinuous bubbles, and led to water-lock. According to the conversion of injection and production rates in this experiment, the continuous flow aperture ratio under slow and fast injection conditions was:where Rsi is the critical pore radius under slow injection conditions, Rfi is the critical pore radius under fast injection conditions, vsi and vfi are the injection speed under slow injection and fast injection conditions, respectively. Therefore, the probability of forming a water lock under fast injection conditions was approximately twice that under slow injection conditions. At the same time, rapid production resulted in a higher gas flow velocity, and the ability of water in the wetting phase to occupy space was relatively weak. Due to the inability to measure the gas flow velocity in the pores, it cannot be quantified here, but the trend remained stable. Therefore, under the dual effect of slow injection and fast production, there are relatively more water-lock points, which was also shown in Figure 9. High-rate and high-pressure gas injection can lead to discontinuous gas paths and the formation of water-locks. However, the pore medium was a wetting phase, and the viscosity of water was much higher than that of gas. Therefore, high-rate gas production caused rapid extraction, while the water phase at the gas-water boundaries occupied the throats. Slow gas production caused the water phase retreated and occupied the throats for a longer period, further led to more water-lock points. Therefore, the influence of injection and production rates on the gas-containing spaces is significant, and the water-lock caused by different injection and production rates was the main factor affecting the gas-containing spaces.

FIGURE 9

The research has proven that slow injection and fast production could effectively improve gas storage spaces. To further clarify the efficiency of the impact of injection and production rates on gas storage boundaries, this study compared the changes in gas-containing pore spaces at different time periods under slow injection and fast production conditions. The gas storage dynamic boundaries, which indicate by the utilization of gas-containing pore spaces, showed rapid increase followed by slow increase (Figure 10). One of the reasons was that the continuous increase of gas-containing spaces and the total amount of compressed gas, the increasing number of binding points at the gas-water boundary, where many capillary valves hinder the expansion of the dynamic boundary was another reason. A more noteworthy pattern was that as the number of injection and production rounds increased, the efficiency of reached the maximum gas content efficiency at different time periods continuously decreased. For example, within the same time period (taking 28–35 min as an example), the first round of injection and production reached the maximum capacity quickly, while the fourth round took longer time. This further proved that as the number of injection and production rounds increased, the number of water-lock points in the gas storage spaces increased, and the increasing number of water-lock points corresponds to more capillary valves, resulted in the subsequent injection of gas playing a more significant role in expanding the edges. The extension time corresponded to slow injection and fast production was relatively long. Due to the presence of the starting pressure of gas-water movement, the available pore spaces at the boundary also increased, resulted in a more significant improvement in the gas storage boundaries. The above analysis indicates that slow injection and fast production have significant benefits in improving the dynamic boundary of gas storages. The injection and production rates are important influencing factors of the dynamic boundaries of gas storages.

FIGURE 10

3.5 Analysis of the dynamic mechanism of gas transportation

For a single pore throat, there are two cases: Case 1: θ-β>90°, under gas-water drainage conditions, the gas-water interface stays at the junction of pores and throats (at the throat entrance), where the resistance is the highest. In this case, the gas enters as resistance and retreats as resistance, resulting in a dilemma (Figure 11A). Case 2: θ-β<90°, under gas-water drainage conditions, the gas-water interface stays at the junction of pores and throats (at the throat entrance), where the resistance is the highest. In this case, the gas retreats as resistance and advances as power, making it difficult to advance (Figure 11B).

FIGURE 11

Under strong hydrophilic conditions, due to the strong adsorption of water on the pore surfaces, the water phase tends to protrude along the corners of the rectangle, forming the corner flow (). This phenomenon is due to the preferential path provided by the corner region for the flow of the water phase, allowing the water phase to pass through these areas more quickly (; ). The formation of corner flow can lead to heterogeneous distributions of the water phase in the pores, which may hinder the effective flow and storage of gas, thus adversely affecting the performance of the gas storage. In practical operations, it is necessary to pay attention to the corner flow phenomenon under strong water wet conditions and take corresponding measures to optimize fluid distribution and improve the injection-production efficiency of the gas storage (Figure 12).

FIGURE 12

For the radius of curvature, there is a formula:where r is radius, h is corner flow thickness, and α is the angle see in Figure 12.

We know the capillary pressure formula (; ; ):where σ is the interfacial tension, θ is the contact angle.

Combine the above two formulas to get:

Therefore, strong hydrophilic was corresponding to small the pore angle with thin water film. In this situation, the capillary effect would become more obvious. While in the hydrophobic situation, the gas is displaced along the center, and the pores that the gas reaches are basically driven to be relatively clean (Figure 13).

FIGURE 13

When there is a simultaneous distribution of gas and water, there are also two cases: Case 1: When the pores are filled with water, there is only one meniscus with additional pressure, and the gas-driven water resistance is small (Figure 14A). The pressure is shown as follows:

FIGURE 14

Case 2: When there are bubbles in the pores, there are three menisci, and the gas enters the pore through the throat with a larger curvature, which results in a smaller curvature. The additional pressure generated by the latter is much greater than that of the former, making it difficult for the gas in the pores to be utilized (Figure 14B). Capillary resistance of gas drive:where R represent pore radius. Since R >> r, the above formula can be rewritten as:

The above analysis showed that the wettability difference is one of the microscopic mechanisms that cause the spatial difference of gas content at different injection and production rates. However, at present, only chip with unique wettability has been studied, and we are developing chips with different wettability by means of surface modification.

3.6 Brine movement behavior assessment in liquid-liquid system

In order to better illustrated the reliability of large-scale chips in the study of the law of gas-water movements, we conducted the liquid-liquid system experiment. According to the two key influencing factor indexes obtained from gas-water tests, namely, injection rates and wettability, we carried out water flooding experiments under different injection rates and wettability conditions. The results showed that the capillary fingering was prominent at low speed, and the viscous displacement was enhanced with the increase of injection rate, but too high displacement rate led to more oil trapped by water, resulted in many beads like residual oil (Supplementary Figure S1). With the decrease of hydrophilic degree, the water phase sweep increased, and the water wet intensity is high, which is easier to form a wall film flow, causing water to rush along the wall, forming a thicker water film, and the oil displacement efficiency in a single hole is higher (Supplementary Figure S2).

At the same time, we also statistically analyzed the quantitative data of oil displacement efficiency under different injection production rates. With the increase of injection rate, oil displacement efficiency increased first and then decreased (Supplementary Figure S3A). This is due to the dispersion of oil droplets caused by high-speed water injection and the low recovery of residual oil. In the range of hydrophilic, the contact angle was negatively correlated with oil displacement efficiency, that is, the more hydrophilic the reservoir was, the higher the oil displacement efficiency was (Supplementary Figure S3A). This is consistent with the understanding of gas-water system: the recovery of non-wetting phase is closely related to the degree of affinity between non-wetting phase and wall. These studies supplement our understanding of gas water systems, and also showed that our large-scale chips can be applied to various fluid systems.

4 Conclusion

  • 1.

    Gas-bearing pore space reached dynamic stability after four injection–production cycles. Microfluidic experiments showed that gas-bearing pore space expanded progressively with increasing cycles and reached dynamic equilibrium after approximately four rounds. Quantitative results indicated that the incremental increase after the fourth cycle was <2%, with subsequent growth <1.1%, satisfying the stable operation criterion for gas storage (<3%).

  • 2.

    Slow injection and fast production enhanced gas-bearing space utilization. Compared with fast injection/slow production, the slow injection (10 μL/min) and fast production (20 μL/min) strategy resulted in larger dynamic gas-bearing boundaries and higher effective pore utilization. Slow injection promoted more uniform displacement, while fast production reduced water retention in pore throats, leading to higher growth rates of gas-bearing space from the second to fourth cycles.

  • 3.

    Fast injection and slow production promoted water locking and bound pore–throat formation. Fast injection (20 μL/min) combined with slow production (10 μL/min) caused rapid early gas invasion but significantly increased water-locked capillary valves. Microscopic observations and statistical analysis showed that the probability of water-lock formation under fast injection was approximately twice that under slow injection, severely restricting subsequent gas migration.

  • 4.

    Capillary number and wettability jointly controlled pore-scale gas–water transport. Capillary number analysis showed that lower values during injection favored stable displacement, whereas higher values during production enhanced non-wetting phase recovery. Under strong water-wet conditions, corner flow increased capillary resistance and water locking, requiring higher breakthrough pressures and reducing gas utilization efficiency during later cycles.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author contributions

DW: Conceptualization, Data curation, Writing – original draft, Formal Analysis. WW: Formal Analysis, Writing – original draft. WA: Methodology, Writing – original draft. JS: Data curation, Validation, Writing – review and editing. JZ: Investigation, Methodology, Writing – original draft. DL: Resources, Supervision, Writing – review and editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. The research is co-funded by National Natural Science Foundation of China (No. 242402146) and the PetroChina Company (No. 2023–4025 and No. 2025–6428).

Acknowledgments

The authors appreciate Dazhong Ren from Xi’an Shiyou University, and Huifang Hu, Deshuo Tao, and Li Tao from Sinopec for valuable advice.

Conflict of interest

Authors DW, WW, and WA were employed by Exploration and Development Research Institute, PetroChina Changqing Oilfield Company.

The remaining author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2026.1669602/full#supplementary-material

SUPPLEMENTARY FIGURE S1

Water flooding experiment results with different wettability. (A) Strong hydrophilic; (B) Medium-strong hydrophilic; (C) Strong-medium hydrophilic; (D) Medium hydrophilic.

References

  • 1

    AkpinarM.YumusakN. (2016). Year ahead demand forecast of city natural gas using seasonal time series methods. Energies9 (9), 727. 10.3390/en9090727

  • 2

    AlAbdullahM. B.DelshadM.NiuF.MohantyK.SepehrnooriK.LiJ.et al (2025). Simulation of low-tension-gas flood in a high-temperature and heterogeneous sandstone reservoir. SPE J.30 (04), 2155–2170. 10.2118/218186-pa

  • 3

    AliM.IsahA.YekeenN.HassanpouryouzbandA.SarmadivalehM.OkoroaforE. R.et al (2025). Recent progress in underground hydrogen storage. Energy and Environ. Sci.18 (12), 5740–5810. 10.1039/d4ee04564e

  • 4

    AlizardF.RobinetJ.-C.RistU. (2010). Sensitivity analysis of a streamwise corner flow. Phys. Fluids22 (1), 014103. 10.1063/1.3292009

  • 5

    Amann-HildenbrandA.GhanizadehA.KroossB. M. (2012). Transport properties of unconventional gas systems. Mar. Petroleum Geol.31 (1), 90–99. 10.1016/j.marpetgeo.2011.11.009

  • 6

    BahramiM.MahaniH.ZivarD.AyatollahiS. (2024). Microfluidic investigation of pore-scale flow behavior and hysteresis in underground hydrogen storage in sandstones. J. Energy Storage98, 112959. 10.1016/j.est.2024.112959

  • 7

    BazmiA. A.ZahediG. (2011). Sustainable energy systems: role of optimization modeling techniques in power generation and supply—A review. Renew. Sustainable Energy Reviews15 (8), 3480–3500. 10.1016/j.rser.2011.05.003

  • 8

    BetancurS.QuevedoL.OlmosC. M. (2024). Microfluidic devices, materials, and recent progress for petroleum applications: a review. Can. J. Chem. Eng.102 (7), 2583–2607. 10.1002/cjce.25214

  • 9

    BlasingameT. A. (2008). “The characteristic flow behavior of low-permeability reservoir systems,” SPE unconventional resources conference/gas technology symposium. SPE.

  • 10

    BluntM. J. (1997). Effects of heterogeneity and wetting on relative permeability using pore level modeling. SPE Journal2 (01), 70–87. 10.2118/36762-pa

  • 11

    BoissonJ. Y.BertrandL.HeitzJ.-F.GolvanY. (2001). In situ and laboratory investigations of fluid flow through an argillaceous formation at different scales of space and time, tournemire tunnel, southern France. Hydrogeology J.9 (1), 108–123. 10.1007/s100400000119

  • 12

    BugajeA. A. B.DiohaM. O.Abraham-DukumaM. C.WakilM. (2022). Rethinking the position of natural gas in a low-carbon energy transition. Energy Res. and Soc. Sci.90, 102604. 10.1016/j.erss.2022.102604

  • 13

    CadalenS.LanceM. (2011). A mechanics approach for wet gas flow metering, theory and application to flow loop tests. Int. Journal Multiphase Flow37 (3), 260–267. 10.1016/j.ijmultiphaseflow.2010.09.006

  • 14

    ChangH.DuoT.LiB.WenK.MengT.BaiC.et al (2025). Optimization of the injection-production process of underground gas storage in depleted gas reservoir. Gas Sci. Eng.138, 205605. 10.1016/j.jgsce.2025.205605

  • 15

    Emami-MeybodiH.MaM.ZhangF.RuiZ.RezaeyanA.GhanizadehA.et al (2024). Cyclic gas injection in low-permeability oil reservoirs: progress in modeling and experiments. SPE J.29 (11), 6217–6250. 10.2118/223116-pa

  • 16

    EnabK.Emami-MeybodiH. (2021). Effects of diffusion, adsorption, and hysteresis on huff-n-puff performance in ultratight reservoirs with different fluid types and injection gases. Energies14 (21), 7379. 10.3390/en14217379

  • 17

    FeldmannF.HagemannB.GanzerL.PanfilovM. (2016). Numerical simulation of hydrodynamic and gas mixing processes in underground hydrogen storages. Environ. Earth Sci.75 (16), 1165. 10.1007/s12665-016-5948-z

  • 18

    GaoJ.KongD.PengY.ZhouY.LiuY.ZhuW. (2023). Pore-scale mechanisms and hysteresis effect during multi-cycle injection and production process in underground hydrogen storage reservoir. Energy283, 129007. 10.1016/j.energy.2023.129007

  • 19

    GaoX.YangJ.HanH.PengH.ChenM.GouX. (2024). Experimental study on reservoir particle migration induced by the injection-production process of underground gas storage. Heliyon10 (15), e34989. 10.1016/j.heliyon.2024.e34989

  • 20

    GrenobleZ.TrabelsiS. (2021). Interfacial tension measured at nitrogen–liquid and liquid–liquid interfaces using model microemulsions at high-pressure and high-temperature conditions. Energy and Fuels35 (16), 13055–13064. 10.1021/acs.energyfuels.1c01358

  • 21

    HeX.WangY.ZhengY.ZhangW.DaiY.ZouH. (2025). Pore-scale gas–water two-phase flow mechanisms for underground hydrogen storage: a mini review of theory, experiment, and simulation. Appl. Sci.15 (10), 5657. 10.3390/app15105657

  • 22

    JuanesR.MacMinnC. W.SzulczewskiM. L. (2010). The footprint of the CO2 plume during carbon dioxide storage in saline aquifers: storage efficiency for capillary trapping at the basin scale. Transp. Porous Media82 (1), 19–30. 10.1007/s11242-009-9420-3

  • 23

    KamaliA.GhassemiA. (2018). Analysis of injection-induced shear slip and fracture propagation in geothermal reservoir stimulation. Geothermics76, 93–105. 10.1016/j.geothermics.2018.07.002

  • 24

    KangZ.LiY.WangL.ZhangY.ZhouY.YangD.et al (2024). Characterization of fracture and percolation processes in oil shale under conduction and convection heating modes based on micro-CT technology. J. Anal. Appl. Pyrolysis181, 106601. 10.1016/j.jaap.2024.106601

  • 25

    KoplikJ.BanavarJ. R. (1995). Corner flow in the sliding plate problem. Phys. Fluids7 (12), 3118–3125. 10.1063/1.868619

  • 26

    KoptyugI. V.MatveevA. V.AltobelliS. (2002). NMR studies of hydrocarbon gas flow and dispersion. Appl. Magn. Reson.22 (2), 187–200. 10.1007/bf03166102

  • 27

    KorsonL.Drost-HansenW.MilleroF. J. (1969). Viscosity of water at various temperatures. J. Phys. Chem.73 (1), 34–39. 10.1021/j100721a006

  • 28

    KuchukF.ZhanL.MaS. M.Al-ShahriA. M.RamakrishnanT.AltundasY. B.et al (2010). Determination of in situ two-phase flow properties through downhole fluid movement monitoring. SPE Reserv. Eval. and Eng.13 (04), 575–587. 10.2118/116068-pa

  • 29

    LarsonR.MorrowN. (1981). Effects of sample size on capillary pressures in porous media. Powder Technology30 (2), 123–138. 10.1016/0032-5910(81)80005-8

  • 30

    LiY.LiuC.LiH.ChenS.LuK.ZhangQ.et al (2022). A review on measurement of the dynamic effect in capillary pressure. J. Petroleum Sci. Eng.208, 109672. 10.1016/j.petrol.2021.109672

  • 31

    LiR.YangY.ZhangY.ZhanW.YangJ.ZhouY. (2024). A novel regional‐minima image segmentation method for fluid transport simulations in unresolved rock images. Water Resour. Res.60 (6), e2023WR036855. 10.1029/2023wr036855

  • 32

    LiuY.BergS.JuY.WeiW.KouJ.CaiJ. (2022). Systematic investigation of corner flow impact in forced imbibition. Water Resour. Res.58 (10), e2022WR032402. 10.1029/2022wr032402

  • 33

    LysyyM.ErslandG.FernøM. (2022). Pore-scale dynamics for underground porous media hydrogen storage. Adv. Water Resour.163, 104167. 10.1016/j.advwatres.2022.104167

  • 34

    LysyyM.LiuN.SolstadC. M.FernøM. A.ErslandG. (2023). Microfluidic hydrogen storage capacity and residual trapping during cyclic injections: implications for underground storage. International Journal Hydrogen Energy48 (80), 31294–31304. 10.1016/j.ijhydene.2023.04.253

  • 35

    MasoodiR.PillaiK. (2012). A general formula for capillary suction-pressure in porous media. J. Porous Media15 (8), 775–783. 10.1615/jpormedia.v15.i8.60

  • 36

    MoridisG. J.CollettT. S.Pooladi-DarvishM.HancockS.SantamarinaC.BoswellR.et al (2011). Challenges, uncertainties, and issues facing gas production from gas-hydrate deposits. SPE Reservoir Evaluation and Engineering14 (01), 76–112. 10.2118/131792-pa

  • 37

    NiazS.ManzoorT.PandithA. H. (2015). Hydrogen storage: materials, methods and perspectives. Renew. Sustain. Energy Rev.50, 457–469. 10.1016/j.rser.2015.05.011

  • 38

    NowrouziI.ManshadA. K.MohammadiA. H. (2019). Effects of ions and dissolved carbon dioxide in brine on wettability alteration, contact angle and oil production in smart water and carbonated smart water injection processes in carbonate oil reservoirs. Fuel235, 1039–1051. 10.1016/j.fuel.2018.08.067

  • 39

    PurswaniP.GuiltinanE. J.ChenY.KangQ.MehanaM. Z.NeilC. W.et al (2024). Pore‐scale modeling of carbon dioxide and hydrogen transport during geologic gas storage. Geophys. Research Letters51 (12), e2024GL109216. 10.1029/2024gl109216

  • 40

    QiaoJ.ZengJ.JiangS.FengS.FengX.GuoZ.et al (2019). Heterogeneity of reservoir quality and gas accumulation in tight sandstone reservoirs revealed by pore structure characterization and physical simulation. Fuel253, 1300–1316. 10.1016/j.fuel.2019.05.112

  • 41

    RazavifarM.MukhametdinovaA.NikooeeE.BurukhinA.RezaeiA.CheremisinA.et al (2021). Rock porous structure characterization: a critical assessment of various state-of-the-art techniques. Transp. Porous Media136 (2), 431–456. 10.1007/s11242-020-01518-6

  • 42

    RebouillatS.LiksonovD. (2010). Fluid–structure interaction in partially filled liquid containers: a comparative review of numerical approaches. Comput. and Fluids39 (5), 739–746. 10.1016/j.compfluid.2009.12.010

  • 43

    RingroseP. S.FurreA.-K.GilfillanS. M.KrevorS.LandrøM.LeslieR.et al (2021). Storage of carbon dioxide in saline aquifers: physicochemical processes, key constraints, and scale-up potential. Annu. Review Chemical Biomolecular Engineering12 (1), 471–494. 10.1146/annurev-chembioeng-093020-091447

  • 44

    SafariA.DasN.LanghelleO.RoyJ.AssadiM. (2019). Natural gas: a transition fuel for sustainable energy system transformation?Energy Sci. and Eng.7 (4), 1075–1094. 10.1002/ese3.380

  • 45

    SeibtD.VogelE.BichE.ButtigD.HasselE. (2006). Viscosity measurements on nitrogen. J. Chem. and Eng. Data51 (2), 526–533. 10.1021/je050399c

  • 46

    ShiH.ZhuQ.ChenZ.LiJ.FengD.ZhangS.et al (2023). Pore-scale modeling of water–gas flow in heterogeneous porous media. Phys. Fluids35 (7), 072114. 10.1063/5.0157655

  • 47

    SongH.ZhouY.XieZ.LaoJ.YueM. (2024). A pore-scale simulation of the effect of heterogeneity on underground hydrogen storage. Water16 (22), 3264. 10.3390/w16223264

  • 48

    SunX.YaoY.LiuD.ZhouY. (2018). Investigations of CO2-water wettability of coal: NMR relaxation method. Int. J. Coal Geol.188, 38–50. 10.1016/j.coal.2018.01.015

  • 49

    TangiralaS.ShengJ. J. (2019). Investigation of oil production and flowback in hydraulically-fractured water-wet formations using the Lab-on-a-Chip method. Fuel254, 115543. 10.1016/j.fuel.2019.05.126

  • 50

    TarkowskiR.Uliasz-MisiakB.TarkowskiP. (2021). Storage of hydrogen, natural gas, and carbon dioxide–geological and legal conditions. Int. J. Hydrogen Energy46 (38), 20010–20022. 10.1016/j.ijhydene.2021.03.131

  • 51

    ToosehE. K.JafariA.TeymouriA. (2018). Gas-water-rock interactions and factors affecting gas storage capacity during natural gas storage in a low permeability aquifer. Petroleum Explor. Dev.45 (6), 1123–1128. 10.1016/s1876-3804(18)30115-0

  • 52

    UsmanM. R. (2022). Hydrogen storage methods: review and current status. Renew. Sustain. Energy Rev.167, 112743. 10.1016/j.rser.2022.112743

  • 53

    van RooijenW.HashemiL.BoonM.FarajzadehR.HajibeygiH. (2022). Microfluidics-based analysis of dynamic contact angles relevant for underground hydrogen storage. Adv. Water Resour.164, 104221. 10.1016/j.advwatres.2022.104221

  • 54

    WangY.ChengH. (2023). Advances in microscopic pore structure characterization of fine-grained mudrocks. Energy and Fuels37 (3), 1495–1510. 10.1021/acs.energyfuels.2c03144

  • 55

    WangJ.LiY.XuS.LiC.WangJ.HeR.et al (2022a). Change mechanism of pore structure and filling efficiency during injection production of sandstone underground gas storage. J. Nat. Gas Sci. Eng.97, 104366. 10.1016/j.jngse.2021.104366

  • 56

    WangJ.WangJ.XuS.WuR.LvJ.LiZ.et al (2022b). A novel mode for “three zones” collaborative reconstruction of underground gas storage and its application to large, low-permeability lithologic gas reservoirs. Energy253, 124148. 10.1016/j.energy.2022.124148

  • 57

    WangL.ZhuL.CaoZ.LiuJ.XueY.WangP.et al (2025). Thermo-mechanical degradation and fracture evolution in low-permeability coal subjected to cyclic heating–cryogenic cooling. Phys. Fluids37 (8), 086617. 10.1063/5.0282266

  • 58

    WangH.ZhouY.KartalM.MaL. (2025). “A critical review on transport and mechanical property of well cement for underground hydrogen storage,” SPE offshore Europe conference and exhibition. SPE.

  • 59

    WangS.PanS.TangY.GaoY.WangK. (2025). Experimental study of hydrogen flow behavior in sandstone at the core scale during cyclic hydrogen and water injections. J. Energy Storage123, 116810. 10.1016/j.est.2025.116810

  • 60

    XiaY.YaoM.LiT.YangH.TangCa (2024). Numerical analysis of hydraulic fracture propagation in deep shale reservoir with different injection strategies. J. Rock Mech. Geotechnical Eng.16 (9), 3558–3574. 10.1016/j.jrmge.2024.02.013

  • 61

    XuD.YangJ.FanJ.ZhouY. (2025). A critical review on deep-water oil and gas development equipment research and application: floating platform, mooring cable and pile foundation. Petroleum Sci.Available online at: https://abdn.elsevierpure.com/en/publications/a-critical-review-on-deep-water-oil-and-gas-development-equipment/.

  • 62

    ZhangJ.YanS.YuanD.AliciG.NguyenN.-T.WarkianiM. E.et al (2016). Fundamentals and applications of inertial microfluidics: a review. Lab a Chip16 (1), 10–34. 10.1039/c5lc01159k

  • 63

    ZhangG.YangS.MoC.ZhangZ. (2022). Experimental research on capacity expansion simulation of multi-cycle injection-production in reconstruction of oil reservoir to underground gas storage. J. Energy Storage54, 105222. 10.1016/j.est.2022.105222

  • 64

    ZhangH.ChengJ.ZhangX.MiaoY.JiangH.YuB.et al (2025). An evaluation of the applicability of the steady-state productivity approach for horizontal Wells in low-permeability heterogeneous gas reservoirs. Processes13 (1), 173. 10.3390/pr13010173

  • 65

    ZhaoY.LiuX.ZhangL.TangH.XiongY.GuoJ.et al (2021). Laws of gas and water flow and mechanism of reservoir drying in tight sandstone gas reservoirs. Nat. Gas. Ind. B8 (2), 195–204. 10.1016/j.ngib.2020.09.018

  • 66

    ZhaoY.-L.FengH.-Y.BuC.-Z.ZhouL.-S.WuJ.-F.ZhangL.-H.et al (2024). Interwell interference model of horizontal Wells in shale gas reservoirs based on multi-connected boundary element method. Petroleum Sci.21 (6), 4278–4297. 10.1016/j.petsci.2024.08.009

  • 67

    ZhouY.HatzignatiouD. G.HellandJ. O. (2017). On the estimation of CO2 capillary entry pressure: implications on geological CO2 storage. Int. Journal Greenhouse Gas Control63, 26–36. 10.1016/j.ijggc.2017.04.013

Summary

Keywords

corner flow, gas injection and production rates, gas-water movement law, low-permeability heterogeneous gas storage, pore-scale

Citation

Wang D, Wang W, An W, Su J, Zhang J and Liu D (2026) The influence of gas injection and production rates on pore-scale gas-water movement in low-permeability heterogeneous gas storage based on microfluidic technique. Front. Earth Sci. 14:1669602. doi: 10.3389/feart.2026.1669602

Received

20 July 2025

Revised

13 December 2025

Accepted

28 January 2026

Published

04 March 2026

Volume

14 - 2026

Edited by

Vahid Niasar, The University of Manchester, United Kingdom

Reviewed by

Zhengzheng Cao, Henan Polytechnic University, China

Jicheng Zhang, Shandong University of Science and Technology, Qingdao, China

Updates

Copyright

*Correspondence: Dengke Liu,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics