Flow Characteristics of Foam in Fractures and Prediction of Preferred Paths

2021 ◽  
Author(s):  
Zhengxiao Xu ◽  
Zhaomin Li ◽  
Binfei Li ◽  
Danqi Chen ◽  
Xianghui Zeng ◽  
...  

Abstract Foam is widely used in fractured reservoirs. The flow characteristics in complex fracture networks are still unclear, and there are few numerical simulations of foam fluid flow in fractures. In this study, a variety of combined visual fracture models were used to observe the flow characteristics of foam in the fracture. Firstly, based on the parallel fracture model, the foam flow characteristics under different fracture depths were explored, and then based on the complex fracture network model, the foam flow path and sweep efficiency are evaluated. Finally, the Dijkstra’s algorithm was used to determine the weighted graph of the fracture network nodes, and the preferred flow paths of the foam were predicted. The results show that when foam flows in parallel fractures with different depths, it preferentially flows in high permeability (100 μm) fractures, and there is gas trapping in low permeability (50 μm) fractures. In the irregular fracture network model, the sweep efficiency of the foam fluid is greatly affected by the foam quality, and the sweep volume is the widest when the foam quality is about 90%. The simulation results based on the Dijkstra’s algorithm can be fitted to the experimental results to a certain extent. By controlling the number of preferred paths and the weight of nodes, the plugging and regulating performance of the foam are characterized. These findings reflect the necessity of considering fractures when foam flows in reservoirs, and provide a certain experimental basis and theoretical guidance for the development of fractured reservoirs.

SPE Journal ◽  
2016 ◽  
Vol 21 (04) ◽  
pp. 1140-1150 ◽  
Author(s):  
M. A. Fernø ◽  
J.. Gauteplass ◽  
M.. Pancharoen ◽  
A.. Haugen ◽  
A.. Graue ◽  
...  

Summary Foam generation for gas mobility reduction in porous media is a well-known method and frequently used in field applications. Application of foam in fractured reservoirs has hitherto not been widely implemented, mainly because foam generation and transport in fractured systems are not clearly understood. In this laboratory work, we experimentally evaluate foam generation in a network of fractures within fractured carbonate slabs. Foam is consistently generated by snap-off in the rough-walled, calcite fracture network during surfactant-alternating-gas (SAG) injection and coinjection of gas and surfactant solution over a range of gas fractional flows. Boundary conditions are systematically changed including gas fractional flow, total flow rate, and liquid rates. Local sweep efficiency is evaluated through visualization of the propagation front and compared for pure gas injection, SAG injection, and coinjection. Foam as a mobility-control agent resulted in significantly improved areal sweep and delayed gas breakthrough. Gas-mobility reduction factors varied from approximately 200 to more than 1,000, consistent with observations of improved areal sweep. A shear-thinning foam flow behavior was observed in the fracture networks over a range of gas fractional flows.


SPE Journal ◽  
2021 ◽  
pp. 1-14
Author(s):  
Bing Wei ◽  
Qingtao Tian ◽  
Shengen Chen ◽  
Xingguang Xu ◽  
Dianlin Wang ◽  
...  

Summary There exist two main issues hampering the wide application and development of carbon dioxide (CO2) foam in conformance improvement and CO2 mobility reduction in fractured systems: (1) instability of foam film under reservoir conditions and (2) uncertainties of foam flow in complex fractures. To address these two issues, we previously developed a series of nanocellulose-strengthened CO2 foam (referred to as NCF-st-CO2 foam), while the primary goal of this work is to thoroughly elucidate generation, propagation, and sweep of NCF-st-CO2 foam in a visual 2D heterogeneous fracture network model. NCF-st-CO2 foam outperformed CO2 foam in reducing gas mobility during either coinjection (COI) or surfactant-alternating-gas (SAG) injection, and the threshold foam quality was approximately 0.67. Foam creation was increased with the total superficial velocity for CO2 foam and almost stayed constant for NCF-st-CO2 foam in fractures during COI. For SAG, large surfactant slug could prevent CO2 from early breakthrough and facilitate foaming in situ. The improved sweep efficiency induced by NCF-st-CO2 foam occurred near the producer for both COI and SAG. Film division and behind mainly led to foam generation in the fracture model. Gravity segregation and override was insignificant during COI but became noticeable during SAG, which caused the sweep efficiency decrease by 3 to 9%. Owing to the enhanced film, NCF-st-CO2 foam enabled mitigation of the gravitational effect, especially around the producer.


SPE Journal ◽  
2017 ◽  
Vol 22 (04) ◽  
pp. 1064-1081 ◽  
Author(s):  
Sanbai Li ◽  
Dongxiao Zhang ◽  
Xiang Li

