Effects of Solvent Properties and Injection Strategies on SolventEnhanced Steam Flooding for Thin Heavy Oil Reservoirs with SemiAnalytical Approach
^{1}
Department of Petroleum Engineering, China University of Petroleum, Beijing, 18 Fuxue road, Changping
102249  PR China
^{2}
The University of Oklahoma, 660 Parrington Oval, Norman, OK
730190390  USA
email: liuhao880501@gmail.com
^{*} Corresponding author
Received:
8
August
2016
Accepted:
24
May
2017
Compared with conventional steam flooding and SteamAssisted Gravity Drainage (SAGD), SolventEnhanced Steam Flooding (SESF) is considered a more effective method for improving heavy oil recovery in thin reservoirs in terms of higher thermal efficiency and oil production rate. However, there remains a deficiency of accurate and efficient methods to evaluate and design an SESF project in the field. A semianalytical model is proposed in this paper to predict the recovery performance of SESF and investigate the effects of solvent properties and injection strategies on the SESF process for thin heavy oil reservoirs. The proposed model provides a simple method to simulate not only single solvent injection but also multisolvent injection by cooperating different values of solvent operating thickness and solvent solubility. To validate the model’s accuracy, comparisons are made between the proposed model results and the numerical simulation results for a specific heavy oil reservoir case. The results indicate that SESF can achieve a considerably higher oil production rate at the early recovery stage than steam flooding. Moreover, the paper also demonstrates that a higher injection rate results in a lower thermal efficiency increment when well spacing is constant. Nevertheless, a high injection rate may also be suitable for longer well spacing owing to the improvement of the viscosity profile beyond the edge of the steam zone caused by longer contact time between the solvent and crude oil.
© H. Liu et al., published by IFP Energies nouvelles, 2017
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Nomenclature
A : Area of heated zone, m^{2}
(t): Solvent diffusion area at time t, m^{2}
c_{P} : Heat capacity at constant pressure within the reservoir, J/(kg °C)
c_{cap} : Heat capacity of understrata and overburden rock
C_{si} : Solvent solubility at the steam front
C_{s} : Solvent concentration at ε
: Average solvent concentration in the solvent penetration depth
D_{t } : Effective diffusion coefficient of solvent at time t, m^{2}/d
D : Solvent molecular diffusion coefficient, m^{2}/s
dr/dt : Radial velocity of steam front, m/s
dx/dt : Longitudinal velocity of steam front, m/s
dδ_{c }/dt : Solvent diffusion rate, m/s
f_{w} : Water fractional flow at specific water saturation
: Slope of tangent of the water fractional flow curve at specific water saturation
h_{w} : Enthalpy of saturated water, kJ/kg
H_{o} : Heat injection rate per unit well length, kJ/md
L_{f} : Horizontal displacement of cold water front in the second stage, m
L_{e} : Effective operating thickness of solvent, m
L_{v} : Latent heat of vaporization of steam, kJ/kg
q_{hloss} : Heat loss rate to understrata and overburden per area, kJ/m^{2}d
q_{heff} : Required heat rate associated with the steam zone expansion per unit well length, kJ/md
S_{or} : Residual oil saturation
S_{wc} : Connate water saturation
S_{wt} : Water saturation at the end of operating thickness
S_{wp} : Water saturation at the production well
T_{r} : Initial reservoir temperature, °C
t_{c} : Time when the first stage ends, d
r : Radius of steam zone area, m
: Mass rate of cold steam condensate per unit well length, kg/md
x : Horizontal displacement of steam front, m
Greek letters
α_{r} : Longitudinal dispersity factor
δ_{c } : Solvent penetration depth, m
ε : Distance beyond the edge of the steam front, m
λ_{cap} : Thermal conduction coefficient of understrata and overburden rock, kJ/md °C
Subscripts
cap: Understrata and overburden rock
Introduction
Heavy oil resources have an important role in crude oil reserve replacement to meet the world’s future energy demands [1]. However, a considerable portion of the heavy oil resides in thin formations and in some countries more than 80% lies in reservoirs with a less than a 10 m pay zone; this has posed a major challenge for efficient development [2]. To date, the predominant in situ recovery methods for reservoirs with thin pay zone are steam injection by horizontal wells. Among these, the steamflooding strategy has been regarded as a significantly superior technology compared to SteamAssisted Gravity Drainage (SAGD) operations in terms of energy efficiency [3] and has been applied widely in thin heavy oil reservoirs in China [46]. However, significant heat losses to the understrata and overburden can render the steamflooding process uneconomic because of the small thickness of the reservoir [7]. SolventEnhanced Steam Flooding (SESF) requires substantially less steam usage and net injected energy and thus is considered as a highly effective method for improving heavy oil recovery in thin reservoirs [7, 8].
Solvent properties and injection strategies are important factors that control the performance of SESF and thus determine the economic viability of the process. Therefore, it is necessary to fully investigate the effects of these factors on production performance before implementing in the field. By conducting 63 experiments on a threedimensional elemental physical model, Redford [9] first demonstrated the effects of solvent properties on the production performance and indicated that naphtha with steam could significantly reduce oil viscosity. Besides, bitumen recovery could be improved further by adding CO_{2} or ethane because they can provide more drive energy for the naphtha to enter the reservoir. Nasr and Pierce [10] conducted a series of experiments on steamCO_{2} injection strategies in a scaled model and revealed that steamCO_{2} continuous injection results in improved performance compared to steamCO_{2} sequential injection. Moreover, oil recovery rates are typically accelerated and improved during the initial stage of the process. Although the effects of the solvent properties and injection strategies of several solvents on steam flooding can be identified by welldesigned experiments, the number of solvents tested is limited whereas the experiment requirements are large. Therefore, the experimental data is difficult to analyse.
Based on laboratory experiments, several attempts using numerical simulation have been undertaken to study the SESF process. Shutler and Boberg [11] first used numerical simulations to delineate the recovery mechanism of the solventsteam process. The results indicate that mediumalkane solvents provide the greatest increase in total oil production, whereas heavyalkanes solvents do not improve recovery. They also proposed the corollary that the production of the process is primarily determined by the placement of the solvent in the reservoir, which in turn is controlled by the steam movement and solvent volatility. Zhao et al. [8] reported reservoir simulation results on the application of SESF for a bitumen reservoir with a net pay of 4 m. The results present optimum steam and steamsolvent flooding injection strategies by comparing the cumulative Energy injected to produced Oil Ratio (cEOR) between different cases. Hence, numerical simulation was proven to be a powerful method to calculate SESF production. However, the calculation of the numerical simulation is not only timeconsuming but also relies heavily on a large amount of experimental data, which limits the application of the numerical simulation method.
For quick calculation and an easy method for the solventsteam process, several analytical and semianalytical methods were developed. Sharma and Gates [12] considered mathematical models to calculate the length scales of the mass and heat transfer beyond the edge of a coinjection chamber. Gupta and Gittins [13] developed a semianalytical approach to estimate oil drainage in the coinjection of steam and solvent. The solvent distribution beyond the steam chamber calculated by their models can be an indicator of different solvents without introducing excessive experimental data. However, their mathematical models are based on Butler’s SAGD model [14, 15] for oil drainage and do not consider the blocking effects of water on the solvent diffusion, which makes the models not applicable for SESF process.
To the best of our knowledge, applying a mathematical model to understand the SESF process for heavy oil reservoirs with thickness less than 10 m has never been undertaken. In this study, to predict SESF production performance and study the effects of solvent properties and injection strategies on the SESF process, a new semianalytical model is proposed based on mass and energy balance principles, diffusion theory, and the theory of BuckleyLeverett. First, the mathematical model is solved in a semianalytical method. Next, comparisons are made between the calculated model and numerical simulations to verify the proposed model. Finally, with the validated model, the effects of the solvent properties and injection strategies on the production performance of SESF are studied in detail.
1 Mathematical Model
A simplified diagram of the SESF process is presented in Figures 1 and 2. To decouple the primary mechanisms from the other complexities that accompany the SESF process, important assumptions and simplifications are provided below:

