Next Article in Journal
Early Fault Detection of Wind Turbines Based on Operational Condition Clustering and Optimized Deep Belief Network Modeling
Next Article in Special Issue
An Experimental Investigation of Flow Regimes in Imbibition and Drainage Using a Microfluidic Platform
Previous Article in Journal
Monitoring Wind Turbine Gearbox with Echo State Network Modeling and Dynamic Threshold Using SCADA Vibration Data
Previous Article in Special Issue
Embedded Discrete Fracture Modeling as a Method to Upscale Permeability for Fractured Reservoirs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulation of Surfactant Oil Recovery Processes and the Role of Phase Behaviour Parameters

by
Pablo Druetta
*,†,‡ and
Francesco Picchioni
†,‡
Department of Chemical Engineering, ENTEG, University of Groningen, 9700 AB Groningen, The Netherlands
*
Author to whom correspondence should be addressed.
Current address: Nijenborgh 4, 9747 AG Groningen, The Netherlands.
These authors contributed equally to this work.
Energies 2019, 12(6), 983; https://doi.org/10.3390/en12060983
Submission received: 28 January 2019 / Revised: 5 March 2019 / Accepted: 8 March 2019 / Published: 13 March 2019
(This article belongs to the Special Issue Enhanced Oil Recovery 2019)

Abstract

:
Chemical Enhanced Oil Recovery (cEOR) processes comprise a number of techniques which modify the rock/fluid properties in order to mobilize the remaining oil. Among these, surfactant flooding is one of the most used and well-known processes; it is mainly used to decrease the interfacial energy between the phases and thus lowering the residual oil saturation. A novel two-dimensional flooding simulator is presented for a four-component (water, petroleum, surfactant, salt), two-phase (aqueous, oleous) model in porous media. The system is then solved using a second-order finite difference method with the IMPEC (IMplicit Pressure and Explicit Concentration) scheme. The oil recovery efficiency evidenced a strong dependency on the chemical component properties and its phase behaviour. In order to accurately model the latter, the simulator uses and improves a simplified ternary diagram, introducing the dependence of the partition coefficient on the salt concentration. Results showed that the surfactant partitioning between the phases is the most important parameter during the EOR process. Moreover, the presence of salt affects this partitioning coefficient, modifying considerably the sweeping efficiency. Therefore, the control of the salinity in the injection water is deemed fundamental for the success of EOR operations with surfactants.

1. Introduction

Oil is and has been for the last century the main energy source, and the economy on a global scale still depends largely on it [1,2,3,4]. After the traditional exploitation techniques, which include natural driven mechanisms (primary recovery) and waterflooding (secondary recovery), up to 50% of the Original Oil in Place (OOIP) might be extracted from the oilfield [3,5,6,7]. Additional production techniques, or Enhanced Oil Recovery processes, target the remaining oil still trapped [8]. The most commonly used EOR methods are [6]: thermal (fire flooding—continuous/cyclic steam injection); chemical (polymers, surfactants and/or caustic); miscible (CO2 or miscible solvent); and others (e.g., microbial EOR). The cEOR techniques are one of the most used for low and medium viscosity crude oils and they can be employed in a wide range of rock formations. Another important advantage of cEOR processes is that the chemicals to be used as sweeping agents can be synthesized or improved based on the desired characteristics for determined crude oils, rock formations and porous media conditions in general, optimizing the overall process and increasing the operational life of the reservoir. The objective in this paper is to present the numerical simulation of surfactant EOR flooding. Thus, a novel numerical simulator is presented to analyze tertiary recovery processes in a 2D oilfield [9,10,11], simulating a two-phase, four-component flow.

1.1. Surfactant Flooding

Surfactant flooding for EOR is not new, and it has been used for more than 40 years [7,12,13,14,15,16,17,18,19]. Surfactants adsorb onto the oil/water interface and thus reduce the interfacial tension (IFT) and capillary forces, which are the main reason for the trapped oil in pores (capillary trapping) [20]. Hence, what makes surfactants useful in EOR is that they can be soluble in both organic compounds and water. The objective in designing surfactants for EOR is to achieve low interfacial tension at low concentrations, and acceptable adsorption levels on the rock formation [21,22,23]. A multicomponent, multiphase system is present in porous media and several authors have already studied how to numerically simulate this kind of flow for both compressible and incompressible systems [24,25,26,27,28,29,30,31,32,33]. This has led to the development of academic and commercial simulators like UTCHEM (University of Texas at Austin), CMG (Stars) and Eclipse (Schlumberger), as well as several others that can be found in the literature [34,35,36,37,38,39,40,41,42,43,44]. Nonetheless, most of the commercial simulators are based on complex models that require a strong knowledge of both the chemicals involved and their interactions with the rock formation. For instance, the phase behaviour in UTCHEM is based on the Hand’s rule, a representation of the ternary phase diagram [45], expressed as a function of the effective salinity. The model used reduces the components present in porous media to a pseudo-ternary component mixture [41,42,46,47,48,49]. UTCHEM is capable to model Winsor Type II(−) (oil-in-water), Type II(+) (water-in-oil), and Type III (three phase microemulsion) systems [50].
This renders a complex system of non-linear partial differential equations which might create numerical stability problems if these parameters are not thoroughly evaluated. On the other hand, simulators reported in the literature [34,35,36] do not perform a comprehensive study of component partitioning and phase behaviour since they are based in an insufficient amount of parameters describing this behaviour. This can be the main source of incorrect results when developing new chemicals for EOR [9]. Moreover, some of these simulators [37,38] were developed for one-dimensional systems (e.g., coreflooding), restricting their application to specific and determined cases, since they do not consider the areal sweeping efficiency. This new model considers a fourth component which represents the salt content in the porous media. Besides, the compressibility of the rock formation is taken into account in the flow and mass equations. A simplified ternary phase diagram is used to determine the value of the volumetric concentration relationships, necessary to solve the compositional model. Numerically speaking, the new model uses a second order discretization scheme with a flux limiter function in order to minimize numerical diffusion and dispersion phenomena. The two-phase system was chosen since this presents greater difficulties in order to determine the component partitioning. In three-phase systems (water-oil-microemulsion), the treatment is relatively straightforward and simulators use the same approach: aqueous and oleous phases are considered pure, and the composition of the microemulsion is assumed to be constant and function only of salinity and temperature. Due to these factors, the simulation of three-phase systems in fields with high salinity or thermal gradients may cause numerical issues (e.g., appearance of spurious oscillations, numerical instabilities), rendering the method ineffective [9,11]. In contrast, in two-phase systems, the components’ volume concentrations in the phases respond to a more complex model. In this case the phase behaviour plays an important role in the amount of oil recovered. It is for this reason that it is considered important to study and modeling two-phase systems. Finally, the influence of the parameters that may alter the model behaviour is also taken into account when choosing the partitioning factors (e.g., the salinity expressed as a function of the concentration of monovalent or divalent cations, temperature) [51].

1.2. Aim of this Work

The improvement of this study resides in the development of a novel simulator for a two-dimensional oilfield, capable of simulating the flow of a two-phase, four-component system for cEOR. This simulator is the improvement of a previous 2D, three-component simulator for surfactant flooding, based on previously reported research [9,37,38,52,53]. The component partitioning is modeled in an accurate, although relatively simple and robust way, solving the inconvenient described in previous simulators with respect to the phase behaviour. The former is assumed to be a function of the total dissolved solids (TDS) present in the water phase, using a piece-wise function based on previous literature. The phase behaviour modeling was done increasing the number of factors used to describe the chemical behaviour, but avoiding creating a system of non-linear partial differential equations which would complicate the simulation process. This leads to a new set of optimum design parameters which can be used during the synthesis of future surfactants.
The combination of these concepts has resulted in a novel simulator, which can be used for the design and screening of new surface active agents to be used in EOR. Nevertheless, the 2-dimensional modelling does not allow to consider the volumetric sweeping efficiency in EOR agents and how the difference in the densities between crude oil and the EOR agents could affect the recovery efficiency, specifically when the vertical permeability cannot be neglected with respect to horizontal ones. Furthermore, the influence of the temperature and energy transfer was not considered in the present model, which mainly affects the phases’ rheological behaviour as well as the influence of the salinity in the surfactant flooding. Besides, the model proposed is limited to two phases, aqueous and oleous, neglecting the presence of the gas phase, which is usually present in many reservoir (e.g., gas cap). The mentioned phenomena should be the objective of further developments, tackling out one of its current disadvantages with respect with commercial reservoir simulators. Simulation of multiphase and multicomponent flow in porous media requires solving a number of coupled, highly non-linear equations dealing with temporal and spatial derivatives of pressure and mass concentrations. It is adopted in this model a second-order accuracy discretization scheme, along with a Total Variation Diminishing (TVD) flux limiter, which reduces the occurrence of numerical diffusion and dispersion mechanisms [11]. The compositional flow is a suitable approach to study the different cEOR techniques, which can be described as the mass transfer of a number of components (e.g., polymers, surfactants, salts, etc.) in a two- or three-phase systems. Generally speaking, commercial simulators present more advantages that academic ones. However, the objective of this paper is to present a novel cEOR simulator. The second-order accuracy in the discretization allows obtaining better results than standard academic simulators and it is comparable to commercial ones. Chemically speaking, the adaptation of the ternary diagram to model the phase behaviour (and the influence of the salinity in the calculation of the partition coefficient) and the phases’ rheological behaviour (e.g., the calculation of aqueous viscosity in a step-wise approach considering the influence of all the components) were never considered in academic simulators and some of them are not considered by commercial ones. Thus, mathematically speaking it represents a novelty with respect to academic simulators and chemically speaking, it represents an advance when compared to academic simulators and presents some minor advantages when compared to commercial simulators. A major disadvantage with respect to these is the number of dimensions analyzed, although it is considered that a 2D approach is appropriate for a first version of the simulator in order to study the behaviour of chemical agents in EOR. Moreover, this approach allows to study the chemical performance in different oilfields in order to determine a set of desired features when synthesizing new products.

