Abstract

Flow characteristics and phase distribution of concentrated brine storage in mining area have been core factors linking to the leakage risk assessment and ecological evaluation. Notably, saturation plays a crucial role in impacting the flow characteristics and distribution of brine, while the existence of oil left by mining machines and original reservoir and gas produced from coal bed gas and air has complicated the issue. In this work, we conducted the microfluidic visualization experiments to reveal the saturation distribution during brine storage in mining area. We applied machine learning model to extract saturation data from experimental images with over 95% accuracy. Eventually, we found that the existence of gas significantly impacts on the saturation distribution in micropores accounting for more than 80% contribution. We clarified that the gas production rate of median 200 μL/min impacts the least on saturation variation. Results in this research are of significance for deeper comprehension on three-phase saturation characteristics of concentrated brine storage in mining area.

1. Introduction

Saturation has been proved to be a crucial parameter impacting flow characteristics of multiphase flow, especially the three-phase porous flow [13]. Notably, the three-phase porous flow widely exists in oil field development, soil sewage invasion, brine storage in mining area, and so on [4, 5]. Due to the fact like this, the three-phase saturation distribution (TSD) is of great significance for a deeper understanding on the three-phase porous flow and better strategies for energy development and ecological conservation [1, 69]. Gas production and flooding, especially during coal mining [1012], have a nonnegligible impact on brine storage in mining area [1316]. In specific, the residual oil in micropores is conducive to plug the predominant pathways for a better brine storage [17, 18], while the gas flooding probably displaces the brine and oil in micropores causing brine transfer, even leakage [4, 1922].

Accordingly, effects of gas flooding and oil residence in micropores during brine storage in mining area require clarification. Since the existence of oil left by mining machines and original reservoir and gas produced from coal bed gas and air has complicated the issue, an effective method is critical. One of the most visual and effective ways to study the mechanisms of TSD by gas flooding is to conduct the visualization flooding experiments. Notably, no matter what type of rocks the reservoir consists of, the inside gas flooding process is invisible [2325]. Exploring and describing the inside pores and fluids via imaging techniques such as CT scanning and nuclear magnetic resonance are one of the methods to achieve the visible development of reservoir by gas flooding. Due to the detail distortion caused by the conversion of photo-electric and magneto-electric signals, overlapping ambiguity, and interference of environmental factors, however, the CT scanning and nuclear magnetic resonance failed to either effectively indicate the microdetails during TSD by gas flooding or clarify the mechanisms of TSD by gas flooding [13, 15, 2628].

The microfluidic visualization experiment is a visible experimental method for studying multiphase porous flow in microscale pores [29]. According to the parameters from actual oil reservoir, the pores in microscale are integrated onto a transparent microfluidic chip in an area of certain square centimeters. Moreover, the visualization experiments of multiphase porous flow are conducted in the microfluidic chip, which could be classified to the application category of reservoir on a chip (ROC). The core advantages of the microfluidic visualization experiment are directly visible, highly efficient, less space occupied, and chip design flexible [3032]. The application of microfluidic visualization experiment is of great significance on clarifying the mechanisms of TSD by gas flooding and further optimizing the TSD by gas flooding technology [3338].

In this research, we conducted the microfluidic visualization experiments according to the reservoir properties of Ling Xin mining area to reveal the saturation distribution during brine storage. We applied machine learning model to extract saturation data from experimental images. We indicated the impact of gas existence on the saturation distribution in micropores. We clarified how the gas production rate impacts on the three-phase saturation variation. Results in this research are of significance for deeper comprehension on three-phase saturation characteristics of concentrated brine storage in mining area.

2. Methodology

Figure 1 displays the route of microfluidic visualization-based mechanism and optimization research on TSD by gas flooding. The whole research was divided into microstructure design, microfluidic chip fabrication, and injection experiment, three parts. First, the target pore structure was formed by reservoir parameters restrained random micropore generation. Based on the target pore structure, the mask pattern was determined. Microfluidic chips with mask pattern were fabricated via photolithography. Afterwards, the gas flooding experiments were established based on the microfluidic chips. Flooding pressure and velocity and experimental images were measured and recorded. TSD by gas flooding was extracted via image recognition algorithm according to experimental images. Mechanisms of TSD by gas flooding were revealed while the flooding optimization was proposed from data and experimental images eventually.