based on the position of the steam front, the entire production process is divided into three stages: before the steam front reaches the upper and lower boundary (steam front rising stage) (Fig. 1b), after the steam front reaches the boundary without steam condensate reaching the production well (steam front spreading stage) (Fig. 2b), after the cold steam condensate reaches the production well (steam condensateproducing stage). The pressure inside the steam zone is constant during all stages [16, 17];

an interface model proposed by van Lookeren [18] has been used widely to describe the shape of the steam front without considering heat loss. It is a reasonable model for the steamflooding process in thick reservoirs; however, it is not applicable for thin reservoirs with considerably more severe heat losses to the understrata and overburden. Thus, according to the reservoir numerical simulation results, the geometries of the steam and water fronts are assumed to be a cosine curve and circular shaped, respectively;

all mobile oil is assumed to be displaced by the steam and hot water; thus, only residual oil is left in the steam and hot water zone;

the gasphase solvent can move freely in the hot water zone owing to the high temperature and low oil saturation [17, 19], whereas the cold steam condensate may block the diffusion of solvent into the heavy oil. By investigating VAPorEXtraction (VAPEX) experiments, Das and Butler [20] proposed that the contact area of the solvent and crude oil was directly related to the rate of the mass transfer, as indicated in Figure 3a. Thus, it is assumed that the diffusion coefficient is proportional to the contact area in a limited distance beyond the edge of the steam front, as illustrated in Figure 3b. The limited distance is defined as the effective operating thickness of the solvent, which is determined by the volatility of solvent;
Figure 1 Illustration of steam front rising stage. a) Threedimensional diagram of the stage; b) twodimensional diagram of the stage.
Figure 2 Illustration of steam front spreading stage. a) Threedimensional diagram of the stage; b) twodimensional diagram of the stage.
Figure 3 Schematic of contact area change at interface. a) Das and Butler’s perspective; b) simplification based on the Das and Butler’s perspective.

medium solvents (C_{6}C_{12}) are considered to provide the greatest increase in total oil production because they can both travel with the steam in a vaporized phase and condense in the cooler regions of the reservoir and diffuse in the crude oil [21]. Light solvents (CO_{2}, C_{1}C_{5}), which have considerably higher volatility, are believed to transport further into the reservoir and thus increase the operating thickness when coinjecting with medium solvents. Therefore, the solvents used in this model are assumed to be medium solvents and a combination of medium and light solvents;

local equilibrium is assumed along the steam front such that the equilibrium state inside the steam zone is not affected by the phase behaviour at the steam front [22]. Other assumptions, such as no mutual solubility between the water and hydrocarbon components and no variation in the rockfluid properties along the horizontal section of the well are identical to those proposed by Keshavarz et al. [23].
The mathematical model is divided into three parts based on the three stages (steam front rising stage, steam front spreading stage, steam condensateproducing stage). These are introduced in detail in the following.
1.1 Steam Front Rising Stage
The process of the steam front rising stage is short because of the small thickness of the reservoir, which makes the generating of steam condensate negligible. Therefore, the velocity of the steam front can be calculated by the energy conservation equation:(1) (2)where H_{o} is the heat injection rate per unit well length; A is the area of the steam front; ρ_{r}, ρ_{o}, ρ_{s}, and ρ_{w} are the densities of the rock, crude oil, steam, and water, respectively; c_{Pr}, c_{Po}, c_{Ps}, and c_{Pw} are the heat capacities of the rock, crude oil, steam, and water, respectively; ϕ is the reservoir porosity; S_{or} is the residual oil saturation; S_{wc} is the connate water saturation; dl is the unit length of well; T_{s} is the steam temperature; T_{r} is the initial reservoir temperature; r is the radius of steam zone area; and U is the velocity of steam front.
The unsteadystate solvent concentration distribution along the steam front can be described by Equation (3) [24]:(3)
The corresponding boundary conditions are:(4)where ε is the distance beyond the edge of the steam front; C_{si} is the solvent solubility at the steam front; C_{s} is the solvent concentration at ε; D is the solvent molecular diffusion coefficient; α_{L} is the longitudinal dispersity factor; and t_{s} is the solvent injection time.
Owing to the nonlinearity of Equation (3), the Heat Integral Method (HIM) is used to determine the concentration profile. This method is a powerful approximation method for solving a variety of diffusion problems [25, 26]. Rabiei Faradonbeh et al. [24] realized that the exponential form of the HIM method provides more accurate predictions of the concentration profile than the polynomial form; therefore, the exponential form of the concentration distribution is used for the calculations of the fluid mixture and oil flow rate in this study.
The concentration distribution in exponential form ahead of the interface can be expressed as:(5)where δ_{c } is the solvent penetration depth. Transforming Equation (3) into integral form in the domain of the solvent penetration depth and substituting Equation (5) for C_{s}, Equation (3) can be simplified into:(6)
According to the Butler’s theory, all mobile oil is displaced by the steam and only residual oil is left in the steam zone. Therefore, based on the mass conservation equation and assumption (6), the oil production rate in this stage can be expressed as:(7) (8)where is the solvent diffusion area at time t and is the average solvent concentration in the solvent penetration depth.
As indicated in Figures 4a and 4b, the steam front rising stage can be divided into two secondary stages according to the position of the solvent boundary. They can be described by Equations (9) and (10), respectively.(9) (10)
Figure 4 Schematic of solvent diffusion process in first stage. a) Before solvent reaches boundary; b) after solvent reaches boundary. 
1.2 Steam Front Spreading Stage
The steam front spreading stage is schematically presented in Figure 2. In this stage, the steam is already touching the understrata and overburden; therefore, the heat injected is divided into two parts based on the energy conservation equation:(11)where q_{hloss} is the heat loss rate to the understrata and overburden per area and q_{heff} is the required heat rate associated with the steam zone expansion per unit well length.
The heat consumption rate used to expand the steam zone is expressed by the energy increase of the area. As described in assumption (2) and Figure 2b, the area of the steam zone is twice that of the hot water zone. Therefore, the required heat rate associated with the steam expansion per unit well length in this stage can be expressed as:(12)
According to assumption (2), the geometry of the water front is assumed to be a diameterconstant circular shape; thus, the change of the heated area can be transformed to the change of distance:(13)
Carslaw and Jaeger [27] built and solved the heat loss model of a semiinfinite plate with a constant temperature boundary. As the hot water mainly contacts the understrata and overburden, the heat loss rate is written as:(14)where T_{s} and T_{w} are the temperature of the steam and hot water, respectively; ρ_{cap} is the density of the understrata and overburden rock; λ_{cap} is the thermal conduction coefficient of the understrata and overburden rock; c_{cap} is the heat capacity of the understrata and overburden rock; x is the horizontal displacement of the steam front; and t_{c } is the time when the first stage terminates.
Rabiei Faradonbeh et al. [28] also employed HIM to determine the temperature profile. Therefore, the temperature of the hot water can be expressed as:(15)
Based on the mass conservation law, the generating rate of the cold steam condensate per unit well length can be given as:(16)where is the mass rate of cold steam condensate per unit well length; X is the steam quality; L_{v} is the latent heat of steam vaporization; h_{w} is the enthalpy of saturated water; and A_{s} and A_{w} are the areas of the steam zone and hot water zone, respectively.
The areas of the steam zone and hot water zone in Equation (16) can be rewritten as Equation (17) and Equation (18) because of the cosine curve and circular shape of the steam and water fronts, respectively:(17) (18)
Shutler and Boberg [11] proposed that pistonlike flow through a thin layer is assumed; therefore, the isothermal twophase flow theory of Buckley and Leverett [29] can be introduced into isothermal zones. Based on Shutler’s theory, the horizontal displacement of the cold water front can be expressed as:(19)where is the slope of the tangent of the water fractional flow curve at the cold water front. Because the solvent penetration depth is extremely thin when the cold water begins to generate, the horizontal displacement of the steam front is assumed to be not affected by the viscosity reduction by the solvent.
Cold water is assumed to block the diffusion of the solvent into the crude oil. Based on assumption (4), the effective diffusion coefficient of the solvent can be given by Figure 5 and Equations (20)(22):(20) (21) (22)where L_{e} is the effective operating thickness of the solvent; is the average water saturation in L_{e}; x_{1} and x_{2} are the coordinates of the steam edge and the effective operating front, respectively; and D_{t } is the effective diffusion coefficient of solvent at time t.
Figure 5 Schematic of solvent operating thickness. 
The solvent diffusion rate and oil production rate can also be calculated by Equations (6) and (7).
1.3 Steam CondensateProducing Stage
The steam condensateproducing stage begins when the steam condensate is produced from the production well. The water saturation profile between the steam edge and production well is determined by the oil viscosity profile, which is in turn influenced by the solvent concentration profile.
The following mixing rule is used for calculating the mixture viscosity of the oil phase [30]:(23)where μ_{o}, μ_{s}, and μ_{m} are the viscosities of the heavy oil, solvent, and mixture of heavy oil and solvent, respectively.
By employing the theory of Buckley and Leverett, the slope of the tangent of the water fractional flow curve at the production well can be written as:(24)
The oil production rate can be rewritten from Equation (7) by considering the water and solvent production:(25)where S_{wp} is the water saturation at the production well, which is determined by and μ_{m}, and is the fractional flow of water at the production well, which can be acquired by considering μ_{m} and S_{wp}.
2 Calculation for Mathematical Model
2.1 Model Treatment
A semianalytical method is used to solve the model due to the complexity of the equations. To solve the semianalytical model conveniently, several equations must first be solved.
The equation describing the movement of the steam front can be obtained by integrating Equations (11)(14).(26)
Equation (26) belongs to the Volterra integral equations of the second kind, which makes it difficult to solve directly. In this article, Equation (26) is solved with the support of a Laplace transform and the displacement of the steam front is given as:(27)where(28) (29)
2.2 Calculation Flow for Mathematical Model
Because the solvent diffusion rate and steam front velocity change with time, time is first discretized into small intervals. As mentioned above, the following steps are used for the calculation of the oil production rate:

for a given heat injection rate and solvent component, Equations (1), (2), (5), and (6) are coupled and solved for the steam front velocity U, radius of the steam zone area r, solvent diffusion rate dδ _{ c }/dt, and solvent penetration depth δ _{ c };

the production rate in the first stage is calculated using Equations (7)(10);

repeat steps 2 and 3 substituting δ _{ c } calculated in step 2 into Equation (6) until the steam front reaches the understrata and overburden, i.e., when the process enters the next stage;

the steam front velocity dx/dt and the steam front displacement x in the second stage are calculated using Equations (15) and (27);

the cold steam condensate generating rate is calculated using Equations (16)(18). Then, Equations (5), (6), and (19)(22) are coupled and solved for the solvent dδ _{ c }/dt, δ _{ c }, and the cold steam condensate distribution beyond the edge of the steam front;

the production rate in the second stage is also calculated using Equations (7)(10);

repeat steps 57 substituting δ _{ c } calculated in step 6 into Equation (6) until the cold steam condensate reaches the production well, when the process enters the third stage;

, dδ _{ c } /dt, and δ _{ c } are calculated in the same manner as step 6. The water fractional flow f _{w}(S _{wo}) in the third stage is calculated using Equations (5), (23), and (24). Then, the oil production rate of the third stage can be calculated using Equation (25);

