Visualization and Analysis of Surfactant Imbibition Into Oil-Wet Fractured Cores

SPE Journal ◽  
2016 ◽  
Vol 21 (01) ◽  
pp. 101-111 ◽  
Author(s):  
Mohammad Mirzaei ◽  
David A. DiCarlo ◽  
Gary A. Pope

Summary Imbibition of surfactant solution into the oil-wet matrix in fractured reservoirs is a complicated process that involves gravity, capillary, viscous, and diffusive forces. The standard experiment for testing imbibition of surfactant solution involves an imbibition cell, in which the core is placed in the surfactant solution and the recovery is measured vs. time. Although these experiments prove the effectiveness of surfactants, little insight into the physics of the problem is achieved. In this study, we performed water and surfactant flooding into oil-wet fractured cores and monitored the imbibition of the surfactant solution by use of computed-tomography (CT) scanning. From the CT images, the surfactant-imbibition dynamics as a function of height along the core was obtained. Although the waterflood only displaced oil from the fracture, the surfactant solution imbibed into the matrix; the imbibition is frontal, with the greatest imbibition rate at the bottom of the core, and the imbibition decreases roughly linearly with height. Experiments with cores of different sizes showed that increase in either the height or the diameter of the core causes decrease in imbibition and fractional oil-recovery rate. We also perform numerical simulations to model the observed imbibition. On the basis of the experimental measurements and numerical-simulation results, we propose a new scaling group that includes both the diameter and the height of the core. We show that the new scaling groups scale the recovery curves better than the traditional scaling group.

SPE Journal ◽  
2018 ◽  
Vol 23 (06) ◽  
pp. 2243-2259 ◽  
Author(s):  
Pengfei Dong ◽  
Maura Puerto ◽  
Guoqing Jian ◽  
Kun Ma ◽  
Khalid Mateen ◽  
...  

Summary Oil recovery in heterogeneous carbonate reservoirs is typically inefficient because of the presence of high-permeability fracture networks and unfavorable capillary forces within the oil-wet matrix. Foam, as a mobility-control agent, has been proposed to mitigate the effect of reservoir heterogeneity by diverting injected fluids from the high-permeability fractured zones into the low-permeability unswept rock matrix, hence improving the sweep efficiency. This paper describes the use of a low-interfacial-tension (low-IFT) foaming formulation to improve oil recovery in highly heterogeneous/fractured oil-wet carbonate reservoirs. This formulation provides both mobility control and oil/water IFT reduction to overcome the unfavorable capillary forces preventing invading fluids from entering an oil-filled matrix. Thus, as expected, the combination of mobility control and low-IFT significantly improves oil recovery compared with either foam or surfactant flooding. A three-component surfactant formulation was tailored using phase-behavior tests with seawater and crude oil from a targeted reservoir. The optimized formulation simultaneously can generate IFT of 10−2 mN/m and strong foam in porous media when oil is present. Foam flooding was investigated in a representative fractured core system, in which a well-defined fracture was created by splitting the core lengthwise and precisely controlling the fracture aperture by applying a specific confining pressure. The foam-flooding experiments reveal that, in an oil-wet fractured Edward Brown dolomite, our low-IFT foaming formulation recovers approximately 72% original oil in place (OOIP), whereas waterflooding recovers only less than 2% OOIP; moreover, the residual oil saturation in the matrix was lowered by more than 20% compared with a foaming formulation lacking a low-IFT property. Coreflood results also indicate that the low-IFT foam diverts primarily the aqueous surfactant solution into the matrix because of (1) mobility reduction caused by foam in the fracture, (2) significantly lower capillary entry pressure for surfactant solution compared with gas, and (3) increasing the water relative permeability in the matrix by decreasing the residual oil. The selective diversion effect of this low-IFT foaming system effectively recovers the trapped oil, which cannot be recovered with single surfactant or high-IFT foaming formulations applied to highly heterogeneous or fractured reservoirs.


SPE Journal ◽  
1900 ◽  
Vol 25 (02) ◽  
pp. 867-882
Author(s):  
Pengfei Dong ◽  
Maura Puerto ◽  
Guoqing Jian ◽  
Kun Ma ◽  
Khalid Mateen ◽  
...  