2.1. Microstructure Design

The pore structures in porous media such as the reservoir are complex and random. In general, the experimental methods only obtain statistical average features such as core porosity and permeability. While actually, to more accurately reflect the multiphase porous flow process of gas flooding in pores, the more detailed description on pore topology, connectivity, and tortuosity is required. By considering the random features of pores in oil reservoir, the QSGS (quartet structure generation set) algorithm was adopted in this research for generating the pore structures with high similarity to oil reservoir [39]. Moreover, the QSGS algorithm was improved and enhanced to satisfy the demands for pore topology, connectivity, and tortuosity in this research.

Random pore generation methods and the improved QSGS algorithm adopted in this research are presented in Figure 2. Generally, there are three main methods of random pore design and generation consisting of digital core, computer assisted drawing (CAD), and optimization algorithm, as shown in Figure 2(a). The optimization algorithm, specifically the improved QSGS algorithm, was adopted in this research. Figure 2(b) presents the pore generation process of the improved QSGS. First, certain number of generation cores were distributed in the region. Afterwards, the generation cores grew in eight directions at given probabilities and occupied the blank region constantly, forming blocks in the region. Compared with the original QSGS, the improved QSGS would judge whether the next step growth resulted in the connection of two blocks or not. If so, the improved QSGS would stop the growth of blocks in the connection direction to ensure the fine connectivity of generated pores. When the porosity reached the setting value, the improved QSGS ended and output the final generated pore structure. Route of pore generation by improved QSGS could be seen in Figure 2(c) [40, 41].

2.2. Microfluidic Device Fabrication

The microfluidic chip was fabricated in laboratory at room temperature 24.3°C and relative humidity 70%. The first step was to fabricate the male mold. First, the silicon wafer was cleaned via ultrasonication at 120 W in 10 minutes. Then, the SU8-3050 negative photoresist was coat spun upon the silicon wafer. Afterwards, the photoresist coated silicon wafer was soft baked at 95°C for 10 minutes. Next, the wafer was put in vacuum chamber, covered by the designed mask and exposed in UV light at 50 mW/cm2 for 10 seconds. The exposed wafer was postexpose baked at initially 65°C for 2 minutes and at 95°C for 8 minutes later. After this, the wafer was put in the PGMEA (propylene glycol methyl ether acetate) solution for 8-minute development. Eventually, the developed wafer was dried at 150°C for 10 minutes to evaporate the remaining SU-8 photoresist and fasten the male mold with the wafer [42].

The second step was to fabricate the microfluidic device. First was to mix the component A (dimethylvinyl siloxane) and component B (curing agent, dimethyl hydrosiloxane) of PDMS at the ratio of 9 : 1. Then, the PDMS was evacuated until small bubbles were completely discharged. Next, the PDMS was poured into a Petri dish placing the male mold. The Petri dish was put in a baker at 75°C to cure the PDMS. Afterwards, the PDMS device was cut out along the edge of male mold and holes were punched at inlets and outlets of the microchannels. The PDMS device was plasma treated together with the glass slide at 100 W for 1 minute, eventually was to bond the PDMS device with the glass slide with slight press and the microfluidic chips were completely fabricated. The chip was heated at 65°C for 30 minutes for stronger bonding [43, 44].

2.3. Injection Experiment Scheme

We conducted the microfluidic injection experiments at a laboratory at room temperature of 24.0°C with 80% relative humidity. The atmospheric pressure of the laboratory site was 99.2 kPa. Detailed properties of fluids used in experiments are shown in Table 1. All the fluid properties were obtained from or according to the Ling Xin mining area, China. The three-phase porous flow conducted in this work could be considered as an incompressible two-dimensional flow [4549].

The schematic of microfluidic injection experiment and data record is presented in Figure 3, referring to the real flowing in the Ling Xin mining area, China. Generally, four modules working for the driving, pressure monitoring, experiment visualization, and data record were set. The equipment and materials contained two constant-speed microinjection pumps for liquid and gas injection, respectively. Medical plastic syringes, capsule pressure gauge, four-way valve, microfluidic chips, optical microscope with electron lens, computer system, waste collection tube, and experimental recording camera were used. The connection between equipment and materials is as shown in Figure 3. The microfluidic chip was initially saturated by the oil. The injection rate was set in the microinjection pump corresponding to the injection rate at the inlet of microfluidic chip [12, 5054].

During the experiment, the fluid pump was switched on to inject the water. In real time, the pressure gauge showed the inlet pressure of the microfluidic chip, while the camera recorded the pressure data and the electronic lens on the microscope recorded the experimental images. Pressure data and experimental images were transmitted to the mobile workstation for storage. After the flow in the microfluidic chip stabilized, the fluid pump was stopped and the gas pump was turned on to inject gas into the microfluidic chip, starting the TSD by gas flooding experiment. The operation processes of each module were consistent with that of the water injection. After the flow in the microfluidic chip, we stabilized again, stopped the gas pump, removed the connecting pipeline between the modules, cleaned up the test table, and collected waste liquid [46, 50, 51, 55].

3. Results and Discussion

3.1. Microstructure and Microfluidic Device

In this research, the improved QSGS (quartet structure generation set) algorithm was adopted for the pore structure design and generation. The pore structure generated by original QSGS algorithm, the pore structure generated by improved QSGS algorithm, the mask pattern for microfluidic chip fabrication, and the fabricated microfluidic chip are shown in Figure 4. The size ratio on pore structure generated by improved QSGS algorithm, mask pattern, and microfluidic chip are 1 : 1 : 1 particularly.

According to Figure 4, it could be found that compared to the pore structure generated by original QSGS, the improved QSGS enhanced the connectivity of the pore structure, which was more fitting in the actual pore structures in rocks. Moreover, the improved QSGS maintained the advantages from the original version consisting of settability, high efficiency, clear physical significance, and so on. The dead zones of the pore structure were eliminated and the edges of pores were smoothed before drawing the mask pattern. Finally, the microfluidic chips with fine connectivity, topology features, and visualization were obtained.

3.2. Saturation Extracted from Experimental Images

Differed from conventional core experiments, saturation of fluids during the flooding process was extracted from experimental images based on image recognition algorithm as demonstrated in Figure 5 instead of the measurement and subtraction of fluid volume between the inlet and outlet. Besides, the gas flooding experiments conducted in this research were classified as the nonsteady state flooding method.

There were two reasons why image recognition-based saturation measurement and nonsteady state flooding method were adopted in this research. On the one hand, the microfluidic visualization experiments focused on the microscale of TSD by gas flooding. Hence, few samples were used during the complete experiments, which was hard to accurately and precisely measure the liquid volume at the microfluidic outlet. On contrast, the experimental images of fluid state in pores at every time step were zoomed in by microscope and recorded by a CCD high definition sensor. Liquid saturation and its change were easily and precisely captured and analyzed by the image recognition algorithm. On the other hand, since that the practical gas flooding in developed oil reservoirs is a nonsteady state process, the nonsteady state flooding method was applied for more appropriate research on TSD by gas flooding.

The TSD at diverse time was calculated based on the oil saturation difference captured from the image recognition algorithm. As shown in Figure 5, first a set of experimental images captured throughout the gas flooding was input into the algorithm. The study area, the pore structure actually, and the different liquid phase were extracted by the algorithm, respectively. Then, the saturation of each type of liquid was counted. Via calculating the ratio of the saturation of each liquid at a certain moment to the initial oil saturation, the TSD at diverse time was obtained eventually.

3.3. Characteristics of Three-Phase Saturation Distribution in Mining Area

The complete process of TSD by gas flooding based on microfluidic visualization experiment is shown in Figure 6. Experiments under diverse injection rate conditions have been conducted while the specific experiment shown in Figure 6 was at the gas injection rate of 200 μL/min. The complete process included three sections. First, the pores in the microfluidic chip were saturated by oil. Next, the water was injected for oil flooding. Only the image of water flooding finished is presented in Figure 6. Last but not least, the gas was injected for oil flooding. Images of gas flooding for 6 s, 8 s, 10 s, 120 s, 1200 s, and 1500 s at which the gas flooding stopped is presented in Figure 6. The complete experiment process was designed and conducted according to the practical gas flooding in oil reservoirs. The water flooding and gas flooding were stopped when the liquid saturation and distribution in pores reached to be stable.

According to the microfluidic visualization experiment shown in Figure 6, it is worth noting that since the materials of the microfluidic chip were hydrophobic and oil wetting, oil in large pores was well displaced by water flooding while that in fine pores was residual. It could be explained that the capillary force on water was resistance which would be stronger in fine pores. Hence, water flooded through the large pores in which the capillary force resisted less.

The same thTSDy could be applied on the gas flooding. According to Figure 6, the gas invaded into the large pores and displaced the water initially before invading into the fine pores. Moreover, it could be found that the gas rapidly formed an escape path initially as shown from gas flooding for 6 s to gas flooding for 10 s in Figure 6. After the escape path was formed, the invasion of gas into pores slowed down, as shown from gas flooding for 120 s to gas flooding for 1200 s in Figure 6. Nevertheless, the TSD was achieved by gas flooding by comparing the image of water flooding finished with that of gas flooding for 1500 s in Figure 6, with a saturation contribution of 82.74%.

To further clarify the mechanism of the TSD achieved by gas flooding on microscale, Figure 7 presents the oil displacement by gas flooding via various ways in details. In Figure 7, the images of water flooding finished and that of gas flooding stopped in Figure 6 were picked for comparison to clearly present the main mechanism of the TSD by gas flooding.

According to Figure 7, two dominant mechanisms of TSD by gas flooding were presented. Mechanism A is that the gas would invade into pores that water failed to permit to enhance oil recovery. It is due to that the density and viscosity of gas are approximately 1% comparing to that of water and oil. Though the process was gradual, the gas achieved to permit into the fine pores at the similar flooding pressure of water flooding. Mechanism B is that gas would strip off the oil film attached to the pore wall to enhance oil recovery, since the injection of gas made up the pressure loss caused by stop of water flooding. It is worth noting that the mechanism A contributed more to the TSD than the mechanism B, as demonstrated in Figure 7.

3.4. Effects of Gas Production Rate on Saturation Distribution in Mining Area

In addition to clarifying the mechanisms of TSD by gas flooding, the optimization of gas flooding was studied either in this research. The core way to optimize the gas flooding is to figure out the effects of injection rate which is a crucial working condition in practical reservoir development. In this research, 100 μL/min, 200 μL/min, and 300 μL/min, three injection rates of gas flooding, were studied.

Figure 8 shows the continuous gas injection time, flooding pressure curve at diverse gas injection rate. Generally, it could be found that all curves reached a peak in certain time because of the compression of the gas at initial flooding period. Specifically, the scheme of 100 μL/min reached the peak most early while the scheme of 200 μL/min and 300 μL/min reached the peak at similar time. It is apparent that the higher the injection rate was, the higher the peak of flooding pressure could reach and the later the gas flooding achieved stability. However, it is worth noting that the relationship between injection rate and flooding pressure not strictly followed the Darcy’s law, revealing a characterization of nonlinear flow. After reaching the peak, flooding pressure started to decline due to the gas decompression caused by most gas escape.

The continuous gas injection time-oil saturation curve at diverse gas injection rate is presented in Figure 9. According to Figure 9, it is surprising that there was also a nonlinear relationship between injection rate and oil saturation. More findings come out with comprehension on Figure 9 combining with Figure 8.