Summary A fully coupled thermal/hydromechanical (THM) model for hydraulic-fracturing treatments is developed in this study. In this model, the mixed finite-volume/finite-element method is used to solve the coupled system, in which the multipoint flux approximation L-method is used to calculate interelement fluid and heat flux. The Gu et al. (2011) crossing criterion is extended to a 3D scenario to delineate the crossing behaviors as hydraulic fractures meet inclined natural fractures. Moreover, the modified Barton et al. (1985) model proposed by Asadollahi et al. (2010) is used to estimate the fracture aperture and model the shear-dilation effect. After being (partially) verified by means of comparison with results from the literature, the developed model is used to investigate complex-fracture-network propagation in naturally fractured reservoirs. Numerical experiments show that the key factors controlling the complexity of the induced-fracture networks include stress anisotropy, injection rate, natural-fracture distribution (fracture-dip angle, strike angle, spacing, density, and length), fracture-filling properties (the degree of cementation and permeability), fracture-surface properties (cohesion and friction angle), and tensile strength of intact rock. It is found that the smaller the stress anisotropy and/or the lower the injection rate, the more complex the fracture network; a high rock tensile strength could increase the possibility of the occurrence of shear fractures; and under conditions of large permeability of fracture filling combined with small cohesive strength and friction coefficient, shear slip could become the dominant mechanism for generating complex-fracture networks. The model developed and the results presented can be used to understand the propagation of complex-fracture networks and aid in the design and optimization of hydraulic-fracturing treatments.


Fractals ◽  
2019 ◽  
Vol 27 (01) ◽  
pp. 1940008 ◽  
Author(s):  
LIMING ZHANG ◽  
CHENYU CUI ◽  
XIAOPENG MA ◽  
ZHIXUE SUN ◽  
FAN LIU ◽  
...  

The distribution of fractures is highly uncertain in naturally fractured reservoirs (NFRs) and may be predicted by using the assisted-history-matching (AHM) that calibrates the reservoir model according to some high-quality static data combined with dynamic production data. A general AHM approach for NFRs is to construct a discrete fracture network (DFN) model and estimate model parameters given the observations. However, the large number of fractures prediction required in the AHM process could pose a high-dimensional optimization problem. This difficulty is particularly challenging when the fractures form a complex multi-scale fracture network. We present in this paper an integrated AHM approach of NFRs to tackle these challenges. Two essential ingredients of the method are (1) a 2D fractal-DFN model constructed as the geological simulation model to describe the complex fracture network, and (2) a mixture of multi-scale parameters, built according to the fractal-DNF model, as an inversion parameter model to alleviate the high-dimensional optimization burden caused by complex fracture networks. A reservoir with a multi-scale fracture network is set up to test the performance of the proposed method. Numerical results demonstrate that by use of the proposed method, the fractures well recognized by assimilating production data.


2017 ◽  
Vol 156 ◽  
pp. 484-496 ◽  
Author(s):  
Lidong Mi ◽  
Bicheng Yan ◽  
Hanqiao Jiang ◽  
Cheng An ◽  
Yuhe Wang ◽  
...  

1999 ◽  
Vol 2 (01) ◽  
pp. 14-24 ◽  
Author(s):  
T.L. Hughes ◽  
F. Friedmann ◽  
D. Johnson ◽  
G.P. Hild ◽  
A. Wilson ◽  
...  