2. Methods

2.1. Physical Model

The physical model represents a reservoir ( Ω ) of known geometric dimensions, which has an absolute permeability (K) and porosity ( ϕ ). These can be constant or represented by a normally distributed fields (Figure 1). The flow is considered isothermal, Newtonian, incompressible and 2-dimensional (it is assumed as in previous models that the vertical permeability is negligible when compared to horizontal permeabilities); it is also considered that the system is in local thermodynamic (phase) equilibrium. The Darcy’s law is valid and gravitational forces are neglected [38].
Surfactant EOR flooding, such as other chemical techniques, involves the flow of fluids in two phases (aqueous and oleous), and several components (water, salt, surfactant and petroleum). It is noteworthy that these components can be mixtures of a number of pure components, since for instance, petroleum is a mixture of many hydrocarbons, water may contain dissolved monovalent and divalent salts, and the surfactant is composed by a number of different molecules (e.g., cosurfactants) [3,8,16]. The recovery process consists in injecting during a first stage an aqueous solution with the surfactant, followed by a water bank (water or brine) in order to drive the chemical slug. This model is represented by a system of strongly non-linear differential equations, complemented by a set of algebraic relationships representing physical properties of the fluid and the rock, which are: component partitioning as a function of the salinity, interfacial tension, residual phase saturations, relative permeabilities, rock wettability, phase viscosities, capillary pressure, adsorption on the formation, and dispersion. It is clear that most of these properties are strongly dependent on the surfactant concentration, which alters the flowability of the aqueous phase. The interfacial tension depends on the former and it affects the residual saturations, which are dependent on the capillary number, a dimensionless parameter based on the IFT and the viscosity.

2.2. Mathematical Model

The compositional approach offers the versatility to model a number of components in porous media affecting the properties of the phases involved. In a previous paper three different components were considered, whilst in this new model the salt is included as a fourth component. This does not modify the equations, but only the number of auxiliary relationships to determine the system. The mathematical relationship, defined by N c o m p · N p h a s e s 1 , can be derived from the compositional approach in order to determine numerically the system.
The numerical simulator presented in this paper is based on previous ones developed also for EOR flooding [38,53]. The 2D model was validated against academic and commercial simulators in a previous paper, discussing also its order of convergence (Figure 2) [9,11]. Since this model it is an extension of a previous one, it is considered that the validation process previously done is valid for this model [54,55,56]. Thus, the Darcy, mass conservation and aqueous pressure equations describing the system are presented in Equations (1)–(3),
u j = K ̲ ̲ · k r j μ j · p j ; j = o , a
ϕ z i t + · j V i j · u j · j D i j ̲ ̲ · · V i j = ϕ A d i t + q i ; i = p , c , w , s
ϕ c r p a t + · ( λ · p a ) = t ϕ · i A d i · λ o · p c + q t

2.2.1. Discretization of the Partial Differential Equations

The differential form of the Darcy’s and mass transport renders a system of non-linear parabolic partial differential equations, discretized using a Finite Difference Method scheme. The aqueous phase pressure is implicitly solved using a centered scheme and a second-order Taylor approximation for the time derivatives. Moreover, the total and aqueous Darcy velocities are also discretized and explicitly solved using a centered difference stencil. Therefore, Equations (1) and (3) are discretized as follows,
u m , n < n + 1 > , [ k + 1 ] = λ x , m , n [ k ] 2 · Δ x · p m + 1 , n a p m 1 , n a [ k + 1 ] λ x , m , n o , [ k ] 2 · Δ x · p c , m + 1 , n a p c , m 1 , n a [ k + 1 ] i n t 1 1 < n + 1 > · i ^ + λ y , m , n [ k ] 2 · Δ y · p m , n + 1 a p m , n 1 a [ k + 1 ] λ y , m , n o , [ k ] 2 · Δ y · p c , m , n + 1 a p c , m , n 1 a [ k + 1 ] i n t 1 1 < n + 1 > · j ^
u m , n a , < n + 1 > , [ k + 1 ] = λ x , m , n a , [ k ] 2 · Δ x · p m + 1 , n a p m 1 , n a [ k + 1 ] < n + 1 > · i ^ + + λ y , m , n a , [ k ] 2 · Δ y · p m , n + 1 a p m , n 1 a [ k + 1 ] < n + 1 > · j ^
c r ϕ + Δ t 2 ϕ t m , n < n + 1 > , [ k ] p m , n a , < n + 1 > p m , n a , < n > Δ t [ k + 1 ] + λ x , m + 1 2 , n < n + 1 > , [ k ] Δ x 2 · p m + 1 , n a p m , n a < n + 1 > , [ k + 1 ] λ x , m 1 2 , n < n + 1 > , [ k ] Δ x 2 · p m , n a p m 1 , n a < n + 1 > , [ k + 1 ] + λ y , m + 1 2 , n < n + 1 > , [ k ] Δ y 2 · p m , n + 1 a p m , n a < n + 1 > , [ k + 1 ] λ y , m 1 2 , n < n + 1 > , [ k ] Δ y 2 · p m , n a p m , n 1 a < n + 1 > , [ k + 1 ] = ϕ + Δ t 2 ϕ t m , n < n + 1 > , [ k ] A d m , n < n + 1 > A d m , n < n > Δ t [ k + 1 ] + + A d + Δ t 2 A d t m , n < n + 1 > , [ k ] ϕ m , n < n + 1 > ϕ m , n < n > Δ t [ k + 1 ] + λ x , m + 1 2 , n o , < n + 1 > , [ k ] Δ x 2 · p m + 1 , n c p m , n c < n + 1 > , [ k + 1 ] λ x , m 1 2 , n o , < n + 1 > , [ k ] Δ x 2 · p m , n c p m 1 , n c < n + 1 > , [ k + 1 ] + λ y , m + 1 2 , n o , < n + 1 > , [ k ] Δ y 2 · p m , n + 1 c p m , n c < n + 1 > , [ k + 1 ] λ y , m 1 2 , n o , < n + 1 > , [ k ] Δ y 2 · p m , n c p m , n 1 c < n + 1 > , [ k + 1 ] + q m , n t , < n + 1 > , [ k ]
Finally, the discretization of the mass conservation is presented. Equation (2) is the advection-diffusion PDE used in hydraulic, mass transport and flow processes in porous media. The advective terms have an hyperbolic nature and the traditional first-order Upwind schemes cause numerical diffusion. In this model a complete second-order explicit discretization scheme is employed based on flux limiting techniques. This increases the numerical accuracy of the simulator and decreases the influence of numerical phenomena on the results [11,25,57,58,59]. Thus, the discretized mass conservation equation is,
C 1 Δ t z i < n + 1 > = C 2 z i < n > + C 3 Δ x · j F L I M , x j , < n + 1 > , [ k + 1 ] u x , m , n j , [ k + 1 ] · V i , m , n j , [ k ] u x , m 1 , n j , [ k + 1 ] · V i , m 1 , n j , [ k ] < n + 1 > + + C 3 Δ y · j F L I M , y j , < n + 1 > , [ k + 1 ] u y , m , n j , [ k + 1 ] · V i , m , n j , [ k ] u y , m , n 1 j , [ k + 1 ] · V i , m , n 1 j , [ k ] < n + 1 > + + 1 Δ x 2 · j S j ϕ d m i j m + 1 2 , n · V i , m + 1 , n j V i , m , n j S j ϕ d m i j m 1 2 , n · V i , m , n j V i , m 1 , n j < n + 1 > , [ k ] + + 1 Δ y 2 · j S j ϕ d m i j m , n + 1 2 · V i , m , n + 1 j V i , m , n j S j ϕ d m i j m , n 1 2 · V i , m , n j V i , m , n 1 j < n + 1 > , [ k ] 1 Δ t ϕ + Δ t ϕ t m , n < n + 1 > , [ k + 1 ] · A d i < n + 1 > A d i < n > m , n [ k ] A d i , m , n < n + 1 > , [ k ] Δ t · ϕ < n + 1 > ϕ < n > m , n [ k ] + + q i , m , n < n + 1 > , [ k + 1 ] + j u x , m , n j , [ k + 1 ] Δ t 2 ϕ m , n < n + 1 > , [ k + 1 ] · A d x ϕ t + A d 2 ϕ t x + ϕ x A d t i , m , n < n + 1 > , [ k ] + + j u y , m , n j , [ k + 1 ] Δ t 2 ϕ m , n < n + 1 > , [ k + 1 ] · A d y ϕ t + A d 2 ϕ t y + ϕ y A d t i , m , n < n + 1 > , [ k ] j u x , m , n j , [ k + 1 ] Δ t 2 2 A d t x + u y , m , n j , [ k + 1 ] Δ t 2 2 A d t y i , m , n < n + 1 > , [ k ] + + j Δ t 2 ϕ m , n < n + 1 > , [ k + 1 ] u x , m , n j , [ k + 1 ] 2 2 V i , m , n j , [ k ] x 2 + u y , m , n j , [ k + 1 ] 2 2 V i , m , n j , [ k ] y 2 + 2 u x , m , n j , [ k + 1 ] u y , m , n j , [ k + 1 ] 2 V i , m , n j , [ k ] x y < n + 1 >
C 1 = ϕ m , n + Δ t ϕ t u t x , m , n Δ t 2 ϕ m , n ϕ x u t y , m , n Δ t 2 ϕ m , n ϕ y < n + 1 > , [ k + 1 ] C 2 = ϕ m , n Δ t u t x , m , n 2 ϕ m , n ϕ x u t y , m , n 2 ϕ m , n ϕ y + u t x , m , n Δ t 2 ϕ m , n 2 ϕ t x + u t y , m , n Δ t 2 ϕ m , n 2 ϕ t y < n + 1 > , [ k + 1 ] C 3 = 1 Δ t 2 ϕ m , n ϕ t < n + 1 > , [ k + 1 ]

