Summary The formation/precipitation and deposition of pipeline solids, such as gas hydrates, asphaltenes, and waxes have long plagued production fields. Given the vast differences in chemistries of these solids, any current prevention or mitigation strategy, particularly for cases in which multiple issues are a concern, is likely to involve an extensive assortment of chemical additives that are both costly and add complexity to the system. Surface treatments (coatings), on the other hand, present a relatively viable option for management strategies. A chemically and physically robust surface treatment with the ability to address deposition issues for multiple pipeline solids could not only decrease the operating expenditures for a field through material cost savings and obviation of downstream separation, but could also simplify produced fluids by eliminating additional chemicals from the mixture. The purpose of this study was to explore the feasibility of a particular surface treatment as part of a solids management strategy. This work used an omniphobic surface treatment to probe its effects on gas hydrate, asphaltene, and wax deposition. Specifically, an interfacial tensiometer (IFT) collected contact angle measurements for wettability studies. High-pressure rocking cells studied gas hydrate nucleation and deposition. A bench-scale flow loop quantified the deposition of oil and asphaltenes after a set time period. Finally, a mechanical shear device measured the adhesion force of wax deposits on untreated/treated surfaces. Static contact angle results showed that the omniphobic surface treatment had reduced surface interaction with water droplets in air, altering from the complete wetting on corroded surfaces to slightly hydrophobic conditions of greater than 100°. In addition, rocking-cell tests indicated that these omniphobic surface treatments may prevent gas hydrate deposition for up to 72 hours. The gas hydrate rocking-cell tests also demonstrated possible increases in induction time and occasional elimination of hydrate nucleation with the surface treatment. Finally, the surface treatment application, which also decreased surface roughness, showed that crude oil and asphaltene particles deposition, as well as the shear required to remove solidified wax deposits, could be reduced by a large factor. Overall, promising results were recorded for all major flow-assurance solids in the presence of the surface treatment.

Summary Nonuniform mixing during low-salinity waterflooding (LSWF) is a function of the pore geometry and flow patterns within the porous system. Salinity-dependent wettability alteration (WA) changes the entry capillary pressure, which may mobilize the trapped oil depending on the flow regime and salt dispersion pattern. The complex interplay between the wettability, capillary number ( N Ca ), and salt dispersion caused by pore-scale heterogeneity on the efficiency of LSWF is not well understood. In this paper, direct numerical simulations in a pore-doublet model (PDM) were carried out with OpenFOAM® (OpenCFD, Berkshire, UK) using the volume-of-fluid (VOF) method. Oil trapping and remobilization were studied at relevant N Ca as low as 10 − 6 under different initial wettability states. Depending on the effective salinity ranges (ESRs) for the low-salinity effect (LSE), three WA models were implemented, and the effects of WA degree and salinity distribution on LSWF flow dynamics were investigated. The slow process of WA by means of thin film phenomena was captured by considering a diffuse interface at the three-phase contact line. Because of the pore structure of the pore doublet, only in nonwater-wet cases, oil is trapped in the narrower side channel (NSC) after high-salinity waterflooding (HSWF) and may be remobilized by LSWF. In strongly oil-wet cases, oil is recovered gradually by LSWF by means of a film-flow mechanism near the outlet. In moderately oil-wet cases, however, the entire trapped oil ganglion can be mobilized, provided that the entry capillary pressure is sufficiently reduced. The degree of WA, ESR, kinetics of WA, and the wettability of pore surface at the outlet are determining factors in the drainage of the trapped oil. The salt dispersion pattern in the flowing region [i.e., wider side channel (WSC)] controls the wettability distribution and the rate and magnitude of oil recovery from the stagnant region (i.e., NSC). The difference between the WA models is more apparent near the outlet, where the salinity profile is more dispersed. The ESR in which WA occurs determines the speed of the entry capillary pressure reduction and, thus, the recovery factor. In cases where WA occurs at a salinity threshold (ST), the highest recovery is obtained, whereas with the full-salinity-range WA model, the oil recovery performance is lowest. From the capillary desaturation perspective, it is found that the LSE becomes more pronounced when N Ca is less than 10 − 5 , and the dispersion regime is in the power-law interval. Because the adverse effect of salt dispersion in the flowing region is delayed, the LSE is intensified. For the simulations to be representative of the actual conditions in the porous medium, much lower N Ca than currently used in many research works must be studied. Otherwise, the simulations may lead to over- or underestimation of the LSE. The synergetic or antagonistic effects caused by the interplay between viscous and capillary forces and dispersion may lead to total recovery or entrapment of oil, regardless of WA. Based on the pore geometry, initial wettability state, and balance of forces, the mobilized oil may flow past the conjunction (favorable) or in the backward direction (unfavorable) to the WSC and get retrapped. Successful drainage of oil from the pore system after WA is essential for observing incremental oil recovery by LSWF.

