With the dramatic increase in the worldwide demand for the crude oil and with the fact that the oil and gas resources are depleting, the enhanced oil recovery process plays an important role to increase the production from the existing hydrocarbon reservoirs. Chemical enhanced oil recovery is one of the most important techniques to unlock significant amount of trapped oil from oil reservoirs. Surface agent materials (Surfactants) are used to lower the interfacial tension (IFT) between water and oil phases to ultralow values and mobilize the trapped oil. When surfactant, water, and oil are mixed together they form a thermodynamically stable phase called microemulsion which can be characterized by ultralow interfacial tension and the ability to solubilize both aqueous and oil compounds. Another characteristic of microemulsion solution is its viscosity which plays an important role in the creation and movement of the oil bank. The microemulsion micro-structure is complex and its viscosity is difficult to predict. Various viscosity models and correlations are presented in the literature to describe microemulsion viscosity behavior, but they fail to represent the rheological behavior of many microemulsion mixtures. Most of these models are valid in the lower and higher ranges of solute where one of the domains is discontinuous. The majority of the models fail to calculate the rheology of microemulsion phase in bicontinuous domains. In this work, we present a systematic study of the rheological behavior of microemulsion systems and the effect of additives such as polymer and co-solvent on rheological properties of microemulsions. Several laboratory experiments were conducted to determine the rheological behavior of surfactant solutions. A new empirical model for the viscosity of microemulsion phase as a function of salinity is introduced. The model consists of three different correlations one for each phase type of Windsor phase behaviors. The proposed model is validated using a number of experimental results presented in this document. The proposed viscosity model is implemented in the UTCHEM simulator and the simulator results are compared with the coreflood experiments. Excellent matches were obtained for the pressure. We further improved the proposed viscosity model to incorporate the effect of polymer and co-solvent on the microemulsion viscosity. / text

Polymer flooding is a commercial enhanced oil recovery (EOR) method used to increase the sweep efficiency of water floods. Hydrolyzed polyacrylamide (HPAM), a synthetic commercial polymer, is widely used in commercial polymer floods and it is also used for mobility control of chemical floods using surfactants such as surfactant-polymer flooding and alkaline-surfactant-polymer flooding. The increase in the degree of hydrolysis of HPAM at elevated temperature or pH with time affects the polymer solution viscosity and its adsorption on rock surfaces. A polymer hydrolysis model based on published laboratory data was implemented in UTCHEM, a chemical EOR simulator, in order to assess the effect of hydrolysis on reservoir performance. Both 1D and 3D simulations were performed to validate the implementation of the model. The simulation results are consistent with the laboratory observations that show an increase in polymer solution viscosity as hydrolysis progresses. The numerical results indicate that hydrolysis occurs very rapidly and impacts the near wellbore region polymer injectivity. / text

Three compositional reservoir simulators have been developed in the Department of Petroleum and Geosystems Engineering at The University of Texas at Austin (UT-Austin): UTCOMP (miscible gas flooding simulator), UTCHEM (chemical flooding simulator), and GPAS (General Purpose Adaptive Simulator). UTCOMP and UTCHEM simulators have been used by various oil companies for solving a variety of field problems. The efficiency and accuracy of each simulator becomes critically important when they are used to solve field problems. In this study, two well-developed solver packages, SAMG and HYPRE, along with existing solvers were compared. Our numerical results showed that SAMG can be an excellent solver for the usage in the three simulators for solving problems with a high accuracy requirement and long simulation times, and BoomerAMG in HYPRE package can also be a good solver for application in the UTCHEM simulator. In order to investigate the flexibility and the efficiency of a partitioned coupling method, the second part of this thesis presents a new implementation using a partition method for a thermal module in an equation-of-state (EOS) compositional simulator, the General Purpose Adaptive Simulator (GPAS) developed at The University of Texas at Austin. The finite difference method (FDM) was used for the solution of governing partial differential equations. Specifically, the new coupled implementation was based on the Schur complement method. For the partition method, two suitable acceleration techniques were constructed. One technique was the optimized choice of preconditioner for the Schur complement; the other was the optimized selection of tolerances for the two solution steps. To validate the implementation, we present simulation examples of hot water injection in an oil reservoir. The numerical comparison between the new implementation and the traditional, fully implicit method showed that the partition method is not only more flexible, but also faster than the classical, fully implicit method for the same test problems without sacrificing accuracy. In conclusion, the new implementation of the partition method is a more flexible and more efficient method for coupling a new module into an existing simulator than the classical, fully implicit method.The third part of this thesis presents another type of coupling method, iterative coupling methods, which has been implemented into GPAS with thermal module, FICM (Fully, Iterative Coupling Method) and GICM (General, Iterative Coupling Method), LICM (Loose, Iterative Coupling Method). The results show that LICM is divergent, and GICM and FICM can work normally. GICM is the fastest among the compared methods, and FICM has a similar efficiency as CFIM (Classic Fully Implicit Method). Although GICM is the fastest method, GICM is less accurate than FICM for in the test cases carried out in this study. / text

