Toward Field-Scale Wettability Modification—The Limitations of Diffusive Transport

2008 ◽  
Vol 11 (03) ◽  
pp. 633-640 ◽  
Author(s):  
Martin Stoll ◽  
Jan Hofman ◽  
Dick J. Ligthelm ◽  
Marinus J. Faber ◽  
Paul van den Hoek

Summary Densely-fractured oil-wet carbonate fields pose a true challenge for oil recovery that traditional primary and secondary processes fail to meet. The difficulty arises from the combination of two unfavorable characteristics: First, the dense fracturing frustrates an efficient waterflood; second, because of the oil-wetness, the water pressure exceeds the oil pressure inside the matrix blocks, thus inhibiting spontaneous imbibition of water. In the past decade, using a new class of surfactants, enhanced oil recovery (EOR) researchers have studied the options to chemically revert the wettability of carbonate rock without drastically decreasing the oil-water interfacial tension. These chemicals, termed "wettability modifiers" (WMs), effectively reverse the sign of capillary pressure at the prevalent saturation. With the oil pressure exceeding the water pressure, the capillary pressure becomes the driving force for oil expulsion from the matrix and into the fracture system. Previous publications on chemical wettability modification focused on the performance of different chemical wettability modifiers for a chosen rock/oil/brine system. In some cases, they demonstrated an almost full oil recovery from core plugs. Little attention, however, has been given to the mechanism underlying the transport of the chemical into the matrix block and to the proper scaling of laboratory results to reservoir size. The present study aims to demonstrate that imbibition after wettability modification is diffusion-limited. To this end, the recovery profiles for spontaneous capillary imbibition, as well as for imbibition after wettability modification, are calculated. The results are then used to compare with the data of Amott cell imbibition experiments. It is confirmed that in both cases, the cumulative recovery is initially proportional to the square root of time. Imbibition after wettability modification, however, takes approximately 1,000 times longer than spontaneous capillary imbibition into a water-wet medium. The slow recovery observed in the case of imbibition after wettability modification is in excellent agreement with the assumption that, in the absence of significant spontaneous imbibition, the WM, to unfold its action, must first diffuse into the porous medium. In any diffusion process, the time scale is linked to the square of the length scale of the medium. Therefore, it would take up to 1,000 times longer (an equivalent of 200 years) before the same recovery is obtained from a meter-scale matrix block as is obtained from a centimeter-scale plug in a laboratory in 100 days. Consequently, unless a significantly faster transport mechanism for the wettability modifier is identified, or unless viscous forces or buoyancy enable forced imbibition, the chemical wettability modification of fractured oil-wet carbonate rock does not provide an economically interesting opportunity. Introduction Rock fractures provide comparatively highly permeable flow paths through oil reservoirs. In a densely fractured reservoir, the permeability contrast between the fracture network and the oil-bearing matrix can be significant. In that case, the viscous pressure differential across individual matrix blocks can be too small to release oil from the blocks under waterflood, thus leading to a poor recovery. Depending on the wetting state of the matrix and its initial water saturation, Swi, capillary action can cause imbibition of water up to a "spontaneous" equilibrium saturation, commonly denoted as Sspw. At this saturation, however, the capillary pressure inside the matrix block coincides with that in the fracture, and the recovery ceases. Experience has shown that carbonate fields often range from intermediate-wet to preferentially oil-wet (Treiber et al. 1972; Chilingar and Yen 1983), which is synonymous with Sspw being close or equal to Swi ; thus, they exhibit very limited recovery during primary and secondary production. Recently, a new EOR technique, designed specifically to tackle the challenges outlined previously, has been suggested by Austad and coworkers (Austad and Milter 1997; Standnes and Austad 2000a, b). In their pioneering work, these authors show that certain chemicals, when dissolved in the surrounding brine, can initiate water imbibition into oil-saturated core plugs and, hence, lead to the recovery of oil. One possible mechanism that explains these observations is the solubilization of adsorbed hydrocarbon components from the pore surface—as demonstrated by an atomic force microscopy study by Kumar et al. (2005), this exposes the intrinsically hydrophilic matrix. Another possibility is the formation of an additional chemical layer covering the adsorbed hydrophobic material. In either case, the pore surface becomes more hydrophilic, and the wettability of the matrix is thus modified. In a capillary rise experiment into parallel plates, Kumar et al. also observed different time scales for different types of wettability-modifying chemicals (2005). Using the cationic wettability modifier dodecyl trimethyl ammonium bromide (DTAB, also known as C12TAB), Standnes and Austad deduced that wettability modification was achieved through the comparatively slow process of partitioning the chemical into the oil phase, followed by desorption and solubilization of anionic hydrocarbon components (2000a, b). Shen et al. (2006) and Rao et al. (2006) measured the effect of surfactants on the relative water/oil permeabilities at different interfacial tensions. Wu et al. (2006) studied the properties and ranked the efficiency of chemical model compounds, based on their chemical structure, to modify the wettability and enhance recoveries. Several groups have taken initiative to model wettability modification in numerical simulators (Adibhatla et al. 2005; Delshad et al. 2006). So far, no significant attention has been given to time dependence and to the subsequent upscaling of the laboratory results to matrix block scale. This subject will be addressed in the present work. The structure of the article is as follows: In the Theory section, the basic results for countercurrent capillary imbibition will be briefly reviewed and compared to Fick's law of molecular diffusion. The oil recovery as a function of time for both capillary imbibition and imbibition after wettability modification will be predicted. The experimental approach to imbibition at different wetting situations will be described in the section Materials and Preparation. The recovery results will then be analyzed using the previously derived equations. Finally, tentative conclusions for the upscaling will be drawn.