repeat steps 9 substituting δ _{ c } calculated in step 9 into Equation (6) until f _{w}(S _{wo}) = 1, when the process ends.
3 Results and Discussion
3.1 Model Validation
It is difficult to address all combinations of different solvents by experiments and numerical simulation, making the validation of the SESF process far from comprehensive. Given this situation, the model with solvent concentration value of zero is validated to ensure that the proposed model is a credible base model for the SESF calculation. This method of validation is considered as reasonable in this article because the major objective of the model is to investigate the effects of the solvent properties on the production of SESF, not that of a specific solvent.
STARS is an effective thermal recovery reservoir numerical simulator developed by the CMG Company. As the results of STARS are exact numerical solutions that can closely reveal the real reservoir production process, it frequently serves as a reference to test and verify the correctness of other models. Hence, to validate the correctness of the proposed model, the results of the model are compared with those of STARS. The parameters used in the proposed model and STARS are all based on the average property parameters of one oil field in China. The well spacing between the production well and injection well is 40 m. The thickness of the reservoir is 7 m. The steam injection rate per unit length is 0.09 t/md, which is low because of the thin thickness of the reservoir. The temperature of the steam is 300 °C and the oil viscosity is 300 mPa s at the reservoir condition. Table 1 lists the reservoir properties and injection parameters used in the calculation of both the new model and STARS.
Average property parameters of X oil field in China.
The steam zone development comparison between the proposed model and STARS is made using Figure 6. The geometries and movements of the steam front and hot water front predicted by the proposed model are similar to the numerical simulation results, which also proves the reasonability of assumption (2). Figure 7 depicts the cold steam condensate front position at the end of the steam front spreading stage. The difference of the condensate front position between the two methods is only 2 m, which also proves the feasibility of the BuckleyLeverett theory in the steamflooding process in a thin layer.
Figure 6 Comparison of steam front movement between proposed model and STARS. a) Steam flooding time = 46 days; b) steam flooding time = 200 days; c) steam flooding time = 400 days. 
Figure 7 Comparison of steam condensate front between proposed model and STARS. 
Figure 8 displays the oil production rate as a function of the steam flooding production time. The oil production rate increases before the steam front reaches the understrata and overburden. In the second stage, the production rate remains virtually stable because the hot water zone between the steam zone and boundaries alleviates the heat loss. After the steam condensate reaches the production well, the oil production rate begins to decrease sharply; in this stage the steam zone and hot water zone stop expanding at the point when the heat loss rate to the overburden and understrata equals the heat injection rate [31, 32].
Figure 8 Comparison of calculated oil production with STARS results. 
RMSRE (RootingMeanSquare of Relative Error) (Eq. 30) and RE (Relative Error) (Eq. 31) are employed to evaluate the accuracy of the proposed model in calculating the cumulative oil production and oil production rate [33]:(30) (31)where NPTS is the number of data points; is the calculated value at point i by the proposed model; and is the calculated value at point i by STARS.
In Figure 9, the cumulative oil calculated by the proposed model and STARS are compared. The RE of the cumulative oil for each point is less than 7%; the RMSRE is only 3.5%. Thus, the proposed model is accurate in calculating the steamflooding process. Figure 10 indicates the RE distribution of the oil production rate with time. It can be clearly observed that the RE is significantly higher in the third stage than in the first and second stages. According to our analyses, the greater difference between the two results in the third stage is largely due to the onedimensional assumption of the BuckleyLeverett theory. The RE and absolute error (difference between the calculated oil production rate by the proposed model and STARS) in the later portion of the third stage is less than that in the middle of the same stage. As is known from a previous study, the solvent mainly affects the early and late stages of the steam injection process [10, 34], which is consistent with the low RE stage in the proposed model. Therefore, the proposed model is sufficiently accurate to calculate the SESF process.
Figure 9 RE distribution along cumulative oil production. 
Figure 10 RE distribution along oil production rate and time. 
After the validation of the proposed model, the oil production rate of the steam flooding is compared with that of SESF, as illustrated in Figure 11. The solubility of the solvent at the steam front is 0.2, the diffusion coefficient of the solvent into the heavy oil is 5.5 × 10^{−2} m^{2}/d, which is one order of magnitude greater than the Das and Butler’s results that pertain to Heleshaw laboratory models of VAPEX [35]. Gupta and Gittins [13] regard the higher diffusion coefficient as reasonable because of the higher temperature and existence of a blanket layer. The longitudinal dispersity factor is 5 × 10^{−3} m and the viscosity of solvent is 4.5 mPa s. It can be seen that oil rate increases rapidly in the first stage with the contact area increasing until the steam front reaches the boundaries when the oil rate reaches its highest point. Then, it begins to decline owing to the lower solvent concentration gradient beyond the steam front and the blocking effects of the steam condensate on the solvent diffusion. Moreover, the production time is shortened because of the improvement of the oil viscosity profile beyond the edge of the steam caused by the dissolution of the solvent into the heavy oil. These are consistent with the working stage of the solvent mentioned in other literatures [10, 34], which in turn verifies the proposed model.
Figure 11 Comparison of calculated oil production rate between steam flooding and SESF. 
3.2 Effects of Solvent Properties on SESF
Figure 12 displays the effects of different effective operating thicknesses of the solvent on both the production rate and reduced steam oil ratio. The reduced steam oil ratio is defined as the steam oil ratio difference between SESF and steam flooding, which represents the thermal efficiency increment by the solvent. We can see that the oil rate decreases slower with higher effective operating thickness, while the production time remains unaffected. The larger the effective operating thickness, the greater the thermal efficiency increment. This means that by pushing the solvent further into the reservoir, the thermal efficiency can be increased without additional consumption of solvent, thus increasing the economic benefits. This is similar to the phenomenon in the experiments of Redford [9].
Figure 12 Effects of solvent operating thickness on production performance of SESF. 
Figure 13 indicates that the solvent solubility in heavy oil also has significant effects on the SESF process. It can be observed that both the oil production rate in the early stage and thermal efficiency increment increase with the solvent solubility. Moreover, according to Equations (23) and (24), for the benefit of a lower mixture viscosity beyond the steam edge, a higher solvent solubility tends to result in a shorter production time without decreasing the cumulative oil production, such that heavy oil is recovered considerably more efficiently.
Figure 13 Effects of solvent solubility on production performance of SESF. 
However, it should be noted that a larger operating thickness always means a lower solvent solubility in terms of a single solvent. Consequently, an optimal solvent must be carefully selected for the best performance of SESF when injecting a single solvent with steam. Multisolvent injection, which is rarely mentioned in other analytical models, is provided with a simple method to simulate in the proposed model using different values of operating thickness and solubility.
3.3 Effects of Injection Strategy on SESF
Figure 14 illustrates the effects of different injection rates on SESF production. It can be concluded that the differences between the largest oil rate and the stable oil rate are virtually the same among different injection rates. The reason is that the diffusion rate is largely determined by the concentration gradient and contact area, not the interface velocity. Further, the generating of water within the operating thickness accelerates with an increasing injection rate, which results in a sharp decline of oil production rate at an early recovery stage for a high injection rate. Moreover, Figure 14 indicates that the thermal efficiency increment increases with a decrease of the injection rate because of the longer contacting time of the solvent with the oil and thus relatively shorter production times compared with steam flooding.
Figure 14 Effects of injection rate on production performance of SESF. 
Figure 15 displays the effects of different well spacing on both the production rate and reduced steam oil ratio. It is clear that the oil rates in the early stage are identical for different well spacing. Moreover, the larger the well spacing, the lower the thermal efficiency increment. This is because the diffusion rates are the same with different well spacing and thus the oil is unaffected by the solvent increases related to the well spacing. Therefore, solvents with higher operating thickness, which means wider working range, are more suitable for longer well spacing.
Figure 15 Effects of well spacing on production performance of SESF. 
From Figures 14 and 15, it can also be concluded that despite the fact that a higher injection rate may result in lower thermal efficiency increment when well spacing is constant, it may be more suitable for longer well spacing owing to the improvement of saturation and oil viscosity profiles beyond the steam edge caused by the longer contact time. Thus, an optimal injection rate must be carefully selected for the best performance of SESF with a given well spacing.
Conclusion
In this study, a semianalytical model for predicting SESF production performance was established. Then the theoretical model was validated using STARS results and error analysis. The production performances of steam flooding and SESF were also compared. Moreover, with the theoretical model, the effects of solvent properties, injection strategies on oil production rate, and thermal efficiency increment were investigated. The findings can be summarized as follows:

