Numerical Modeling of Waterflooding Experiments in Artificially Fractured and Gel Treated Core Plugs by Embedded Discrete Fracture Model of a Reservoir Simulation Toolbox

© 2022, Society of Petroleum Engineers.The fluid flow dynamics of the matrix and fractures are significantly different from each other. Fractures are high-permeability flow channels that serve as the main flow units. On the other hand, the Matrix takes up the majority of the reservoir volume and is generally regarded as the main storage unit. The primary goal of this research is to investigate numerically the effects of fractures and polymer gel treatment on oil recovery during waterflooding of artificially fractured core plugs. In this study, the MATLAB Reservoir Simulation Toolbox (MRST) was used for the numerical solution. Different numerical models were developed using MRST to describe three main cases: non-fractured core plug, fractured core plug, and polymer gel treated core plug. Following the creation of the physical models, 2 PV water was introduced into all core plugs. Oil recovery and water saturation profiles vs. time plots were obtained. The standard Buckley-Leveret solution is utilized to evaluate the numerical model, and the fractures are modeled using the Embedded Discrete Fracture Model (EDFM). The results of the simulations were compared with the results of the experiments. In the experiments, results were recorded after 2 PV water injections. For the polymer gel treated core plugs, 2 PV more water was injected after the polymer gel operation. same injection volumes as used in the MRST model. For an artificially fractured core sample, initial oil recovery was measured as 28.57% experimentally and 28.87% with MRST. Then polymer gel was applied to the core plug, increasing the oil recovery to 42.85% experimentally and to 40.83% with MRST. Similarly, before and after polymer gel operation, mean water saturation was measured as 58.34% and 66.5%, respectively. MRST results showed mean water saturation of 58.38% and 65.45%. It is clear from both numerical and experimental models that the existence of fractures decreases the overall hydrocarbon recovery. Polymer gel treatment decreases fracture permeability, resulting in a more uniform sweep and increased overall recovery. Additional oil recovery was observed after polymer gel treatment. Besides, polymer gel treatment of the matrix is also efficient for increasing the recovery and leads to the same results. Moreover, the effects of the fracture aperture and fracture permeability on the recovery were also investigated. Fracture aperture directly impacts the recovery of the low aperture values when the permeability is constant. Similarly, permeability directly affects recovery for high values when the aperture is constant. Finally, the results showed that experimental and numerical findings are significantly close to each other for all non-fractured, fractured, and polymer gel-treated cases.
2022 SPE International Conference and Exhibition on Formation Damage Control, FD 2022


Kaya, Onur Alp; Durgut, İsmail; Department of Petroleum and Natural Gas Engineering (2021-9)
The dynamics of fluid flow in the matrix and in the fractures are significantly different from each other. Fractures are highly permeable flow pathways and are generally considered as the main flow unit. On the other hand, the matrix occupies most of the porous medium’s volume and is considered as the main storage unit. The main objective of this thesis is numerical modeling of water flooding experiments in artificially fractured and gel-treated core plugs. MATLAB Reservoir Simulation Toolbox (MRST) is used...
Numerical simulation of fluid flow and heat transfer in a trapezoidal microchannel with COMSOL multiphysics: A case study
Turgay, Metin Bilgehan; Güvenç Yazıcıoğlu, Almıla (2018-01-01)
In this study, fluid flow and heat transfer in a trapezoidal microchannel are numerically investigated. For this purpose, a reference study with experimental and numerical solutions is adopted from the literature and solved with COMSOL multiphysics. Good agreement with the results of the reference work is obtained. In addition, effects of stabilization methods and element discretization options that are offered by the program on the results are investigated and discussed with examples. In addition, two diff...
Numerical analysis of a projection-based stabilization method for the natural convection problems
Çıbık, Aytekin Bayram; Kaya Merdan, Songül; Department of Mathematics (2011)
In this thesis, we consider a projection-based stabilization method for solving buoyancy driven flows (natural convection problems). The method consists of adding global stabilization for all scales and then anti-diffusing these effects on the large scales defined by projections into appropriate function spaces. In this way, stabilization acts only on the small scales. We consider two different variations of buoyancy driven flows based on the projection-based stabilization. First, we focus on the steady-sta...
Sensitivity Analysis of Major Drilling Parameters on Cuttings Transport during Drilling Highly-inclined Wells
Ozbayoglu, E. M.; Miska, S. Z.; Takach, N.; Reed, T. (Informa UK Limited, 2009-01-01)
In this study, a layered cuttings transport model is developed for high-angle and horizontal wells, which can be used for incompressible non-Newtonian fluids as well as compressible non-Newtonian fluids (i.e., foams). The effects of major drilling parameters, such as flow rate, rate of penetration, fluid density, viscosity, gas ratio, cuttings size, cuttings density, wellbore inclination and eccentricity of the drillsting on cuttings transport efficiency are analyzed. The major findings from this study are,...
Numerical solution of buoyancy MHD flow with magnetic potential
Pekmen, B.; Tezer, Münevver (2014-04-01)
In this study, dual reciprocity boundary element method (DRBEM) is applied for solving the unsteady flow of a viscous, incompressible, electrically conducting fluid in channels under the effect of an externally applied magnetic field and buoyancy force. Magnetohydrodynamics (MHD) equations are coupled with the energy equation due to the heat transfer by means of the Boussinessq approximation. Then, the 20 non-dimensional full MHD equations in terms of stream function, temperature, magnetic potential, curren...
Citation Formats
O. A. Kaya, İ. Durgut, and S. Canbolat, “Numerical Modeling of Waterflooding Experiments in Artificially Fractured and Gel Treated Core Plugs by Embedded Discrete Fracture Model of a Reservoir Simulation Toolbox,” Louisiana, Amerika Birleşik Devletleri, 2022, vol. 2022-February, Accessed: 00, 2022. [Online]. Available: