Search

Summary Polymer flooding is an enhanced oil recovery (EOR) method which improves the mobility ratio and sweep efficiency of a waterflood. In theory, the high viscosity of the polymer reduces its injectivity compared to water. However, field studies have documented much higher polymer injectivity than predicted by theoretical models. There are various reasons for high polymer injectivity. The objective of this work is to predict polymer injectivity in granular media accounting for fluid-induced fractures, water quality, polymer rheology, and undissolved polymers. We perform grain-scale, coupled fluid dynamics and granular mechanics modeling. Fluid-particle interactions are modeled by coupling computational fluid dynamics (CFD) and the discrete element method (DEM). Simulation results show that polymer injection can create fractures in the granular media along the direction perpendicular to the minimum principal stress, thereby reducing wellbore pressure buildup at a constant polymer injection rate. The polymer tends to flow in the direction of fracture propagation in granular media, so the direction of the fracture affects the swept area of the polymer. Polymer rheology, water quality, and undissolved polymer also affect the polymer injectivity. Suspended solid particles may plug pores and reduce the injectivity of polymer by ~25%. Mechanically trapped undissolved polymers can greatly reduce polymer injectivity in low-permeability granular media. This work shows for the first time initiation of polymer-driven fractures in a granular model and demonstrates its implications on polymer injectivity.

Summary An efficient tertiary chemical flood involves the chemicals displacing the oil in a stable manner, so the initially disconnected trapped oil ganglia coalesce into a bank as the flood progresses. However, the location of surfactants and polymers in the core during a chemical flood is difficult or impossible to determine at very fine time and space resolutions. In this work, we bridge this gap by visualizing saturation profiles in a coreflood-on-a-chip, a micromodel that is the same length as a core, and monitoring fluorescent aqueous injectant as displacements occur. We visualize, for the first time, surfactants in and around oil banks at the pore and Darcy scales during low-tension displacements. We present 12 chemical floods in the coreflood-on-a-chip, imaged at the centimeter scale, including six surfactant and six alkali-surfactant-polymer (ASP) floods at low, intermediate, and high initial water saturations. Additionally, we present a micron-scale visualization of an ASP flood under UV light to observe the distribution of surfactant around an oil bank at the pore scale. We found that oil banks formed even during very unfavorable displacements, when surfactant solution (1 cp) displaced oil (80 cp) in a micromodel at intermediate and residual oil saturations. Although our results show saturation profiles are mostly well described by fractional flow theory, the distribution of injected aqueous phase is not. Significant aqueous injectant, including surfactant, flows within and ahead of oil banks.

We propose a simple microfluidic approach: Dissolution-After-Precipitation (DAP), to investigate regimes of carbonate rock dissolution and multiphase reactive transport. In this method, a carbonate porous medium is created in a glass microchannelivia/icalcium carbonate precipitation, after which an acid is injected into the channel to dissolve the precipitated porous medium. Utilizing the DAP method, for the first time we realized all five classical single-phase carbonate rock dissolution regimes (uniform, compact, conical, wormhole, ramified wormholes) in a microfluidic chip. The results are validated against the established theoretical dissolution diagram, which shows good agreement. Detailed analysis of these single-phase dissolutions suggests that the heterogeneity of the porous medium may significantly impact how the dissolution patterns evolve over time. Furthermore, DAP is utilized to investigate multiphase dissolution. As examples we tested the cases of an oleic phase (tetradecane) and a gaseous phase (COsub2/sub). Results show that the presence of a nonaqueous phase in pore spaces induces the formation of wormholes despite weak advection, and these wormholes ultimately become pathways for nonaqueous phase transport. However, the transport of tetradecane in the wormhole is very slow, causing acid breakthrough into neighboring regions. This mechanism enhances lateral connectivity between wormholes and may lead to a wormhole network. In contrast, COsub2/submoves rapidly and continuously seeks to enter a widening wormhole from a narrower wormhole or the porous regions, generating phenomena such as ganglia redistribution and counterflow (flow of gas opposite to acid flow). Extensive independent experiments are conducted to verify the reproducibility of the observed phenomena/mechanisms and further analyze them. Real-time monitoring of fluid pressure drop during dissolution is implemented to complement microscopy image analysis. Our method can be implemented repeatedly on the same chip, which offers a convenient and inexpensive option to study pore-scale reactive transport mechanisms.

SummaryNew experiments using polyethylene oxide (PEO) polymer were performed to evaluate its potential for enhanced oil recovery (EOR) applications in low-permeability reservoirs. This is the first time that high molecular weight PEO solutions have been shown to have favorable transport in low-permeability (~20 md) carbonate cores and the first time PEO has been shown to improve oil recovery in a fractured carbonate core. Rheology measurements in synthetic seawater show the higher viscosity of PEO solutions compares favorably to the viscosity of acrylamide–sodium acrylate (AM-AA) copolymers of similar molecular weight because PEO is less sensitive to hardness and high salinity. Filtration experiments using 0.45 μm cellulose filters show very favorable filtration ratios of PEO with a molecular weight of 4 million g/mol, which is consistent with its favorable transport in low-permeability cores. Four coreflood experiments in Texas Cream Limestone (TC Limestone) cores demonstrate the viability of PEO for EOR in low-permeability carbonate rocks. Single-phase experiments show 4 million g/mol PEO solutions transported through 18 and 28 md TC Limestone cores. Oil recovery experiments show 4 million g/mol PEO solutions transported and was more efficient than waterflooding in aged TC Limestone with favorable retention of 40 µg/g rock. An oil recovery experiment in an artificially fractured TC Limestone core improved oil recovery by a remarkable 15% considering the very large fracture-matrix permeability contrast (>7,000). These experimental results as well as other favorable properties of PEO reported in the literature indicate PEO should be considered for some EOR applications, especially in low-permeability reservoirs.

