Skauge, T. (CIPR Uni Research) | Skauge, A. (CIPR Uni Research) | Salmo, I. C. (CIPR Uni Research) | Ormehaug, P. A. (CIPR Uni Research) | Al-Azri, N. (PDO) | Wassing, L. M. (Shell Global Solutions International BV) | Glasbergen, G. (Shell Global Solutions International BV) | Van Wunnik, J. N. (Shell Global Solutions International BV) | Masalmeh, S. K. (Shell Global Solutions International BV)
Polymer injectivity is a critical parameter for implementation of polymer flood projects. An improved understanding of polymer injectivity is important in order to facilitate an increase in polymer EOR implementation. Typically, injectivity studies are performed using linear core floods. Here we demonstrate that polymer flow in radial and linear models may be significantly different and discuss the concept in theoretical and experimental terms.
Linear core floods using partially hydrolyzed polyacrylamides (HPAM) were performed at various rates to determine in-situ viscosity and polymer injectivity. Radial polymer floods were performed on Bentheimer discs (30 cm diameter, 2-3 cm thickness) with pressure taps distributed between a central injector and the perimeter production well. The in-situ rheological data are also compared to bulk rheology. The experimental set up allowed a detailed analysis of pressure changes from well injection to production line in the radial models and using internal pressure taps in linear cores.
Linear core floods show degradation of polymer at high flow rates and a severe degree of shear thickening leading to presumably high injection pressures. This is in agreement with current literature. However, the radial injectivity experiments show a significant reduction in differential pressure compared to the linear core floods. Onset of shear thickening occurs at significantly higher flow velocities than for linear core floods. These data confirm that polymer flow is significantly different in linear and radial flow. This is partly explained by the fact that linear floods are being performed at steady state conditions, while radial injections go through transient (unsteady state) and semi-transient pressure regimes.
History matching of polymer injectivity was performed for radial injection experiments. Differences in polymer injectivity are discussed in the framework of theoretical and experimental considerations. The results may have impact on evaluation of polymer flood projects as polymer injectivity is a key risk factor for implementation.
Enhanced oil displacement in a reservoir is highly affected by wettability alterations in conjunction with the lowering of viscosities during steam assisted gravity drainage (SAGD) for bitumen extraction. The impartation of energy in the form of heat to the fluid by injecting steam triggers an alteration to a more water-wet state during SAGD. However, the presence of three distinct phases in the reservoir has implications for the effective modeling of the complex fluid dynamics. Dependency of the relative permeability endpoints on the temperature realized as a function of the introduction of steam is difficult to model. Optimization of any steam process requires simulation in order to adequately characterize years of flow and so a model that is capable of representing three phase flow is necessary. To obtain this a pseudo-two phase relative permeability is proposed that assumes fractional flow theory is valid and treats the experiments as a waterflood.
In this study, experimental recovery data for two SAGD experiments and one hot water flood are empirically matched by manipulating relative permeabilities. The analytical approach implemented allows for the representation of fluid flow in the reservoir by achieving a pseudo-two phase relative permeability that results in comparable performance to the experiments. Waterflooding techniques were utilized which allowed for the negation of the steam phase in the model and so two-phase flow was established.
The sensitivity of the relative permeability curves to temperature change results in the inability to formulate a generic three-phase curve and so the pseudo-two phase curve is valuable for the purpose of simulation. The methodology presented enables the formulation of a simplified relative permeability that is unique to each process used and in that specific location. The model that was established was validated and proven credible by the good match with the experimentally obtained values.
Luo, Haishan (The University of Texas at Austin) | Mohanty, Kishore K. (The University of Texas at Austin) | Delshad, Mojdeh (The University of Texas at Austin) | Pope, Gary A. (The University of Texas at Austin)
Upscaling of unstable immiscible flow remains an unsolved challenge for the oil industry. The absence of a reliable upscaling approach greatly hinders the effective reservoir simulation and optimization of heavy oil recoveries using waterflood, polymer flood and other chemical floods, which are inherently unstable processes. The difficulty in upscaling unstable flow lies in estimating the propagation of fingers smaller than the gridblock size. Using classical relative permeabilities obtained from stable flow analysis can lead to incorrect oil recovery and pressure drop in reservoir simulations.
In a recent study based on abundant experimental data, it is found that the heavy-oil recovery by waterfloods and polymer floods has a power-law correlation with a dimensionless number (named viscous finger number in this paper), which is a combination of viscosity ratio, capillary number, permeability, and the cross-section area of the core. Based upon this important finding as well as the features of unstable immiscible floods, an effective-finger model is developed in this paper. A porous medium domain is dynamically identified as three effective zones, which are two-phase flow zone, oil single-phase flow zone, and bypassed oil (isolated oil island) zone, respectively. Flow functions are derived according to effective flows in these zones. This new model is capable of history-matching a set of heavy-oil waterflood corefloods under different viscosity ratios and injection rates. Model parameters obtained from the history match also have a power-law correlation with the viscous finger number.
The build-up of this correlation contains reasonable physical meanings to quantitatively characterize the upscaled behavior of viscous fingering effects. Having such a correlation enables the estimation of model parameters in any gridblock of the reservoir by knowing the local viscous finger number in reservoir simulations. The model is applied to several heavy-oil field cases with waterfloods and polymer floods with different heterogeneities. Oil recovery in water flooding of viscous oils is overpredicted by classical simulation methods which do not incorporate viscous fingering properly. Simulation results indicate that the new model reasonably differentiates the oil recoveries at different viscous finger numbers, e.g., lower injection rate leads to higher oil recovery. In contrast, classical simulations obtain close oil recoveries under different injection rates or degrees of polymer shear-thinning, which is apparently incorrect for unstable floods. Moreover, coarse-grid simulations using the new model are able to obtain consistent saturation and pressure maps with fine-grid simulations when the correlation lengths are not smaller than the coarse gridblock size. Furthermore, it is well captured by the model that the shear-shinning polymer solution can strengthen the fingering in high-permeability regions due to increased capillary number and viscosity ratio, which is not observed in waterflood. As a whole, the new model shows encouraging capability to simulate unstable water and polymer floods in heavy oil reservoirs, and hence can facilitate the optimization of heavy-oil EOR projects.
Skrettingland, K. (Statoil ASA) | Ulland, E. N. (Statoil ASA) | Ravndal, O. (Statoil ASA) | Tangen, M. (Statoil ASA) | Kristoffersen, J. B. (Statoil ASA) | Stenerud, V. R. (Statoil ASA) | Dalen, V. (Statoil ASA) | Standnes, D. C. (Statoil ASA) | Fevang, Ø. (Statoil ASA) | Mevik, K. M. (Knutsen Subsea Solutions) | McIntosh, N. (Knutsen Subsea Solutions) | Mebratu, A. (Halliburton) | Melien, I. (Halliburton) | Stavland, A. (Intl Research Inst of Stavanger)
Declining oil production and increasing water cut in mature fields highlight the need for improved conformance control. Here we report on a successful in-depth water diversion treatment using sodium silicate to increase oil recovery at the Snorre field, offshore Norway, utilizing a new operational concept of using a stimulation vessel as a platform for the large-scale injection into a subsea well. A custom modified 35,000 DWT shuttle tanker was employed for the field pilot. This paper describes the vessel preparations and the large-scale interwell silicate injection operation. The operational aspects of the large-scale interwell silicate injection include; identification of injection vessel requirements, major vessel modifications, chemical logistic, general logistics to site, major equipment set-up on vessel, subsea connection, mixing and pumping schedules, onsite QC, and real time monitoring. Experience from these operations and lessons learned are included in this paper.
After the injection of approximately 400,000 Sm3 (113,000 Sm3 preflush, followed by 240,000 Sm3 of sodium silicate gelant and 49,000 Sm3 of postflush fluid) at injection rates up to 4,000 Sm3/d, the injection from the vessel was stopped and the well was put on regular seawater injection. Following more than two years of regular production, transient pressure measurements, tracer testing and water cut data are presented from the ongoing comprehensive data acquisition program. These results demonstrate clearly the achieved in-depth flow diversion through a delayed breakthrough of injected tracers and lower water cut in the relevant production well.
Although geochemical reactions are the fundamental basis of the alkaline/surfactant/polymer (ASP) flooding, their importance is commonly overlooked and not fully assessed. Common assumptions made when modeling geochemical reactions in ASP floods include: 1) ideal solution (i.e., using molalities rather than ion activities) for the water and aqueous geochemical species 2) limiting the number of reactions (i.e., oil/alkali and alkali consumptions) rather than including the entire thermodynamically-equilibrated system 3) ignoring the effect of temperature and pressure on reactions 4) local equilibrium ignoring the kinetics. To the best of our knowledge, the significance of these assumptions has never been discussed in the literature. In this paper we investigate the importance of geochemical reactions during alkaline/surfactant/polymer floods using a comprehensive tool in the sense of surfactant/soap phase behavior as well as geochemistry.
We coupled the United States Geological Survey (USGS) state-of-the-art geochemical tool, with 3D flow and transport chemical flooding module of UTCHEM. This geochemical module includes several thermodynamic databases with various geochemical reactions, such as ion speciation by applying several ion-association aqueous models, mineral, solid-solution, surface-complexation, and ion-exchange reaction. It has capabilities of saturation index calculation, reversible and irreversible reactions, kinetic reaction, mixing solutions, inverse modeling and includes impacts of temperature and pressure on reaction constants and solubility products. The chemical flood simulator has a three phase (water, oil, microemulsion) phase behavior package for the mixture of surfactant/soap, oil, and water as a function of surfactant/soap, salinity, temperature, and co-solvent concentration. Hence, the coupled software package provides a comprehensive tool to assess the significance of geochemical assumptions typically imposed in modeling ASP floods. Moreover, this integrated tool enables modeling of variations in mineralogy present in reservoir rocks. We parallelized the geochemistry module of this coupled simulator for large-scale reservoir simulations.
Our simulation results show that the assumption of ideal solution overestimates ASP oil recovery. Assuming only a subset of reactions for a coupled system is not recommended, particularly when a large number of geochemical species is involved, as is the case in realistic applications of ASP. Reservoir pressure has a negligible effect but temperature has a significant impact on geochemical calculations. Although mineral reaction kinetics is largely a function of the temperature and in-situ water composition, some general conclusions can be drawn as follows: to a good approximation, minerals with slow rate kinetic reaction (e.g., quartz) can be excluded when modeling ASP laboratory floods. However, minerals with fast rate kinetic reactions (e.g., calcite) must be included when modeling lab results. On the other hand, in modeling field-scale applications, local equilibrium assumption (LEA) can be applied for fast rate kinetic minerals, whereas kinetics should be used for slow rate kinetic minerals.
Preliminary studies have been done to characterize rock-fluid properties, and flow mechanisms in the shale reservoirs. Most of these studies, through modifying methods used for conventional reservoirs, fail to capture dynamic features of shale rock and fluids in confined nano-pore space. In unconventional reservoirs, interactions between the wall of shale and the contained fluid significantly affect phase and flow behaviors. The inability to model capillarity with the consideration of pore size distribution characteristics using commercial software may lead to an inaccurate oil production performance in Bakken. This paper presents a novel formulation that consistently evaluates capillary force and adsorption using pore size distribution (PSD) directly from core measurements. The new findings could better address differences in flow mechanisms in unconventional reservoirs, and thus lead to an optimized IOR practice.
Wang, D. (University of North Dakota) | Dawson, M. (Statoil Gulf Services LLC) | Butler, R. (University of North Dakota) | Li, H. (Statoil Gulf Services LLC) | Zhang, J. (University of North Dakota) | Olatunji, K. (University of North Dakota)
With the recent dramatic drop in oil price, production from ultra-tight resources, like the Bakken formation, may drop substantially. Since expenditures for drilling, completion, and fracking have already been made, existing wells will continue to flow, but oil rates will decline—rapidly in many cases. In a low oil-price environment, what can be done to sustain oil production from these tight formations?
We are testing a surfactant imbibition process to recovery oil from shales. We measured surfactant imbibition rates and oil recovery values in laboratory cores from the Bakken shale. After optimizing surfactant formulations at reservoir conditions, we observed oil recovery values up to 10–20% OOIP incremental over brine imbibition. However, whether or not surfactant imbibition will be a viable recovery process depends on achieving sufficiently high oil production rates in a field setting—which requires that we identify conditions that will maximize imbibition rate, as well as total oil recovery. In this paper, we describe laboratory evaluations of oil recovery using different core plugs. These recovery studies involved
(1) surfactant formulation optimization on concentration, salinity and pH, (2) characterization of phase behavior, (3) spontaneous imbibition, and (4) forced imbibition (flooding) with gravity drainage assistance.
In preserved cores, we observed: (1) Formulations using 0.1% surfactant concentration at 4% TDS salinity showed favorable oil recoveries (up to 40% OOIP). (2) Generally, surfactant formulations at optimal concentration and salinity were stable at high temperature (115°C). (3) Injectivity/permeability enhancements up to 75 percent occurred after acidification using acetic acid or HCl. (4) Wettability alteration is the dominant mechanism for surfactant imbibition. Of course, actions that increase fracture width will aid gravity drainage and oil recovery. This information is being used to design and implement a field application of the surfactant imbibition process in an ultra-tight resource.
Once a shale gas condensate reservoir is produced, the reservoir pressure falls below the dew-point pressure, and the condensate liquid will be formed in the pore space; the condensate can then accumulate near the wellbore. This condensate blockage would reduce the gas relative permeability and decrease the gas production. The condensate is formed by the heavy components in reservoir fluid, and these heavy componenets are very valuable economically in the industry. Therefore, operators are seeking ways to maximize condensate recovery from gas-condensate reservoirs.
Huff-n-puff gas injection is an effective approach for enhancing condensate recovery in shale gas condensate reservoirs, as shown by our earlier papers (
In this experimental study, a binary gas condensate mixture was used to investigate the dominant mechanism. The core was saturated with a gas condensate mixture at 2200 psi to simulate the initial reservoir condition. Then, the pressure was depleted to 1500 psi, which was lower than the dew point pressure. During the depletion, the produced gas was collected in a vacuumed gas sample bag. After depletion, the huff-n-puff method was applied. After every cycle of huff-n-puff, the produced gas was collected. GC was used to analyze the compositions of the different gas samples. Also, a CT scanner was used to determine the condensate saturation in the core. From the GC analysis, by comparing the gas sample after primary depletion with the gas sample after the first cycle, it was found that the heavy component-butane increased significantly. This means that most of the heavy components of condensate were revaporized and flew out with the dry gas. This proves the revaporization mechanism of the huff-n-puff gas injection.
Our experiment results show that huff-n-puff was an effective way to enhance the condensate recovery, and revaporization was the main mechanism of huff-n-puff. When the pressure was increased in the huff process, the heavy components were revaporized and flowed out with gas in the puff process. Though in the gas flooding method, the reservoir pressure was also increased, but the near-wellbore pressure was not increased very much in the shale gas reservoirs; thus, heavy components would still be formed near the wellbore. However, in the huff-n-puff method, because of the same well, the pressure near the wellbore would be higher than the dew point pressure in the beginning of production process. Therefore, the heavy component would be recovered with gas. Our GC results visually showed the revaporization mechanism of huff-n-puff in the shale gas condensate.
Achieving maximum oil recovery utilizing CO2 has limitations when operating at, or very close, to the Minimum Miscibility Pressure (MMP) of the CO2 in the oil. A modular source of CO2 would allow Enhanced Oil Recovery (EOR) flooding of "stranded" and shallow reservoirs. Unfortunately, modular sources of CO2 production often include CO and N2 mixed with the CO2. Thus, testing for EOR application of a mixed gas-containing CO2, N2, and CO was initiated.
Bench scale testing using Rising Bubble Apparatus (RBA), Slim Tubes, and linear core flood have been conducted on oils ranging from 16-42° gravities having viscosities of 0.5-280 cp. All tests were conducted at reservoir temperatures and pressures. CO, being a strong reducing agent, was further tested on reservoir rock containing swelling clays with hydrated ferric hydroxides. Due to the apparent reduction of the ferric hydroxide, and the liberation of its water of hydration, an increase in matrix permeability and clay stabilization, was observed.
For most oils tested, the CO2/CO mixture increased rate of oil recovery by 2-3X, using only 50-60% as much gas/bo as compared to pure CO2. Recovery factors of 80%, at immiscible pressures 30-40% below CO2 MMP, were achieved. Addition of 15% N2 (v/v) to the CO2/CO mixture did not impair oil recovery. Interfacial testing (IFT) of oils, using pure CO, demonstrated a lowering of the IFT. RBA testing of asphaltine-rich heavy oils has shown that a mixture of CO2/CO dissolves into the oil at a far faster rate than either CO2 or CO individually and faster than the sum of both individual gases. A similar test using non-asphaltine type oils did not display this unique characteristic. Slim tube testing suggests that CO facilitates the mobilization of asphaltine-rich heavy oils and lowers viscosity. A linear corefloods of a reservoir containing 5% smectite + illite/smectite + and chlorite demonstrated a 275% increase in matrix permeability. Packed column tests, containing quartz sand and bentonite, demonstrated up to 300-900% increase in permeability in the presence of CO.
Thus a method to recover oil faster, from stranded reservoirs, at pressures below MMP, using significantly less gas, appears possible. In addition the use of CO, either alone or in combination with CO2 and/or N2, has been shown to increase matrix permeability. Such a gas mixture may be beneficial to achieving low pressure EOR from shallow, "stranded" reservoirs, non-conventional shale oil reservoirs, and viscous, heavy oil reservoirs at low temperatures. Incorporation of CO, or CO2/CO, in a frac fluid, or alternately as a post frac cleanup for shale oil and gas applications appears to warrant investigation.
Waterflood implementation accounts for more than half of the oil production worldwide. Despite the observations and extensive research from a large number of floods and thousands of simulation studies, managing waterfloods and Enhanced Oil Recovery (EOR) floods is still a technical challenge. A major contributor to this challenge are waterflood induced fractures (WIF). Managing waterfloods is a multivariable problem although WIF are one aspect, it is by no means the only controlling factor.
The best evidence that WIF are one of the main factors controlling flow in reservoirs is the insensitivity of injection pressure to injection rates. With our experience, in hundreds of waterfloods, we have frequently observed this phenomenon in the field data. If fluid flow depended on diffusive Darcy flow alone, we would expect higher injection rates with higher injection pressures. However, it is common to observed relatively constant injection pressures over a wide range of water injection rates. Rapid well communication and changes in water cuts that vary with injection rates also support an interpretation of high permeability induced fractures between injector and producer. In some reservoirs, interwell tracer data can be used to determine the influence of induced fracture features. The interwell tracers usually show very fast water movement.
Induced fractures in waterfloods and EOR projects can be caused by a number of mechanisms such as but not limited to, pressure depletion, changing pressure regimes, thermal effects, or plugging effects. These fractures can either be beneficial to the reservoir performance or effect performance negatively. Benefits include improved injectivity and increased throughput of the displacing fluid. Negative effects can come in the form of reduced volumetric sweep efficiency, impaired ultimate recovery or injected fluid losses out of zone.
Case studies, theory, and available literature from Western Canada will be reviewed in order to suggest and improve reservoir management strategies for waterfloods. We have completed hundreds of waterflood feasibility, waterflood management and EOR flood studies worldwide and continue to be amazed and humbled by the complexity that many waterfloods and EOR floods exhibit due to induced fracturing. WIF and EOR induced fractures (EIF) are common and should be analysed to optimize production. Growth of the WIF, response to waterflood with the presence of WIF, implication of WIF and reservoir management are the main areas which will be addressed.