In fluence of formation heterogeneity on foam flooding performance using 2D and 3D models: an experimental study

Ling?Zhi Hu · Lin Sun · Jin?Zhou Zhao · Peng Wei · Wan?Fen Pu

Abstract The formation heterogeneity is considered as one of the major factors limiting the application of foam flooding. In this paper, in fluences of formation properties, such as permeability, permeability distribution, interlayer, sedimentary rhythm and 3D heterogeneity, on the mobility control capability and oil displacement efficiency of foam flooding, were systematically investigated using 2D homogeneous and 2D/3D heterogeneous models under 120 °C and salinity of 20 × 104 mg/L. The flow resistance of foam was promoted as the permeability increased, which thus resulted in a considerable oil recovery behavior.In the scenario of the vertical heterogeneous formations, it was observed that the permeability of the high-permeable layer was crucial to foam mobility control, and the positive rhythm appeared favorable to improve the foam flooding performance.The additional oil recovery increased to about 40%. The interlayer was favorable for the increases in mobility reduction factor and oil recovery of foam flooding when the low permeability ratio was involved. For the 3D heterogeneous formations, foam could efficiently adjust the areal and vertical heterogeneity through mobility control and gravity segregation,and thus enhancing the oil recovery to 11%–14%. The results derived from this work may provide some insight for the field test designs of foam flooding.

Keywords Foam flooding · Permeability ratio · 2D/3D models · Heterogeneity · Enhanced oil recovery

1 Introduction

The technical feasibility of foam flooding for enhanced oil recovery (EOR) has been approved previously (Wei et al.2018; Zhang et al. 2010). It was recognized that the high flow resistance created by foam in high-permeability zones could promote fluid diversion to the previously unswept zones, and the gravity difference appeared in the presence of foam could help to sweep the upper part of reservoirs(Ma et al. 2013; Pang et al. 2018). Meanwhile, owing to the viscoelastic liquid films, the superposed gas bubbles could squeeze and pull oil droplets out of the dead end pores(Arnold and Wit 1990; Ashoori et al. 2011). Furthermore,the foam agents, which are surfactant-based materials, help in reducing the interfacial tension between oil and water and altering rocky wettability (Kamal 2016; Kumar et al. 2008;Pu et al. 2016). Therefore, a foam flooding process is not only capable of improving the sweep efficiency but also able to improve the microscopic displacement efficiency, which therefore results in a signi ficantly high oil recovery factor(Farajzadeh et al. 2009a, 2012; Shokrollahi et al. 2014).

It was well-documented that the sweep efficiency of foam flooding could be noticeably improved through screening foam agents and optimizing injected parameters (Hou et al.2013; Li et al. 2011; Liu et al. 2014; Wang et al. 2015). In addition to the aforementioned, the physical properties of the formation also play a crucial role in foam flooding behavior.As some researchers claimed, the mobility reduction factor (MRF) that represents the mobility control capacity of foam in the porous media is signi ficantly improved with the increases in core permeability and core length (Farajzadeh et al. 2009b; Wang et al. 2012; Zhao et al. 2012). In terms of heterogeneous formations, foam creates a high sweep efficiency by diverting fluids to relatively low-permeable zones (Liu et al. 2008), especially for the positive rhythm cases (Liu et al. 2010; Ma et al. 2013). Siddiqui et al. (1997)stated that fluid diversion can be enhanced by increasing the permeability ratio within a speci fic range. While placing an impermeable layer in the multilayer core, the cross flow and fluid diversion are both hindered in foam flooding. Although pioneering works have made great efforts to address the factors that in fluence the foam sweep efficiency, there are very limited studies of the impact of formation heterogeneity,which is believed to be extremely important for design of foam flooding pilot in oil fields.

Given this issue, the primary objective of this work is to systematically investigate the EOR performance of foam flooding in different types of formations, and figure out the in fluence of formation structure on the foam flooding effi-ciency. To accomplish this research objective, long sand packs, which had different permeabilities, were initially set up to examine the flow behavior of foam including mobility control and flow resistance in the subsequent water flooding. Thereafter, 2D homogeneous/heterogeneous cores were used to study the effects of permeability, permeability distribution, interlayer and sedimentary rhythm on the foam mobility control and oil recovery. Eventually, foam flooding tests were carried out in the 3D heterogeneous models to comprehensively investigate the displacement behavior in complex heterogeneous conditions.

2 Experimental

2.1 Materials