2.2.2. Boundary Conditions

At the beginning of the EOR operations, the oil saturation in the oilfield is assumed to be equal to values after a primary recovery, or the saturation after a waterflooding which reached a threshold fractional flow at the producing well, rendering the exploitation no longer profitable. The pressure is assumed to be constant and moreover there is no chemical agents present.
t = 0 ; ( x , y ) Ω : z c = 0 ; z p = S o r H ; p a = p i
The EOR process starts with the injection of the chemical agent (surfactant) with constant concentration. After this period, a water-bank follows in order to sweep the remaining oil. The boundary condition for the oilfield consists in a ‘no flow’ on the contour ( Γ ), since it is assumed that the porous formation is surrounded by an impermeable rock layer. Regarding the advection, this boundary condition is satisfied by imposing the mobilities to be zero.
Injecting   Well 0 t t i n : z c = z i n t > t i n : z i n = z w , z c = 0
Boundaries λ m , n j = 0 z i n ˇ Γ = 0 ; i = p , c ; t ( m , n ) Γ

2.3. Physical Properties

2.3.1. Chemical Component Partition

The most important part of a numerical simulation using a compositional model is to understand how the components distribute into the phase, what it is called phase behaviour of the system. A surfactant flooding can be reasonably well modeled in a ternary phase diagram, wherein the chemical component is located in the apex, and water and petroleum occupy the lower vertices. The composition of a mixture is then determined by any point inside the triangle, which is known with two of the total concentrations of the components [33,37,60]. The numerical simulation of this model comprises a two-phase, four pseudo-components system, which can be modeled using this ternary phase diagram. If the concentration of the surfactant increases, oil and water become miscible, hence the triangle can be divided into two zones: a miscible in the upper part and the immiscible in the bottom. The curves delimiting these regions are determined by volumetric concentration ratios calculated as,
Solubilization Coefficient = L p c a = V p a V c a
Swelling Coefficient = L w c o = V w o V c o
Partition Coefficient = k c = V c o V c a
Depending on the value of k c , two different behaviours can be observed: Type II(−) (for k c < 1 ), and Type II(+) systems (for k c > 1 ). These parameters affect greatly the efficiency of the EOR process. They should be correctly measured to predict recovery efficiencies with existing surfactants, or be the cardinal values when the objective is to design new chemical products. The partition coefficient depends on the composition of the surfactant as well as on the water characteristics, such as temperature and salinity. Sheng [51] described this two-phase approximation to surfactant flooding and the dependence on the salinity of the partition coefficient. He presented an empirical, piece-wise function (Equation (15) and Figure 3) to describe this, which also depends on the optimum volumetric salt concentration in the aqueous phase ( V s , o p t a ) [49],
k c = 10 2 V s a V s , o p t a 1 if V s a > V s , o p t a 10 2 1 V s a V s , o p t a if V s a < V s , o p t a
In this case is assumed that all the salinity present in the reservoir is composed by monovalent cations and that the salt component is only present in the aqueous phase ( V s o = 0 ). With these extra relationships the physical system becomes, mathematically speaking, determined with a unique solution, and the parameters previously introduced can now be calculated for each representative elementary volume (REV). Figure 4 shows the original and simplified ternary diagrams used for this simulation.

2.3.2. Interfacial Tension

The IFT of the system depends on the presence and the concentration of the surfactant as well as the emulsion type present in the reservoir. In this simulator, a simplified IFT correlation is used [32,37,38,52,53]. For Type II(−) systems (oil/water emulsion):
log ( σ ) = log ( F ) + ( 1 L p c a ) · log ( σ H ) + G 1 1 + G 2 · L p c a ; L p c a < 1 log ( σ ) = log ( F ) + G 1 1 + L p c a · G 2 ; L p c a 1
For Type II(+) systems (oil emulsion/water):
log ( σ ) = log ( F ) + ( 1 L w c o ) · log ( σ H ) + G 1 1 + G 2 · L w c o ; L w c o < 1 log ( σ ) = log ( F ) + G 1 1 + L w c o · G 2 ; L w c o 1
where G 1 and G 2 are input parameters, while the term F is the correction factor introduced by Hirasaki [32],
F = 1 e i = p , w , c ( V i o V i a ) 2 1 e 2
In cEOR processes, the surfactant decreases the interfacial energy, allowing the mobilization of the oil trapped in the reservoir. Thus, the residual saturations of the phases depend on the IFT (Figure 5). The IFT of the water-oil system (with no surfactant present) is assumed to be constant throughout the simulation.

2.3.3. Residual Saturation

The residual saturations represent an important value in conventional and enhanced oil recovery processes. They establish a limit to how much oil can be mobilized by the sweeping agents. If these values can be reduced, this will increase the efficiency of the whole process. As explained before, they depend on the IFT in the water-oil system. The presence of the surfactant in an EOR flooding can modify the residual saturations in porous media. This relationship is ruled by a dimensionless group, the capillary number,
N v c = u · K λ · σ
The functionality between this group and the residual saturation is described by the following model [38]:
S j r S j r H = 1 if N v c < 10 1 T 1 j T 2 j T 1 j · [ log ( N v c ) + T 2 j ] if 10 1 T 1 j T 2 j N v c 10 T 2 j 0 if N v c > 10 T 2 j
The piecewise function is defined by a set of constant parameters which depend on the fluids and the rock formation. The relationship between the residual saturation after EOR and waterflooding processes is known as normalized residual saturation of phase j. The Equation (20) determines what is known as capillary desaturation curves (Figure 6). At low capillary numbers, the behaviour is similar to a secondary recovery and the normalized residual saturation does not decrease. As the IFT decreases and/or the viscosity increases, the capillary number raises to values several orders of magnitude higher than those of secondary recovery. It is for this reason that in areas of high speeds (i.e., nearby the wells) oil saturations lower than those attained in a waterflooding can be achieved. As can be seen, the aqueous phase requires much higher values of N v c to achieve a full desaturation [8].

2.3.4. Relative Permeabilities

The relative permeabilities influence the phase velocities in the Darcy’s equation, and thus the efficiency of oil recovery depends greatly on them. These are functions of the residual saturations calculated in the previous point. The model adopted to calculate these permeabilities is taken from Camilleri [41,42], which is generally used for most chemical flooding processes. Knowing beforehand the phase saturations, the relative permeabilities are calculated according to,
k r j = k r j 0 · S j S j r 1 S j r S j r e j ; j = o , a ; j j
where k r j 0 and e j represent the end point and the curvature of the function k r j ( S j ) . These values are calculated by the following equations:
k r j 0 = k r j 0 H + ( 1 k r j 0 H ) · 1 S j r S j r H ; j = o , a ; j j
e j = e j H + ( 1 e j H ) · 1 S j r S j r H ; j = o , a ; j j
where k r j 0 H and e j H are the endpoint values of curvature and relative permeability function system for water-oil without the presence of chemical agents, respectively.

2.3.5. Phase Viscosities

The viscosity of each phase depends on its composition as a function of the volumetric concentration of each component, according to the following function [38,41,42]:
μ j = V w j · μ a H · e α 1 · ( V p j + V c j ) + V p j · μ o H · e α 1 · ( V w j + V c j ) + V c j · α 3 · e α 2 · ( V w j + V o j ) ; j = o , a
where α k are constants, and μ a H and μ o H are the viscosities in the water-oil system without surfactant, respectively. Nevertheless, in this work a modification is proposed to the formulation presented in simulators (e.g., UTCHEM) and the model developed by Porcelli [37], introducing the dependence on the salt in the pure water/brine viscosity and replacing the latter in Equation (24), according to the following expression [61],
μ b r i n e = μ a · 1 + A s a l V s a + B s a l V s a 2
where A s a l and B s a l are constants based on rheology experiments. For the oil phase, as did before, a Newtonian behaviour is adopted. According to the literature [62], light and medium oil cuts exhibit Newtonian behaviour while heavy oil might present a slight shear-thinning region.

3. Results and Discussion

3.1. Introduction

The tertiary recovery process starts with an initial period in which an aqueous surfactant solution is injected, followed by a water bank used to sweep the surfactant slug and oil bank towards the producing wells. The chemical concentration profile has a shape similar to a wave propagation in a two-dimensional field, travelling at a speed equal to Darcy velocities in each block. For the first part of simulations, the influence of the salt will be neglected in order to study the phase behaviour (three-component system). Subsequently, a set of simulations is analyzed comprising the influence of the fourth component (salt) on the oil recovery mechanisms during surfactant flooding, establishing guidelines for the synthesis of new surfactants.
The initial conditions are set to consider that the surfactant EOR process starts after a waterflooding (or secondary recovery). This is aimed at studying how the surfactant can mobilize the remaining oil which is trapped and can no longer be swept by conventional recovery techniques. Thus, this initial oil saturation allows to study how the surfactant modifies the oil mobility.

Data

The simulation conditions and physical properties are established aimed at emulating an existing oilfield apt for a EOR recovery process with surfactants (Table 1 and Table 2). These parameters were used during the simulator validation process as well as in previous papers in which different oil recovery processes were analyzed [9,11,37,38].

3.2. Three-Component System

This part of the study was done considering the salt is not present in the system, with no influence on the partitioning component and phase behaviour. The goal is to determine the influence of these parameters on the oil recovery process. Subsequently, this fourth component was added to the system in order to determine its influence on the abovementioned parameters. The latter will be used as reference case for the analysis when salt is considered in the mathematical model.