compared with the STARS results for one case based on average reservoir parameters, the proposed model proved to be sufficiently accurate in calculating the SESF production performance;

SESF can achieve a considerably higher oil production rate at the early recovery stage than steam flooding. Further, the production time is shortened owing to the improvement of the viscosity profile beyond the edge of the steam, which is caused by the solvent dissolution into the heavy oil;

the operating thickness and solubility of a single solvent are proportional to the oil production increment of SESF; however, they cannot be increased simultaneously. Thus, an optimal solvent must be carefully selected for the best performance of SESF when injecting a single solvent with steam. The proposed model also provided a simple method to simulate multisolvent injection using different values of operating thickness and solubility;

a higher injection rate can result in a lower thermal efficiency increment when well spacing is constant. However, it may be suitable for longer well spacing owing to the improvement of the viscosity profile beyond the steam edge, which is caused by longer contact time between the solvent and crude oil.
Acknowledgments
The authors wish to thank the Research Institute of Petroleum Exploration & Development, Liaohe Oilfield Company, PetroChina (2013JS9399). This work was also supported in part by a grant from the National Science and Technology Major Projects of China (2016ZX05012005004) for financial support. Discussions with colleagues of an oil reservoir numerical simulation group provided great insight for our paper.
References
 Besson C. (2005) Resources to Reserves: Oil & Gas Technologies for the Energy Markets of the Future, OECD Publishing, Paris. [Google Scholar]
 Adams D. (1982) Experiences with waterflooding Lloydminster heavyoil reservoirs, J. Petro. Technol. 34, 8, 1643–1650. [CrossRef] [Google Scholar]
 Stalder J.L. (2007) Unlocking bitumen in thin and/or lower pressure pay using cross SAGD (XSAGD). Proceedings of the Canadian International Petroleum Conference, F. Petroleum Society of Canada. [Google Scholar]
 Fu C. (2007) The research and implementation of Qi40 steamflood techniques, Petrol. Geol. Eng. 21, 2, 39–42. [Google Scholar]
 Zhao S. (2012) Research and application on key technique of efficient development of heavy oil resources in ultrathin layer in south region of Chenjiazhuang oilfield, Petrol. Geol. Reco. Effic. 19, 3, 98–100. [Google Scholar]
 Song Y. (2013) Optimization design of steam flooding for horizontal well in thin heavy oil reservoir, FaultBlock Oil & Gas Field 2, 239–241. [Google Scholar]
 Zhao D.W., Wang J., Gates I.D. (2014) Thermal recovery strategies for thin heavy oil reservoirs, Fuel 117, 431–441. [CrossRef] [Google Scholar]
 Zhao D.W., Wang J., Gates I.D. (2013) Optimized solventaided steamflooding strategy for recovery of thin heavy oil reservoirs, Fuel 112, 50–59. [CrossRef] [Google Scholar]
 Redford D. (1982) The use of solvents and gases with steam in the recovery of bitumen from oil sands, J. Can. Pet. Technol. 21, 1, 45–53. [CrossRef] [Google Scholar]
 Nasr T., Pierce G. (1995) SteamCO_{2} recovery processes for bottom water oil reservoirs, J. Can. Pet. Technol. 34, 7, 42–49. [CrossRef] [Google Scholar]
 Shutler N., Boberg T. (1972) A onedimensional analytical technique for predicting oil recovery by steamflooding, SPE J. 12, 6, 489–498. [CrossRef] [Google Scholar]
 Sharma J., Gates I.D. (2010) Steamsolvent coupling at the chamber edge in an in situ bitumen recovery process, Proceedings of the SPE Oil & Gas India Conference & Exhibition, Mumbai, F. Society of Petroleum Engineers. [Google Scholar]
 Gupta S.C., Gittins S. (2012) An investigation into optimal solvent use and the nature of vapor/liquid interface in solventaided SAGD process with a semianalytical approach, SPE J. 17, 4, 1255–1264. [CrossRef] [Google Scholar]
 Butler R.M. (1985) A new approach to the modelling of steamassisted gravity drainage, J. Can. Pet. Technol. 24, 3, 42–51. [CrossRef] [Google Scholar]
 Butler R.M. (1994) Steamassisted gravity drainage: concept, development, performance and future, J. Can. Pet. Technol. 33, 2, 44–50. [CrossRef] [Google Scholar]
 Thomas S. (2008) Enhanced oil recovery – an overview, Oil Gas Sci. Technol. – Rev. IFP 63, 1, 9–19. [CrossRef] [EDP Sciences] [Google Scholar]
 He C., Mu L., Fan Z., Xu A., Zeng B., Ji Z., Han H. (2017) An improved steam injection model with the consideration of steam override, Oil Gas Sci. Technol. – Rev. IFP 72, 6, 1–14. [CrossRef] [EDP Sciences] [Google Scholar]
 van Lookeren J. (1983) Calculation methods for linear and radial steam flow in oil reservoirs, SPE J. 23, 3, 427–439. [CrossRef] [Google Scholar]
 de Haas T.W., Fadaei H., Guerrero U., Sinton D. (2013) Steamonachip for oil recovery: the role of alkaline additives in steam assisted gravity drainage, Lab Chip 13, 19, 3832–3839. [CrossRef] [PubMed] [Google Scholar]
 Das S.K., Butler R.M. (1998) Mechanism of the vapor extraction process for heavy oil and bitumen, J. Pet. Sci. Eng. 21, 1, 43–59. [CrossRef] [Google Scholar]
 Shu W., Hartman K. (1988) Effect of solvent on steam recovery of heavy oil, SPE Reserv. Eng. 3, 2, 457–465. [CrossRef] [Google Scholar]
 Edmunds N.R. (2013) Observations on the mechanisms of solventadditive SAGD processes, Proceedings of the SPE Heavy Oil ConferenceCanada, F. Society of Petroleum Engineers. [Google Scholar]
 Keshavarz M., Okuno R., Babadagli T. (2015) A semianalytical solution to optimize singlecomponent solvent coinjection with steam during SAGD, Fuel 144, 400–414. [CrossRef] [Google Scholar]
 Rabiei Faradonbeh M., Harding T.G., Abedi J. (2012) Semianalytical modeling of steamsolvent gravity drainage of heavy oil and bitumen, part 1: enhanced flow rate at mobile zone. Proceedings of the SPE Annual Technical Conference and Exhibition, F. Society of Petroleum Engineers. [Google Scholar]
 Goodman T.R. (1958) The heatbalance integral and its application to problems involving a change of phase, Trans. ASME 80, 2, 335–342. [Google Scholar]
 Closmann P.J., Smith R.A. (1983) Temperature observations and steamzone rise in the vicinity of a steamheated fracture, SPE J. 23, 4, 575–586. [CrossRef] [Google Scholar]
 Carslaw H.S., Jaeger J.C. (1988) Conduction of heat in solids, 2nd edn., Oxford University Press, New York, USA. [Google Scholar]
 Rabiei Faradonbeh M., Harding T.G., Abedi J. (2014) Semianalytical modeling of steamsolvent gravity drainage of heavy oil and bitumen. Proceedings of the SPE Heavy Oil ConferenceCanada, F. Society of Petroleum Engineers. [Google Scholar]
 Buckley S.E., Leverett M. (1942) Mechanism of fluid displacement in sands, Trans. ASME 146, 1, 107–116. [Google Scholar]
 Poling B.E., Prausnitz J.M., O’connell J.P. (2001) The properties of gases and liquids, McGrawHill, New York. [Google Scholar]
 Mandl G., Volek C. (1969) Heat and mass transport in steamdrive processes, SPE J. 9, 1, 59–79. [CrossRef] [Google Scholar]
 Chandra S. (2006) Improved steamflood analytical model, Texas A&M University. [Google Scholar]
 Wang X., Cai W., Lu J., Sun Y., Zhao L. (2015) Modelbased optimization strategy of chiller driven liquid desiccant dehumidifier with genetic algorithm, Energy 82, 939–948. [CrossRef] [Google Scholar]
 Hedden R., Verlaan M., Lastovka V. (2014) Solvent enhanced steam drive. Proceedings of the SPE Improved Oil Recovery Symposium, F. Society of Petroleum Engineers. [Google Scholar]
 Das S.K., Butler R.M. (1996) Diffusion coefficients of propane and butane in peace river bitumen, Can. J. Chem. Eng. 74, 6, 985–992. [CrossRef] [Google Scholar]