Summary Nanoparticles have great potential to mobilize trapped oil in reservoirs because of their chemical, thermal, and interfacial properties. However, the direct application of magnetic forces on superparamagnetic nanoparticles in reservoir engineering applications has not been extensively investigated. We demonstrate the enhanced oil recovery (EOR) potential of hydrophilic superparamagnetic nanoparticles in oil production by direct observation using microfluidics. We studied the mobilization of oil blobs by a ferrofluid (a suspension of hydrophilic superparamagnetic nanoparticles in water) both in a converging/diverging micromodel channel and in a foot-long pore network micromodel, both with varying depth (so-called 2.5D micromodels). The water-based ferrofluid in all cases was the wetting fluid. Initial ferrofluid flooding experiments in single channels were performed without and then with a static magnetic field. This magnetic field caused oil droplet deformation, dynamic breakup into smaller droplets, and subsequent residual oil saturation reduction. During the flooding, after the magnetic field was applied, significant oil displacement was observed within 2 hours [6 pore volumes injected (PVI)], and 86.2% of the oil that was not mobilized without a magnetic field was mobilized within 64 hours (192 PVI). Then, in experiments in the micromodel and in a Hele-Shaw cell without flooding, we observed self-assembly of oil droplets, indicating the formation of the hydrophilic magnetic nanoparticle microstructures (chains under the magnetic field) and their interaction with the oil blobs. Further ferrofluid flooding experiments were performed in a foot-long micromodel under a rotating magnetic field. The oil saturation was reduced from 44.6 to 33.3% after 17 hours (8.5 PVI) of ferrofluid flooding after the rotating magnetic field was applied. Finally, a discussion of field application of ferrofluid flooding is presented.

This study presents a workflow to optimize the location of CO2 injectors in order to maximize stored volume and prevent fault reactivation due to increases of pore pressure. We combine coupled reservoir flow and geomechanics simulations with neural network surrogate models and find the best injector location based on net present environmental value (NPEV). The surrogate models can complement the numerical reservoir simulations and efficiently optimize the well location and injection rate via the steepest ascent algorithm. We apply the workflow to the Brugge field, a closed synthetic reservoir that incorporates multiple faults, as an example. The surrogate model predicts a NPEV of approximately 0.19 billion carbon credits for the Brugge field, which is consistent with the corresponding reservoir simulation. NPEV will change under different safety factors that serve as an additional injection control to enhance the safety of injection. This proposed workflow could be readily applied to other similar reservoir models.

Geometric straining of particles in porous media is of critical importance in a broad range of natural and industrial settings, such as contaminant transport in aquifers and the permeability decline due to pore plugging in oil reservoirs. Pore network modeling is a computationally efficient approach to simulate transport problems in porous media that has been used to simulate particle deposition and size exclusion. However, existing network models are unable to simulate particle jamming due to the simplification of geometry and the lack of capability for simulating particle–particle interactions. Here, we develop a novel pore network model for particle jamming in porous media. The jamming in pore throats is predicted by the probability of jamming, which is a function of pore/particle size ratio, particle concentration, and coefficient of friction (COF). A unified probability model is developed based on Discrete Element Method (DEM) simulations of particle jamming in porous media consisting of a single layer of spherical grains. The geometry is based on two packing extremes, those with grains arranged in a triangle and a square. The model is then implemented into the pore network model to predict jamming in porous media. This framework combines the efficiency of network modeling and the accuracy of direct simulations on particle jamming. We verify the model against CFD-DEM simulations. The results of network simulations show that the probability of jamming increases with COF. At low particle concentrations (C = 5%), the probability of jamming for large pore throats is small. There is a critical particle concentration (C = 9%) when the grain/particle size ratio is 10, above which jamming is self-reinforcing and the inlet face of the porous medium will be completely blocked (C = 11%). The effect of velocity on the probability of jamming is insignificant compared with the effects of particle concentration and pore/particle size ratio.

Countercurrent spontaneous imbibition (SI) is an important flow mechanism for oil recovery in fractured reservoirs during waterflooding. SI plays a key role in the mobilization of oil in the matrix because it facilitates water infiltration by capillarity even when the matrix permeability is low, which limits fluid transport by advection. However, the modeling of SI in fractured media under dynamic conditions has been insufficiently studied. Most imbibition models assume conventional exponential functions and empirical constants based on experimental results of oil recovery under static conditions. Thus, the modeling of water distribution in the fracture and the matrix has been ignored which may lead to incorrect estimates of the efficiency of countercurrent SI to recover oil. Using the classic fractional flow equation, we present a semi-analytical solution to model countercurrent SI in a water-wet fractured medium by including a transfer function to account for continuous fluid exchange between the fracture and the matrix. The model is numerically solved using finite differences by including an effective imbibition time as a function of water advance in the fracture, which overcomes the difficulty of solving iterative numerical summations as shown in other imbibition models. The novelty of the presented solution is that it enables the modeling of water infiltration in the matrix under dynamic conditions. We verified the semi-analytical model against 2D numerical simulations and validated against experimental data to demonstrate that the model accurately predicts oil recovery and water infiltration in the matrix.

ABSTRACT: Fluid injection into sand may initiate fluid-driven fractures. The fluid-sand interaction patterns depend on the injection velocity, fluid viscosity, rheology and stress boundary conditions. However, the extension to non-Newtonian fluids and anisotropic stress conditions, such as in the case of injection of a polymer in the subsurface, is largely unexplored. In this paper, we study the injection of a Newtonian fluid (glycerol aqueous solution) and a shear-thinning fluid (scleroglucan) into a sand-filled Hele-Shaw cell. The experimental results demonstrate complex fracture patterns depending on injection velocity, fluid viscosity, boundary conditions and fluid rheology. The presence of stress-free boundaries causes the fracture opening to be almost perpendicular to the minimum principal stress. The peak injection pressure is also higher compared to the fixed-displacement boundaries. Compared to Newtonian fluid (glycerol), shear thinning polymer (scleroglucan) causes thinner fracture openings and more radial infiltration morphology in sand. The fracturing mechanism in sand is associated with the shear dilation at the tip of fractures and fracture offshoots. This study helps to clarify the basic flow and failure mechanisms of sand when subjected to injection of a non-Newtonian fluid. The insights obtained from these experiments can aid in the prediction of fracture morphology and sweep efficiency for on-site polymer injection. 1. INTRODUCTION Hydro-mechanical interactions between fluid flow and solids are common in many industries, such as oil and gas exploration, formation damage caused by particle jamming during oil and gas production (Li et al., 2021a), hydraulic fracturing (Zheng et al., 2020; Zheng and Sharma, 2021, Li et al., 2020), spray coating processes in food processing (Hilton et al., 2013) and tablet production in pharmaceuticals (Hemamanjushree and Tippavajhala, 2020). Different fluid infiltration and fracture opening morphologies caused by the injection of Newtonian fluids affect the fluid injectivity (Li et al., 2021b). However, extensions to non-Newtonian fluids and anisotropic stress conditions, such as in the case of injection of polymers in geological formations, are largely unexplored.

Many non-Newtonian fluids, including polymers, exhibit both shear-thinning and viscoelastic rheological properties. A lattice Boltzmann (LB) model is developed for simulation of the flow of thinning–elastic fluids through porous media. This model applies the Oldroyd-B constitutive equation and the Carreau model, respectively, to account for the viscoelastic and shear-thinning behaviors of the thinning–elastic fluid in porous media. Both rheological features are captured well by this model and are verified against analytical solutions. The thinning-then-thickening viscosity curve of the thinning–elastic fluid observed in experiments is reproduced by the present pore-scale simulations. In addition to the traditional extensional theory, we propose other important mechanisms for the increase in apparent viscosity of viscoelastic fluids at higher shear rates. The mechanisms proposed include the reduction in conductivity due to stagnant fluid, the compressed effective flow region, and larger energy dissipations caused by the viscoelastic instability. We find that the viscoelastic thickening effect is more prominent in porous geometries with a large pore–throat ratio.

Summary Surfactant flooding is an effective enhanced oil recovery method in which the oil/water interfacial tension (IFT) is reduced to ultralow values (10 mN/m). However, the microscopic displacement dynamics can be significantly different when the IFT is ultralow because the dynamic contact angle increases with the increase of the capillary number. In this study, surfactant flooding was performed and visualized in micromodels to investigate the dynamics of multiphase displacement at ultralow IFT. Although the micromodels used were strongly water-wet, the displacements of oil by surfactant solutions at ultralow IFT appeared as drainage. Furthermore, a macroscopic oil film was left behind on the surface, which indicates that a contact line instability occurred during displacements. The shape of the oil/water meniscus was determined by the balance between viscous forces and capillary forces. The meniscus can be significantly distorted by viscous forces at ultralow IFT. Therefore, the water-wet micromodel exhibits an oil-wet behavior at ultralow IFT, and the displacements of oil by surfactant solutions at ultralow IFT manifested as drainage rather than imbibition. The flow behavior is further complicated by the spontaneous formation of microemulsion during displacement. The microemulsion is mainly formed from the residual oil. The formation of a microemulsion bank made the surfactant solution discontinuous, with transport in the form of droplets in the microemulsion bank and displacement front. The novelty of this work is to reveal the effects of dynamic contact angle on the ultralow IFT displacement.

Summary The reduction of pore pressure caused by depletion can induce significant reservoir compaction and loss of permeability, especially in unconsolidated reservoirs. In this paper, we develop a numerical approach on the basis of computer-based simulations of unconsolidated rock samples subjected to mechanical tests that replicate the one-dimensional (1D) strain depletion path and allow for a prediction of permeability loss. The 1D strain stress path is a good approximation for long and thin conventional reservoirs with a compliant caprock. The numerical sample consists of crushable stiff and soft grains (proxies for sand and shale) simulated with the discrete element method (DEM) coupled with the bonded-particle model (BPM). Model parameters are calibrated through numerical single-grain-crushing tests which reproduce the experimentally measured sand strength. Grain crushing induced by the uniaxial strain stress path results in a pronounced reduction of porosity and permeability, which manifests more readily for samples with large grain size. The change of particle-size distribution indicates that high effective stresses cause grain crushing and production of a significant amount of fines that lower permeability. Simulation results indicate that the presence of soft grains and inclusions (e.g., shale fragments) facilitates grain crushing. Reservoir simulations—incorporating the change of porosity and permeability as a compaction table—show that the compaction can enhance cumulative production due to compaction drive but also reduces production rate by impairing the reservoir permeability. This multiscale numerical workflow bridges grain-scale compaction behavior and field-scale reservoir production.

Permanent trapping of an oscillating, nonwetting droplet is observed in a converging-diverging microchannel when aqueous, viscoelastic fluids are injected. Classical theories based on the balance between capillary and viscous forces suggest that the droplet size should decrease with increasing flow rates of a displacing Newtonian fluid, and the droplet should be completely displaced at high injection rates. However, droplets in viscoelastic fluids cannot be removed by increasing flow rates due to the oscillation. The oscillation amplitude linearly increases with the Deborah number (De), which further inhibits the droplet's passing through the constriction, "permanently." Our microfluidic experiments show that the onset of oscillation is determined by a critical De, which is near 1. We derive a linear relationship for the trapped droplet length with Ec^{1/3}, where Ec is the elastocapillary number, by introducing the elastic force into the force balance in addition to the capillary and viscous forces.

Harsh conditions, such as high temperature (>100 oC) and high salinity (>50,000 ppm TDS), can make the application of chemical enhanced oil recovery (EOR) challenging by causing many surfactants and polymers to degrade. Carbonate reservoirs also tend to have higher concentrations of divalent cations as well as positive surface charges that contribute to chemical degradation and surfactant adsorption. The objective of this work is to develop a surfactant-polymer (SP) formulation that can be injected with available hard brine, achieve ultra-low IFT in these harsh conditions, and yield low surfactant retention. Phase behavior experiments were performed to identify effective SP formulations. A combination of anionic and zwitterionic surfactants, cosolvents, brine, and oil was implemented in these tests. High molecular weight polymer was used in conjunction with the surfactant to provide a high viscosity and stable displacement during the chemical flood. Effective surfactant formulations were determined and five chemical floods were performed to test the oil recovery potential. The first two floods were performed using sandpacks from ground Indiana limestone while the other three floods used Indiana limestone cores. The sandpack experiments showed high oil recovery proving the effectiveness of the formulations, but the oil recovery was lower in the cores due to complex pore structure. The surfactant retention was high in the sandpacks, but it was lower in Indiana Limestone cores (0.29-0.39 mg/gm of rock). About 0.4 PV of surfactant slug was enough to achieve the oil recovery. A preflush of sodium polyacrylate improved the oil recovery.