3.2.1. Influence of the Component Partitioning Parameters

The first part of the analysis, and the one deemed as the most important in surfactant flooding, is the study of the influence of the phase behaviour parameters on the three-component system. This importance has already been discussed previously [33,41,42,44,52,63,64,65,66,67,68], and it plays an important role in the recovery process, because it determines the fractions of chemical, oil and water in each of the phases. It is noteworthy that, while these parameters are dependent on the salt content present in porous media, in this section this dependence will not be considered since this component is not taken into account. The objective in following this approach is to study firstly the influence of these parameters and determine, if possible, a trend and relationship between the values of these and the oil recovery. Subsequently, the salt is included as fourth component and also its functionality with the partition coefficient.
The partition coefficient was considered by several authors as the most important parameter to study the phase behaviour and hence it was chosen as the first to be evaluated. This parameter determines the type of emulsion present: lower values ( k c < 1 ) render oil-in-water emulsions, known as Type II(−), whilst for higher values it renders water-in-oil emulsions, also known as Type II(+). Several recovery processes were simulated keeping constant the other two parameters ( L p c a = L w c o = 1 ), obtaining the results presented in Table 3 and Figure 7.
Table 3 and Figure 7 show a clear relationship between oil recovered and the partition coefficient, a trend that agrees with the results published by Porcelli [52] and Bidner [38]. Numerical simulations were stopped after the chemical flowrate decrease to negligible values, which virtually stops any further oil recovery. For partition coefficients larger than the unity ( k c > 1 ), complete desaturation profiles are almost achieved in the area around the injection well, which also coincides with the results of the 1D model published by Porcelli (Figure 8 and Figure 9). Due to increased sweeping efficiency, the time employed for the chemical breakthrough to take place is also proportional to the partition coefficient. As more chemical component goes into the oleous phase, the chemical bank is less dispersed and it rises above its injecting concentration [38]. This was also observed during the simulations, with maximum concentration values of surfactant (obtained at the injection node) proportional to k c . The oleous phase swells up with water and chemical, and a lower amount of oil is then a left behind after the chemical slug, leading to higher oil recovery efficiencies is as k c increases.
In turn, Figure 7c shows the pressure drop observed during the whole simulation. Firstly, it is important to note that the maximum values obtained with surfactant are significantly lower than those with normal polymer EOR techniques due to its rheological characteristics. In addition, at the end of flooding process, there is no difference between waterflooding and surfactant EOR flooding because in this the case the smaller size of its molecules do not cause the disproportionate permeability reduction (DPR) phenomenon, commonly found with polymers. Furthermore, there is no significant difference between the different partition coefficients with respect to the final pressure. This is due to the fact that, when no chemical is present, the oil relative permeability becomes negligible and all the pressure drop comes from the aqueous phase, which is nearly at the saturation when its relative permeability reaches its maximum. This renders that pressure drops are approximately the same at the end of the flooding. With respect to the time evolution, there is an initial decrease during the injection period related to the fact that surfactant in the injection well region starts mobilizing the trapped oil. Afterwards, as the oil bank is formed and displaced towards the production well, there is a gradual increase of the pressure drop. When the oil bank (along with the surfactant) finally reaches the production well and the total saturation in the reservoir decreases, there is a sudden decrease in the pressure drop. The minimum in the latter is also a function of the partition coefficient: as the k c is increased, the oil recovery is improved and therefore the oil saturation in the reservoir decreases, thus the pressure difference between injector and producer is also diminished. Finally, where the surfactant concentration is lower than the critical micelle concentration (CMC), the capillary number decreases due to an increase in the IFT, and thus the relative permeability values become similar to those at the beginning of the EOR process. Moreover, the pressure drop returns to the initial values. Figure 7b presents the flows for each component as a function of time. As mentioned above, the chemical bank suffers a delaying as k c increases due to increased sweep efficiency. It is also noteworthy in this graph that the total injected chemical differs in all the cases. When the simulations were performed, the injection time was selected so as to have approximately the same maximum concentration of the chemical slug in the injection well region in the reservoir. Nonetheless, this trend was verified by simulating these scenarios injecting the same total amount of chemical and the results, although differ to those presented in this paper, show the same behaviour: higher partition coefficients yield higher recovery efficiencies. Figure 8 and Figure 9 depict the oil saturation profile, surfactant concentration, IFT and pressure during simulations for different partition coefficient cases.
Figure 8 and Figure 9 (c and d) show the oil saturation profiles during the simulations carried out for two different cases of surfactant emulsions. In the case of oil-in-water emulsions, the full desaturation region in porous media has reached a greater reservoir area. In the top left (a) of these Figures, the surfactant concentration profiles are plotted for the abovementioned cases of k c . These show what it was described above: for cases of k c > 1 , the maximum values are higher and the surfactant bank is less dispersed than in the case with k c < 1 . The rheology behaviour of the emulsion is considered independent of shear rate, but only a function of the concentrations of the different components. Therefore, the viscosity does not vary significantly due to partition coefficient, nor it reached similar values to those obtained in polymer flooding processes. It is also due to this reason that surfactants must be used in petroleum reservoirs with low/medium viscosity, so as not to diminish the sweep efficiency because of increased mobility ratios. Finally, top right (b) plots depict the IFT and pressure profiles, respectively. In the absence of surfactant, the system yields the water-oil IFT, while ultra-low values are achieved in the chemical bank area, increasing the Capillary number and reducing the residual oil saturation. The sweeping efficiency in the vertexes with no wells is reduced due to the negligible pressure drop in that region, which renders the lowest velocities in porous media.
Following the study of the influence of the partition coefficient, the attention is set on the other two parameters: the solubilization and swelling coefficients. The partial miscibility of both aqueous and oleous phases due to the presence of the surfactant is quantified by these parameters. The former is generated either by the solubilization of the petroleum component into the aqueous phase or by the swelling of the oleous phase with water and chemical components [38,52]. Taking into consideration this goal, the following simulations were performed for each of the previous partition coefficients (Table 4), varying the remaining two parameters in order to estimate their degree of influence on the oil recovery process (Figure 10).
The summary of influence of the phase behaviour parameters is presented in Table 4. As mentioned previously, the partition coefficient is the most important factor in terms of recovery efficiency. This determines the type of emulsion (water-in-oil or oil-in-water), affecting the phases’ flow behavior and the amount of oil that is mobilized by the presence of the surfactant. The solubilization and swelling coefficients also affect the amount of oil recovered, though to a lesser extent. From these two, the results showed a bigger sensitivity to the solubilization coefficient, whilst the importance of the swelling coefficient was more appreciable in Type II(+) emulsions.
The results of the sensitivity to the swelling and solubilization parameters show similar trends to those reported by Bidner [38]. Both define the area and shape of the interior triangle (two-phase zone) in the ternary diagrams (Figure 4). The position of the apex of the triangle, its area, and the distance between the injection point and the binodal interior triangle along the extended tie-line which passes through the injection determine the sweeping effectiveness. Furthermore, the values of these parameters affect the IFT: low values of L p c a and L w c o increase the IFT, which reduces the capillary number and therefore increases the relative desaturation of the system.
As in the study of the partition coefficient, the oil recovery resulted proportional to the value of these parameters (Table 4 and Figure 10). Nevertheless, in the case of the swelling coefficient, this influence was less significant. In the case of Type II(−) emulsions, the values of oil recovery in the case of L w c o were substantially similar, and showed greater sensitivity in Type II(+) emulsions due to the reasons mentioned above in this paragraph (Figure 11 and Figure 12). In the case of the solubilization coefficient ( L p c a ), the values obtained showed greater sensitivity. This is because the distance between the injection point and the binodal interior triangle along the extended tie-line, which passes through the injection, varies with the values of L p c a , which does not occur with L w c o [38]. This decreases the influence of the latter in the recovery process.
As conclusion of this analysis, the trends and sensitivities shown were confirmed for the three parameters defining the phase behaviour. The values found in the 2D model are consistent with those reported [9,38,53,69]. As mentioned in the introduction, the phase behaviour plays a fundamental role in surfactant flooding and therefore the parameters must be carefully studied and chosen according to the chemicals being used. However, it is important to mention that these values depend on the salt content (mono- and divalent) present in the reservoir. This was not considered in previous simulations using this simplified model.
All in all, based on previous results [9,38,69], it is considered that the main difference between them and this simulator is the treatment of the last two parameters: the simplified model treats them as constants whilst it is evident from Figure 4 that their values change as function of the concentration. This is because in the model presented there is a direct proportionality between the recovery efficiency and these coefficients. In order to match the results of the simulations with the observed values, it is necessary to modify this and establish the functionality with the concentration and the partition coefficient. Then a set of functions would be obtained defining all the phase behaviour parameters based on the chemical species’ properties, thus they can be adjusted in order to achieve a matching between the simulation results and the laboratory tests. This is deemed a critical point in the model to be addressed in future research, in order to achieve a full dependability of the parameters on the system conditions. The partition coefficient will take the form presented by Sheng [51], and the solubilization and swelling coefficients will be represented by functions of the concentrations and the partition coefficient. Such a system would allow, for example, obtain maximum oil recovery values for salt content near the critical value, as it was published by several authors [49].

3.2.2. Influence of the Partition Coefficient in Random Permeability Fields

The second part of the study verifies and extends the previous simulations to a case with a heterogeneous absolute permeability field, as found in real oil reservoirs. Based on a previous work [11], a Gaussian random permeability field was adopted to model the heterogeneities in the oilfield. However, it is well known that permeability field usually follow a log-normal distribution. The objective of this simulator was not focused on the geostatistical reservoir-modeling but on analyze the influence of heterogeneities on the EOR process. The upgrade and modification of the permeability field distribution to account for more real oilfields will be considered in future studies. The influence of this in the performance of the EOR agent can be studied using different distribution laws, as presented in this paper. For this study it was generated an artificial permeability field according to a Gaussian random distribution (Figure 13).
The operating conditions and physical properties of the phases involved were kept constant from the previous simulations. An additional difference, besides the random permeability field, is the number of elements used to discretize the domain, since the simulation in this study was done in a field with refined grid ( 40 × 40 elements). Table 5 and Figure 14, Figure 15 and Figure 16 present the results of this analysis.
As expected, the heterogeneity in the absolute permeability affected severely both the recovery and the sweeping efficiency, rendering a major part of the oilfield to not be contacted by the EOR agent (Figure 15). It is well-known the effect of such heterogeneity in secondary and tertiary recovery techniques, causing that large areas of the oilfield may remain unaffected by the recovery process. Figure 15a–d show the mentioned phenomenon in the results from the simulator. Nevertheless, the field permeabilities did not alter the proportionality between the partition coefficient and the oil recovery efficiency. The variations in permeability primarily affect the surfactant concentration reaching every cell of the reservoir, which can significantly alter the results with respect to a constant permeability field. The heterogeneity in the permeability field did not affect however the trends regarding pressure drop, chemical breakthrough or residual saturations, which are similar to those presented above, although it did affected the absolute values of the variables mentioned. Although the permeability does not depend on the characteristics of the surfactants used, it should be taken into account during previous operations to determine an optimum well-pattern distribution, because a non-optimal arrangement in the wells location can cause a major by-pass by the chemical slug of an important part of the OOIP, reducing the effectiveness of the EOR method employed.

3.3. Four-Component System

The second part of this study comprises the analysis and discussion of an EOR process considering a fourth component, namely the salt. On this new simulator no distinctions are made between mono- or divalent salt ions. The salt changes the phase behaviour in systems with surfactants [49,64,65,66,70], so this is crucial to be taken into account during the simulation and laboratory or field experiments, in order to understand the actual behaviour of the system. The partition coefficient is affected by the latter according to Figure 3 [51]. To model this in the simulator, the last curve was adapted to a function (Equation (15)), dependent on the salt content (and the critical salt content, when k c = 1 ). During this paper an initial constant salt content was assumed in porous media.
The scope of the study is then to analyze the influence of salt content on the two phenomena flow, specially in the phase behaviour and the adsorption of chemical species onto the rock formation. The simulations performed in this section considered porous media with several initial salt contents, both below and above the critical salt concentration. This describes the possible two-phase systems in porous media, Type II(+) and II(−), together with different contents of salt in the injection water (brine), in order to study the influence of salinity with the injection process.

Influence of the Salt Component on the Phase Behaviour

This analysis and its subsequent discussion is divided into two different streams: how the salt present in the reservoir affects the recovery (with no salt content in the water being injected), and how the salinity in the water injection influences the sweeping process. In the mathematical model the salt is only present in the aqueous phase, so the transport phenomenon of this component depends only on the properties of the abovementioned phase. A number of simulations were performed with different initial salt contents (constant) in the reservoir, which, according to Equation (15), generate different phase behaviours and types of emulsions (Table 6 and Figure 17). The injection water contains only surfactant and the water itself.
The results demonstrate that, although the salt content has influenced the recovery efficiency, the difference in oil recovered between the different initial salt contents is not as significant as in the previous cases, for different partition coefficients. This is mainly due to what was exposed above: the salt component is present in the aqueous phase and travels faster than the oleous one, so that the salt content in the surfactant bank decreases rapidly, and thus it does the partition coefficient. However, the salt content increases the sweeping efficiency, which agrees with previous results [64,65,66]. This is shown in Figure 18 and Figure 19, where the partition coefficient and the oil saturation are plotted. These simulations show that not only the content of salt in the reservoir is important for the effectiveness of the EOR process, but also the salinity in the injection water influences the performance of the whole process. Therefore, the second part of this four-component study is devoted to present and discuss the results of simulations when brine is being injected and how that can significantly impact in the oil recovery efficiency.
The influence of salinity in the water injection has proved to be more significant than the initial one in the reservoir. This is because a constant input of salt in the water phase, along with the salt already present in the formation, modifies the properties of phase behaviour throughout the reservoir. Different configurations and concentrations of brine being injected were simulated in order to study the influence on the process efficiency. The results are presented in Table 7 and Figure 20.
Regarding the influence of salt in the reservoir, for the same concentration in the injection water, the same trend is observed as at the beginning of this section: the salt content influences the recovery efficiency, but it is not as significantly as the other parameters. This is valid for both Type II(−) and II(+) emulsions. However, the most important difference occurs with the injection water. Small variations in the salinity affect the oil recovered, and this is accentuated for Type II(+) emulsions, where the partition coefficient is more sensitive to the salt content (Equation (15)). Slightly increasing the salt content causes a substantial increase in the partition coefficient and the oil recovered (Figure 21 and Figure 22). This is clearly appreciated in Table 7: the case with a salt concentration in the injection of 1.5% incremented substantially the recovery efficiency due to the formation of a type II(+) emulsion as well as the time required to sweep the reservoir (compared to the other cases). The reason is that the injected salt keeps the partition coefficient in non-zero values. In the previous case the salt, being only present in the aqueous phase and moving faster than the oil, was transported faster than the surfactant bank and that generated low or negligible partition coefficient values. This was avoided with continuous injection of low salt concentration brine. The effect of the salt injected is clearly appreciated in cases with initial salt content equal to 1.5 % . A 3-fold increase in the salt content caused an oil recovery growth of 2-fold, due to the significant increase in the partition coefficient (Table 7). Therefore, it is advisable in oil recovery processes with surfactants the injection of low salt concentration brine in order to improve sweeping efficiency, which agrees the results reported by several authors [47,64,65,66,71].

4. Conclusions

This study aimed at designing a new two-dimensional simulator for a four pseudo-components, two phase flow, applied to a surfactant enhanced oil recovery process. The problems involved, both physical and economic, were formulated and discussed, concluding that new methods of exploitation are necessary to improve the current oil production, until new energy sources are developed. The physical model is described by a system of non-linear partial differential equations, solved by a Finite Difference Method, using an algorithm implemented in MathWorks MATLAB. The higher-order numerical scheme coupled with the TVD flux limiter functions allowed reducing the numerical errors. The compositional approach proved to be a versatile tool performing both secondary and tertiary recovery processes without complications. The phase behaviour is represented by three parameters: the partition, solubilization and swelling coefficients. Along with these three parameters, there is a set of ancillary values which must be carefully defined so as to match the experimental results with the simulation. These values depend on the surfactant chosen (e.g., CMC), so experiments must be carried out beforehand to determine these properties accurately. Besides, further research is deemed necessary to determine a set of mathematical functions describing the phase behaviour parameters in order to match the theoretical results with the laboratory or field tests.
The simplification of the components ternary diagram proved to be an useful and valid approach to model the surfactant properties and phase behaviour, which evidenced to play a significant role in the EOR method. The partition coefficient showed to be the most important parameter in terms of oil recovered, followed by swelling and solubilization ratios. When the surfactant formed Type II(+) (water-in-oil) emulsions, the highest recovery efficiencies were obtained. Furthermore, the influence of the salt only proved to be relevant when it was injected (as brine) in the oilfield, since only in these cases its value remained above the critical salinity in the area where the surfactant slug was present, increasing the sweeping efficiency. Surfactant EOR simulations presented the potential of cEOR methods to sweep the residual oil by means of reducing the IFT. The simplification of the ternary diagram may provide a simple way to design chemical surfactants with the desired properties before carrying out experimental or field tests. Nevertheless, even though the numerical validation has been done, future work is deemed necessary before a full field-scale application of this model. It is considered that lab-scale experiments should be done in 2D or 3D models with surfactants in order to demonstrate the applicability of this numerical simulator. Core flooding is also recommended, though it is advisable the usage of models in which the area and volumetric sweeping efficiencies can be also studied. This will allow determining the parameters in an accurate way before a field-scale application is pursued in order to design novel EOR agents.

Author Contributions

Methodology, Formal Analysis, Validation, Investigation, Resources, Writing-Original Draft Preparation, P.D.; Writing—Review & Editing, Supervision, Project Administration, F.P.

Funding

This research received no external funding.

Acknowledgments

P.D. acknowledges the support of the Roberto Rocca Fellowship Program and the Erasmus Mundus EURICA scholarship program (Program Number 2013-2587/001-001-EMA2).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
cEORChemical Enhanced Oil Recovery
EOREnhanced Oil Recovery
IMPECImplicit Pressure, Explicit Concentration
OOIPOriginal Oil in Place
IFTInterfacial Tension
TDSTotal Dissolved Solids
TVDTotal Variation Diminishing
PDEPartial Differential Equation
REVRepresentative Elementary Volume
DPRDisproportionate Permeability Reduction
CMCCritical Micelle Concentration

Nomenclature

A d Component Adsorption [day–1] μ Absolute Viscosity [Pa·s]
c r Rock Compressibility [Pa–1] σ Interfacial Tension [mN/m]
DDispersion Tensor ϕ Formation Porosity
d m Molecular Diffusion [m2/s] Ω Reservoir Domain
d l Longitudinal Dispersion [m2/s]
d t Transversal Dispersion [m2/s]Superscripts
F L I M Flux Limiting FunctionaAqueous Phase
KAbsolute Permeability [mD]cCapillary
k r Relative PermeabilityHWater-Oil System (no Chemical)
pReservoir Pressure [Pa]jPhase
p w f Bottomhole Pressure [Pa] < n > Time-Step
qFlowrate [m3/day]oOleous Phase
r w Well Radius [m]rResidual
SPhase Saturation
sWell Skin FactorSubscripts
uDarcy Velocity [m/day]cSurfactant Component
VVolumetric ConcentrationiComponent
zOverall Concentration i n Injection
m , n Spatial Grid Blocks
Greek LetterspPetroleum Component
Γ Domain BoundarysSalt Component
δ i j Kronecker DeltatTotal
λ Phase Mobility [m2/(Pa·s)]wWater Component

References

  1. Owen, N.A.; Inderwildi, O.R.; King, D.A. The status of conventional world oil reserves-Hype or cause for concern? Energy Policy 2010, 38, 4743–4749. [Google Scholar] [CrossRef]
  2. BP. Statistical Review of World Energy; Annual Report 63; BP p.l.c.: London, UK, 2014. [Google Scholar]
  3. Dake, L.P. Fundamentals of Reservoir Engineering; Elsevier: Amsterdam, The Netherlands, 1978; ISBN 0-444-41830-X. [Google Scholar]
  4. Bidner, M.S. Propiedades de la Roca y los Fluidos en Reservorios de Petroleo; EUDEBA: Buenos Aires, Argentina, 2001; ISBN 9-502-31150-7. [Google Scholar]
  5. Satter, A.; Iqbal, G.M.; Buchwalter, J.L. Practical Enhanced Reservoir Engineering; PennWell Books: Tulsa, OK, USA, 2008; ISBN 978-1-59370-056-0. [Google Scholar]
  6. Schmidt, R.L.; Venuto, P.B.; Lake, L.W. A Niche for Enhanced Oil Recovery in the 1990’s. Oilfield Rev. 1992, 4, 55–61. [Google Scholar]
  7. Hirasaki, G.J.; Miller, C.A.; Puerto, M. Recent Advances in Surfactant EOR. SPE J. 2011, 16, 889–907. [Google Scholar] [CrossRef]
  8. Lake, L.W. Enhanced Oil Recovery; Prentice-Hall Inc.: Englewood Cliffs, NJ, USA, 1989; ISBN 0-13-281601-6. [Google Scholar]
  9. Druetta, P.; Yue, J.; Tesi, P.; Persis, C.D.; Picchioni, F. Numerical modeling of a compositional flow for chemical EOR and its stability analysis. Appl. Math. Model. 2017, 47, 141–159. [Google Scholar] [CrossRef]
  10. Druetta, P.; Tesi, P.; Persis, C.D.; Picchioni, F. Methods in Oil Recovery Processes and Reservoir Simulation. Adv. Chem. Eng. Sci. 2016, 6, 39. [Google Scholar] [CrossRef]
  11. Druetta, P.; Picchioni, F. Numerical Modeling and Validation of a Novel 2D Compositional Flooding Simulator Using a Second-Order TVD Scheme. Energies 2018, 11, 2280. [Google Scholar] [CrossRef]
  12. Garrett, H.E. Surface Active Chemicals; Pergamon Press: New York, NY, USA, 1972; ISBN 0-080-16422-6. [Google Scholar]
  13. Farrell, H.H.; Gregory, M.D.; Borah, M.T. Progress Report: Big Muddy Field Low-Tension Flood Demonstration Project With Emphasis on Injectivity and Mobility. In Proceedings of the SPE Enhanced Oil Recovery Symposium, Tulsa, OK, USA, 15–18 April 1984. [Google Scholar] [CrossRef]
  14. Reppert, T.R.; Bragg, J.R.; Wilkinson, J.R.; Snow, T.M.; Maer, N.K., Jr.; Gale, W.W. Second Ripley Surfactant Flood Pilot Test. In Proceedings of the SPE/DOE Enhanced Oil Recovery Symposium, Tulsa, OK, USA, 22–25 April 1990. [Google Scholar] [CrossRef]
  15. Maerker, J.M.; Gale, W.W. Surfactant flood process design for Loudon. SPE Reserv. Eng. 1992, 7, 36–44. [Google Scholar] [CrossRef]
  16. Green, D.W.; Willhite, G.P. Enhanced Oil Recovery; Society of Petroleum Engineers: Richardson, TX, USA, 1998; ISBN 978-1-55563-077-5. [Google Scholar]
  17. Iglauer, S.; Wu, Y.; Shuler, P.; Tang, Y.; Goddard, W.A., III. New surfactant classes for enhanced oil recovery and their tertiary oil recovery potential. J. Petrol. Sci. Engi. 2010, 71, 23–29. [Google Scholar] [CrossRef]
  18. Sheng, J.J. Status of surfactant EOR technology. Petroleum 2015, 1, 97–105. [Google Scholar] [CrossRef]
  19. Larson, R.G.; Davis, H.T.; Scriven, L.E. Elementary Mechanisms of Oil-Recovery by Chemical Methods. J. Petrol. Technol. 1982, 34, 243–258. [Google Scholar] [CrossRef]
  20. Pashley, R.M.; Karaman, M.E. Surfactants and Self-Assembly; Applied Colloid and Surface Chemistry; John Wiley & Sons, Ltd.: Chichester, UK, 2004; pp. 61–77. ISBN 978-0-47001-470-7. [Google Scholar] [CrossRef]
  21. Adibhatla, B.; Mohanty, K.K. Parametric analysis of surfactant-aided imbibition in fractured carbonates. J. Colloid Interface Sci. 2008, 317, 513–522. [Google Scholar] [CrossRef]
  22. Akin, S.; Kovscek, A.R. Computed tomography in petroleum engineering research. Geol. Soc. Lond. Spec. Publ. 2003, 215, 23–38. [Google Scholar] [CrossRef]
  23. Yu, Q.; Jiang, H.; Zhao, C. Study of Interfacial Tension Between Oil and Surfactant Polymer Flooding. Petrol. Sci. Technol. 2010, 28, 1846–1854. [Google Scholar] [CrossRef]
  24. Kou, J.; Sun, S.; Wang, X. Linearly Decoupled Energy-Stable Numerical Methods for Multicomponent Two-Phase Compressible Flow. SIAM J. Numer. Anal. 2018, 56, 3219–3248. [Google Scholar] [CrossRef]
  25. El-Amin, M.; Kou, J.; Sun, S.; Salama, A. Adaptive Time-Splitting Scheme for Two-Phase Flow in Heterogeneous Porous Media. Adv. Geo-Energy Res. 2017, 1, 182–189. [Google Scholar] [CrossRef]
  26. Lake, L.; Pope, G.; Carey, G.; Sepehrnoori, K. Isothermal, Multiphase, Multicomponent Fluid-Flow in Permeable Media. 1. Description and Mathematical Formulation. In Situ 1984, 8, 1–40. [Google Scholar]
  27. Pope, G.A.; Carey, G.F.; Sepehrnoori, K. Isothermal, Multiphase, Multicomponent Fluid-Flow in Permeable Media. 2. Numerical Techniques and Solution. In Situ 1984, 8, 41–97. [Google Scholar]
  28. Helfferich, F. Theory of Multicomponent, Multiphase Displacement in Porous-Media. Soc. Petrol. Eng. J. 1981, 21, 51–62. [Google Scholar] [CrossRef]
  29. Fleming, P.; Thomas, C.; Winter, W. Formulation of a General Multiphase, Multicomponent Chemical Flood Model. Soc. Petrol. Eng. J. 1981, 21, 63–76. [Google Scholar] [CrossRef]
  30. Thomas, C.; Fleming, P.; Winter, W. A Ternary, 2-Phase, Mathematical-Model of Oil-Recovery with Surfactant Systems. Soc. Petrol. Eng. J. 1984, 24, 606–616. [Google Scholar] [CrossRef]
  31. Larson, R.G. Analysis of the Physical Mechanisms in Surfactant Flooding. Soc. Petrol. Eng. J. 1978, 18, 42–58. [Google Scholar] [CrossRef]
  32. Hirasaki, G.J. Application of the Theory of Multicomponent, Multiphase Displacement to 3-Component, 2-Phase Surfactant Flooding. Soc. Petrol. Eng. J. 1981, 21, 191–204. [Google Scholar] [CrossRef]
  33. Larson, R.G. The Influence of Phase Behavior on Surfactant Flooding. Soc. Petrol. Eng. J. 1979, 19, 411–422. [Google Scholar] [CrossRef]
  34. Morshedi, S.; Foroughi, S.; Beiranvand, M.S. Numerical Simulation of Surfactant Flooding in Darcy Scale Flow. Petrol. Sci. Technol. 2014, 32, 1365–1374. [Google Scholar] [CrossRef]
  35. Yin, D.-Y.; Pu, H. Numerical Simulation Study on Surfactant Flooding for Low Permeability Oilfield in the Condition of Threshold Pressure. J. Hydrodyn. 2008, 20, 492–498. [Google Scholar] [CrossRef]
  36. Yin, D.; Zhou, Y. Numerical Simulation of Surfactant Flooding and Its Application in Low Permeability Reservoir. Int. J. Control Autom. 2013, 6, 451–472. [Google Scholar]
  37. Porcelli, P.; Bidner, M. Simulation and Transport Phenomena of a Ternary 2-Phase Flow. Transp. Porous Media 1994, 14, 101–122. [Google Scholar] [CrossRef]
  38. Bidner, M.S.; Savioli, G.B. On the Numerical Modeling for Surfactant Flooding of Oil Reservoirs. Mec. Comput. 2002, XXI, 566–585. [Google Scholar]
  39. Rai, S.; Bera, A.; Mandal, A. Modeling of surfactant and surfactant-polymer flooding for enhanced oil recovery using STARS (CMG) software. J. Petrol. Explor. Prod. Technol. 2015, 5, 1–11. [Google Scholar] [CrossRef]
  40. Sulaiman, W.; Lee, E.S. Simulation of surfactant based enhanced oil recovery. J. Petrol. Gas Eng. 2012, 3, 58–73. [Google Scholar]
  41. Camilleri, D.; Engelson, S.; Lake, L.W.; Lin, E.C.; Ohnos, T.; Pope, G.; Sepehrnoori, K. Description of an Improved Compositional Micellar/Polymer Simulator. SPE Reserv. Eng. 1987, 2, 427–432. [Google Scholar] [CrossRef]
  42. Camilleri, D.; Fil, A.; Pope, G.A.; Rouse, B.A.; Sepehrnoori, K. Comparison of an Improved Compositional Micellar/Polymer Simulator With Laboratory Corefloods. SPE Reserv. Eng. 1987, 2, 441–451. [Google Scholar] [CrossRef]
  43. Delshad, M.; Najafabadi, N.F.; Anderson, G.; Pope, G.A.; Sepehrnoori, K. Modeling Wettability Alteration By Surfactants in Naturally Fractured Reservoirs. SPE Reserv. Eval. Eng. 2009, 12, 361–370. [Google Scholar] [CrossRef]
  44. Patacchini, L.; de Loubens, R.; Moncorge, A.; Trouillaud, A. Four-Fluid-Phase, Fully Implicit Simulation of Surfactant Flooding. SPE Reserv. Eval. Eng. 2014, 17, 271–285. [Google Scholar] [CrossRef]
  45. Hand, D.B. Dineric Distribution. J. Phys. Chem. 1929, 34, 1961–2000. [Google Scholar] [CrossRef]
  46. Prouvost, L.; Pope, G.; Rouse, B. Microemulsion Phase-Behavior—A Thermodynamic Modeling of the Phase Partitioning of Amphiphilic Species. Soc. Petrol. Eng. J. 1985, 25, 693–703. [Google Scholar] [CrossRef]
  47. Prouvost, L.P.; Satoh, T.; Sepehrnoori, K.; Pope, G.A. A New Micellar Phase-Behavior Model for Simulating Systems with Up to Three Amphiphilic Species. In Proceedings of the SPE Annual Technical Conference and Exhibition, Houston, TX, USA, 16–19 September 1984. [Google Scholar] [CrossRef]
  48. Saad, N.; Pope, G.A.; Sepehrnoorl, K. Application of higher-order methods in compositional simulation. SPE Reserv. Eng. 1990, 5, 623–630. [Google Scholar] [CrossRef]
  49. Delshad, M.; Pope, G.; Sepehrnoori, K. UTCHEM Version 9.0 Technical Documentation; Center for Petroleum and Geosystems Engineering, The University of Texas at Austin: Austin, TX, USA, 2000; Volume 78751. [Google Scholar]
  50. Yuan, C. Commercial Scale Simulations of Surfactant/Polymer Flooding. Ph.D. Thesis, The University of Texas at Austin, Austin, TX, USA, 2012. [Google Scholar]
  51. Sheng, J. Modern Chemical Enhanced Oil Recovery; Elsevier: Amsterdam, The Netherlands, 2011; ISBN 978-0-08096-163-7. [Google Scholar]
  52. Bidner, M.; Porcelli, P. Influence of phase behavior on chemical flood transport phenomena. Transp. Porous Media 1996, 24, 247–273. [Google Scholar] [CrossRef]
  53. Bidner, M.; Porcelli, P. Influence of capillary pressure, adsorption and dispersion on chemical flood transport phenomena. Transp. Porous Media 1996, 24, 275–296. [Google Scholar] [CrossRef]
  54. Bear, J. Dynamics of Fluids In Porous Media; American Elsevier Publishing Company: New York, NY, USA, 1972; ISBN 978-0-44400-114-6. [Google Scholar]
  55. Chen, Z.; Huan, G.; Ma, Y. Computational Methods for Multiphase Flows in Porous Media; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2006; ISBN 978-0-89871-606-1. [Google Scholar] [CrossRef]
  56. Kuzmin, D. A Guide to Numerical Methods for Transport Equations; University Erlangen-Nuremberg: Erlangen, Germany, 2010. [Google Scholar]
  57. Liu, J.; Delshad, M.; Pope, G.A.; Sepehrnoori, K. Application of higher-order flux-limited methods in compositional simulation. Transp. Porous Media 1994, 16, 1–29. [Google Scholar] [CrossRef]
  58. Barrett, R.; Berry, M.; Chan, T.; Demmel, J.; Donato, J.; Dongarra, J.; Eijkhout, V.; Pozo, R.; Romine, C.; Van der Vorst, H. Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1994; ISBN 978-0-89871-328-2. [Google Scholar] [CrossRef]
  59. Kamalyar, K.; Kharrat, R.; Nikbakht, M. Numerical Aspects of the Convection-Dispersion Equation. Petrol. Sci. Technol. 2014, 32, 1729–1762. [Google Scholar] [CrossRef]
  60. Attwood, D.; Florence, A.T. Surfactant Systems—Their Chemistry, Pharmacy and Biology; Chapman and Hall: London, UK, 1983; ISBN 978-9-40095-775-6. [Google Scholar]
  61. El-Dessouky, H.; Ettouney, H.M. Fundamentals of Salt Water Desalination; Elsevier Science: Amsterdam, The Netherlands, 2002; ISBN 978-0-08053-212-7. [Google Scholar]
  62. Ghannam, M.T.; Hasan, S.W.; Abu-Jdayil, B.; Esmail, N. Rheological properties of heavy & light crude oil mixtures for improving flowability. J. Petrol. Sci. Eng. 2012, 81, 122–128. [Google Scholar] [CrossRef]
  63. Bhuyan, D.; Lake, L.W.; Pope, G.A. Mathematical Modeling of High-pH Chemical Flooding. SPE Reserv. Eng. 1990, 5, 213–220. [Google Scholar] [CrossRef]
  64. Delshad, M.; Han, C.; Sepehrnoori, K.; Najafabadi, N.F. Development of a Three Phase, Fully Implicit, Parallel Chemical Flood Simulator. In Proceedings of the SPE Reservoir Simulation Symposium, The Woodlands, TX, USA, 2–4 February 2009. [Google Scholar] [CrossRef]
  65. Han, C.; Delshad, M.; Pope, G.A.; Sepehrnoori, K. Coupling EOS Compositional and Surfactant Models in a Fully Implicit Parallel Reservoir Simulator Using EACN Concept. In Proceedings of the SPE Annual Technical Conference and Exhibition, San Antonio, TX, USA, 24–27 September 2006. [Google Scholar] [CrossRef]
  66. Lotfollahi, M.; Varavei, A.; Delshad, M.; Farajzadeh, R.; Pope, G.A. A Four-Phase Flow Model to Simulate Chemical EOR with Gas. In Proceedings of the SPE Reservoir Simulation Symposium, Houston, TX, USA, 23–25 February 2015. [Google Scholar] [CrossRef]
  67. Nelson, R.C.; Pope, G.A. Phase Relationships in Chemical Flooding. Soc. Petrol. Eng. J. 1978, 18, 325–338. [Google Scholar] [CrossRef]
  68. Pope, G.A.; Nelson, R.C. A Chemical Flooding Compositional Simulator. Soc. Petrol. Eng. J. 1978, 18, 339–354. [Google Scholar] [CrossRef]
  69. Aulia, A.; El-Khatib, N.A. Mathematical Description of the Implementation of the Adaptive Newton-Raphson Method in Compositional Porous Media Flow. Int. J. Basic Appl. Sci. 2010, 10, 108–115. [Google Scholar]
  70. Fadili, A.; Kristensen, M.R.; Moreno, J. Smart Integrated Chemical EOR Simulation. In Proceedings of the International Petroleum Technology Conference, Doha, Qatar, 7–9 December 2009. [Google Scholar] [CrossRef]
  71. Liu, S.; Miller, C.A.; Li, R.F.; Hirasaki, G. Alkaline/Surfactant/Polymer Processes: Wide Range of Conditions for Good Recovery. SPE J. 2010, 15, 282–293. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of the quarter 5-spot oilfield used for the EOR simulations [9].
