Mechanisms Involved in Microbially Enhanced Oil Recovery

Microbial enhanced oil recovery (MEOR) represents a possible cost-effective tertiary oil recovery method. Although the idea of MEOR has been around for more than 75 years, even now little is known of the mechanisms involved. In this study, Draugen and Ekofisk enrichment cultures, along with Pseudomonas spp. were utilized to study the selected MEOR mechanisms. Substrates which could potentially stimulate the microorganisms were examined, and l-fructose, d-galacturonic acid, turnose, pyruvic acid and pyruvic acid methyl ester were found to be the best utilized by the Ekofisk fermentative enrichment culture. Modelling results indicated that a mechanism likely to be important for enhanced oil recovery is biofilm formation, as it required a lower in situ cell concentration compared with some of the other MEOR mechanisms. The bacterial cells themselves were found to play an important role in the formation of emulsions. Bulk coreflood and flow cell experiments were performed to examine MEOR mechanisms, and microbial growth was found to lead to possible alterations in wettability. This was observed as a change in wettability from oil wet (contact angle 154°) to water wet (0°) due to the formation of biofilms on the polycarbonate coupons.     

Microscopic Mechanisms of Oil Recovery By Near-Miscible Gas Injection

As gas flooding becomes a more viable means of enhanced oil recovery, it is important to identify and understand the pore-scale flow mechanisms, both for the development of improved gas flooding applications and for the predicting phase mobilisation under secondary and tertiary gas flooding. The purpose of this study was to visually investigate the pore-level mechanisms of oil recovery by near-miscible secondary and tertiary gas floods. High-pressure glass micromodels and model fluids representing a near-miscible fluid system were used for the flow experiments. A new pore-scale recovery mechanism was identified which significantly contributed to oil recovery through enhanced flow and cross-flow between the bypassed pores and the injected gas. This mechanism is strongly related to a very low gas/oil interfacial tension (IFT), perfect wetting conditions and simultaneous flow of gas and oil in the same pore, all of which occur as the gas/oil critical point is approached. The results of this study helps us to better understand the pore-scale mechanisms of oil recovery in very low-IFT (near-miscible) systems. In particular we show that in near-miscible gas floods, behind the main gas front, the recovery of the oil continues by cross-flow from the bypassed pores into the main flow stream and as a result almost all of the oil, which has been contacted by the gas, could be recovered. Our observations in high-pressure micromodel experiments have demonstrated that this mechanism can only occur in near-miscible processes (as opposed to immiscible and completely miscible processes), which makes oil displacement by near-miscible gas floods a very effective process.     

Mechanisms of Microscopic Displacement During Enhanced Oil Recovery in Mixed-Wet Rocks Revealed Using Direct Numerical Simulation

Transport in Porous Media - We demonstrate how to use numerical simulation models directly on micro-CT images to understand the impact of several enhanced oil recovery (EOR) methods on microscopic...     

Fundamental Study of Pore Scale Mechanisms in Microbial Improved Oil Recovery Processes

A fundamental study of microscopic mechanisms and pore-level phenomena in the Microbial Improved Oil Recovery method has been investigated. Understanding active mechanisms to increase oil recovery is the key to predict and plan MIOR projects successfully. This article presents the results of visualization experiments carried out in a transparent pore network model. In order to study the pore scale behavior of bacteria, dodecane and an alkane oxidizing bacterium, Rhodococcus sp. 094, suspended in brine, are examined for evaluating the performance of bacterial flooding in the glass micromodel. The observations show the effects of bacteria on remaining oil saturation, allowing us to get better insight on the mechanisms. Bacterial mass composed of bacteria and bioproducts growth in the fluid interfaces and pore walls have been recorded and are presented. No gas is observed throughout any of the experiments. The biomass blocks some pores and pore-throats, and thereby changing the flow pattern. As a consequent, the flow pattern change together with the previously proposed mechanisms, including the interfacial tension reduction and wettability changes are recognized as active mechanisms in the MIOR process.     

Quantitative In Situ Enhanced Oil Recovery Monitoring Using Nuclear Magnetic Resonance