2021 ◽  
Author(s):  
Aymen AlRamadhan ◽  
Yildiray Cinar ◽  
Arshad Hussain ◽  
Nader BuKhamseen

Abstract This paper presents a numerical study to examine how the interplay between the matrix imbibition capillary pressure (Pci) and matrix-fracture transfer affects oil recovery from naturally-fractured reservoirs under waterflooding. We use a dual-porosity, dual-permeability (DPDP) finite difference simulator to investigate the impact of uncertainties in Pci on the waterflood recovery behavior and matrix-fracture transfer. A comprehensive assessment of the factors that control the matrix-fracture transfer, namely Pci, gravity forces, shape factor and fracture-matrix permeabilities is presented. We examine how the use of Pci curves in reservoir simulation can affect the recovery assessment. We present two conceptual scenarios to demonstrate the impact of spontaneous and forced imbibition on the flood-front movement, waterflood recovery processes, and ultimate recovery in the DPDP reservoir systems of varying reservoir quality. The results demonstrate that the inclusion of Pci in reservoir simulation delays the breakthrough time due to a higher displacement efficiency. The study reveals that the matrix-fracture transfer is mainly controlled by the fracture surface area, fracture permeability, shape factor, and the uncertainty in Pci. We underline a discrepancy among various shape factors proposed in the literature due to three main factors: (1) the variations in matrix-block geometries considered, (2) how the physics of imbibition forces that control the multiphase fluid transfer is captured, and (3) how the assumption of pseudo steady-state flow is addressed.


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.


2021 ◽  
Author(s):  
Sarah Abdullatif Alruwayi ◽  
Ozan Uzun ◽  
Hossein Kazemi

Abstract In this paper, we will show that it is highly beneficial to model dual-porosity reservoirs using matrix refinement (similar to the multiple interacting continua, MINC, of Preuss, 1985) for water displacing oil. Two practical situations are considered. The first is the effect of matrix refinement on the unsteady-state pressure solution, and the second situation is modeling water-oil, Buckley-Leverett (BL) displacement in waterflooding a fracture-dominated flow domain. The usefulness of matrix refinement will be illustrated using a three-node refinement of individual matrix blocks. Furthermore, this model was modified to account for matrix block size variability within each grid cell (in other words, statistical distribution of matrix size within each grid cell) using a discrete matrix-block-size distribution function. The paper will include two mathematical models, one unsteady-state pressure solution of the pressure diffusivity equation for use in rate transient analysis, and a second model, the Buckley-Leverett model to track saturation changes both in the reservoir fractures and within individual matrix blocks. To illustrate the effect of matrix heterogeneity on modeling results, we used three matrix bock sizes within each computation grid and one level of grid refinement for the individual matrix blocks. A critical issue in dual-porosity modeling is that much of the fluid interactions occur at the fracture-matrix interface. Therefore, refining the matrix block helps capture a more accurate transport of the fluid in-and-out of the matrix blocks. Our numerical results indicate that the none-refined matrix models provide only a poor approximation to saturation distribution within individual matrices. In other words, the saturation distribution is numerically dispersed; that is, no matrix refinement causes unwarranted large numerical dispersion in saturation distribution. Furthermore, matrix block size-distribution is more representative of fractured reservoirs.


2021 ◽  
pp. 1-18
Author(s):  
Takaaki Uetani ◽  
Hiromi Kaido ◽  
Hideharu Yonebayashi