Figure 1. Schematic representation of the quarter 5-spot oilfield used for the EOR simulations [9].
Energies 12 00983 g001
Figure 2. Oil recovery during a waterflooding, comparing the results with UTCHEM and GPAS [11].
Figure 2. Oil recovery during a waterflooding, comparing the results with UTCHEM and GPAS [11].
Energies 12 00983 g002
Figure 3. Partition coefficient dependence on the salinity. In this example, according to Equation (15), the optimum volumetric concentration of salt in water ( V s , o p t a ) was adopted from Sheng [51] and it is equal to 2%.
Figure 3. Partition coefficient dependence on the salinity. In this example, according to Equation (15), the optimum volumetric concentration of salt in water ( V s , o p t a ) was adopted from Sheng [51] and it is equal to 2%.
Energies 12 00983 g003
Figure 4. Ternary phase diagrams for Type II(−) (left) and II(+) (right) systems (top) and their simplified representations (bottom) [9].
Figure 4. Ternary phase diagrams for Type II(−) (left) and II(+) (right) systems (top) and their simplified representations (bottom) [9].
Energies 12 00983 g004
Figure 5. Interfacial Tension [mN/m] as a function of the solubilization coefficient ( L p c a ) for a Type II(−) emulsion with a partition coefficient k c = 0.5 .
Figure 5. Interfacial Tension [mN/m] as a function of the solubilization coefficient ( L p c a ) for a Type II(−) emulsion with a partition coefficient k c = 0.5 .
Energies 12 00983 g005
Figure 6. Capillary desaturation curves for non-wetting (oleous) and wetting (aqueous) phases used for this simulation [9].
Figure 6. Capillary desaturation curves for non-wetting (oleous) and wetting (aqueous) phases used for this simulation [9].
Energies 12 00983 g006
Figure 7. Oil recovered (a), total surfactant flowrates (b), and pressure drop (c) as a function of time for different partition coefficients in a surfactant flooding.
Figure 7. Oil recovered (a), total surfactant flowrates (b), and pressure drop (c) as a function of time for different partition coefficients in a surfactant flooding.
Energies 12 00983 g007
Figure 8. Chemical (a) and IFT [mN/m] profiles (b) after 1000 days of simulation. Oil saturation profiles at 1000 days (c) and at the end of the EOR process (d) for a Type II(−) emulsion with k c = 0.5 .
Figure 8. Chemical (a) and IFT [mN/m] profiles (b) after 1000 days of simulation. Oil saturation profiles at 1000 days (c) and at the end of the EOR process (d) for a Type II(−) emulsion with k c = 0.5 .
Energies 12 00983 g008
Figure 9. Chemical (a) and pressure profiles (b) after 3000 days of simulation. Oil saturation profiles at 3000 days (c) and at the end of the EOR process (d) for a Type II(+) emulsion with k c = 2 .
Figure 9. Chemical (a) and pressure profiles (b) after 3000 days of simulation. Oil saturation profiles at 3000 days (c) and at the end of the EOR process (d) for a Type II(+) emulsion with k c = 2 .
Energies 12 00983 g009
Figure 10. Oil recovered as a function of time for two different phase behaviour parameters (II(−)—a and II(+)—b) in a surfactant flooding.
Figure 10. Oil recovered as a function of time for two different phase behaviour parameters (II(−)—a and II(+)—b) in a surfactant flooding.
Energies 12 00983 g010
Figure 11. Surfactant concentration (after 1000 days—a) and oil saturation (after 4000 days—b) profiles for a cEOR flooding with k c = 0.25 , L p c a = 0.5 and L w c o = 1 .
Figure 11. Surfactant concentration (after 1000 days—a) and oil saturation (after 4000 days—b) profiles for a cEOR flooding with k c = 0.25 , L p c a = 0.5 and L w c o = 1 .
Energies 12 00983 g011
Figure 12. Surfactant concentration (after 5000 days—a), oil saturation (after 15,000 days—b), aqueous phase viscosity [mPa·s] (after 5000 days—c), and auqeous phase saturation (after 5000 days—d) profiles for a cEOR flooding with k c = 2 , L p c a = 2 and L w c o = 1 .
Figure 12. Surfactant concentration (after 5000 days—a), oil saturation (after 15,000 days—b), aqueous phase viscosity [mPa·s] (after 5000 days—c), and auqeous phase saturation (after 5000 days—d) profiles for a cEOR flooding with k c = 2 , L p c a = 2 and L w c o = 1 .
Energies 12 00983 g012
Figure 13. Permeability random fields [mD] with normal distribution in X (left) and Y (right) axes.
Figure 13. Permeability random fields [mD] with normal distribution in X (left) and Y (right) axes.
Energies 12 00983 g013
Figure 14. Oil recovered in a random permeability field as a function of time.
Figure 14. Oil recovered in a random permeability field as a function of time.
Energies 12 00983 g014
Figure 15. Oil recovered in a random permeability field for different partition coefficients: k c = 0.25 (a), k c = 0.50 (b), k c = 2 (c), and k c = 4 (d).
Figure 15. Oil recovered in a random permeability field for different partition coefficients: k c = 0.25 (a), k c = 0.50 (b), k c = 2 (c), and k c = 4 (d).
Energies 12 00983 g015
Figure 16. Surfactant concentration after 500 days (a), 1000 days (b), and 2000 days (c), and pressure profile (d) at the end of an EOR flooding with k c = 0.50 in a random permeability reservoir.Random Various
Figure 16. Surfactant concentration after 500 days (a), 1000 days (b), and 2000 days (c), and pressure profile (d) at the end of an EOR flooding with k c = 0.50 in a random permeability reservoir.Random Various
Energies 12 00983 g016aEnergies 12 00983 g016b
Figure 17. Oil recovered as a function of time for different initial salt concentrations in the reservoir.
Figure 17. Oil recovered as a function of time for different initial salt concentrations in the reservoir.
Energies 12 00983 g017
Figure 18. Partition coefficient after 1000 days with initial salt content of 0.5% (a) and 2.5% (b).
Figure 18. Partition coefficient after 1000 days with initial salt content of 0.5% (a) and 2.5% (b).
Energies 12 00983 g018
Figure 19. Oil saturation at the end of the EOR process considering an initial salt content of 0.5% (a) and 2.5% (b).
Figure 19. Oil saturation at the end of the EOR process considering an initial salt content of 0.5% (a) and 2.5% (b).
Energies 12 00983 g019
Figure 20. Oil recovered as a function of time for different sets of initial and injected salt concentrations in the reservoir.
Figure 20. Oil recovered as a function of time for different sets of initial and injected salt concentrations in the reservoir.
Energies 12 00983 g020
Figure 21. Oil saturation at the end of the EOR process under two different configurations of salt (initial/injected): 1.5/0.5% (a) and 1.5/1.5% (b).
Figure 21. Oil saturation at the end of the EOR process under two different configurations of salt (initial/injected): 1.5/0.5% (a) and 1.5/1.5% (b).
Energies 12 00983 g021
Figure 22. Partition Coefficient at the end of the EOR process under two different configurations of salt (initial/injected): 1.5/0.5% (a) and 1.5/1.5% (b).
Figure 22. Partition Coefficient at the end of the EOR process under two different configurations of salt (initial/injected): 1.5/0.5% (a) and 1.5/1.5% (b).
Energies 12 00983 g022
Table 1. General and initial conditions used for the simulations.
Table 1. General and initial conditions used for the simulations.
Geometrical Data of the Reservoir
Length in axis X500 mLength in axis Y500 mLayer thickness5 m
n x elements25/40 n y elements25/40
Rock Properties
Porosity0.25 k x x 200 mD k y y 200 mD
Initial Conditions
S o 0.36 S o r (EOR)0.35 S a r H = S o r H 0.15
Simulation Data
Total time6000 daysChem. injection period200 days z c I N 0.025
Phases’ Physical Data
μ a H 1 cP μ o H 5 cPOil density850 kg/m 3
Water density1020 kg/m 3 IFT50 mN/m
Table 2. Wells operational conditions.
Table 2. Wells operational conditions.
Physical Data
Number of wells2Well radius0.25 mSkin factor0
Operating Conditions
Total production flowrate1650 STB/dayBottomhole pressure55,160 kPa
Table 3. Results of the recovery process as function of the partition coefficient.
Table 3. Results of the recovery process as function of the partition coefficient.
Partition Coefficient k c Solubilization/Swelling L pc a = L wc o Oil Recovered
--m 3 %OOIP
0.25132,06828.5
0.5142,16037.5
2175,93167.5
4192,18181.9
Table 4. Influence of the swelling and solubilization parameters on the oil recovery.
Table 4. Influence of the swelling and solubilization parameters on the oil recovery.
k c L pc a L wc o Oil Recovered k c L pc a L wc o Oil Recovered
---m 3 %OOIP---m 3 %OOIP
0.250.5119,12517.020.5166,18358.8
0.2510.531,75228.2210.540,37935.9
0.251132,06828.521175,93167.5
0.252137,57233.422183,09273.9
0.251232,34828.821276,70068.2
0.50.5123,32220.740.5185,48876.0
0.510.541,96137.3410.545,64140.6
0.51142,16037.541192,18181.9
0.52152,38546.642196,47885.8
0.51242,62937.941292,82282.5
Table 5. Oil recovery efficiency in a random permeability field as a function of the partition coefficient.
Table 5. Oil recovery efficiency in a random permeability field as a function of the partition coefficient.
Partition Coefficient k c Solubilization/Swelling L pc a = L wc o Oil Recovered
--m 3 %OOIP
0.25137,07533.0
0.5144,60539.7
2162,24655.3
4171,59663.6
Table 6. Oil recovery efficiency as a function of the initial salt concentration in the reservoir.
Table 6. Oil recovery efficiency as a function of the initial salt concentration in the reservoir.
Initial Salt ContentCritical Salt ContentOil Recovered
%%m 3 %OOIP
0.5229,15325.9
1.5231,86628.3
2.5234,40530.6
4235,30231.4
Table 7. Oil recovery efficiency as a function of the initial and injected salt concentrations in the reservoir.
Table 7. Oil recovery efficiency as a function of the initial and injected salt concentrations in the reservoir.
Initial Salt ContentCritical Salt ContentSalt InjectedOil Recovered
%%%m 3 %OOIP
0.52152,46046.6
1.520.540,92636.4
1.52152,99747.1
1.521.580,75871.8
42157,15250.8

Share and Cite

MDPI and ACS Style

Druetta, P.; Picchioni, F. Simulation of Surfactant Oil Recovery Processes and the Role of Phase Behaviour Parameters. Energies 2019, 12, 983. https://doi.org/10.3390/en12060983

AMA Style

Druetta P, Picchioni F. Simulation of Surfactant Oil Recovery Processes and the Role of Phase Behaviour Parameters. Energies. 2019; 12(6):983. https://doi.org/10.3390/en12060983

Chicago/Turabian Style

Druetta, Pablo, and Francesco Picchioni. 2019. "Simulation of Surfactant Oil Recovery Processes and the Role of Phase Behaviour Parameters" Energies 12, no. 6: 983. https://doi.org/10.3390/en12060983

APA Style

Druetta, P., & Picchioni, F. (2019). Simulation of Surfactant Oil Recovery Processes and the Role of Phase Behaviour Parameters. Energies, 12(6), 983. https://doi.org/10.3390/en12060983

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop