Sharif Digital Repository

The effectiveness of secondary recovery methods in reservoir development studies depends on the knowledge about how fluid-carrying regions (i.e. good-quality rock types) are connected between injection and production wells. To estimate reservoir performance uncertainty, comprehensive simulations on many reservoir model realisations are necessary, which is very CPU consuming and time demanding. Alternatively, we can use much simpler and physically based methods such as percolation approach. Classic percolation assumes connectivity between opposite 2-D faces of a 3-D system; whereas, hydrocarbon production is achieved through active wells that are one-dimensional lines (e.g. vertical,... 

In this work, hybrid of silica nanoparticles (NPs) with sulfonated copolymer has been prepared in order to improve the copolymer properties for Enhanced Oil Recovery (EOR). Some tests are done to find the effectiveness of injecting fluid on Recovery Factor (RF), such as Contact Angle (CA), interfacial tension (IFT), inherent viscosity, and eventually Micro-Model (MM) flooding. In CA test, wettability alteration from Oil-Wet (OW) to Water-Wet (WW) is reached by sedimentation and adsorption of NPs on the rock slice. In addition, IFT reduction is obtained by increasing the NaCl concentration. The viscosity change is investigated for the hybrid and copolymer under simulated high-temperature oil... 

Relative permeability curves of two-phase flow in a fracture have been a subject of study in recent years. The importance of these curves have been widely observed in multidisciplines, such as water subsurface resources, geothermal energy and underground hydrocarbon resources, especially fractured oil and gas reservoirs. Extensive experimental studies have been cited alongside the numerical studies in this area. However, simple analytical and practical solutions are still attractive. In the current study, wettability effects and phase interference explicitly were tried to be implemented in a simple analytical formula. The wettability effects are represented by residual saturations which... 

In the optimization of naturally fractured reservoirs, it is required to take into account their complex flow behavior due to high conductivity fractures. In this regard, the possible effects of fractures must be included in the optimization procedure. In a water-flooding project, fast water breakthrough from injection to production wells may be occurred because of high permeability fractures. To consider the effect of the fracture system, a multi-criterion optimization procedure is proposed in this work. For this contribution, Non-dominated Sorting Genetic Algorithm version II (NSGA-II) is implemented for the optimization purposes. Considering the effect of the fracture system on the flow... 

Stress sensitivity is a phenomenon that affects reservoir rock properties, such as permeability and therefore changes the well pressure transient behavior. This paper aims to study these behaviors in stress sensitive reservoirs and evaluate the pressure loss in such reservoirs during the process of hydrocarbon production. A power model is used to correlate the changes in permeability with pore pressure. A novel semi-implicit three-dimensional finite element method has been employed to numerically solve the flow problem. The numerical results have been validated by analytical results obtained in a non-sensitive reservoir. Pressure drawdown test for different scenarios has been studied. The... 

A sequential implicit numerical method based on discrete-fracture model and the Galerkin Finite Element method, for time-dependent coupled fluid flow and geomechanics problems in fractured subsurface formations is presented. Discrete-fracture model has been used to explicitly represent the fracture network inside porous media. The Galerkin Finite Element method with adaptive unstructured gridding is implemented to numerically solve the spatially three-dimensional and time-dependent problem. The presented method is validated with previously obtained solutions. Two problems are numerically solved by applying the presented methodology in a three-dimensional fractured petroleum reservoir under... 

In this study, impact of initial ensembles on posterior distribution of ensemble-based assimilation methods is statistically analyzed. Along with, sampling performance as well as uncertainty quantification of these methods are compared in terms of their ability to accurately and consistently evaluate unknown reservoir model parameters and reservoir future performance. For this purpose, a synthetic test problem, which is a small but highly nonlinear reservoir model under two-phase flow, is utilized. Subsequently, different initial ensemble sets are considered and are updated through the assimilation process performed on the test problem using ensemble-based assimilation methods. Afterwards,... 

Colloidal gas aphron (CGA) based fluids have recently been introduced to develop depleted hydrocarbon reservoirs due to their ability in controlling fluid losses. Bubbles size of CGAs plays an important role in pore blockage ability as the controlling mechanism in fluid invasion reduction. However, fundamental understanding of how bubble size distribution is controlled by polymer and surfactant concentrations is not well discussed in the available literature. Almost all reported experiences on CGAs sizing were conducted on single bubble behavior, and little attention has been given to the variation of bubble size distribution (BSD) of CGAs at different levels of polymer/surfactant... 

Concerns about climate change have persuaded the researchers to examine CO2 injection in the form of carbonated water (CW) into oil reservoir as a safe and effective CO2 storage and enhanced oil recovery process. Although interfacial tension (IFT) between crude oil and injected fluid has a vital role on the displacement of fluids in porous media, the effect of CO2 and crude oil type on the dynamic IFT of crude oil under different operational conditions is not well understood. Accordingly, this study was carried out to assess the effects of temperature, pressure, crude oil type and CO2 on the dynamic IFT of crude oil/CW. To achieve this goal, two types of crude oil were provided from southern... 

Ensemble-based assimilation methods are the most promising tools for dynamic characterization of reservoir models. However, because of inherent assumption of Gaussianity, these methods are not directly applicable to channelized reservoirs wherein the distribution of petrophysical properties is multimodal. Transformation of facies field to level set functions have been proposed to alleviate the problem of multimodality. Level set representation ensures that the estimated fields are facies realizations as well as no modification of the assimilation method is required. Moreover, due to the complexity of the history matching problem in the channelized reservoirs, most researchers conventionally... 

Accurate determination of relative permeability (kr) curves and their hysteresis is needed for reliable prediction of the performance of oil and gas reservoirs. A few options (e.g., Carlson, Killough and Jargon models) are available in commercial reservoir simulators to account for hysteresis in kr curves under two-phase flow. Two-phase kr curves are also needed for estimating kr hysteresis under three-phase flow during water-alternating-gas (WAG) injection. Although, most oil reservoirs are mixed-wet, the existing hysteresis predictive approaches have been developed based on water-wet conditions. Experimentally measured data are needed to assess the performance of these methodologies under... 

Scaling imbibition curves are of great importance in the characterization and simulation of oil production from naturally fractured reservoirs. Different parameters such as matrix porosity and permeability, oil and water viscosities, matrix dimensions, and oil/water interfacial tensions have an effective on the imbibition process. Studies on the scaling imbibition curves along with the consideration of different assumptions have resulted in various scaling equations. In this work, using an artificial neural network (ANN) method, a novel technique is presented for scaling imbibition recovery curves, which can be used for scaling the experimental and field-scale imbibition cases. The... 

Purpose - This paper aims to perform a comparative study between capabilities of two numerical schemes from two main branches of numerical methods for solving hyperbolic conservation equations. Design/methodology/ approach - The accuracy and performance of a newly developed high-resolution central scheme vs a higher-order Godunov-based method are evaluated in the context of black-oil reservoir simulations. Both methods are modified enabling study of applications that are not strictly hyperbolic and exhibit local linear degeneracies in their wave structure. Findings - The numerical computations show that while both schemes produce results with virtually the same accuracy, the Godunov method... 

This paper presents a new combined method for accurate upscaling of two-phase displacements in highly heterogeneous reservoirs. The method has the capability to retain its high performance for various flow regimes, from viscous to gravity dominant displacements, without the need for further modifications and computational steps. Two different grids are incorporated for simulation. The grid on fine scale is used to recognize the complicated physics of flow which depends on dominated driving forces and their interaction with heterogeneity. However, to achieve a fast simulation, the global flow calculation is performed on the coarse scale grid using upscaled equivalent properties. The... 

The Steam Assisted Gravity Drainage (SAGD) process, a developed Enhanced Oil Recovery (EOR) process to recover oil and bitumen, has been studied theoretically and experimentally in conventional reservoirs and models and is found a promising EOR method for certain heavy oil reservoirs. In this work simulation studies of the SAGD process were made on different fractured models consisting of fractures in both Near Well Region (NWR) and Above Well Region (AWR) and even in the presence of networked fractures. At early stage of the SAGD process in fractured system, steam moves through the fractures first and then the matrix blocks are heated primarily by conduction and possibly some steam... 

Nitrogen has emerged as a suitable alternative to carbon dioxide for injection into hydrocarbon reservoirs worldwide to enhance the recovery of subsurface energy. Nitrogen typically costs less than CO2 and natural gas, and has the added benefit of being widely available and non-corrosive. However, the underlying mechanisms of recovery following N2 injection into fractured reservoirs that make up a large portion of the world's oil and gas reserves are not well understood. Here we present the laboratory results of N2 injection into carbonate rocks acquired from a newly developed oil reservoir in Iran with a huge N2-containing natural gas reservoir nearby. We investigate the effectiveness of N2... 

Wettability alteration of rock surfaces toward gas wetting have been recognized as a practical approach for maximizing the production from the gas condensate reservoirs. Most of the reported work in this area applied the so called sessile drop contact angle measurement technique to examine the change in wetting state of a surface. However, the size-dependent wetting behavior of drop which could affect the exact determination of wettability and wettability changes was not well discussed in the previous studies. Therefore, in this work, the size dependency of contact angle for four different liquid-solid-gas systems; i.e., water-calcite-air, water-treated calcite-air (nanofluid treated... 

Acidizing is one of the most effective techniques to remove the formation damage and restore (or even increase) the permeability of the near wellbore region. Although this technique is widely used to resolve the skin problem, it could initiate new damages itself that hinder the fluid flow and decrease the well production, consequently. Acid-oil emulsion and sludge formation are known as two major induced formation damage and the main reasons for the oil well acid treatment failures. Despite its critical effects, no comprehensive study has been addressed this specific type of formation damages in details. In this study, the acid-oil emulsion and asphaltic sludge formation were evaluated... 

Capillary imbibition is an important recovery mechanism in naturally fractured reservoirs when water-filled fractures surround water-wet matrix blocks. A large amount of studies of imbibition process is simply total or partial immersion of nonwetting phase saturated rock in aqueous wetting phase. However, water advance in fractures during water flooding or water encroachment from an active aquifer introduces time dependent boundary conditions where invariant exposure of rock surface to water is not representative. In this work, a laboratory simulated matrix-fracture system was used to investigate different aspects of imbibition in the presence of fracture fluid flow (namely dynamic... 

Gas condensate production using technology of multi-stage hydraulically fracturing in shale gas condensate reservoirs' horizontal wells is a new topic of unconventional resources studies. Thus, shale gas condensate as a new source of energy can be considered as an important issue for development and further studies. In this work, a semi-analytical solution of gas and oil two-phase flow is presented for pressure transient analysis (PTA) and rate transient analysis (RTA) of a shale gas condensate reservoir's production data. Fluid flow assumption here is flow in a pseudo triple-porosity porous media, which are matrix, natural fractures and adsorbed gas. Adsorbed gas is a form of gas in porous...