Summary Controlled-ionic-composition waterflooding is an economic and effective method to improve oil recovery in carbonate oil reservoirs. Recent studies show controlling the salinity and ionic composition of injection water can alter the wettability of carbonate mineral surfaces. The pore-scale oil connectivity and displacement by controlled-ionic-composition waterflooding in heterogeneous carbonate reservoirs, especially at the early stage, is still unclear. The goal of this study is to examine the role of ion concentrations and types in the oil displacement efficiency and investigate the impact of the waterflooding on the pore-scale oil displacement using the national synchrotron facility. A carbonate rock sample was flooded with synthetic high-salinity water and other water solutions with different sulfate concentrations. The waterflooding processes were visualized with synchrotron X-ray microtomography to follow the evolution of pore-scale oil/brine interactions at typical field flow rates. Experimental results show that the water with lower sulfate concentration and higher salinity did not change the wettability of the pore surfaces. Higher sulfate ion concentrations in the water, in contrast, altered the wettability of carbonate pore surfaces from oil-wet to neutral-wet within the first few minutes of waterflooding. Novel insight was gained on the ability of water with high-sulfate concentration to displace oil in the small pores and through abundant oil channels, which could consequently lead to higher oil recovery from the carbonate rock.

Summary In this study, to investigate how gravity affects foam in open vertical fractures, we report foam experiments in three 1-m-long, 15-cm-wide glass-model fractures. Each fracture has a smooth wall and a roughened wall. Between the two walls is a slit-like channel representing a single geological fracture. Three model fractures (Models A, B, and C) share the same roughness and have different hydraulic apertures of 78, 98, and 128 µm, respectively. We conduct foam experiments by horizontal injection in the three model fractures placed horizontally and sideways (i.e., with the model fractures turned on their long side), and in Model A placed vertically with injection upward or downward. Direct imaging of the foam inside the model fracture is facilitated using a high-speed camera. We find that foam reaches local equilibrium (LE; where the rate of bubble generation equals that of bubble destruction) in horizontal-flow experiments in all three model fractures and in vertical-flow experiments in Model A. In fractures with a larger hydraulic aperture, foam is coarser because of less in-situ foam generation. In the vertical-flow experiments in Model A, we find that the properties of the foam are different in upward and downward flow. Compared with downward flooding, upward flooding creates a finer-texture foam, as sections near the inlet of this experiment are in a wetter state, which benefits in-situ foam generation. Moreover, less gas is trapped during upward flooding, as gravitational potential helps overcome the capillarity and moves bubbles upward. In the sideways-flow experiments, gravity segregation takes place. As a result, drier foam propagates along the top of the fractures and wetter foam along the bottom. The segregation is more significant in fractures with a larger hydraulic aperture. At foam quality 0.8, gas saturation is 27.7% greater at the top than the bottom for Model C, and 19.3% and 10.8% for Models B and A, respectively. Despite the gravity segregation in all three model fractures, water and gas are not completely segregated. All three model fractures thus represent a capillary transition zone, with greater segregation with increasing aperture. Our results suggest that the propagation of foam in vertical natural fractures meters tall and tens of meters long, with an aperture of hundreds of microns or greater, is problematic. Gravity segregation in foam would weaken its capacity in the field to maintain uniform flow and divert gas in a tall fracture over large distances.

Summary The expansion of the steam chamber is very important for the recovery performance of steamflooding. In this paper, we discuss 1D and 2D sandpack experiments to performed analyze the effect of flue gas on steam chamber expansion and displacement efficiency in steamflooding. In addition, we examine the effect of flue gas acting on the steam condensation characteristics. The results show that within a certain range of injection rate, flue gas can significantly enlarge the swept volume and oil displacement efficiency of steam. However, when the flue gas injection rate is excessively high (the ratio of gas injection rate to steam injection rate exceeds 4), gas channels may form, resulting in a decline of oil recovery from steamflooding. The results of the 2D visualization experiments reveal that the swept volume of the steam chamber during steamflooding was small, and the remaining oil saturation in the reservoir was high, so the recovery was only 28%. The swept volume of the steam chamber for flue-gas-assisted steamflooding was obviously larger than that of steamflooding, and the recovery of flue-gas-assisted steamflooding in 2D experiments could reach 40.35%. The results of the steam condensation experiment indicate that flue gas could reduce the growth and coalescence rates of steam-condensed droplets on the cooling wall and increase the shedding period of the droplets. Macroscopically, flue gas could reduce the heat exchange rate between the steam and the reservoir and inhibit the rapid condensation and heat exchange of the steam near the injection well. As a result, flue gas could expand the steam chamber into the reservoir for heating and displacing oil.

Summary In this paper, we propose a fully implicit space-time multiscale scheme to improve computational efficiency in solving nonlinear multiphase flow in porous media. Here, error estimators are used for adaptively changing the spatial-temporal mesh. This algorithm applied to the black-oil model is compared to a standard control volume approach using a fine time and spatial mesh. Error estimators are introduced to determine subdomains of the reservoir in which high nonlinearity hinders Newtonian convergence. This is followed by applying local fine timesteps to these marked regions, whereas the remaining regions retain the coarse time scale. Once a temporal discretization is determined for different parts of the reservoir, the spatial mesh is refined for treating saturation fronts. The nonmatching interfaces arising from different temporal and spatial scales are resolved by the enhanced velocity method, which enforces strongly the continuity of fluxes. This whole system is solved monolithically. Results from a three-phaseblack-oil model are described. The multiscale solution is compared to a uniformly fine spatial mesh and fine timestepping solution to confirm accuracy. Solutions from both Gaussian and channelized permeability fields are presented. In the multiscale solution, we observe temporal refinements being applied to the water and gas saturation fronts by reducing the timestep size to guarantee Newtonian convergence. We also observe that the extent of refinements is balanced between the two saturation fronts based on their respective nonlinearity. For Gaussian permeability fields, the algorithm only treats saturation fronts spatially, whereas for channelized permeability fields, the geological features are also considered. Production profiles for the three phases match well between the two solutions under specific refinement criteria. We also investigate the improvement in computational efficiency, as well as the algorithm scalability in regard to increasing problem sizes. We observe a speedup of approximately 10 for solving the linear system, and such speedup increases as the problem size expands. In reservoir simulation, schemes that attempt to decouple the diffusion and advection process become less robust as the physics becomes more complex. Our approach, which is focused on handling the high nonlinearities, can be treated fully implicitly and provide a reduction in simulation costs. This approach can also be applied to model order reduction in providing accurate snapshot solutions. NOTE: This paper is published as part of the 2021 Reservoir Simulation Conference Special Issue.

Summary This work studies 1D steady-state flow of gas from compressible shale matrix subject to water blocking toward a neighboring fracture. Water blocking is a capillary end effect causing wetting phase (e.g., water) to accumulate near the transition from matrix to fracture. Hydraulic fracturing is essential for economical shale gas production. Water is frequently used as fracturing fluid, but its accumulation in the matrix can reduce gas mobility and production rate. Gas transport is considered at a defined pressure drop. The model accounts for apparent permeability (slip), compressibility of gas and shale, permeability reduction, saturation tortuosity (reduced relative permeability upon compaction), and multiphase flow parameters like relative permeability and capillary pressure, which depend on wettability. The behavior of gas flow rate and distributions of gas saturation, pressure, and permeability subject to different conditions and the stated mechanisms is explored. Water blockage reduces gas relative permeability over a large zone and reduces the gas flow rate. Despite gas flowing, strong capillary forces sustain mobile water over the entire system. Reducing drawdown gave lower driving force and higher resistance (by water blockage) for gas flow. The results show that 75% reduction of drawdown made the gas flow rate a couple orders of magnitude lower compared to if there was no blockage. The impact was most severe in more water-wetsystems. The blockage caused most of the pressure drop to occur near the outlet. High pressure in the rest of the system reduced effects from gas decompression, matrix compression, and slip-enhanced permeability, whereas rapid gradients in all these effects occurred near the outlet. Gas decompression resulted in an approximately 10 times higher Darcy velocity and pressure gradient near the outlet compared to inlet, which contributed to removing blockage, but the added resistance reduced the gas production rate. Similarly, higher gas Corey exponent associated gas flow with higher pressure drop. The result was less blockage but lower gas production. Slip increased permeability, especially toward the outlet, and contributed to increase in gas production by 16%. Significant matrix compression was associated with permeability reduction and increased Corey exponent in some examples. These effects reduced production and shifted more of the pressure drop toward the outlet. Upstream pressure was more uniform, and less compression and permeability reduction were seen overall compared to a system without water blockage.

Summary Seawater breakthrough percentage monitoring is critical for offshore oil reservoirs because seawater fraction is an important parameter for estimating the severity of many flow assurance issues caused by seawater injection and further developing effective strategies to mitigate the impact of those issues on production. The validation of using natural ions as a tracer to calculate the seawater fraction was investigated systematically by studying the natural chemical composition evolution in porous media using coreflood tests and static bottle tests. The applicable range of ions was discussed based on the interaction between ion and rock. The barium sulfate reactive model was improved by integrating interaction between ions and rock as well as fluid flow effect. The results indicate that chloride and sodium interact with rock, but the influence of the interaction can be minimized to a negligible level because of the high concentrations of chloride and sodium. Thus, chloride and sodium can be used as conservative tracers during the seawater flooding process. However, adsorption/desorption may have a large influence on chloride and sodium concentrations under the scenario that both injection water and formation water have low chloride and sodium content. Bromide shows negligible interaction with rock even at low concentrations and can be regarded as being conservative. The application of a barium and sulfate reaction model in coreflood tests does not work as well as in bottle tests because fluid flow in porous media and ion interaction with rock is not taken into account. Although sulfate and barium adsorption on clay is small, it should not be neglected. The barium sulfate reaction model was improved based on the simulation of ion transport in porous media. Cations (magnesium, calcium, and potassium) are involved in the complicated cation-exchange process, which causes large deviation. Therefore, magnesium, calcium, and potassium are not recommended to calculate seawater fraction. Boron, which exists as anions in formation water and is used as a conservative tracer, has significant interactions with core matrix, and using boron in an ion tracking method directly can significantly underestimate the seawater fraction. The results give guidelines on selecting suitable ions as tracers to determine seawater breakthrough percentages under different production scenarios.

Summary Oil/water interfacial tension (IFT) is an important parameter in petroleum engineering, especially for enhanced-oil-recovery (EOR) techniques. Surfactant and low-salinity EOR target IFT reduction to improve oil recovery. IFT values can be determined by empirical correlation, but widely used thermodynamic-based correlations do not account for the surface-activities characteristic of the polar/nonpolar interactions caused by naturally existing components in the crude oil. In addition, most crude oils included in these correlations come from conventional reservoirs, which are often dissimilar to the low-asphaltene crude oils produced from shale reservoirs. This study presents a novel oil-composition-based IFT correlation that can be applied to shale-crude-oil samples. The correlation is dependent on the saturates/aromatics/resins/asphaltenes (SARA) analysis of the oil samples. We show that the crude oil produced from most unconventional reservoirs contains little to no asphaltic material. In addition, a more thorough investigation of the effect of oil components, salinity, temperature, and their interactions on the oil/water IFT is provided and explained using the mutual polarity/solubility concept. Fifteen crude-oil samples from prominent US shale plays (i.e., Eagle Ford, Middle Bakken, and Wolfcamp) are included in this study. IFT was measured in systems with salinity from 0 to 24% and temperatures up to 195°F.

Summary In this paper, we present a solvent-assistedsteam-assisted gravity drainage (SA-SAGD) experiment with multicomponent solvent (i.e., condensate) using a large physical model. The sandpack for the experiment had a porosity of 0.33 and a permeability of 5.6 darcies in the cylindrical pressure vessel that was 1.22 m in length and 0.425 m in internal diameter. The sandpack was initially saturated with 93% Athabasca bitumen and 7% deionized water. The main objective of this research was to study the in-situ thermal/compositional flow and produced bitumen properties in SA-SAGD with condensate. After the preheating of the sandpack for 24 hours, SA-SAGD with 2.8-mol% condensate was performed at 50 cm 3 /min (cold-water equivalent) at 3500 kPa for 3 days. The experimental data of production, injection, and temperature distribution were recorded. Also, 10 samples of produced oil were taken and analyzed for density and asphaltene content. The sandpack was excavated after the experiment to analyze the asphaltene content in the remaining oil at different locations. A numerical simulation model was calibrated based on the data of material balance and temperature distribution, and it was validated with properties of the produced and excavated samples. The simulation model used fluid models based on experimental data of viscosities, densities, and bubblepoints for four condensate/bitumen mixtures. Results showed that SA-SAGD was efficient in bitumen recovery with a cumulative steam-to-oil ratio (SOR) that was two to three times smaller than that in SAGD using the same physical model. Detailed analysis of the calibrated simulation model indicated that SA-SAGD enabled the steam chamber to expand more efficiently with a smaller amount of water throughput than SAGD. Volatile solvent components tended to remain in the chamber, and the condensed solvent components acted as a miscible carrier for bitumen components. The analysis further showed that the more efficient oil recovery in SA-SAGD occurred with predominantly cocurrent flow of oil and water near the chamber edge. SA-SAGD recovered a larger amount of asphaltene components (i.e., less in-situ upgrading) than SAGD likely because of its lower chamber temperature, shorter production period, and enhanced local displacement efficiency.

Summary The limiting capillary pressure of foam ( P c * ) and foam trapping in porous media are pore-scale foam properties that affect foam transport in porous media. They are strongly influenced by the characteristics of rock pores and throats. Because of experimental limitations, these foam properties are difficult to measure at core scale. As a result, our understanding of their relationship with different pore characteristics is limited. In this paper, novel coreflood and graphical analysis techniques were used to measure P c * and the foam-trapping coefficient (FTC) at core scale. FTC is a new parameter synonymous to Land’s (1968) trapping coefficient, which describes foam-trapping behavior across an entire range of saturation as opposed to a single endpoint trapped saturation. The scalability of these two foam properties with permeability and other pore characteristics [average pore size (PS), average throat size (TS), average aspect ratio (AR), coordination number (CN), surface area/volume ratio, and reservoir-quality index (RQI)] were also investigated. Pore characteristics of 12 different rock samples were measured from 3D pore-network models generated from high-resolution X-raycomputed-microtomography images. The heterogeneity of the rock samples were quantified by the Dykstra-Parsons index ( Dysktra and Parsons 1950 ), while the RQI and J -function methods were used to classify them according to their storage and flow properties. Each of the measured pore characteristics and their combination [combined pore character (CPC)] were then correlated with P c * and FTC to understand their respective roles. Furthermore, the data points obtained from the graphical analysis of the coreflood data provided the required input data for a mechanistic foam model for relative permeability of foamed gas ( Kovscek and Radke 1994 ). The estimated relative permeability of foamed gas was then used to study foam mobility in the different pore geometries. The overall results showed the following: P c * has strong negative correlations with all pore characteristics except AR, which has a weak positive correlation. P c * has the strongest correlation with RQI, CPC, and permeability; a moderate correlation with CN and TS; and a very weak correlation with PS. Foam trapping has positive correlations with all pore characteristics except AR, which has a negative correlation. Low AR appears to be responsible for significant trapping of foam in high-permeabilityrocks. Low AR favors more foam trapping, while high AR favors trapping of oil and gas during water imbibition in water-wetrocks. Foam trapping appears to have the dominant control on foam mobility.

Summary Lateral propagation of foam in heterogeneous reservoirs, where pore geometries vary laterally, depends on the roles of pore geometries on the foam properties. In this paper, the pore attributes of 12 different rock samples were characterized in terms of porosity, permeability, pore shape, pore size, throat size, aspect ratio, coordination number, and log mean of surface relaxation times ( T 2LM ). These were measured from gas porosimeter and permeameter, X-ray microcomputed tomography (CT)-basedpore-network models, thin-section photomicrographs, and nuclear magnetic resonance (NMR) surface relaxometry. The samples have a wide range of porosity: 12 to 29%; permeability: 1 to 5,000 md; average pore size: 3.7 to 9 µm; average throat size: 2.4 to 8 µm; average aspect ratio: 1 to 1.7; average coordination number: 2.6 to 5.2; and T 2LM : 9.4 to 740 ms. Nitrogen foam flow experiments (without oil) were then conducted on each rock sample using a specialized coreflood apparatus. A graphical analysis of the coreflood data was used to estimate the total saturation of trapped foam (10 to 66%), flowing foam (3 to 14%), and apparent viscosity of foam (3.2 to 73 cp). Trapped foam saturation and apparent viscosity values were then correlated with each of the measured pore attributes. The results revealed that all pore attributes, except aspect ratio, have positive correlations with foam trapping and apparent viscosity. The best correlation with trapped foam saturation was obtained when the most influential pore attributes (pore size, throat size, aspect ratio, and coordination number) were combined into a single mathematical function. Foam apparent viscosity also appears to be mostly influenced by trapped foam saturation, permeability, and coordination number of pore systems. Trapping is also likely enhanced by the presence of fenestral or channel pores. Furthermore, the shape and angularity of pores seem to facilitate snap-off and trapping of foam, because rock samples with angular pores trapped the highest foam saturation compared with other samples with rounded and subrounded pores. It was also shown that the correlation between trapped foam saturation (and foam apparent viscosity) and the absolute permeability of porous media may reverse at some high-permeability values (greater than several darcies), when one or both of the following conditions exist: (1) The aspect ratio of a lower-permeability porous medium is lower than that of a higher-permeability porous medium, and (2) the coordination number of a lower-permeability porous medium is higher than that of a higher-permeability porous medium. Finally, T 2LM showed a good correlation with foam trapping, making NMR logging a prospective tool for pre-evaluating foam performance in targeted reservoir sections.

Summary Numerical simulation of coupled multiphase multicomponent flow and transport in porous media is a crucial tool for understanding and forecasting of complex industrial applications related to the subsurface. The discretized governing equations are highly nonlinear and usually need to be solved with Newton’s method, which corresponds with high computational cost and complexity. With the presence of capillary and gravity forces, the nonlinearity of the problem is amplified even further, which usually leads to a higher numerical cost. A recently proposed operator-based linearization (OBL) approach effectively improves the performance of complex physical modeling by transforming the discretized nonlinear conservation equations into a quasilinear form according to state-dependent operators. These operators are approximated by means of a discrete representation on a uniform mesh in physical parameter space. Continuous representation is achieved through the multilinear interpolation. This approach provides a unique framework for the multifidelity representation of physics in general-purpose simulation. The applicability of the OBL approach was demonstrated for various energy subsurface applications with multiphase flow of mass and heat in the presence of buoyancy and diffusive forces. In this work, the OBL approach is extended for multiphase multicomponent systems with capillarity. Through the comparisons with a legacy commercial simulator using a set of benchmark tests, we demonstrate that the extended OBL scheme significantly improves the computational efficiency with the controlled accuracy of approximation and converges to the results of the conventional continuous approach with an increased resolution of parametrization.

Summary Phase behavior and physical properties including saturation pressures, swelling factors (SFs), phase volumes, dimethyl ether (DME) partition coefficients, and DME solubility for heavy-oil mixtures containing polar substances have been experimentally and theoretically determined. Experimentally, novel phase behavior experiments of DME/water/heavy-oil mixtures spanning a wide range of pressures and temperatures have been conducted. More specifically, a total of five pressure/volume/temperature (PVT) experiments consisting of two tests of DME/heavy-oil mixtures and three tests of DME/water/heavy-oil mixtures have been performed to measure saturation pressures, phase volumes, and SFs. Theoretically, the modified Peng-Robinson equation of state (EOS) (PR EOS) together with the Huron-Vidal mixing rule, as well as the Péneloux et al. (1982) volume-translation strategy, is adopted to perform phase-equilibrium calculations. The binary-interaction parameter (BIP) between the DME/heavy-oil pair, which is obtained by matching the measured saturation pressures of DME/heavy-oil mixtures, works well for DME/heavy-oil mixtures in the presence and absence of water. The new model developed in this work is capable of accurately reproducing the experimentally measured multiphase boundaries, phase volumes, and SFs for the aforementioned mixtures with the root-mean-squared relative error (RMSRE) of 3.92, 9.40, and 0.92%, respectively, while it can also be used to determine DME partition coefficients and DME solubility for DME/water/heavy-oil systems.

Summary Enhanced oil recovery (EOR) in fractured carbonate reservoirs is challenging because of the heterogeneous and oil-wet nature. In this work, a new application of using polymer nanospheres (PNSs) and diluted microemulsion (DME) is presented to plug fractures and enhance water imbibition to recover oil from the tight, naturally fractured carbonate reservoirs. DME with different electric charges is compared through contact-angle and core-imbibition tests to evaluate their performances on EOR. The cationic DME is chosen because it has the fastest wettability-alteration rate and thus the highest oil recovery rate. Migration and plugging efficiency tests are conducted to identify the screened particle sizes of PNSs for the target reservoir cores. PNSs with a particle size of 300 nm are demonstrated to have the best performance of in-depth propagation before swelling and plugging after swelling within the naturally fractured cores are used in this study. Then coreflooding experiments are conducted to evaluate the EOR performance when PNSs and DME are used together, and results indicate that the oil recovery rate is increased by 24.3 and 44.1% compared to using PNSs or DME alone. In the end, a microfluidic experiment is carried out to reveal how DME works with PNSs.

Summary The objectives of this study are to perform a fundamental analysis of the mutual interactions between crude oil components, water, hydrocarbon solvents, and clays, and to determine the optimum hydrocarbon solvent in solvent steamflooding for a particular reservoir type. The water/oil emulsion formation mechanism in the obtained oil for steam and solvent steamflooding processes has been studied via intermolecular associations between asphaltenes, water, and migrated clay particles. A series of 21 steam and solvent-steamflooding experiments has been conducted, first without any clays in the oil/sand packing, and then using two different clay types in the reservoir rock: Clay 1, which is kaolinite, and Clay 2, which is a mixture of kaolinite and illite. Paraffinic (propane, n -butane, n -pentane, n -hexane, n -heptane) and aromatic (toluene) solvents are coinjected with steam. Cumulative oil recovery is found to decrease in the following order: no clay, Clay 1, Clay 2. Based on the obtained produced oil analyses, Clay 1 and Clay 2 are found to have an affinity with the water and oil phases, respectively. Moreover, the biwettable nature of Clay 2 makes it dispersed in the oil phase toward the oil/water interface, stabilizing the water/oil emulsions. Paraffinic solvent n -hexane is found to be an optimum coinjector for solvent steamflooding in bitumen recovery.

Summary Surfactant-polymer (SP) flooding has been regarded as an efficient technique for enhanced oil recovery in the development of mature oil fields, especially for those with heterogeneous conglomerate reservoirs. However, people are still unclear about the optimal SP flooding initiation timing (OSPT) that is expected to contribute to the maximum ultimate recovery factor in the case with a limited amount of SP solution injection. Accordingly, this study aims to investigate OSPT through conducting a series of experiments, including nuclear magnetic resonance (NMR) online monitoring, full-diameter coreflooding, and microfluidic study. The fractional-flow curve is used to identify OSPT, of which the effect on the oil recovery is analyzed. OSPT is demonstrated to be dependent on the amount of injected SP solution. An earlier-started SP flooding is favorable for achieving higher oil recovery factors under the premise of sufficiently high SP solution injection [more than 1.5 pore volume (PV)]. With the commonly used 0.65 PV of SP solution in the reservoir scale, OSPT is suggested to be at the moment when a water cut of 80 to 90% is reached. The formation of dense emulsions in the early-started SP flooding affects the performance of the post-waterflooding, which eventually decreases the ultimate oil recoveries because of inadequacy of SP solution. An earlier-started SP flooding contributes to a larger swept volume, but the initial efficiency of the SP flooding is lower than that of the waterflooding when the injection pressure is constant. OSPT is proposed through analyzing the fractional-flow curve in the case of 0.65 PV of SP injection, and the determined OSPT is validated by coreflooding experiments and field data. Moreover, OSPT for the conglomerate reservoir is suggested to be earlier than that for the relatively homogenous sandstone reservoir.