Summary Large-volume foam-gel treatments can provide a cost-effective method to achieve in-depth conformance improvement in fractured reservoirs. The applicability and cost effectiveness of the approach depends on the availability of a cheap source of gas, the efficiency with which the foam can be placed into the high permeability thief zone(s), and the effectiveness of the gelled foam barrier in diverting reservoir drive fluids to improve oil recovery. This paper reviews progress in the application of large-volume CO2-foam-gel treatments to improve conformance in the Rangely Weber Sand Unit (RWSU), Colorado. During the period November 1996-November 1997 three large-volume foam-gel treatments were successfully placed into the Rangely reservoir. The first 36?400 bbl treatment, implemented November 1996, increased the pattern oil rate from 260 barrels of oil per day (BOPD) in March 1997 to ±330 BOPD in August 1998; a conservative estimate of incremental oil recovery is ±40?000 bbl by the end of August 1998. The second 43?450 bbl treatment, implemented August-September 1997, increased the pattern oil rate from ±430 BOPD in March 1998 to ±470 BOPD in August 1998; post-treatment, the pattern oil rate data is described by a linear regression with slope, +56 BOPD but it is too early to make a firm estimate of incremental oil recovery. The third 44?700 bbl treatment, implemented October-November 1997, increased the pattern oil rate from ±330 BOPD in May 1998 to ±375 BOPD in July-August 1998; a linear regression of the post-treatment data gives a positive slope but again it is too early to estimate incremental oil recovery. Some general features in the pattern production response given by the three foam-gel treatments were observed. First, each of the treatments induces a stabilization in the pattern oil rate which, for treatments I and II, is accompanied by a decrease in the pattern gas rate. Second, the first positive oil rate response given by each of the treatments is observed 6-8 months after treatment execution and is dominated by the response at producer wells lying to the west/southwest and/or east/southeast of the treated injector well. For a given treatment volume, the cost of a foam-gel treatment at Rangely is 40%-50% below the average cost of polymer gel treatments. As the foam is injected at a higher rate, the total pump time required for a 40?000 bbl foam-gel treatment is similar to a 20?000 bbl polymer gel treatment. Early during pumping treatments II and III, we attempted to increase the CO2 content of the foam from 80 to 85 vol?%; this resulted in a wellhead pressure which was too close to the CO2 pressure limit necessitating a decrease in foam injection rate. Thus, in optimizing foam-gel treatment cost, there is a balance between maximizing the content of the inexpensive CO2 phase and minimizing total pump time. For Treatments II and III, the cost of the liquid phase formulation was reduced by decreasing the concentrations of surfactant and buffer. The implementation and evaluation of three large-volume foam-gel treatments at Rangely indicates that the foam-gel approach provides a cost-effective method to achieve in-depth conformance improvement in fractured reservoirs. Introduction A recent survey1 indicated that the proportion of U.S. EOR production recovered by gas injection has increased from 18% to 41% during the period 1986-1996. A major contribution to this trend has been the strong increase in the number of miscible carbon dioxide (CO2) projects which now account for > 70% of the total number of ongoing gas injection projects in the U.S. The Rangely CO2 flood began in 1986; currently, there are 372 active producer wells and 300 active injector wells, 259 of which are injecting CO2 using the water-alternating-gas (WAG) process. In the application of gas injection to heterogeneous reservoirs, oil recovery efficiency can be limited by poor conformance as an increasing proportion of the injected gas flows through higher permeability thief zones and/or fractures. The importance of conformance improvement has long been recognized at Rangely. The main problem being addressed is poor CO2 conformance due to preferential flow through the natural fracture network leading to premature gas breakthrough at the associated producers. This process increases operating costs and reduces oil recovery. The objective of the Rangely Conformance Improvement Team (CIT) is to improve conformance in order to reduce operating costs and increase the oil recovery to >1 billion bbl (>50% OOIP) compared to the current 815 million bbl (43% OOIP). A number of mechanical methods and chemical treatments have been employed to improve conformance at Rangely. While dual injection strings and selective injection equipment (SIE) have been used for improved injection profile control, chemical treatments using polymer gels2 and CO2 foam3 have been used to improve volumetric sweep efficiency and oil recovery. During the period 1994-1997, 49 injector wells were treated by placing a MARCIT™ gel4 into the fracture network.5 While these treatments have improved local sweep efficiency and oil recovery, economics limit the maximum treatment volume per injector well to 15?000-20?000 bbl. Certain regions of the Rangely reservoir require considerably larger treatment volumes to reduce the permeability of a larger volume of the fracture network and improve conformance in a larger volume of the well pattern.


2017 ◽  
Vol 10 (1) ◽  
pp. 180-186
Author(s):  
Siddhartha Biswas

In this paper the author introduces the notion of Z-weighted graph or Z-graph in Graph Theory, considers the Shortest Path Problem (SPP) in a Z-graph. The classical Dijkstra’s algorithm to find the shortest path in graphs is not applicable to Z-graphs. Consequently the author proposes a new algorithm called by Z-Dijkstra's Algorithm with the philosophy of the classical Dijkstra's Algorithm to solve the SPP in a Z-graph.


2018 ◽  
Vol 140 (5) ◽  
Author(s):  
Miao Zhang ◽  
Luis F. Ayala

Modeling fractured reservoirs, especially those with complex, nonorthogonal fracture network, can prove to be a challenging task. This work proposes a general integral solution applicable to two-dimensional (2D) fluid flow analysis in fractured reservoirs that reduces the original 2D problem to equivalent integral equation problem written along boundary and fracture domains. The integral formulation is analytically derived from the governing partial differential equations written for the fluid flow problem in reservoirs with complex fracture geometries, and the solution is obtained via solving system of equations that combines contributions from both boundary and fracture domains. Compared to more generally used numerical simulation methods for discrete fracture modeling such as finite volume and finite element methods, this work only requires discretization along the boundary and fractures, resulting in much fewer discretized elements. The validity of proposed solution is verified using several case studies through comparison with available analytical solutions (for simplified, single-fracture cases) and finite difference/finite volume finely gridded numerical simulators (for multiple, complex, and nonorthogonal fracture network cases).


Sign in / Sign up

Export Citation Format

Share Document