We observed the highest flooding pressure and shortest stabilization time, but the most remaining oil saturation at the injection flow rate of 300 μL/min. It is dominantly because the high injection rate results in serious gas escape failing to enhance oil recovery. Point A and Point B in Figure 9 respectively represent the oil saturation equivalent point of 300 μL/min with 200 μL/min and 300 μL/min with 100 μL/min, showing that the injection rate of 300 μL/min achieved the highest oil saturation initially while left behind by 200 μL/min and 100 μL/min gradually. Besides, it indicates that the scheme of 100 μL/min and 200 μL/min achieved the similar highest oil saturation where the oil saturation by gas flooding at 200 μL/min was 0.23% higher than that at 100 μL/min from Figure 9. According to Figure 8 and the mechanism A as mentioned previously, the scheme of 100 μL/min failed to offer adequate flooding pressure for gas to permit into the fine pores displacing the residual oil. Furthermore, the scheme of 200 μL/min became stable much earlier than 100 μL/min, which means it is more economical in practical reservoir development. Hence, the 200 μL/min is considered as the optimal injection rate in this research rather than 100 μL/min or 300 μL/min, inspiring the choice of injection rate in practical flooding engineering.

Furthermore, to figure out the effects of gas on oil during the flooding process, the relationship between TSD and gas saturation and variation of TSD and gas saturation at diverse gas injection time of various gas injection rate were considered and studied in this research, as specifically revealed by Figures 10 and 11.

According to Figure 10, especially Figure 11, it could be found that at the initial gas flooding period, especially the first 300 s, the variety of brine saturation (Δ brine saturation) fluctuated intensively. Meanwhile, the variety of oil saturation (Δ oil saturation) fluctuated along with the Δ brine saturation showing diverse features at different injection rate, however. Particularly, the Δ brine saturation at 300 μL/min reached the peak ahead of the Δ oil saturation, which means that the gas displaced the water first and then performed the oil flooding. The same could be seen on the injection rate of 100 μL/min in Figure 11 though the peak reaching moment was 100 s behind that of 300 μL/min. It is worth noting that at 200 μL/min, the raising and decreasing periods of Δ brine saturation and Δ oil saturation were coincident demonstrating that the gas displaced the oil and water simultaneously, which is valuable in practical reservoir development. It also indicates that the injection rate of 200 μL/min in this research was the optimal injection rate from another point of view.

In the middle and later periods of gas flooding, especially after 500 s injection, both the Δ brine saturation and Δ oil saturation declined at a similar rate no matter what injection rate it was according to Figure 11. Hence, in the middle and later periods, water had been mostly displaced by the gas and the gas started to flood the oil directly based on both mechanism A and mechanism B. Nonetheless, the most crucial period of oil saturation by gas flooding was the initial injection, comprehensively based on Figures 10 and 11.

4. Conclusion

The flow characteristics and phase distribution of concentrated brine storage in mining area have been core factors linking to the leakage risk assessment and ecological evaluation. Notably, saturation plays a crucial role in impacting the flow characteristics and distribution of brine, while the existence of oil left by mining machines and original reservoir and gas produced from coal bed gas and air has complicated the issue.

In this work, we conducted the microfluidic visualization experiments to reveal the saturation distribution during brine storage in mining area. We applied machine learning model to extract saturation data from experimental images with over 95% accuracy. Eventually, we found that the existence of gas significantly impacts on the saturation distribution in micropores accounting for more than 80% contribution, while the oil contributed nearly 15%. We clarified that the gas production rate of median 200 μL/min from the scheme of 100 μL/min, 200 μL/min, and 300 μL/min impacts the least on saturation variation. In general, advantages including visualization and high efficiency of microfluidic technology have boosted the enhancement and optimization of three-phase saturation distribution researches in mining area. Results in this research are of significance for deeper comprehension on three-phase saturation characteristics of concentrated brine storage in mining area.

Data Availability

Data are available on request. Please contact the corresponding author for the underlying data supporting the results of the research.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

The authors are grateful for financial support from the National Energy Group Ningxia Coal Industry Company Lingxin Coal Industry Demonstration Project (Environmental Risk Research of Underground Brine Storage Reservoir in Mining Area), Open Fund of Chinese State Key Laboratory of Water Resource Protection and Utilization in Coal Mining with grant number SHJT-16-30.3.