Quantitative in situ monitoring of oil recovery from sedimentary rock is demonstrated for the first time using advanced two-dimensional (2D) nuclear magnetic resonance (NMR) correlation measurements on a low field spectrometer. The laboratory-scale NMR system was chosen to provide a common physics of measurement with NMR well-logging tools. The NMR protocols are used to monitor recovery of a heavy Middle East crude oil from high permeability sandstone plugs using a brine (water) flood followed by chemical enhanced oil recovery agents: polymer and alkaline?Csurfactant?Cpolymer solutions. 2D correlations between relaxation time (T ₁, T ₂) and apparent self-diffusion coefficient D _app are used to obtain simultaneously a volumetric determination of the oil and aqueous fluid-phase saturations present in the porous material. The T ₁ ? T ₂ and D _app ? T ₂ correlations are bulk measurements of the entire rock core-plug; excellent agreement is shown between the measures of remaining oil (from NMR) and recovered oil (from gravimetric assay of the effluent). Furthermore, we introduce the capability to measure spatially resolved T ₂ distributions on a low field spectrometer using a rapid frequency-encoded y ? T ₂ map. A non-uniform distribution of remaining oil is observed due to viscous instabilities in the flowing liquids; the final oil saturation ranges from ${S_{\rm o}^{\rm{(final)}} \approx 0}$ to 20?% along the direction of flow. These results highlight the quantitative nature of the NMR data obtainable in low field NMR core analysis and also the importance of spatially resolved measurements when studying short core-plugs.     

Turbulence Modelling of Unsteady Turbulent Flows Using the Stress Strain Lag Model

This paper reports the application of a recently developed turbulence modelling scheme known as the C _as model. This model was specifically developed to capture the effects of stress-strain misalignment observed in turbulent flows with mean unsteadiness. Earlier work has reported the approach applied within a linear k-ε modelling framework, and some initial testing of it within the k-ω SST model of Menter (AIAA J 32:1598–1605, 1994). The resulting k-ε-C _as and SST-C _as models have been shown to result in some of the advantages of a full Reynolds Stress transport Model (RSM), whilst retaining the computational efficiency and stability benefits of a eddy viscosity model (EVM). Here, the development of the the high-Reynolds-number version of the C _as model is outlined, with some example applications to steady and unsteady homogeneous shear flows. The SST-C _as form of the model is then applied to further, more challenging cases of 2-D flow around a NACA0012 aerofoil beyond stall and the 3-D flow around a circular cylinder in a square duct, both being flows which exhibit large, unsteady, separated flow regions. The predictions returned by a range of other common turbulence modelling schemes are included for comparison and the SST-C _as scheme is shown to return generally good results, comparable in some respects to those obtainable from far more complex schemes, for only moderate computing resource requirements.     

Simulation of Three-Phase Displacement Mechanisms Using a 2D Lattice-Boltzmann Model

Using a numerical technique, known as the lattice-Boltzmann method, we study immiscible three-phase flow at the pore scale. An important phenomenon at this scale is the spreading of oil onto the gas–water interface. In this paper, we recognize from first principles how injected gas remobilizes initially trapped oil blobs. The two main flow mechanisms which account for this type of remobilization are simulated. These are the double-drainage mechanism and (countercurrent) film flow of oil. The simulations agree qualitatively with experimental findings in the literature. We also simulate steady-state three-phase flow (fixed and equal saturations) in a small segment of a waterwet porous medium under both spreading and nonspreading conditions. The difference between the two conditions with respect to the coefficients in the generalized law of Darcy (which also includes viscous coupling) is investigated.     

Visualisation of Residual Oil Recovery by Near-miscible Gas and SWAG Injection Using High-pressure Micromodels

We present results of high-pressure micromodel visualizations of pore-scale fluid distribution and displacement mechanisms during the recovery of residual oil by near-miscible hydrocarbon gas and SWAG (simultaneous water and gas) injection under conditions of very low gas–oil IFT (interfacial tension), negligible gravity forces and water-wet porous medium. We demonstrate that a significant amount of residual oil left behind after waterflooding can be recovered by both near-miscible gas and SWAG injection. In particular, we show that in both processes, the recovery of the contacted residual oil continues behind the main gas front and ultimately all of the oil that can be contacted by the gas will be recovered. This oil is recovered by a microscopic mechanism, which is strongly linked to the low IFT between the oil and gas and to the perfect spreading of the oil over water, both of which occur as the critical point of the gas–oil system is approached. Ultimate oil recovery by near-miscible SWAG injection was as high as near-miscible gas injection with SWAG injection using much less gas compared to gas injection. Comparison of the results of SWAG experiments with two different gas fractional flow values (SWAG ratio) of 0.5 and 0.2 shows that fractional flow of the near-miscible gas injected simultaneously with water is not a crucial factor for ultimate oil recovery. This makes SWAG injection an attractive IOR (improved oil recovery) process especially for reservoirs, where continuous and high-rate gas injection is not possible (e.g. due to supply constraint).     

A Model for Spontaneous Imbibition as a Mechanism for Oil Recovery in Fractured Reservoirs