Summary Low-salinity water (LSW) flooding is an attractive enhanced oil recovery (EOR) option, but its mechanism leading to EOR is poorly understood, especially in carbonate rock. In this paper, we investigate the main reason behind two tertiary LSW coreflood tests that failed to demonstrate promising EOR response in reservoir carbonate rock; additional oil recovery factors by the LSW injection were only +2% and +4% oil initially in place. We suspected either the oil composition (lack of acid content) or the recovery mode (tertiary mode) was inappropriate. Therefore, we repeated the experiments using an acid-enriched oil sample and injected LSW in the secondary mode. The result showed that the low-salinity effect was substantially enhanced; the additional oil recovery factor by the tertiary LSW injection jumped to +23%. Moreover, it was also found that the secondary LSW injection was more efficient than the tertiary LSW injection, especially in the acid-enriched oil reservoir. In summary, it was concluded that the total acid number (TAN) and the recovery mode appear to be the key successful factors for LSW in our carbonate system. To support the conclusion, we also performed contact angle measurement and spontaneous imbibition tests to investigate the influence of acid enrichment on wettability, and moreover, LSW injection on wettability alteration.


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.


1962 ◽  
Vol 2 (02) ◽  
pp. 177-184 ◽  
Author(s):  
C.C. Mattax ◽  
J.R. Kyte

Abstract Previous workers have developed differential equations to describe oil displacement by water imbibition, but have not explicitly defined the relationship between recovery behavior for a single reservoir matrix block and its size. In the present work, imbibition theory is extended to show that the time required to recover a given fraction of the oil from a matrix block is proportional to the square of the distance between fractures. Using this relationship, recovery behavior for a large reservoir matrix block is predicted from an imbibition test on a small reservoir core sample. The prediction is then extended to analyze recovery behavior for fractured - matrix, water - drive reservoirs in which imbibition is the dominant oil-producing mechanism. Experimental data are presented to support the basic imbibition theory relating matrix block size, fluid viscosity level and permeability to recovery behavior. Introduction Imbibition has long been recognized as an important factor in recovering oil from water - wet, fractured-matrix reservoirs subjected to water flood or water drive. Recently, two approaches have been published which might be used to predict imbibition oil-recovery behavior for reservoir - sized matrix blocks. Graham and Richardson used a synthetic model to scale a single element of a fractured-matrix reservoir. Blair, on the other hand, used numerical techniques to solve the differential equations describing imbibition in linear and radial systems. This latter method requires auxiliary experimental data in the form of capillary pressure and relative permeability functions. These two approaches, i.e., synthetic models and numerical techniques, have been used to study a variety of reservoir fluid-flow problems. One purpose of this work is to present, with experimental verification, a third method for predicting imbibition oil recovery for large reservoir matrix blocks. This method uses scaled imbibition tests on small reservoir core samples to predict field performance. The imbibition tests are easier to perform than the capillary pressure and relative permeability tests required to apply the numerical method. Furthermore, when suitably preserved reservoir-rock samples are used, the properties of the laboratory system are the same as those of the field. This offers an important advantage over the use of synthetic models because there is usually some question as to how accurately reservoir-rock properties can be duplicated in such models. Based on the recovery behavior for a unit matrix block, an analysis is presented to predict oil recovery for a fractured, water-drive reservoir made up of many such unit blocks. In the analysis, it is assumed that the flow resistance and the volume of the reservoir fracture system are negligible compared with that of the porous matrix. These assumptions are generally consistent with observed characteristics of many fractured-matrix reservoirs, and have been employed in previous studies. It is further assumed that the effect of gravity on flow in the matrix blocks is negligible. On first thought, the latter assumption might appear to seriously limit application of the method. However, in a fractured reservoir, the effect of gravity on flow in a matrix block will be restricted by the height of the block. Furthermore, matrix permeabilities are often very low (10 md or less) in such reservoirs. This means that capillary or imbibition forces will be large, thus tending to minimize the relative importance of gravity. For these reasons, imbibition should be the dominant oil-recovery mechanism in many fractured-matrix, water-drive reservoirs. The predictive method presented in this report is applicable to such reservoirs. SPEJ P. 177^


SPE Journal ◽  
2018 ◽  
Vol 23 (04) ◽  
pp. 1343-1358 ◽  
Author(s):  
Somayeh Karimi ◽  
Hossein Kazemi