Summary The high formation heterogeneity in naturally fractured limestone reservoirs requires mobility control agents to improve sweep efficiency and boost oil recovery. However, typical mobility control agents, such as polymers and gels, are impractical in tight sub-10-md formations due to potential plugging issues. The objective of this study is to demonstrate the feasibility of a low-interfacial-tension (low-IFT) foam process in fractured low-permeability limestone reservoirs and to investigate relevant geochemical interactions. The low-IFT foam process was investigated through coreflood experiments in homogeneous and fractured oil-wet cores with sub-10-md matrix permeability. The performance of a low-IFT foaming formulation and a well-known standard foamer [alpha olefin sulfonate (AOS) C14-16] were compared in terms of the efficiency of oil recovery. The effluent ionic concentrations were measured to understand how the geochemical properties of limestone influenced the low-IFT foam process. Aqueous stability and phase behavior tests with crushed core materials and brines containing various divalent ion concentrations were conducted to interpret the observations in the coreflood experiments. Low-IFT foam process can achieve significant incremental oil recovery in fractured oil-wet limestone reservoirs with sub-10-md matrix permeability. Low-IFT foam flooding in a fractured oil-wet limestone core with 5-md matrix permeability achieved 64% incremental oil recovery compared to waterflooding. In this process, because of the significantly lower capillary entry pressure for surfactant solution compared to gas, the foam primarily diverted surfactant solution from the fracture into the matrix. This selective diversion effect resulted in surfactant or weak foam flooding in the tight matrix and hence improved the invading fluid flow in the matrix. Meanwhile, the low-IFT property of the foaming formulation mobilized the remaining oil in the matrix. This oil mobilization effect of the low-IFT formulation achieved lower remaining oil saturation in the swept zones compared with the formulation lacking low-IFT property with oil. The limestone geochemical instability caused additional challenges for the low-IFT foam process in limestone reservoirs compared to dolomite reservoirs. The reactions of calcite with injected fluids—such as mineral dissolution and the exchange of calcium and magnesium—were found to increase the Ca2+ concentration in the produced fluids. Because the low-IFT foam process is sensitive to brine salinity, the additional Ca2+ may cause potential surfactant precipitation and unfavorable over-optimum conditions. It, therefore, may cause injectivity and phase-trapping issues especially in the homogeneous limestone. Results in this work demonstrated that despite the challenges associated with limestone dissolution, the low-IFT foam process can remarkably extend chemical enhanced oil recovery (EOR) in fractured oil-wet tight reservoirs with matrix permeability as low as 5 md.


SPE Journal ◽  
2019 ◽  
Vol 24 (03) ◽  
pp. 1179-1191 ◽  
Author(s):  
Qingbang Meng ◽  
Jianchao Cai ◽  
Jing Wang

Summary Scaling of imbibition data is of essential importance in predicting oil recovery from fractured reservoirs. In this work, oil recovery by countercurrent spontaneous imbibition from 2D matrix blocks with different boundary conditions was studied using numerical calculations. The numerical results show that the shape of imbibition-recovery curves changes with different boundary conditions. Therefore, the imbibition curves could not be closely correlated with a constant parameter. A modified characteristic length was proposed by a combination of Ma et al. (1997) and theoretical characteristic length. The modified characteristic length is a function of imbibition time, and the shape of imbibition curves could be changed using the modified characteristic length. The overall imbibition curves were closely correlated using the modified characteristic length. Finally, the modified characteristic length was verified by experimental data for imbibition with different boundary conditions.


Fluids ◽  
2018 ◽  
Vol 3 (4) ◽  
pp. 70 ◽  
Author(s):  
Ahmad Zareidarmiyan ◽  
Hossein Salarirad ◽  
Victor Vilarrasa ◽  
Silvia De Simone ◽  
Sebastia Olivella

Geologic carbon storage will most likely be feasible only if carbon dioxide (CO2) is utilized for improved oil recovery (IOR). The majority of carbonate reservoirs that bear hydrocarbons are fractured. Thus, the geomechanical response of the reservoir and caprock to IOR operations is controlled by pre-existing fractures. However, given the complexity of including fractures in numerical models, they are usually neglected and incorporated into an equivalent porous media. In this paper, we perform fully coupled thermo-hydro-mechanical numerical simulations of fluid injection and production into a naturally fractured carbonate reservoir. Simulation results show that fluid pressure propagates through the fractures much faster than the reservoir matrix as a result of their permeability contrast. Nevertheless, pressure diffusion propagates through the matrix blocks within days, reaching equilibrium with the fluid pressure in the fractures. In contrast, the cooling front remains within the fractures because it advances much faster by advection through the fractures than by conduction towards the matrix blocks. Moreover, the total stresses change proportionally to pressure changes and inversely proportional to temperature changes, with the maximum change occurring in the longitudinal direction of the fracture and the minimum in the direction normal to it. We find that shear failure is more likely to occur in the fractures and reservoir matrix that undergo cooling than in the region that is only affected by pressure changes. We also find that stability changes in the caprock are small and its integrity is maintained. We conclude that explicitly including fractures into numerical models permits identifying fracture instability that may be otherwise neglected.


2020 ◽  
Vol 2020 ◽  
pp. 1-10
Author(s):  
Qiang Liu ◽  
Jianjun Liu ◽  
Guihong Pei ◽  
Zhengwen Zhu ◽  
Yun Lei

The fracture-control matrix unit (F-CMU) is a special body present in low-permeability fractured reservoirs that can be distinguished by a fracture system and a matrix system. The imbibition phenomenon of the F-CMU provides the possibility for secondary development of low-permeability fractured reservoirs because of the driving force including capillary force and gravity. However, the F-CMU is difficult to obtain during the field core drilling, which has limited the development for laboratory dynamic imbibition tests. Therefore, a new F-CMU reconstruction method is proposed in this study. According to the geometry and parameters, combining laser engraving technology, the fracture system is designed and engraved. Then, the F-CMU is established using a three-dimensional (3D) printed material called polyvinyl alcohol (PVA) as fracture support material which has a faster dissolution rate and causes less damage to the core due to water being the solvent. Finally, the porosity, permeability, and wettability of the matrix system and the T2 spectra from nuclear magnetic resonance (NMR) before and after reconstruction are measured. In addition, numerical simulation calculation of F-CMU permeability is performed. The results show that the characteristic parameters of the matrix system hardly change, indicating low damage to the core. The reconstructed fracture system is found on the T2 spectra, and the fracture permeability is consistent by comparing with the experimental and numerical simulation results. The permeability of the fracture system is about 104 orders of magnitude of the matrix system, which is closer to real core and meets the requirements needed for dynamic permeability experiments.


SPE Journal ◽  
2010 ◽  
Vol 16 (02) ◽  
pp. 411-428 ◽  
Author(s):  
Hamidreza Salimi ◽  
Johannes Bruining

Summary We use upscaling through homogenization to predict oil recovery from fractured reservoirs consisting of matrix columns, also called vertically fractured reservoirs (VFRs), for a variety of conditions. The upscaled VFR model overcomes limitations of the dual-porosity model, including the use of a shape factor. The purpose of this paper is to investigate three main physical aspects of multiphase flow in fractured reservoirs: reservoir wettability, viscosity ratio, and heterogeneity in rock/fluid properties. The main characteristic that determines reservoir behavior is the Péclet number that expresses the ratio of the average imbibition time divided by the residence time of the fluids in the fractures. The second characteristic dimensionless number is the gravity number. Upscaled VFR simulations, aimed at studying the mentioned features, add new insights. First, we discuss the results at low Péclet numbers. For only small gravity numbers, the effect of contact angle, delay time for the nonequilibrium capillary effect, the heterogeneity of the matrix-column size, and the matrix permeability can be ignored without appreciable loss of accuracy. The ultimate oil recovery for mixed-wet VFRs is approximately equal to the Amott index, and the oil production does not depend on the absolute value of the phase viscosity but on viscosity ratio. However, large gravity numbers enhance underriding, aggravated by large contact angles, longer delay times, and higher viscosity ratios. Layering can lead to an improvement or deterioration, depending on the fracture aperture and permeability distribution. At low Péclet numbers, the fractured reservoir behaves very similarly to a conventional reservoir and depends largely on the viscosity ratio and the gravity number. At high Péclet numbers, after water breakthrough, the oil recovery appears to be proportional to the cosine of the contact angle and inversely proportional to the sum of the oil and water viscosity. In addition, the mixed-wetting effect is more pronounced; there are significant influences of delay time (nonequilibrium effects), matrix permeability, matrix-column size, and the column-size distribution on oil recovery. At low gravity numbers and an effective length/thickness ratio larger than 10, the oil recovery is independent of the vertical-fracture-aperture distribution. For the same amount of injected water, the recovery at low Péclet numbers is larger than the recovery at high Péclet numbers.


Author(s):  
Amin Abolhasanzadeh ◽  
Ali Reza Khaz’ali ◽  
Rohallah Hashemi ◽  
Mohammadhadi Jazini

Without Enhanced Oil Recovery (EOR) operations, the final recovery factor of most hydrocarbon reservoirs would be limited. However, EOR can be an expensive task, especially for methods involving gas injection. On the other hand, aqueous injection in fractured reservoirs with small oil-wet or mixed-wet matrices will not be beneficial if the rock wettability is not changed effectively. In the current research, an unpracticed fabrication method was implemented to build natively oil-wet, fractured micromodels. Then, the efficiency of microbial flooding in the micromodels, as a low-cost EOR method, is investigated using a new-found bacteria, Bacillus persicus. Bacillus persicus improves the sweep efficiency via reduction of water/oil IFT and oil viscosity, in-situ gas production, and wettability alteration mechanisms. In our experiments, the microbial flooding technique extracted 65% of matrix oil, while no oil was produced from the matrix system by water or surfactant flooding.


2010 ◽  
Vol 13 (01) ◽  
pp. 109-117 ◽  
Author(s):  
V.. Er ◽  
T.. Babadagli