Scale deposition, either in the formation or inside the tubing, is a serious problem that can affect the productivity of oil fields. Production sustainability depends on the successful implementation of scale management strategies prior to developing new fields. Such strategies should involve tools capable of addressing the risks of developing scales during the production stage as well as determining the outcomes of possible remediation jobs in the future. UTCHEM, a multi-compositional flow model, was used in this work to present a comprehensive study that includes both precipitation and remediation scenarios. Although there are different mechanisms prompting the deposition of mineral scales, barite and calcite were selected primarily to simulate the effect of mixing incompatible water compositions; an issue that is usually associated with seawater injection. Equilibrium state calculations were carried out using a geochemical model (EQBATCH) to verify the incompatibility of the injection water with the formation water. In this work, we show the evolution, distribution, and remediation of solids over time for several hypothetical cases. The quantity of deposits in the near-wellbore region was found to be less at a highly heterogeneous reservoir model in contrast to the amount precipitated in homogenous reservoirs. This could be critical to wells productivity in the long-run since much of the drop in reservoir pressure occurs near the wellbore. The predictive ability of UTCHEM was extended to include simulating the removal of carbonate scales using a chelating chemical. The optimization of the injected treatment can be achieved mechanically through adjusting the well spacing (during the initial stages of field development) or through adjusting the concentrations of active components in the remediation fluids. The model provides a valuable tool that helps planners to predict scaling-related issues ahead of time, and subsequently to determine the economic viability of the project. This work serves as an opportunity to re-assess this simulator and allows for further work to enhance its capabilities. / text

Surfactants both change the wettability and lower the interfacial tension by various degrees depending on the type of surfactant and how it interacts with the specific oil. Ultra low IFT means almost zero capillary pressure, which in turn indicates little oil should be produced from capillary imbibition when the surfactant reduces the IFT in naturally fractured oil reservoirs that are mixed-wet or oil-wet. What is the transport mechanism for the surfactant to get far into the matrix and how does it scale? Molecular diffusion and capillary pressure are much too slow to explain the experimental data. Recent dynamic laboratory data suggest that the process is faster when a pressure gradient is applied compared to static tests. A mechanistic chemical compositional simulator was used to study the effect of pressure gradient on chemical oil recovery from naturally fractured oil reservoirs for several different chemical processes (polymer, surfactant, surfactant-polymer, alkali-surfactant-polymer flooding). The fractures were simulated explicitly by using small gridblocks with fracture properties. Both homogeneous and heterogeneous matrix blocks were simulated. Microemulsion phase behavior and related chemistry and physics were modeled in a manner similar to single porosity reservoirs. The simulations indicate that even very small pressure gradients (transverse to the flow in the fractures) are highly significant in terms of the chemical transport into the matrix and that increasing the injected fluid viscosity greatly improves the oil recovery. Field scale simulations show that the transverse pressure gradients promote transport of the surfactant into the matrix at a feasible rate even when there is a high contrast between the permeability of the fractures and the matrix. These simulations indicate that injecting a chemical solution that is viscous (because of polymer or foam or microemulsion) and lowers the IFT as well as alters the wettability from mixed-wet to water-wet, produces more oil and produces it faster than static chemical processes. These findings have significant implications for enhanced oil recovery from naturally fractured oil reservoirs and how these processes should be optimized and scaled up from the laboratory to the field. / text

The overall objective of this research was to gain an insight into the challenges encountered during chemical flooding under high hardness conditions. Different aspects of this problem were studied using a combination of laboratory experiments and simulation studies. Chemical Flooding is an important Enhanced Oil Recovery process. One of the major components of the operational expenses of any chemical flooding project, especially Alkali Surfactant Polymer (ASP) flooding is the cost of softening the injection brine to prevent the precipitation of the carbonates of the calcium and magnesium ions which are invariably present in the formation brine. Novel hardness tolerant alkalis like sodium metaborate have been shown to perform well with brines of high salinity and hardness, thereby eliminating the need to soften the injection brine. The first part of this research was aimed at designing an optimal chemical flooding formulation for a reservoir having hard formation brine. Sodium metaborate was used as the alkali in the formulation with the hard brine. Under the experimental conditions, sodium metaborate was found to be inadequate in preventing precipitation in the ASP slug. Factors affecting the ability of sodium metaborate to sequester divalent ions, including its potential limitations under the experimental conditions were studied. The second part of this research studied the factors affecting the ability of novel alkali and chelating agents like sodium metaborate and tetrasodium EDTA to sequester divalent ions. Recent studies have shown that both these chemicals showed good performance in sequestering divalent ions under high hardness conditions. A study of the geochemical species in solution under different conditions was done using the computer program PHREEQC. Sensitivity studies about the effect of the presence of different solution species on the performance of these alkalis were done. The third part of this research focused on field scale mechanistic simulation studies of geochemical scaling during ASP flooding. This is one of the major challenges faced by the oil and gas industry and has been found to occur when sodium carbonate is used as the alkali and the formation brine present in situ has a sufficiently high hardness content. The multicomponent and multiphase compositional chemical flooding simulator, UTCHEM was used to determine the quantity and composition of the scales formed in the reservoir as well as the injection and production wells. Reactions occurring between the injected fluids, in situ fluids and the reservoir rocks were taken into consideration for this study. Sensitivity studies of the effect of key reservoir and process parameters like the physical dispersion and the alkali concentration on the extent of scaling were also done as a part of this study. / text