Crude oil and formation water were kindly provided by an oil field located in western China. The viscosity of the crude oil was about 8.5 mPa s at 120 °C. The composition of formation water is given in Table 1. N2had purity of 99.9% (Xinju Co., Ltd). The foaming agent HTS was synthesized in our laboratory, which was a mixture of glycoside surfactant and zwitterionic surfactant. Based on the Waring blender method, the foam volume and the half-life time of a 0.2% HTS solution were 580 mL and 185 s at 120 °C. The interfacial tension between the 0.2% HTS solution and the crude oil was 1.5 mN/m. Quartz sand was supplied by Pixian Co., Ltd (Chengdu, China). The cores used in this work were arti ficial cores mainly composed of quartz sand, and supplied by the Northeast Petroleum University (Daqing, China).

2.2 Experimental apparatus and procedures

The sand pack and 2D model tests were conducted at 20 MPa and 120 °C. The 3D model tests were carried out at 5 MPa and 120 °C. In all experiments, foam flooding was carried out by co-injecting the nitrogen and foam solution (0.2% HTS, foam quality was 0.67).

2.2.1 Sand pack tests

Figure 1 shows a schematic diagram of the apparatus for sand pack tests. Four pressure monitoring ports were evenly distributed along the sand pack, which was filled with quartz sand (70–140 mesh). The length and diameter of the sand packs were 100 cm and 2.5 cm, respectively.The physical parameters of the sand packs are listed in Table 2.

The formation water was first injected at 1.0 mL/min followed by foam injection. At each stage, the differential pressure ?Pwiof water and ?Pfiof foam at each section were recorded (Fig. 2) when four pressures changed slightly. Thereafter, the corresponding mobility reduction factor (MRF) was calculated from Eq. (1). Formation water injection was resumed at the same rate (1.0 mL/min), and the dynamic differential pressure ?Pswiwas measured in order to obtain the residual resistance factor(RRF) from Eq. (2).

2.2.2 2D model tests

The apparatus for 2D model tests is analogous to Fig. 1with the exception of the porous media. To simulate different permeability properties, three types of cuboid cores were prepared as shown in Fig. 3. For the multilayer core,the high permeability layer was placed at the bottom of the core holder to simulate the positive rhythm if not speci fied.The physical properties of the cuboid cores are given in Tables 3 and 4.

Table 1 Composition of formation water

Fig. 1 Schematic diagram of the apparatus for sand pack tests

Table 2 Physical parameters of the sand packs used

Fig. 2 Diagram of pressure monitoring for sand pack tests. The red arrow points to the flow direction

Fig. 3 Images of different cuboid cores. a Homogeneous core. b Two-layered heterogeneous core. c Two-layered heterogeneous core with an interlayer involved

Table 3 Physical parameters of the homogeneous cores for oil displacement experiments

Table 4 Physical parameters of the heterogeneous cores

In mobility control experiments, formation water was first injected until the steady state was reached; afterward, the foam was injected, during which the dynamic differential pressure (?Pセ) was recorded in order to calculate MRF.

In the case of oil displacement experiments, water flooding was first conducted until the water cut reached 98%.After that, 0.5 pore volume (PV) foam was injected followed by water flooding. The injection rate of 1.2 mL/min was used in the 2D model tests.

2.2.3 3D model tests

A schematic diagram of the apparatus for 3D model tests is shown in Fig. 4. The 3D model had dimensions of 14.4 cm × 14.4 cm × 10.2 cm (length, width, and height) and was equipped with six valves. Valves 1 and 2 were regarded as an injection well and a production well, respectively.Valves a, b, c and d were the spare switches used to operate the water and oil saturation process. Three types of sands(i.e, 70–140 mesh, fine; 40–70 mesh, medium; 20–40 mesh,rough) were used to prepare the 3D models with the three layers for the same thickness, which are particularly shown in Fig. 5 and Table 5.

The procedures of oil displacement experiments conducted on the 3D model were similar to that of the 2D model tests. The only difference was that the injection rate was 4.0 mL/min. Moreover, images of the three sand layers were captured with a high-resolution camera, and then software ImageJ was used to analyze the residual oil distribution.

3 Results and discussion

3.1 Mobility control behavior in sand packs of different permeabilities

3.1.1 Mobility control