Summary CO2 injection has been applied in naturally fractured reservoirs (NFRs) for the purpose of enhanced oil recovery (i.e., the Wey-burn and Midale fields, Canada; the Wasson and Slaughter fields, USA; and the Bati Raman field, Turkey). The matrix part of these types of reservoirs could potentially be a good storage medium as well. Understanding the matrix/fracture interaction during this process and the dynamics of the flow in this dual-porosity system requires visual analyses. We mimicked fully miscible CO2 injection in NFRs using 2D models with a single fracture and oil (solute)/hydrocarbon solvent pairs. The focus was on the visual pore-scale analysis of miscibility interaction, breakthrough of solvent through fracture, transfer between matrix and fracture, and the dynamics of miscible displacement inside the matrix. First, matrix/fracture interaction was studied intensively using 2D glass-bead models experimentally. The model was prepared using acrylic sheets and glass beads saturated with oil as a porous medium while a narrow gap of 1-mm size containing filter paper served as a fracture. The first contact miscible solvent (pentane) was injected into the fracture, and the flow distribution was monitored using an image-acquisition and -processing system. The produced solvent and solute were continuously analyzed for compositional study. The effects of several parameters, such as flow rate, viscosity ratio (oil/solvent), and gravity, were studied. Next, the process was modeled numerically using a commercial compositional simulator, and the saturation distribution in the matrix was matched to experimental data. The key parameters in the matching process were the effective diffusion coefficients and the longitudinal and the transverse dispersivities. The diffusion coefficients were specified for each fluid, and dispersivities were assigned into gridblocks separately for the fracture and the matrix.


2009 ◽  
Vol 12 (02) ◽  
pp. 200-210 ◽  
Author(s):  
Benjamin Ramirez ◽  
Hossein Kazemi ◽  
Mohammed Al-kobaisi ◽  
Erdal Ozkan ◽  
Safian Atan

Summary Accurate calculation of multiphase-fluid transfer between the fracture and matrix in naturally fractured reservoirs is a crucial issue. In this paper, we will present the viability of the use of simple transfer functions to account accurately for fluid exchange resulting from capillary, gravity, and diffusion mass transfer for immiscible flow between fracture and matrix in dual-porosity numerical models. The transfer functions are designed for sugar-cube or match-stick idealizations of matrix blocks. The study relies on numerical experiments involving fine-grid simulation of oil recovery from a typical matrix block by water or gas in an adjacent fracture. The fine-grid results for water/oil and gas/oil systems were compared with results obtained with transfer functions. In both water and gas injection, the simulations emphasize the interaction of capillary and gravity forces to produce oil, depending on the wettability of the matrix. In gas injection, the thermodynamic phase equilibrium, aided by gravity/capillary interaction and, to a lesser extent, by molecular diffusion, is a major contributor to interphase mass transfer. For miscible flow, the fracture/matrix mass transfer is less complicated because there are no capillary forces associated with solvent and oil; nevertheless, gravity contrast between solvent in the fracture and oil in the matrix creates convective mass transfer and drainage of oil. Using the transfer functions presented in this paper, fracture- and matrix-flow calculations can be decoupled and solved sequentially--reducing the complexity of the computation. Furthermore, the transfer-function equations can be used independently to calculate oil recovery from a matrix block.


Energies ◽  
2019 ◽  
Vol 12 (19) ◽  
pp. 3699 ◽  
Author(s):  
Faisal Awad Aljuboori ◽  
Jang Hyun Lee ◽  
Khaled A. Elraies ◽  
Karl D. Stephen

Gravity drainage is one of the essential recovery mechanisms in naturally fractured reservoirs. Several mathematical formulas have been proposed to simulate the drainage process using the dual-porosity model. Nevertheless, they were varied in their abilities to capture the real saturation profiles and recovery speed in the reservoir. Therefore, understanding each mathematical model can help in deciding the best gravity model that suits each reservoir case. Real field data from a naturally fractured carbonate reservoir from the Middle East have used to examine the performance of various gravity equations. The reservoir represents a gas–oil system and has four decades of production history, which provided the required mean to evaluate the performance of each gravity model. The simulation outcomes demonstrated remarkable differences in the oil and gas saturation profile and in the oil recovery speed from the matrix blocks, which attributed to a different definition of the flow potential in the vertical direction. Moreover, a sensitivity study showed that some matrix parameters such as block height and vertical permeability exhibited a different behavior and effectiveness in each gravity model, which highlighted the associated uncertainty to the possible range that often used in the simulation. These parameters should be modelled accurately to avoid overestimation of the oil recovery from the matrix blocks, recovery speed, and to capture the advanced gas front in the oil zone.


Sign in / Sign up

Export Citation Format

Share Document