Summary To understand the flow and transport mechanisms in shale reservoirs, we needed reliable core-measured data that were not available to us. Thus, in 2014, we conducted a series of diverse experiments to characterize pores and determine the flow properties of 12 Middle Bakken cores that served as representatives for unconventional low-permeability reservoirs. The experiments included centrifuge, mercury-intrusion capillary pressure (MICP), nitrogen adsorption, nuclear magnetic resonance (NMR), and resistivity. From the centrifuge measurements, we determined the mobile-fluid-saturation range for water displacing oil and gas displacing oil in addition to irreducible fluid saturations. From MICP and nitrogen adsorption, we determined pore-size distribution (PSD). Finally, from resistivity measurements, we determined tortuosity. In addition to flow characterization, these data provided key parameters that shed light on the mechanisms involved in primary production and the enhanced-oil-recovery (EOR) technique. The cores were in three conditions: clean, preserved, and uncleaned. The hydrocarbon included Bakken dead oil and decane, and the brine included Bakken produced water and synthetic brine. After saturating the cores with brine or oil, a set of drainage and imbibition experiments was performed. NMR measurements were conducted before and after each saturation/desaturation step. After cleaning, PSD was determined for four cores using MICP and nitrogen-adsorption tests. Finally, resistivity was measured for five of the brine-saturated cores. The most significant results include the following: Centrifuge capillary pressure in Bakken cores was on the order of hundreds of psi, both in positive and negative range. Mobile-oil-saturation range for water displacing oil was very narrow [approximately 12% pore volume (PV)] and much wider (approximately 40% PV) for gas displacing oil. In Bakken cores, oil production by spontaneous imbibition of high-salinity brine was small unless low-salinity brine was used for spontaneous imbibition. Resistivity measurements yielded unexpectedly large tortuosity values (12 to 19), indicating that molecules and bulk fluids have great difficulty to travel from one point to another in shale reservoirs.


2021 ◽  
pp. 014459872098420
Author(s):  
Qi Zhang ◽  
Xinyue Wu ◽  
Yingfu He ◽  
Qingbang Meng

Spontaneous imbibition is an important mechanism of oil recovery from fractured reservoirs and unconventional reservoirs. Oil is produced by combining co- and counter-current imbibition when the matrix blocks was partially covered by water. In this paper, we focused on the effect of viscosity ratios on oil production by spontaneous imbibition and established the numerical model for one-dimensional linear imbibition with TEO-OW boundary conditions, which was validated by the experimental data. The effect of viscosity ratio on co- and counter-current imbibition is investigated and scaling result of the imbibition recovery curve for wide range of viscosity ratio using the conventional scaling equation was tested, which indicates that the close correlation was achieved only when oil-water viscosity ratios are higher. Then, a modified scaling equation was developed based on the piston-like assumption for one-dimensional co-current imbibition and close correlation of imbibition recovery curves was achieved when viscosity ratios are lower. Finally, correlation of imbibition recovery curves was improved for wide range of viscosity ratios by combining conventional and modified scaling equation. Results show that since the shape of imbibition recovery curves is not similar for different viscosity ratios, it is difficult to obtain the perfect correlation using the constant viscosity term.


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

Summary Flow modeling in fractured reservoirs is largely confined to the so-called sugar-cube model. Here, we consider a situation where matrix blocks are connected to neighboring blocks so that part of the global flow occurs only in the matrix domain. We call this a partially fractured reservoir (PFR). As opposed to the sugar-cube model, global flow in the matrix blocks plays an important role in the PFR when the interconnections between the matrix blocks are sufficiently large. We apply homogenization to derive an upscaled model for PFRs that combines dual-porosity and dual-permeability concepts simultaneously. We formulate a well-posed fully implicit 3D upscaled numerical model and investigate oil-recovery mechanisms for different dimensionless characteristic numbers. As we found previously for the sugar-cube model, the Péclet number, defined here as the ratio of the capillary diffusion time in the matrix to the residence time of the fluids in the fracture, plays a crucial role. The gravity number and specific fracture/matrix-interface area play a secondary role. For low Péclet numbers and high gravity numbers, the results are sensitive to gravity and water-injection rates, but relatively insensitive to the specific fracture/matrix-interface area, matrix-block size, and reservoir geometry (i.e., sugar cube vs. PFR). At low Péclet numbers and high gravity numbers, ECLIPSE simulations using the Barenblatt or Warren and Root (BWR) approach give poor predictions and overestimate the oil recovery, but, at short injection times, show good agreement with the solution of the PFR model at intermediate Péclet numbers. At high Péclet numbers, the results are relatively insensitive to gravity, but sensitive to the other conditions mentioned. In particular, when the specific fracture/matrix-interface area is large, it enhances the imbibition and, consequently, leads to a higher oil production. If this specific interface area is small, it leads to a considerable retardation of the imbibition process, which leads to an earlier water breakthrough and lower oil recovery. The BWR (commercial simulator) simulations and the sugar-cube model result in inaccurate predictions of the oil-production rate at high Péclet numbers. This can be inferred from the discrepancy with respect to the PFR model for which we assert that it accurately predicts the oil recovery. We conclude that, at low Péclet numbers and large gravity numbers, it is advantageous to increase the water-injection rate to improve the net present value. However, at high Péclet numbers, increasing the flow rate may lead to uneconomical water cuts.


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.


Sign in / Sign up

Export Citation Format

Share Document