The flow of oil and water in naturally fractured reservoirs (NFR) can be highly complex and a simplified model is presented to illustrate some main features of this flow system. NFRs typically consist of low-permeable matrix rock containing a high-permeable fracture network. The effect of this network is that the advective flow bypasses the main portions of the reservoir where the oil is contained. Instead capillary forces and gravity forces are important for recovering the oil from these sections. We consider a linear fracture which is symmetrically surrounded by porous matrix. Advective flow occurs only along the fracture, while capillary driven flow occurs only along the axis of the matrix normal to the fracture. For a given set of relative permeability and capillary pressure curves, the behavior of the system is completely determined by the choice of two dimensionless parameters: (i) the ratio of time scales for advective flow in fracture to capillary flow in matrix $\alpha =\tau ^f/\tau ^m$ ; (ii) the ratio of pore volumes in matrix and fracture $\beta =V^m/V^f$ . A characteristic property of the flow in the coupled fracture–matrix medium is the linear recovery curve (before water breakthrough) which has been referred to as the “filling fracture” regime Rangel-German and Kovscek (J Pet Sci Eng 36:45–60, 2002), followed by a nonlinear period, referred to as the “instantly filled” regime, where the rate is approximately linear with the square root of time. We derive an analytical solution for the limiting case where the time scale $\tau ^{m}$ of the matrix imbibition becomes small relative to the time scale $\tau ^{f}$ of the fracture flow (i.e., $\alpha \rightarrow \infty $ ), and verify by numerical experiments that the model will converge to this limit as $\alpha $ becomes large. The model provides insight into the role played by parameters like saturation functions, injection rate, volume of fractures versus volume of matrix, different viscosity relations, and strength of capillary forces versus injection rate. Especially, a scaling number $\omega $ is suggested that seems to incorporate variations in these parameters. An interesting observation is that at $\omega =1$ there is little to gain in efficiency by reducing the injection rate. The model can be used as a tool for interpretation of laboratory experiments involving fracture–matrix flow as well as a tool for testing different transfer functions that have been suggested to use in reservoir simulators.     

Modelling of a Premixed Swirl-stabilized Flame Using a Turbulent Flame Speed Closure Model in LES

This paper proposes a combustion model based on a turbulent flame speed closure (TFC) technique for large eddy simulation (LES) of premixed flames. The model was originally developed for the RANS (Reynolds Averaged Navier Stokes equations) approach and was extended here to LES. The turbulent quantities needed for calculation of the turbulent flame speed are obtained at the sub grid level. This model was at first experienced via an test case and then applied to a typical industrial combustor with a swirl stabilized flame. The paper shows that the model is easy to apply and that the results are promising. Even typical frequencies of arising combustion instabilities can be captured. But, the use of compressible LES may also lead to unphysical pressure waves which have their origin in the numerical treatment of the boundary conditions.     

An Analytical Model for Determination of the Solvent Convective Dispersion Coefficient in the Vapor Extraction Heavy Oil Recovery Process

In this article, a new model is developed to determine the solvent convective dispersion coefficient in a solvent vapor extraction (VAPEX) heavy oil recovery process. It is assumed that solvent mass transfer by convective dispersion takes place along the transition zone between the solvent chamber and untouched heavy oil, whereas solvent mass transfer by molecular diffusion occurs in the direction normal to the transition zone. It is also assumed that the solvent-diluted heavy oil gravity drainage through the transition zone has a linear or quadratic velocity profile in order to obtain analytical solutions of the solvent convective dispersion coefficients for the solvent chamber spreading and falling phases. As a result, this analytical model correlates the solvent convective dispersion coefficient to the maximum apparent oil gravity drainage velocity at the interface between the solvent chamber and transition zone, solvent molecular diffusion coefficient, transition-zone thickness, and porosity of the porous medium. To determine the solvent convective dispersion coefficient, the maximum apparent oil gravity drainage velocity is calculated by using Darcy’s law and the transition-zone thickness is obtained either from a previous study or by using a time similarity between the solvent molecular diffusion and oil gravity drainage. It is found that such a determined solvent convective dispersion coefficient is two to five orders larger than the solvent molecular diffusion coefficient, depending on the detailed experimental conditions of a specific VAPEX test.     

Near-Wall Modification of an Explicit Algebraic Reynolds Stress Model Using Elliptic Blending

The elliptic blending approach is used in order to modify an Explicit Algebraic Reynolds Stress Model so as to reproduce the correct near wall behaviour of the turbulent stresses. The anisotropy stress tensor is expressed as a linear combination of tensor bases whose coefficients are sensitised to the non-local wall-blocking effect through the elliptic blending parameter γ. This parameter is obtained from a separate elliptic equation. The model does not use the distance from the wall thus it can be easily applied to complex geometries. It is validated against detailed DNS data for mean and turbulence quantities for the case of flow and heat transfer between parallel flat plates at three Reynolds numbers as well as against experimental data for the flow in a backward facing step at Re _H = 28,000. The comparison with DNS results or experiments is quite satisfactory and shows the validity of the approach.     

Modeling of Temperature and Strain-Rate Effects in Metals Using an Internal Variable Model

The thermo-viscoplastic behavior of three metals is characterized in a large range of loading conditions by using a new phenomenological constitutive model. The flow stress is decomposed into the sum of an effective stress with an internal stress depending upon an internal parameter which describes the strain hardening effect. The evolution of the internal stress is sensitive to the history of strain-rate and temperature. A systematic method is used for determining the model’s parameters. The model predictions show a good correlation with experimental data. Temperature history effects are especially analyzed.     

A Numerical Study of Instability Control for the Design of an Optimal Policy of Enhanced Oil Recovery by Tertiary Displacement Processes

In this paper, we consider the problem of control of hydrodynamic instability arising in the displacement processes during enhanced oil recovery by SP-flooding (Surfactant?CPolymer). In particular, we consider a flooding process involving displacement of a viscous fluid in porous media by a less viscous fluid containing polymer and surfactant over a finite length which in turn is displaced by a even less viscous fluid such as water. The maximum stabilization capacities of several monotonic and non-monotonic viscous profiles created by non-uniform polymer concentration are studied in the presence of interfacial tensions created by surfactants. The study has been carried out numerically to determine and characterize the most optimal viscous profiles of each family. Similarities in optimal monotonic viscous profiles of this constant-time injection policy and other injection policies by previous workers are noted. The presence of interfacial instability (due to viscosity jump) and layer instability (due to viscosity gradient) in appropriate proportions has been numerically demonstrated to be a necessary condition for monotonic as well as optimal non-monotonic profiles except in the limiting case of infinite time injection in which case maximum stabilization appears to result from pure layer instability. It has also been demonstrated numerically that the optimal non-monotonic viscous profiles can have better stabilization potential than the optimal monotonic profiles. Many other new features of this injection policy which have not been recognized before have been discussed.     

The PROSPER General Circulation Model,an adapted form of the Sandia Ocean Modelling method: A numerical description

The PROSPER General Circulation Model (PGCM) is a three-dimensional model based on the incompressible Navier-Stokes equations, an equation of state and the heat equation. The hydrostatic approximation and the rigid lid approximation are used. The system of equations is converted into an equivalent form in which the surface pressure is more directly expressed in terms of a two-dimensional Poisson equation. The finite difference method is described and analysed. In particular, the iteration method within every time step to determine the new surface pressure and velocity components, and numerical diffusion aspects due to the use of the staggered Arakawa-C grid are looked at. Since part of the development of the PGCM code is a result of studying the Sandia Ocean Modelling System (SOMS), a comparison is made with respect to the concepts used in both models.     

Identification of Material Damage Model Parameters: an Inverse Approach Using Digital Image Processing

A reliable prediction of ductile failure in metals is still a wide-open matter of research. Several models are available in the literature, ranging from empirical criteria, porosity-based models and continuum damage mechanics (CDM). One major issue is the accurate identification of parameters which describe material behavior. For some damage models, parameter identification is more or less straightforward, being possible to perform experiments for their evaluation. For the others, direct calibration from laboratory tests is not possible, so that the approach of inverse methods is required for a proper identification. In material model calibration, the inverse approach consists in a non-linear iterative fitting of a parameter-dependent load–displacement curve (coming from a FEM simulation) on the experimental specimen response. The test is usually a tensile test on a round-notched cylindrical bar. The present paper shows a novel inverse procedure aimed to estimate the material parameters of the Gurson–Tvergaard–Needleman (GTN) porosity-based plastic damage model by means of experimental data collected using image analysis. The use of digital image processing allows to substitute the load–displacement curve with other global quantities resulting from the measuring of specimen profile during loading. The advantage of this analysis is that more data are available for calibration thus allowing a greater level of confidence and accuracy in model parameter evaluation.     

The Onset of Convection in an Anisotropic Porous Layer Using a Thermal Non-Equilibrium Model

The stability of a horizontal fluid saturated anisotropic porous layer heated from below and cooled from above is examined analytically when the solid and fluid phases are not in local thermal equilibrium. Darcy model with anisotropic permeability is employed to describe the flow and a two-field model is used for energy equation each representing the solid and fluid phases separately. The linear stability theory is implemented to compute the critical Rayleigh number and the corresponding wavenumber for the onset of convective motion. The effect of thermal non-equilibrium and anisotropy in both mechanical and thermal properties of the porous medium on the onset of convection is discussed. Besides, asymptotic analysis for both very small and large values of the interphase heat transfer coefficient is also presented. An excellent agreement is found between the exact and asymptotic solutions. Some known results, which correspond to thermal equilibrium and isotropic porous medium, are recovered in limiting cases.     

Fracture Prediction for an Advanced High-Strength Steel Sheet Using the Fully Coupled Elastoplastic Damage Model with Stress-State Dependence

In this study,the numerical simulations of sheet metal forming processes are performed based on a fully coupled elastoplastic damage model.The effects of stress...