Fluid injection into porous media can initiate fluid-driven fractures. Unlike cohesive and homogeneous media, fluid-driven fractures in cohesionless granular media have complex fracture and leak-off morphologies. In this study, we investigate the mechanisms behind fluid-driven fractures in cohesionless granular media by coupling the discrete element method with computational fluid dynamics at the pore/grain scale. We perform numerical simulations of the injection of a fluid into a granular medium that is saturated with a second, immiscible fluid. Results show that fluid injection can initiate a fracture orthogonal to the minimum principal stress, as expected, through viscous drag forces. However, the displacement field around the fracture is much more complex than what is expected for a linear elastic solid, showing that (1) fracture propagation requires small leak-off (a few times the fracture width) and (2) localized shear deformation ahead of the fracture tip affects fracture morphology. The simulations support a transition from a leakoff-dominated fluid flow regime to a fracture-dominated regime, with the increase in characteristic time t1 (ratio of permeability and injection rate per unit length) and decrease in t2 (ratio of injected fluid viscosity and minimum total principal stress). We develop a correlation of the dimensionless initiation pressure PD,peak and dimensionless stress SD (upper bound: $$P_{D,peak} = 5.647S_{D}^{0.86}$$ , lower bound: $$P_{D,peak} = 0.567S_{D}^{0.83}$$ ) based on our grain-scale simulations, experimental data and field data from the literature, which can be used to predict the fracture initiation pressure in cohesionless granular media.

SummaryWe present a study of the low–salinity effect during oil recovery using microfluidics experiments in an attempt to narrow the gap between pore–scale observations and porous–media–flow mechanisms, and to explain one type of low–salinity effect with delayed oil recovery and without the presence of clay. A microfluidic toolbox is used, including single–pore–scale microchannels, a pore–network–scale (approximately 102 pores) micromodel, and a reservoir–on–a–chip model (approximately 104 pores with heterogeneity), all with 2D connectivity. Experiments at the single–pore scale reveal a time–dependent oil dewetting and swelling behavior when a crude–oil droplet is in contact with low–salinity water. An interplay between water chemical potential and oil–phase polar compounds explains this pore–scale observation well. Experiments at the pore–network scale illustrate that the dewetting and swelling of residual oil in the swept region increase the water–flow resistance, modifying the flow field and thus redirecting the flooding liquid into unswept regions. This pore–network–scale effect is re–expressed into a macroscale model as a sweep–efficiency improvement derived from the change of relative permeabilities, which requires time to develop. Finally, experiments on our “reservoir–on–a–chip” model show significant incremental oil recovery during tertiary low–salinity waterflooding and confirm that late–time sweep–efficiency improvement contributes to most of the incremental oil recovery. On the basis of this microfluidic framework, we emphasize the following three findings: Low–salinity tertiary waterflooding can improve oil recovery by an improvement of sweep efficiency, which is a consequence of residual–oil dewetting and swelling.The low–salinity effect can occur without the existence of clay.The wettability alteration and oil swelling are time–dependent processes and should be expressed as a function of oil/water contact time rather than dimensionless time [pore volume (PV)], which explains some observations from previous coreflood experiments.

Wellbore integrity is a critical subject in oil and gas production and CO2 storage. Successful subsurface deposition of various fluids, such as CO2, depends on the integrity of the storage site. In a storage site, injection wells and pre-existing wells might leak due to over-pressurization, mechanical/chemical degradation, and/or a poor cement job, thus reducing the sealing capacity of the site. Wells that leak due to microannuli or cement fractures on the order of microns are difficult to seal with typical workover techniques. We tested a novel polymer gelant, originally developed for near borehole isolation, in a pilot experiment at Mont Terri, Switzerland to evaluate its performance in the aforementioned scenario. The polymer gel sealant was injected to seal a leaky wellbore drilled in the Opalinus Clay as a pilot test. The success of the pH-triggered polymer gel (sealant) in sealing cement fractures was previously demonstrated in laboratory coreflood experiments [1-3]. pH-sensitive microgels viscosify upon neutralization in contact with alkaline cement to become highly swollen gels with substantial yield stress that can block fluid flow. The leaky wellbore setup was prepared by heating-cooling cycles to induce leakage pathways in the cased and cemented wellbore. The leakage pathways are a combination of fractures in the cement and microannuli at the cement-formation interface. The exact nature of these leakage pathways can be determined by over-coring at the end of the experiment life. We used polyacrylic acid polymer (sealant) to seal these intervals. The process comprises of three stages: (1) injection of a chelating agent as the preflush to ensure a favorable environment for the polymer gel, (2) injection of polymer solution, and (3) shut-in for the polymer gelation. Then, we evaluated the short-/long-term performance of the sealant in withholding the injected fluids (formation brine and CO2 gas). The novel sealant was successfully deployed to seal the small aperture pathways of the borehole at the pilot test. We conducted performance tests using formation brine and CO2 gas to put differential pressure on the polymer gel seal. Pressure and flow rate at the specific interval were monitored during and after injection of brine and CO2. Results of performance tests after polymer injection were compared against those in the absence of the sealant. Several short-term (4 min) constant-pressure tests at different pressure levels were performed using formation brine, and no significant injection flow rate (rates were below 0.3 ml/min) was observed. The result shows more than a ten-fold drop in the injection rate compared to the case without the sealant. The polymer gel showed compressible behavior at the beginning of the short-term performance tests. Our long-term (1-week) test shows even less injectivity (~ 0.15 ml/min) after polymer gelation. The CO2 performance test shows only 3 bar pressure dissipation overnight after injection compared to an abrupt loss of CO2 pressure in the absence of polymer gel. Sealant shows good performance even in the presence of CO2 gas with high diffusivity and acidity. Pilot test of our novel sealant proves its competency to mitigate wellbore leakage through fractured cement or debonded microannuli, where other remedy techniques are seldom effective. The effectiveness of the sealing process was successfully tested in the high-alkaline wellbore environment of formation brine in contact with cement. The results to date are encouraging and will be further analyzed once over-coring of the wellbore containing the cemented annulus occurs. The results are useful to understand the complexities of the cement/wellbore interface and adjust the sealant/process to sustain the dynamic geochemical environment of the wellbore.

Preferential flow that leads to non-uniform displacement, especially in heterogeneous porous media, is usually unwelcome in most practical processes. We propose a self-adaptive preferential flow control mechanism by using dispersed polymers, which is supported strongly by experimental and numerical evidence. Our experiments are performed on a microchip with heterogeneous porous structures where oil is displaced by dispersed polymer microsphere particles. Even though the size of the particles is much smaller than the pore-throat size, the diversion effect by the dispersed microspheres is still proved. Therefore, the plugging effect is not the major mechanism for preferential flow control by dispersed polymers. The mechanisms are further investigated by pore-scale modelling, which indicates that the dispersed polymers exhibit an adaption ability to pressure and resistance in the porous flow field. In such an intelligent way, the displacing fluid with dispersed polymers smartly controls the preferential flow by inducing pressure fluctuations, and demonstrates better performance in both efficiency and economic aspects than the traditional method by simply increasing the viscosity. These insights can be applied to improve techniques in the field, such as enhanced oil recovery and soil wetting.

Reservoir simulators often predict poor injectivity in polymer enhanced oil recovery (EOR) projects because of the high polymer viscosity, which is a deterrent for the project. However, field studies have shown much higher injectivity than predicted by the models. The objective of this work is to perform grain-scale, coupled fluid dynamics and geomechanics modeling to predict the injectivity of viscous, non-Newtonian polymers in wellbores. Fluid-rock interactions are modeled by coupling computational fluid dynamics (CFD) and the discrete element method (DEM). Fluid flow is determined using an open-source CFD software that solves the volume-averaged Navier-Stokes equation using the finite volume method on Eulerian grids. Grain-scale geomechanics (DEM) is used to explicitly solve the particle trajectories in a Lagrangian reference system. The simulation results confirm the hypothesis of fracture initiation and sand failure near the injector. The results show that the polymer-driven fracture initiation is associated with sand shear failure, while the fracture geometry is the result of the localization of sand shear failure and fluidization of unconsolidated sand at the fracture tip. The injection of a viscous fluid can create fractures in the direction perpendicular to the applied minimum principal stress. The presence of fractures increases the injectivity. The peak injection pressure is more than 3 times greater than the applied minimum principal stress. The viscosity increase of polymers promotes the initiation of fractures and results in a greater fracture aperture. The injection of polymer can promote the initiation of fractures, and therefore, increase the polymer injectivity. This work, for the first time, uses a grain-scale model to predict polymer-driven fracture initiation and demonstrates the improved injectivity observed in the field.

We performed coreflood experiments to determine the effect of polymer viscoelasticity on the residual saturation of low viscosity oils by varying the rock type, wettability, polymer rheology, and flow conditions. Several authors have shown that polymers, if viscoelastic, can recover a significant amount of capillary-trapped oil, beyond inelastic floods, in water-wet sandstones. We demonstrate that the polymer viscoelastic effect works for several different rock types, both water-wet and oil-wet media, and with low viscosity oils which broadens the applicability found in previous works. In polymer field studies, the rheology can be tailored to maximize the effect based on our findings. Eleven corefloods were performed in Bentheimer, Boise, and Berea cores. Two of the cores were made oil-wet using dichlorodiphenylsilane (DCDPS) and chlorotrimethylsilane (CTMS). The cores were initially saturated with brine and then displaced by low viscosity (4-10 cp) crude oil until steady state was reached. Brine was then injected to displace the oil until steady state and residual oil saturation was reached. In most experiments, HPAM polymer was then injected at ∼1 ft/day; the polymer rheology was tailored via the salinity, concentration, molecular weight, and degree of hydrolysis. Once steady state was reached, a second polymer solution (with the same viscosity but higher salinity) was injected until a final oil saturation was reached. A reduction in oil saturation was observed in all corefloods in which the polymer was viscoelastic. On average, the oil saturation was reduced by 5.6% in the first polymer flood and 4.0% in the second, high salinity polymer flood. Higher recoveries were found for some experiments in which the first polymer was more elastic (higher dimensionless Deborah number). Final oil saturations as low as 5% were achieved using only polymer (and no surfactant). Results were not dependent on the rock type (e.g. Bentheimer versus Boise) but the effect was more pronounced in water-wet cores than those changed to oil-wet.

Simulations and theory reveal oscillations of droplets when driven by viscoelastic fluids because of their elastic memory effect. The viscoelastic oscillation helps release trapped droplets from hard-to-displace positions.

Surfactant-polymer enhanced oil recovery (EOR) in carbonates is challenging at high salinity/hardness and high temperature conditions due to poor aqueous stability at these conditions. This work finds that a surfactant blend of zwitterionic surfactant (lauryl betaine) and anionic surfactants (carboxylate & internal olefin sulfonate – IOS) that is aqueous stable at 4–6 wt% salinity and 100 °C and has ultralow water/oil interfacial tension (IFT) to allow near-miscible water/oil displacement to enhance oil recovery. The surfactant-polymer aqueous stability is greatly improved after mixing 2 wt% oil to create a homogeneous oil-in-water microemulsion. Flow experiments are conducted at 100 °C in Indiana limestone cores with vugs and microporosity. The zwitterionic-anionic surfactant blend reduces oil saturation to 10% and drains oil in pores ≥200 nm. In comparison, the anionic blend (carboxylate + IOS) reduces oil saturation to 23% and drains oil in pores ≥1 μm. The adsorption of both surfactant blends on carbonate is less than 0.4 mg/g. Particularly, the X-ray computed tomography (CT) images show a negative correlation between pixel-wise oil saturation and porosity during oil displacements by water and by surfactant/polymer. This correlation is likely due to capillary pressure and the mixed wettability of the Indiana limestone. These results provide new surfactant formulations for EOR at extreme conditions and insights for modeling EOR in rocks with porosity heterogeneity and mixed-wettability.

Depletion-induced formation compaction is a concern for production in hydrocarbon reservoirs. In this study, numerical mechanical tests and digital samples are proposed as an alternative to uniaxial strain experiments for investigating the effects of compaction. We propose a grain crushing model with two-sized subgrains based on the discrete element method and bonded-particle model (DEM-BPM) and perform numerical single grain crushing tests to calibrate the model parameters against the experimental sand grain characteristic strength. Numerical uniaxial strain tests are used to predict porosity and permeability reduction caused by reservoir compaction. Grain crushing incurs a significant reduction of porosity under the uniaxial strain stress path, which becomes more evident for digital samples consisting of larger grains. The stress-dependent changes in porosity and permeability obtained from the DEM simulations can be used as inputs for a reservoir-scale correlation to predict the productivity loss due to reservoir compaction.

CO2 geological storage in saline aquifers results in acidification of resident brine. Chemical reactions between acidified brine and rock minerals lead to dissolution and precipitation of minerals at various time scales. Mineral dissolution and precipitation are often neglected in assessing the mechanical integrity of target storage formations, yet, changes in rock strength and deformational behavior can impact trapping mechanisms. This study shows the impact of exposure to CO2-charged brine on shear strength and stiffness of various outcrop rocks evaluated through triaxial testing. The tested rocks were exposed to CO2-charged brine over geological time at a naturally occurring near-surface seepage along the Little Grand Wash Fault and Salt Wash Grabens, which include the Crystal Geyser site near the town of Green River, Utah. Prior work suggests that this site provides a near-surface structural analog for possible fault-controlled CO2 leakage over time scales that exceed expected injection time scales (10–100 years). Results show mechanical alteration in various aspects: (1) CO2-charged brine alteration at near-surface conditions results in mineral dissolution/precipitation and reduction of shear strength and brittleness of Entrada sandstone and Summerville siltstone samples, and (2) carbonate precipitation in fractured Mancos shale leads to matrix stiffening and fracture mineralization resulting in overall stiffer and likely tighter shale. Additional discrete element simulations coupled with a bonded-particle-model confirm the role of cement bond size alteration as one of the main controls for rock chemo-mechanical alteration in sandstones. The chemo-mechanical alteration path that mimics cement dissolution (under stressed subsurface conditions) results in vertical compaction and lateral stress relaxation. Overall, results show that rock exposure to CO2-charged brine can impart distinct petrophysical and geomechanical changes according to rock lithology and location with respect to major CO2 conduits. While mineral dissolution in the storage rock may result in undesired reservoir strains and changes of stresses, mineral precipitation downstream from a leakage path can help seal potentially induced fractures.

The potential leakage of hydrocarbon fluids or carbon dioxide from subsurface formations is a primary concern in wellbore integrity, oil and gas production, and CO2 storage. Leaky wells with fractured cement or debonded microannuli are common sources of subsurface fluid leakage. The hydrocarbon fluid or CO2 can migrate through such pathways to shallower formations and ultimately to surface. Cement fractures may have apertures on the order of microns, which are difficult to seal with typical workover techniques. A material that provides low viscosity during the injection but much higher viscosity after injection, with a minimum pressure gradient to yield flow at the target zone, is a potentially effective approach to seal the leakage pathways through cement fractures. pH-triggered polymers are such a material: aqueous solutions with low viscosity at low pH, containing pH-sensitive microgels which viscosify upon neutralization to become highly swollen gels with substantial yield stress that can block fluid flow. For the wellbore leakage application, the large alkalinity of wellbore cement provides the required neutralization. Our coreflood and rheological experiments show that pH-triggered polymer sealants such as polyacrylic acid polymer provide a robust seal if the process is properly designed; however, its long-term applicability depends on the dynamic geochemical environment of the wellbore. The process comprises three stages: (1) injection of a chelating agent as the preflush to ensure a favorable environment for the polymer gel, (2) injection of polymer solution, and (3) shut-in for the polymer gelation. A systematic study was done to understand the conditions under which the polymer gel remains stable and effectively seals the leakage pathways. A numerical model, based on polymer rheological properties and governing mechanisms observed in the laboratory experiments, was developed to simulate the reactive flow and transport of pH-triggered polymers in narrow fractures. Comparison with experiments shows a generally good agreement, despite the relative simplicity of the model. The numerical model was used to investigate further the underlying mechanisms of the process. The results can be used to design effectively the remediation process for a known fracture aperture size of the target zone.

Many engineering and natural materials, as well as earth systems, exhibit a combination of elastic-plastic and viscous behaviors, but precisely evaluating their rheology, and the micromechanics and chemical-mechanical feedbacks governing such rheology, can be a challenge. The cross-linked polymer Carbopol has long been used to explore fundamental rheological behaviors, most recently including fracturing during viscous creep in nature. Here, through rheometer experiments we establish that Carbopol 940 is, to a first order, a Herschel-Bulkley material. However, we further establish that the yield stress and viscosity are affected by chemically-sensitive micromechanical controls, namely pH and concentration of the polymer mixture. We explore these effects via the novel use of cryogenic scanning electron microscopy (SEM). Through the SEM imaging we show that there is a semi-quantitative relationship between pH, porosity at the >10-μm scale, and yield stress, a result of the ionic repulsion between polymer links at the molecular scale. We appeal to a model wherein the yield stress is a direct function of jamming expressed at the SEM scale, similar to those described in granular systems. As the pH of polymer dispersion increases, and the porosity decreases, the yield strength increases as a result of the increasingly jammed system. The initial viscosity is thus controlled by the yield stress, but after failure evolves with increasing shear rate due to characteristic unjammed flow of the material. The different controls on the yield stress versus viscous flow rates has implications for borehole engineering efforts (carbon capture and storage) employing Carbopol, and could prove instructive for modeling of natural viscoplastic deformation.

Countercurrent spontaneous imbibition (SI) is an important flow mechanism to recover oil in water-wet porous media through capillarity. Experimental imbibition tests conducted in core samples show that oil recovery as a function of time occurs in a characteristic S-shaped curve, which describes the infinite-acting and the boundary-dominated regime. Although several dimensionless time groups have been proposed to model countercurrent SI in porous media, they fail to properly scale the results, causing inaccurate estimates of oil recovery. In this work, we present a new approach, using a hybrid solution, to model countercurrent SI for oil-water systems dominated by capillary forces. The infinite-acting regime is modeled by an early-time solution using a new dimensionless time group, which enables the correct scaling of imbibition results onto a universal single curve. To model the boundary-dominated regime, we introduce a late-time solution derived as a function of a characteristic distribution of saturation to estimate the flow behavior after the imbibing fluid reaches the no-flow boundary. The novelty of the model is that it enables the accurate estimation of fluid imbibition under the boundary-dominated regime, a flow condition critical to evaluate the true potential of countercurrent SI driven by capillarity for oil recovery in porous media. We verified the model against numerical simulations under a wide range of flow conditions relevant in water-oil systems, and used experimental data reported in the literature to validate our model. The solution presented provides an accurate approach to model countercurrent SI, which could be extended for dynamic conditions and the modeling of the flow of fluids in fractured media.

Spontaneous liquid imbibition is a dominant mechanism for moving fluids in confinements with extremely high hydrodynamic resistance; i.e. nanopores. We demonstrate the impact of hydrophilicity and liquid polarity on nanoscale imbibition with dynamic measurements of the uptake of water-isopropanol (polar) and heptane (nonpolar) within 2D glass-silica nanochannels exposed to varied drying conditions and rehydration (rehydroxylation) cycles. The Lucas-Washburn equation, which does not consider interfacial fluidity effects, predicts that water-IPA and heptane should imbibe at similar speeds. However, we observed stymied and hysteretic water-IPA imbibition trends explained by extremely large contact line friction and increased effective viscosity, both surface chemistry-dependent, whereupon the results match a modified version of the Lucas-Washburn equation that accounts for dynamic wetting. In contrast, heptane imbibition, though still slower than the Lucas-Washburn equation prediction, was fairly insensitive to drying history and an order of magnitude faster than the polar mixture. The imbibition of aqueous solutions in analogous siliceous nanoporous materials and structures, ubiquitous in earth science and nanotechnology, may also be subject to hysteretic and large energy dissipation at contact lines and interfaces on account of hydrophilicity.

The injection of CO2 into geological formations leads to geochemical re-equilibrium between the pore fluid and rock minerals. Mineral–brine–CO2 reactions can induce alteration of mechanical properties and affect the structural integrity of the storage formation. The location of alterable mineral phases within the rock skeleton is important to assess the potential effects of mineral dissolution on bulk geomechanical properties. Hence, although often disregarded, the understanding of particle-scale mechanisms responsible for alterations is necessary to predict the extent of geomechanical alteration as a function of dissolved mineral amounts. This study investigates the CO2-related rock chemo-mechanical alteration through numerical modeling and matching of naturally altered rocks probed with micro-scratch tests. We use a model that couples the discrete element method (DEM) and the bonded particle model (BPM) to perform simulations of micro-scratch tests on synthetic rocks that mimic Entrada sandstone. Experimental results serve to calibrate numerical scratch tests with DEM–BPM parameters. Sensitivity analyses indicate that the cement size and bond shear strength are the most sensitive microscopic parameters that govern the CO2-induced alteration in Entrada sandstone. Reductions in cement size lead to decrease in scratch toughness and an increase in ductility in the rock samples. This work demonstrates how small variations of microscopic bond properties in cemented sandstone can lead to significant changes in macroscopic large-strain mechanical properties.

The infiltration of solid particles (fines) into a reservoir during drilling or water reinjection is a long-standing problem in the petroleum industry because of the associated production decline or injectivity loss. Solid particle infiltration is also an active area of research in other fields such as waste water treatment and fines migration in clay-rich reservoirs. In this work, we perform experiments to study particle infiltration into sintered glass bead core plugs and measure the changes in porosity and permeability using computed tomography (CT) scanning and pressure transducers, respectively. Suspensions of prescribed bimodal particle sizes, concentrations, and flowrates are flooded through the core plugs. Permeability changes are measured continuously over time, and the porosity change is measured after flowing a fixed mass of invaded particles. A dual pore-scale numerical model approach (a combination of a direct pore-scale discrete element method (DEM) and a pore-scale network model) is used to predict permeability reduction by particle filtration in porous media. Large particle deposition is modeled using the DEM in a disordered sphere pack geometry; the result is a prediction of deposition as a function of large particle concentration. Small particle deposition is modeled using the network model; the result is a prediction of deposition as a function of small particle concentration. Crucially, the geometry of the network model depends on the amount of large particles deposited. In this way, the deposition of large and small sized particles is coupled together. The permeability and porosity reduction of the network due to the deposition of small and large particles are then calculated from the functions above. We compare the experimental results with simulation predictions and find that the dual pore-scale model is capable of predicting the permeability of the invaded core in the regions away from the injection face. Permeability prediction in the region adjacent to the injection face is improved by incorporating the influence of the external filter cake.