It was known that dynamic equilibrium process between generation and destruction usually occurs during foam propagating in porous media (Blázquez et al. 2016).Hence, if the foam generation exceeds rupture, MRF would increase and thus strengthen the mobility control capability and vice versa (Almajid and Kovscek 2016). The rates of foam generation and destruction vary with permeability, which would alter the mobility control characteristics during foam flooding. As shown in Fig. 6, when the water permeability of the sand pack was less than or equal to 1105 mD (Kw≤ 1105 mD), the MRF progressively reduced from the front to the tail. Moreover, as the sand pack permeability increased, the position where the largest value of MRF occurred shifted from the front to the tail, which revealed that the effective distance of foam mobility control was extended. Regarding the same part of different permeability sand packs, the high permeability resulted in the large MRF. Nevertheless, when the permeability was increased to 1592 mD, the MRF in the front turned out to be smaller than that of the 753-mD one; meanwhile, the generated MRF gradually increased from the front to the tail. With the exception of the MRF in the front, the MRF in the other parts was larger than all the evaluated sand packs. This interesting result could be attributed to the shear forces and entrance effects in the porous media (Ashoori et al. 2011).The weaker shear and more obvious entrance effect were found at a higher permeability porous medium (Parlar et al.1995). As a result, the weak foam texture was produced in the front of the sand pack (1592 mD), and the low shear rate led to high apparent viscosity of foam, which bene fited foam mobility control in the other parts.

Fig. 4 Schematic diagram of the apparatus for 3D model tests

3.1.2 Flow resistance in the subsequent water flooding

The RRF produced during water injection can illustrate the flow resistance in the subsequent water flooding in the presence of foam. As shown in Fig. 7, with an increase in permeability, the shape of RFF curves was signi ficantly changed.This result is believed to be associated with foam destruction and fluid diversion. As a result of the augmented flow resistance in large pores, the subsequent water could sweep both the large and small pores. The water in the large pores accelerated foam destruction and thus caused the RRF to reduce.On the contrary, the water diverted to small pores would increase the flow differential pressure and also RRF, especially in high permeability sand packs with strong mobility control potential. Consequently, when the permeability was high enough, the RRF would be promoted substantially at the early stage of water injection. Furthermore, the equilibrium RRF, the ratio of the maximum RFF to the initial RFF and also the corresponding water-injected volume were all directly proportional to permeability. The data demonstrated that the effective time of foam mobility control prolonged with permeability.

3.2 Foam flooding performance in 2D homogeneous models

Five homogeneous cores with average gas permeabilityKgof 200–2000 mD were used to study the effect of permeability on the performance of foam in EOR. The experimental results are given in Table 6, and the production dynamics are shown in Fig. 8.

Fig. 5 Patterns of sand packing in the 3D models (i.e., Model A, Model B, Model C). a Top layer, b middle layer and c bottom layer of the model. Open circle: injection well and closed circle: production well

Table 5 Permeability properties of the 3D models

Fig. 6 MRF in different parts of sand packs at the steady state

As given in Table 6, with an increase in permeability,the water flooding recovery and the additional oil recovery of foam flooding and subsequent water flooding were both enhanced. It was also observed that when the permeabilityKg≥ 700 mD, the water cut dropped precipitously to 47%–63% due to foam injection, and remained at a low level for quite a long duration of time (Fig. 8). This result can be interpreted as follows: (1) For a higher permeability core, the oil recovery was in a larger enhancement as foam flooding was conducted, and the smaller-sized pores in the oil-bearing zones were supposed to be substantially swept resulting from the mobility control and water flow resistance in the presence of foam. (2) The stable foam in the high permeability cores could improve the oil displacement efficiency (Sun et al. 2014). Moreover, the low residual oil saturation in the high-permeability cores after water flooding was in favor of mitigating the de-foaming effect of oil to foam, which has been previously veri fied(Blázquez et al. 2016; Sun et al.2014).

Fig. 7 Residual resistance factor (RRF) in different permeability sand packs as a function of injection volume. a 295 mD. b 753 mD. c 1105 mD. d 1592 mD

Table 6 Summary of oil displacement tests

Fig. 8 a Water cut and b oil recovery as a function of injection volume for 2D models of different permeabilities

Fig. 9 MRF as a function of injection volume for 2D heterogeneous models of different permeability distributions

3.3 Foam flooding performance in 2D heterogeneous models

In some research, several independent cores were parallelized to simulate heterogeneous formations in the laboratory(Shamekhi 2013; Xiao et al. 2014; Vernáez et al. 2016).However, the system without capillary communication can hardly represent the state of the real flow. Therefore, in our study, two types of heterogeneous cores, i.e., two-layered heterogeneous cores and two-layered heterogeneous cores with interlayer involved were assembled.

3.3.1 Effect of permeability distribution

The heterogeneous cores were set up to examine the effect of permeability distribution on foam flooding behavior.

Table 7 Oil recovery for different permeability distributions

Figure 9 plots the MRF as a function of the volume of foam injected. The data indicated that the maximum permeabilityKmaxwas the dominant factor for foam mobility control in the heterogeneous cores, while the permeability ratio only determined the increasing rate of MRF. At the very early stage of foam injection, foam was prone to be generated and propagated in the high permeability layer(Thompson and Gdanski 1993). Therefore, the increasing rate of MRF was almost identical in Cores 22 and 23, which had the sameKmax. Whereas once the foam was diverted to the low permeability layer, a much higher increasing rate of MRF was found in Core 22 which had a higherKmin. After a sufficient volume of foam had been injected, the high permeability layer would be saturated with strong foam, which subsequently formed a considerable MRF that represented the entire core. The phenomena could be understood by the same equilibrium MRFs of Cores 22 and 23, and the discrepant equilibrium MRFs of Cores 21 and 22.

Fig. 10 a Water cut and b oil recovery as a function of injection volume for 2D heterogeneous models of different permeability distributions

Fig. 11 MRF as a function of injection volume in 2D heterogeneous models with or without interlayer. a Heterogeneous cores of a low permeability ratio (300 mD/700 mD). b Heterogeneous cores of a high permeability ratio (300 mD/1700 mD)

The production dynamics in the heterogeneous cores with different permeability distributions are shown in Table 7 and Fig. 10. It can be seen that both theKmaxand permeability ratio signi ficantly in fluenced the oil recovery of foam flooding. At the given permeability ratio (Cores 24 and 26),the additional oil recovery increased with an increase inKmax. Similarly, Pu et al. (2017) found that the oil recovery enhanced by foam increased with the permeability as foam flooding was conducted in the heterogeneous cores. Moreover, the water cut quickly declined to the lowest level and remained low for a relatively long production time. With an increase inKmax, the foam could not only displace the oil from the high permeability layer, but also form a stable front for the subsequent water flooding. At the givenKmax(Cores 26 and 27), similar tendencies of oil recovery and water cut were also observed with an increase in the permeability ratio. The low water flooding recovery of Core 27 provided a high potentiality of oil recovery for the foam flooding.In addition, a comparison of Figs. 8 and 10 veri fied that foam flooding would be more efficient for heterogeneous reservoirs than that for homogenous reservoirs, which was reported by Liu et al. (2008).

3.3.2 Effect of interlayer

The effect of interlayer on foam flooding behavior was investigated in heterogeneous cores with the same permeability ratio. The details of foam flooding are listed in Table 8.

As shown in Fig. 11, the effect of interlayer varied with permeability ratio mostly embodied two aspects: the larger MRF was generated in the case of interlayered cores with the low permeability ratio (2.3–2.4); On the contrary, the MRF became smaller in the interlayered cores with the high permeability ratio (5.7). Nguyen et al. (2005) also con firmed that the interlayer played a negative role in foam performance when there was a strong heterogeneity. It was known that the foam-triggered cross flow usually occurs in heterogeneous cores with capillary communicating (Li et al. 2010;Bertin et al. 1998). The low density allowed the foam in the bottom layer of high permeability to permeate into the top layer of low permeability, while the capillary pressure of the low permeability layer gave rise to the reverse flow to the high permeability layer. Therefore, when the permeability ratio was low, the foam tended to propagate to the low permeability layer decreasing the MRF. However, the existence of interlayer could hinder the foam cross flow creating strong mobility control. Nevertheless, when the permeability ratio was high (i.e., a larger permeability of the high permeability layer), the prevention of downward foam flow by the impermeable interlayer would diminish the foam mobility control.

Based on foam flooding experiments, Li et al. (2010)demonstrated that cross flow was critical in promoting the sweep efficiency of foam flooding. However, we found that when the permeability ratio was high (5.7), the cross flow favored the EOR performance of foam flooding. If the permeability ratio was low (2.3), the cross flow imposed a negative effect on foam flooding (Figs. 12, 13). The in fluence of interlayer on oil recovery was similar to that on mobility control, whereas the positive response of oil recovery to the interlayer seemed less sensitive than that of mobility control, which was related to the distribution of residual oil and the oil de-foaming action.

Table 8 Oil recovery in heterogeneous cores with or without interlayer

Fig. 12 a Water cut and b oil recovery as a function of injection volume for cores of a low permeability ratio (300 mD/700 mD) with or without interlayer

3.3.3 Effect of sedimentary rhythm

In this section, the heterogeneous cores without interlayer were placed in opposite directions to simulate different sedimentary rhythms. The details of foam flooding are given in Table 9.