Cite this article as: H. Liu, L. Cheng, H. Xiong and S. Huang (2017). Effects of Solvent Properties and Injection Strategies on SolventEnhanced Steam Flooding for Thin Heavy Oil Reservoirs with SemiAnalytical Approach, Oil Gas Sci. Technol 72, 20.
All Tables
All Figures
Figure 1 Illustration of steam front rising stage. a) Threedimensional diagram of the stage; b) twodimensional diagram of the stage. 

In the text 
Figure 2 Illustration of steam front spreading stage. a) Threedimensional diagram of the stage; b) twodimensional diagram of the stage. 

In the text 
Figure 3 Schematic of contact area change at interface. a) Das and Butler’s perspective; b) simplification based on the Das and Butler’s perspective. 

In the text 
Figure 4 Schematic of solvent diffusion process in first stage. a) Before solvent reaches boundary; b) after solvent reaches boundary. 

In the text 
Figure 5 Schematic of solvent operating thickness. 

In the text 
Figure 6 Comparison of steam front movement between proposed model and STARS. a) Steam flooding time = 46 days; b) steam flooding time = 200 days; c) steam flooding time = 400 days. 

In the text 
Figure 7 Comparison of steam condensate front between proposed model and STARS. 

In the text 
Figure 8 Comparison of calculated oil production with STARS results. 

In the text 
Figure 9 RE distribution along cumulative oil production. 

In the text 
Figure 10 RE distribution along oil production rate and time. 

In the text 
Figure 11 Comparison of calculated oil production rate between steam flooding and SESF. 

In the text 
Figure 12 Effects of solvent operating thickness on production performance of SESF. 

In the text 
Figure 13 Effects of solvent solubility on production performance of SESF. 

In the text 
Figure 14 Effects of injection rate on production performance of SESF. 

In the text 
Figure 15 Effects of well spacing on production performance of SESF. 

In the text 