Shale oil production from plays such as the Bakken and Eagle Ford, driven by advances in horizontal well drilling and hydraulic fracturing technology, has helped improve US domestic oil production. However, due to the low permeability of shale oil reservoirs, primary depletion only produces 5–10% of the original oil in place and there is a huge potential for improved recovery. Chemical and gas huff-n-puff have been piloted with varying degrees of success, but no breakthrough recovery method has been discovered. This study examines improved oil recovery from shale reservoirs by thermal stimulation, consisting of primary depletion (early in a well's life), followed by conversion of the well to a heat injector to elevate reservoir temperature, and finally followed by a secondary depletion. The proposed method was tested using a compositional, thermal reservoir simulator (CMG STARS) and key parameters affecting recovery during thermal stimulation were investigated. This study found that 1000 days of thermal stimulation with a 700 °F heater has the potential to economically increase oil recovery from about 7% to more than 11.5%, with potential for even greater recovery if heat injection time and temperature are optimized. The study also found that thermal pressurization of oil is the primary mechanism for the improved recovery. Kerogen decomposition into oil and gas results in a significant increase of hydrocarbons in place but is only a minor contribution to the additional recovery because production is limited by a lack of flow capacity. Furthermore, a two-fold increase in permeability is observed as a result of kerogen decomposition, but its contribution to recovery is also minor because the heated region is not well connected with the fractured region. The heating scheme may be improved in the future to better connect the heated region with the fractured region.

The ultra-low permeability of shale makes injection from a well to other wells difficult. A novel scheme was proposed in a recent patent (Dombrowski et al., 2015) where gas is injected into a hydraulic fracture along a horizontal well and production occurs from an adjacent fracture, intersecting the same well. Compositional reservoir modeling was performed to investigate the effectiveness of the proposed gas injection scheme. The computational domain consists of two hydrofrac half-stages along a horizontal well. The results show 15.7% and 12.5% OOIP incremental recovery over 5000 days of CO 2 injection for the base models with the matrix permeability of 10 μD and 1 μD, respectively, demonstrating that the gas injection scheme has the potential to vastly improve oil recovery in oil-rich shale formations. The effects of reservoir properties and injection conditions on oil recovery were investigated by changing the injection pressure, reservoir heterogeneity, distribution of natural fractures, hydrofrac spacing, size of pore space, and compositions of the injection gas. Most of them affect the oil recovery significantly. Recovery by miscible hydrocarbon gas injection is comparable to CO 2 ; so it should be considered as an alternative.

We investigate the mechanisms of opening-mode fracture initiation in granular media. The study is based on a simulation of grain-scale fluid-grain interactions through a coupled numerical approach in which the discrete element method is used to solve for the mechanics of a solid granular medium, and computational fluid dynamics is used to model fluid flow and drag forces. We present benchmark problems with analytical solutions and validate this numerical model against experiments on a viscous-drag-driven cavity in the literature. Additional simulation results show fracture initiation mechanisms in a random granular packing subjected to constant boundary stresses and to fluid injection with a localized source. The dimensionless variable F_{s}/F_{sk} (ratio of seepage force F_{s} and skeletal force F_{sk}) incorporates the impacts of physical properties and injection parameters including fluid viscosity, injection velocity, grain size, and effective stresses, and it has been used as a criterion separating regimes of fluid invasion and drag-driven fracture opening. Our simulation results show that F_{s}/F_{sk} in combination with τ_{1} (ratio of diffusion time from hydromechanical coupling and injection time) serves as a prediction of fracture opening within granular packing. We suggest a simple criterion (F_{s}/F_{sk}>1 or τ_{1}>0.17) that is valid for various types of granular media and injection conditions to determine if fracture opening will occur. Among other applications, this study is useful to predict the initiation and propagation of fractures in natural sediments.

Long-term integrity and practical storage of CO2 is contingent upon its seal performance and the dynamic sealing capacity of faults for the CO2 storage site. Faults are prone to reactivation with reservoir pressurization caused by CO2 injection. The goal of this study is to create and verify a reservoir elasto-plastic model capable of capturing short-term evolution of fault reactivation and the resulting change of permeability. This model is then used to explore the effects of coupling geomechanics with reservoir fluid flow on the reactivation of faults. In this paper, we introduce a workflow for modeling of fault reactivation with fault elements as gridblocks instead of surfaces. Reservoir simulation, with coupled fluid flow and geomechanics, was used for this purpose. The simulation models utilize a geomechanical module to capture elasto-plasticity and a compositional numerical scheme based on an equation of state (EOS) to calculate CO2-brine interaction. The geomechanical module used in this study is based on Hierarchical Single Surface (HISS) model that captures strain softening and hardening, and therefore post-yield plastic deformations related to fault reactivation. The compositional numerical scheme based on EOS calculates the amount of CO2 solubilization in brine as well as the density and viscosity of the CO2- and aqueous-rich phase. In this approach, the flow properties, i.e. permeability and porosity, dynamically change in response to geomechanical effects. The dynamic change was captured through a volumetric strain-permeability law. Our simulation results show that the model is capable of capturing short-term evolution of fault reactivation and the resulting change of permeability along the fault. The dynamic changes of fault properties control the extent of fault reactivation, the pressure relief during injection, and the fault sealing capacity.

Wellbore integrity is a critical subject in oil and gas production, and CO2 storage. Successful subsurface deposition of various fluids, such as CO2, depends on the integrity of the storage site. In a storage site, injection wells and pre-existing wells might leak due to over-pressurization, mechanical/chemical degradation, and/or a poor cement job, thus reducing the sealing capacity of the site. Wells that leak due to microannuli or cement fractures on the order of microns are difficult to seal with typical workover techniques. We tested a novel polymer gelant, originally developed for near borehole isolation, in a pilot experiment at Mont Terri, Switzerland to evaluate its performance in the aforementioned scenario.The polymer gel sealant was injected to seal a leaky wellbore drilled in the Opalinus Clay as a pilot test. The success of the pH-triggered polymer gel (sealant) in sealing cement fractures was previously demonstrated in laboratory coreflood experiments (Ho et al. 2016, Tavassoli et al. 2018). pH-sensitive microgels viscosify upon neutralization in contact with alkaline cement to become highly swollen gels with substantial yield stress that can block fluid flow. The leaky wellbore setup was prepared by heating-cooling cycles to induce leakage pathways in the cased and cemented wellbore. The leakage pathways are a combination of fractures in the cement and microannuli at the cement-formation interface. The exact nature of these leakage pathways can be determined by over-coring at the end of the experiment life. We used polyacrylic acid polymer (sealant) to seal these intervals. The process comprises of three stages: (1) injection of a chelating agent as the preflush to ensure a favorable environment for the polymer gel, (2) injection of polymer solution, and (3) shut-in for the polymer gelation. Then, we evaluated the short-/long-term performance of the sealant in withholding the injected fluids (formation brine and CO2 gas).The novel sealant was successfully deployed to seal the small aperture pathways of the borehole at the pilot test. We conducted performance tests using formation brine and CO2 gas to put differential pressure on the polymer gel seal. Pressure and flow rate at the specific interval were monitored during and after injection of brine and CO2. Results of performance tests after polymer injection were compared against those in the absence of the sealant.Several short-term (4 min) constant-pressure tests at different pressure levels were performed using formation brine, and no significant injection flow rate (rates were below 0.3 ml/min) was observed. The result shows more than a ten-fold drop in the injection rate compared to the case without the sealant. The polymer gel showed compressible behavior at the beginning of the short-term performance tests. Our long-term (1-week) test shows even less injectivity (~0.15 ml/min) after polymer gelation. The CO2 performance test shows only 3 bar pressure dissipation overnight after injection compared to abrupt loss of CO2 pressure in the absence of polymer gel. Sealant shows good performance even in the presence of CO2 gas with high diffusivity and acidity.Pilot test of our novel sealant proves its competency to mitigate wellbore leakage through fractured cement or debonded microannuli, where other remedy techniques are seldom effective. The effectiveness of the sealing process was successfully tested in the high-alkaline wellbore environment of formation brine in contact with cement. The results to date are encouraging and will be further analyzed once over-coring of the wellbore containing the cemented annulus occurs. The results are useful to understand the complexities of cement/wellbore interface and adjust the sealant/process to sustain the dynamic geochemical environment of the wellbore.