Summary Foam injection is a promising enhanced-oil-recovery (EOR) technology that significantly improves the sweep efficiency of gas injection. Simulation of foam/oil displacement in reservoirs is an expensive process for conventional simulation because of the strongly nonlinear physics, such as multiphase flow and transport with oil/foam interactions. In this work, an operator-based linearization (OBL) approach, combined with the representation of foam by an implicit-texture (IT) model with two flow regimes, is extended for the simulation of the foam EOR process. The OBL approach improves the efficiency of the highly nonlinear foam-simulation problem by transforming the discretized nonlinear conservation equations into a quasilinear form using state-dependent operators. The state-dependent operators are approximated by discrete representation on a uniform mesh in parameter space. The numerical-simulation results are validated by using three-phase fractional-flow theory for foam/oil flow. Starting with an initial guess depending on the fitting of steady-state experimental data with oil, the OBL foam model is regressed to experimental observations using a gradient-optimization technique. A series of numerical validation studies is performed to investigate the accuracy of the proposed approach. The numerical model shows good agreement with analytical solutions at different conditions and with different foam parameters. With finer grids, the resolution of the simulation is better, but at the cost of more expensive computations. The foam-quality scan is accurately fitted to steady-state experimental data, except in the low-quality regime. In this regime, the used IT foam model cannot capture the upward-tilting pressure gradient (or apparent viscosity) contours. 1D and 3D simulation results clearly demonstrate two stages of foam propagation from inlet to outlet, as seen in the computed-tomography (CT) coreflood experiments: weak foam displaces most of the oil, followed by a propagation of stronger foam at lower oil saturation. OBL is a direct method to reduce nonlinearity in complex physical problems, which can significantly improve computational performance. Taking its accuracy and efficiency into account, the data-driven OBL-based approach could serve as a platform for efficient numerical upscaling to field-scale applications.

Summary Spontaneous imbibition (SI) is regarded as an effective method to improve the oil recovery in a tight sandstone reservoir, which leads to a significant change in fracturing design and flowback treatment. However, a longtime shut-in period would aggravate the retention of fracturing fluid, which is in contradiction with high production in the field. It is imperative to understand how SI works during shut-in time, so as to maximize the effect of imbibition in oil recovery enhancement. In this study, a series of experiments were conducted to simulate the status of residual oil saturation so that the inner mechanism of imbibition on oil recovery can be investigated. Low-field nuclear magnetic resonance (LF-NMR) was used to provide direct observation of phase changes in different pore sizes. The experimental results show a positive effect of imbibition on residual oil reduction. This phenomenon further elucidates the observations made during the well shut-in, soaking period, and low flowback efficiency. This study aims to understand the mechanism of SI behavior and help to improve the accuracy of production prediction.

Summary As the crucial step in closed-loop reservoir management, robust life-cycle production optimization is defined as maximizing/minimizing the expected value of a predefined objective (cost) function over geological uncertainties (i.e., uncertainties in the reservoir permeability, porosity, endpoint relative permeability, etc.). However, with robust optimization, there is no control over downside risk defined as the minimum net present value (NPV) among the individual NPVs of the different reservoir models. Yet, field operators generally wish to keep this minimum NPV reasonably large to try to ensure that the reservoir is commercially viable. In addition, the field operator may desire to maximize the NPV of production over a much shorter time period than the life of the reservoir under the limitation of surface facilities (e.g., field liquid and water production rates). Thus, it is important to consider multiobjective robust production optimization with nonlinear constraints and when geological uncertainties are incorporated. The three objectives considered in this paper are; to maximize the average life-cycle NPV, to maximize the average short-term NPV, and to maximize the minimum NPV of the set of realizations. Generally, these objectives are in conflict; for example, the well controls that give a global maximum for robust life-cycle production optimization do not usually correspond to the controls that maximize the short-term average NPV of production. Moreover, handling the nonlinear state constraints (e.g., field liquid production rates and field water production rates for the bottom-hole pressure controlled producers in the robust production optimization) is also a challenge because those nonlinear constraints should be satisfied at each control steps for each geological realization. To provide potential solutions to the multiobjective robust optimization problem with state constraints, we developed a modified lexicographic method with a minimizing-maximum scheme to attempt to obtain a set of Pareto optimal solutions and to satisfy all nonlinear constraints. We apply the sequential quadratic programming filter with modified stochastic gradients to solve a sequence of optimization problems, where each solution is designed to generate a single point on the Pareto front. In the modified lexicographic method, the objective is always considered to be the primary objective, and the other objectives are considered by specifying bounds on them to convert them to state constraints. The temporal damping and truncation schemes are applied to improve the quality of the stochastic gradient on nonlinear constraints, and the minimizing–maximum procedure is applied to enforce constraints on the normal state constraints. The main advantage that the modified lexicographic method has over the standard lexicographic method is that it allows for the generation of potential Pareto optimal points, which are uniformly spaced in the values of the second and/or third objective that one wishes to improve by multiobjective optimization.