As shown in Table 9, the reverse rhythm was favorable for water flooding, whereas the positive rhythm bene fited the foam flooding performance. This effect was even more evident at high permeability ratio. During the stage of water flooding, due to the density difference of oil and water, and the flow resistance difference of high- and lowpermeability layers, the waterline is relatively uniform in the reverse rhythm. A comparison of the two rhythms indicated that the foam was more stable in the high permeability layer because of the low residual oil saturation and was more abundant in the low permeability layer owing to the gravity segregation, which eventually promoted the effi-ciency of oil recovery enhancement and water cut reduction under the condition of positive rhythm (Figs. 14, 15).

Fig. 13 a Water cut and b oil recovery as a function of injection volume for cores of a high permeability ratio (300 mD/1700 mD) with or without interlayer

Table 9 Oil recovery in heterogeneous cores with different sedimentary rhythms

Fig. 14 a Water cut and b oil recovery as a function of injection volume for cores of a low permeability ratio (300 mD/700 mD) with different sedimentary rhythms

3.4 Foam flooding performance in 3D heterogeneous models

The 3D heterogeneous models were established to investigate the effect of areal and vertical heterogeneity on foam flooding performance. The oil recovery data are listed in Table 10. The production details are presented in Figs. 16 and 17.

Fig. 15 a Water cut and b oil recovery as a function of injection volume for cores of a high permeability ratio (300 mD/1700 mD) with different sedimentary rhythms

Table 10 Oil recovery in the 3D model

Fig. 16 a Water cut and b oil recovery as a function of injection volume in 3D heterogeneous models

The experimental results turned out that the maximum waterflooding recovery and additional oil recovery occurred in the model A, which was only vertically heterogeneous (Fig. 16). For the models that were both areally and vertically heterogeneous, when the wells were assigned in the low permeability zone with vertical homogeneity (model C), the oil recovery was found to be higher(Table 10). In order to further analyze foam flooding in the 3D model, we peeled off the sand model layer by layer after foam flooding and carried out image analysis in terms of each layer. As shown in Fig. 17, uniform swept regions (the light area) were observed in all the middle and bottom layers (Fig. 17b, c), although the residual oil saturation was lower for the less heterogeneous layer. This confirmed that the foam would stabilize the displacement front through improving the mobility control. Meanwhile,a stable displacement and a low residual oil saturation were observed around the production well in the top layer. This interesting phenomenon proved that the gravity segregation of foam took place along the displacement direction, which could be considered as a type of deep fluid diversion. The gravity segregation made a great contribution to oil displacement efficiency in the top layer,but also caused the residual oil mainly distributed in the bottom layer. Therefore, it needs a further research on the optimization of the gravity segregation of foam for achieving a higher sweep efficiency.

Fig. 17 Distributions of residual oil in different layers of 3D models. a Top, b middle and c bottom layers. The dark spots represent high oil concentrations, while the lighter spots represent low oil concentrations (open circle: injection well and closed circle: production well)

4 Conclusions

Foam flooding performance in different types of formations was comprehensively investigated to figure out the various EOR mechanisms. Based on the experimental results, some main conclusions can be stated as follows:

1. With an increase in permeability, the mobility control and sequent water seepage resistance of foam were signi ficantly improved. A high flow resistance of foam could be produced in a long distance displacement process, which accordingly resulted in prominent oil recovery efficiency.

2. When foam flooding conducted in the vertical heterogeneous formations, the maximum permeabilityKmaxwas mainly responsible for the equilibrium MRF, and the permeability ratio only altered the increasing rate of MRF. Meanwhile, the highKmaxand large permeability ratio bene fited the oil recovery enhancement and water cut reduction. Furthermore, within appropriate permeability ratio range, foam could produce a higher additional oil recovery in a vertical heterogeneous core instead of a homogeneous core.

3. The effect of interlayer on foam flooding varied with permeability ratio. When the permeability ratio was low, the interlayer was favorable for the increase in MRF and oil recovery, while a high permeability ratio did the opposite.

4. In comparison with the reverse rhythm, the positive rhythm was more suitable for foam flooding, especially under a high permeability ratio.

5. Foam was able to tune the areal and vertical heterogeneities through mobility control and gravity segregation,and thus enhanced the oil recovery. However, the effect was not as signi ficant as that in the vertical heterogeneous formations.

AcknowledgementsThis work was financially supported by the Scienti fic Research Startup Foundation of Xinjiang University (No.620312377) and the National Science and Technology Major Project of China (No. 2016ZX05053-013). Special thanks to Sinopec Northwest Company of China for supplying the geological materials and crude oil sample.

Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source,provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons.org/licen ses/by/4.0/.