On the optimum operating temperature for steam floods

The technical, environmental and economic performances of a steam flood are partly influenced by the operating temperature (pressure). However, the definition and procedure for determining the optimum operating temperature are still debatable. Employing a combination of analytic modelling and numerical simulations, this paper investigates the existence (or otherwise) of an optimum injection temperature Topt for saturated-steam floods. Considering the maximization of productivity and thermal efficiency as objective, an analytic procedure, which explores the effects of temperature on injectivity, total steam enthalpy, oil viscosity and relative permeabilities, shows that the operating temperature (pressure) of a steam flood should not exceed 515 K (3.5 MPa). A simple closed-form expression is proposed for Topt as a function of basic rock and fluid properties. For an example three-dimensional reservoir model comprising an 8-m oil shale unit sandwiched between two sandy units each 15 m thick, numerical simulations show sensitivity to temperature (and viscosity effect) in the range 350–450 K, but becomes increasingly insensitive in the band 500–650 K. It is established that ~500–550 K is the optimum band when the optimization objective is to maximize both discounted oil recovery and cumulative oil-steam ratio. These results agree with an optimum injection temperature of ~501 K estimated from the proposed analytical model in this case. Therefore, based on the results of the analytical model, thermal simulations and other considerations, it is concluded that the optimum steam-injection temperature is project and system specific. The insights gained should find relevance in the design and management of steam floods, as well as other steam-based recovery processes.


Introduction
A combination of increasing oil demand and declining supply from conventional sources continues to fuel the interests of energy developers in unconventional oil resources such as natural bitumen. Apart from their occurrence in vast quantities [1,2], a common feature of these resources is that they are characterized by relatively high viscosities and largely immobile under reservoir conditions. As a result, they are typically associated with subeconomic natural production rates [3,4].
Owing to the ready availability of water and the high heat-carrying capacity of steam, water remains a preferred thermal agent, explaining why steam-based processes are the most popular for the in-situ development of these vast resources [5,6]. In principle, the injected heat raises the reservoir temperature, reducing in-situ oil viscosity, thus improving oil mobility and productivity. Obviously, the net heat available to the reservoir, the attainable viscosity reduction and upgrading of the in-situ oil are functions of the operating injection temperature and pressure [7].
The intensity of CO 2 emissions, non-productive thermal losses as well as fuel and other operating costs of a steam flood can partly be related to the steaming temperature [3,4,8]. Therefore, to improve the recovery, thermal and environmental efficiencies as well as the economics of a steam-based development projects, it is imperative to define and implement the optimum steam injection temperature (pressure).
From the standpoints of operational integrity as well as safety and environmental performances, shallow heavyoil and bituminous reservoirs are particularly vulnerable to high-pressure steam injection. Such reservoirs, which typically have limited or no cap rock, are not uncommon [9][10][11]. Based on thermal reservoir simulations, Al-Turki et al. [9] emphasized the criticality of pressure control when employing steam-based processes to exploit this class of reservoirs.
Leveraging on technical and practical evidence, Edmunds and Chhina [12] investigated the most economical optimum operating pressure for steam-assisted gravity drainage (SAGD) projects in northeast Alberta, Canada. They recommended a general economic optimum pressure range of 0.3-0.9 MPa (407-448 K) for SAGD projects.
In a later study, Nguyen et al. [13] confirmed that SAGD efficiency was controlled by both reservoir properties and operating conditions. Specifically, applying D-optimal design and response surface methodology while focusing on the maximization of net-present value as the objective function, they investigated the impacts of several operating variables, which included preheating, injector/ producer spacing, steam injection pressure and rate. They concluded that 5 MPa (537 K) was the optimum steam injection pressure (temperature) for SAGD in Athabasca, Canada. It is noteworthy that their conclusion was at variance with that of Edmunds and Chhina [12].
Contrary to views supporting low-temperature (pressure) injection [12], Collins [14,15] argued in favor of high-temperature operations. Although low-temperature minimizes the amount of heat injected, high-temperature minimizes the net amount of heat injected because more heat is recoverable from the produced fluids and retained in the formation. As noted by Collins [15], the beneficial effects of geomechanics are more pronounced at high temperatures (pressures). However, for a balanced view, one must be mindful of the interplay between latent and sensible heats of steam versus temperature. As the critical temperature is approached, latent heat decreases dramatically, making high-temperature steaming less optimal. Hence, for every scenario, one would expect some optimum operating temperature.
In searching for the optimum steam temperature to heat a bitumen-bearing reservoir, some workers focused on the viscosity-versus-temperature characteristics of the bitumen in question [16]. As a guide, he suggested that the optimum temperature is that at which the viscosity-versus-temperature curve goes asymptotic. This asymptotic temperature was found by taking successive gradients of the viscosity-temperature curve and establishing the temperature at which the change in gradient starts becoming negligible between successive temperatures. An obvious deduction from this technique is that the optimum temperature is fluid specific, hence should be evaluated on a case-by-case basis. Another shortcoming of this technique is the impression that the reduction of oil viscosity is the only determinant of the optimum operating temperature (pressure), downplaying the effects of other factors such as rock properties, steam flow behavior, completion type and thermal efficiency [3-5, 13, 17].
O'Dell [18] investigated the optimum steam zone pressure in steam flood projects. The author noted the detrimental effect of a low steam-zone pressure (temperature) on oil production rates, as it implies lower drawdown and limited viscosity reduction. Conversely, it was observed that operating costs reduce with lower steam pressure. Based on this work, the competing interplay of productivity and operating costs suggests the existence of an economic optimum steaming pressure and temperature. However, it is instructive that the author did not conclude on a universal optimum point but recommended that a cost-benefit analysis should always be conducted for any steam flood project.
Using artificial neural network, Kam et al. [19] studied the optimal strategy to establish the operating pressures (temperatures) that would maximize the economic value of solvent-aided thermal recovery projects. The objective function was set up as a multi-variable problem, with the optimization variables including bitumen production, steam injection, solvent retention and oil price. However, application of their work to Athabasca oil sands did not reveal the existence of an optimum operating pressure, but rather proposed an optimal schedule of different operating pressures.
Recently, Tian et al. [6] published their findings from a set of three-dimensional (3D) laboratory experiments on SAGD feasibility in high-pressure environments. Considering oil recovery factor, cumulative oil-steam-ratio and the quality of produced oil as key performance indicators, they concluded that SAGD was not feasible in high-pressure environments. They recommended that overpressured heavy-oil reservoirs should be depressurized prior to executing a SAGD process.
In this paper, we investigate the optimum injection temperature for steam floods. Taking advantage of fundamental principles and considering multiple factors, we interrogate the existence (or otherwise) of an optimum steaming temperature. From these theoretical investigations, a simple analytical model is proposed for estimating optimum temperature for saturated steam floods. Afterward, the robustness of this model is examined using 3D, multiphase numerical simulations of a SAGD process.
We are not unaware of recent efforts to achieve dramatic improvement in the recovery efficiency of steam floods by injecting supercritical steam (i.e. above 647.5 K and 22.1 MPa) for the exploitation of heavy oil and bitumenbearing reservoirs [7,20]. However, to keep the problem tractable and limited to the main recovery mechanisms of viscosity reduction and alteration of relative permeabilities, this paper focuses on traditional steam floods, in which sub-critical steam is utilized. This notwithstanding, where supercritical temperature (650 K) is considered, this study ignores complex recovery mechanisms such as solvent action, decoking and pyrolysis, that have been linked to such conditions [7,21].
This paper is organized into three major sections. The theoretical framework and formulation of analytical models are presented. This is followed by discussions on the parametric tests conducted on the new analytical model and the results obtained. Insights from model validations with numerical thermal simulations are discussed. The final section presents the key conclusions and limitations of this work, while offering future directions for research that could advance the current body of knowledge on this subject.

Theoretical framework
In this section, we examine some rock and fluid properties that are known to influence the performance of a steam flood. Leveraging on the temperature-dependency of such properties, we seek a functional relationship between a given performance indicator and the specific rock and fluid property of interest. The same concept of analyzing the temperature dependency is extended to the total steam enthalpy, which influences the thermal efficiency of a steam flood.
Effect of oil viscosity only: In this case, we consider the basic principle of a thermal flood, which is to accomplish significant reduction of in-situ oil viscosity by elevating the temperature. Accordingly, we explore the functional relationship between oil viscosity and temperature to find a possible optimum steam injection (operating) temperature.
Considering the common two-parameter Andrade viscosity model, oil dynamic viscosity and temperature are related as follows [22].
where μ oi and a are empirical constants measured in Pa.s and K, respectively. T is absolute temperature (K). μ o is the dynamic viscosity (Pa.s) at the prevailing temperature T.
In seeking the optimum injection temperature from the standpoint of oil viscosity only, we are looking for that injection temperature that always minimizes oil viscosity in the reservoir. Mathematically, this search requires setting the temperature derivative of Eq. 1 to zero and solving for T i.e.
Solving Eq. 2 yields the solution T = ∞ . This result suggests the non-existence of a unique and physically realistic temperature that guarantees minimum in-situ oil viscosity. Therefore, from the viewpoint of oil-viscosity effect only, it does not appear that there is a realistic optimum injection temperature for steam floods.
Effect of steam injectivity only: This approach considers the impact of temperature on steam injectivity. Given the simple case of a radial incompressible two-phase flow in a porous medium, the injectivity of steam at steady-state conditions is expressed as follows.
To simplify our search for an optimum temperature from a steam-injectivity viewpoint, we assume that only the quantity μ s is temperature-dependent in Eq. 3. As a result, we use the following correlation of Tortike and Farouq-Ali [23] to describe saturated steam viscosity (Pa s) versus temperature (K). This correlation is valid within the range 273.15≤T≤645K.
Following the application of chain rule to Eqs. 3 and 4 i.e. Again, solving Eq. 5 does not yield a unique and physically realistic temperature as the optimum. Therefore, from the standpoint of steam injectivity, there is no unique optimum temperature.
Effect of total steam enthalpy: By summing expressions for specific latent and sensible enthalpies, the specific total enthalpy (in J kg −1 ) of a saturated steam in the temperature range 273.15≤T≤645 K, is given by the following correlation [23].
In the search for a temperature that optimizes the total enthalpy carried by the injected steam, we derive the equation dE ∕ dT = 0 and solve same for T. Fig. 1 displays the resulting dE ∕ dT vs. temperature curve for saturated steam. From this plot, it is deduced that ~515 K is the temperature that yields the optimum specific total enthalpy of injected steam. Therefore, strictly from the standpoint of thermal efficiency, it can be concluded that the operating temperature (pressure) of a steam flood should not exceed 515 K (3.5 MPa). At higher temperatures, there is an increasing risk of operating a suboptimal steam flood, whose thermal efficiency may not be competitive. This increased risk can be attributed to two effects vis-à-vis (i) lower injected total enthalpy at higher temperatures, and (ii) higher thermal losses due to larger temperature difference between injected steam and the surroundings, including along surface injection and production networks, wellbores and adjacent formations [18,24,25].
Combined effects of oil viscosity and relative permeability: In this case, we examine the impacts of temperature on oil productivity (mobility), as influenced by both in-situ viscosity and relative permeability. Given the simple case of a radial incompressible two-phase flow in a (6) porous medium, the productivity of oil at steady-state conditions is expressed as follows.
where J o = oil productivity index (m 3 s −1 Pa −1 ), k ro = oil relative permeability (dimensionless), and B o = oil volume factor (dimensionless). For convenience, we restrict the current treatment of temperature effects on J o to the examination of the mobility terms k ro and μ o only. The temperature-derivative of J o requires applying the following chain rule to Eq. 7.
While the resolution of the second terms on the righthand side of Eq. 8 are straight-forward, the quantity dk ro∕ dT in the first term is more demanding.
For mathematical convenience, we assume an arbitrary relative-permeability model for a two-phase oil-water system.
in which k rom = end-point oil relative permeability at connate water saturation (dimensionless), S w = water saturation (dimensionless), S iw = irreducible-water saturation (dimensionless), and n o = oil relative-permeability exponent (dimensionless). Note that the quantities k ro , k rom , S w and S iw are physically bounded to the range 0-1.Based on the efforts of different generations of researchers, the sensitivities of relative permeability to temperature have been established for both sandstone and carbonate formations [26][27][28]. More important, though to varied degrees, several researchers reported and applied temperature effects on the end-point saturations (either residual oil or connate water), hence oil relative permeability [26][27][28][29].
Based on the foregoing experimental and theoretical studies, we employ the following semi-empirical relationship between S iw and temperature [28]. By deductive argument, the same expression can readily be used to derive a model of k ro vs. temperature, which is our current interest (see Eq. 8). Fig. 1 dE ∕ dT vs. temperature plot for saturated steam, suggesting that the heat-carrying capacity of a steam flood becomes sub-optimal above 515 K (3.5 MPa) in which the quantities α and β are empirical constants. They are measured in K -1 and dimensionless units, respectively. Considering Eqs. 9 and 10 while assuming that k rom , n o and S w are not sensitive to temperature changes, we apply the following chain rule. It is worthy of note that the validity of S w∕ T = 0 has been demonstrated by Nakornthap and Evans [28].
Combining Eqs. 1, 7 and 9-11, then substiting in Eq. 8 yields the temperature-derivative of J o . Setting this derivative to zero, we have Solving Eq. 12 for T yields the following expression for the optimum injection (operating) temperature (in K) which, in principle, maximizes oil productivity under the synergistic effects of reduced oil viscosity and increased oil relative permeability.
For practical applications, the realistic values from Eq. 13 should be limited to the range 273.15≤T opt ≤645 K. This is the practical range of temperature of a saturated steam.
To ensure that Eq. 13 always returns real solutions for T opt , we impose additional constraints on the discriminant i.e.

Parametric tests
For the purpose of illustration, Figs. 2, 3, 4 and 5 show the results of some parametric tests conducted with Eq. 13. From these plots, the optimum temperature is seen to be a function of the rock and fluid properties. It is worthy of note that these results suggest that the optimum temperature depends on the characteristics of the system in question. This deduction is in contrast with the conclusions of some workers [12,13], who recommended 537 K (5 MPa) and 407-448 K (0.3-0.9 MPa), respectively as the (10) S iw = T + , general optimum operating conditions for steam-flooding projects.

Model validation
A 3-D semi-homogeneous multiphase reservoir model is constructed with STARS, which is a commercial thermal simulator developed by the Computer Modelling Group. STARS is a standard industry software for numerical modelling of enhanced recovery processes involving displacing with fluid systems such as steam, air, solvents or chemicals. To evaluate the robustness of our proposed analytical model (Eq. 13), we employed the 3D numerical model to investigate the existence of an optimum injection temperature for steam floods. As an example, a sector of the heavy-oil and bitumen belt in the Eastern Dahomey Basin, southwestern part of Nigeria is considered [16,30]. A layer-cake stratigraphy is assumed with an 8-m oil shale unit sandwiched between two sand-rich units, each with a thickness of 15 m. A uniform average porosity of 30% is considered while oil saturations are 68 and 80% in the oil shale and sandy units, respectively. The oil shale and sandy units are characterized by average horizontal permeabilities of 2000 and 400 mD respectively, while all the units have a uniform vertical anisotropic ratio of 0.5 [16]. There is no initial free gas and aquifer. The reservoir is considered water-wet and hydrostatic, while relativepermeability functions are of Corey type. For completeness, the relative-permeability models account for temperature dependency. The reservoir model, characterized by a dimension of 200 m x 500 m x 38 m, is described by a total of 7600 simulation cells (Fig. 6). While fine scale models are required for full-field simulation studies, relatively coarse layer-cake models are often suitable for preliminary screening of steam flooding projects, assessment of uncertainties as well as history-matching experimental and actual steam-injection projects [31,32]. Each simulation cell is of dimension 40 m x 12.5 m x 1 m.
To keep the problem simple and tractable, the fluid system is characterized by pseudo-components. These are a series of hypothetical close-cut fractions spliced from the boiling range distribution of crude oil fractions above carbon number 6 as a substitute to full compositional analysis. Its viscosity-temperature relationship is described with the two-parameter Andrade model (Eq. 1), assuming μ oi = 7.0 x 10 −9 Pa.s. While mass and heat exchanges occur through the layers overlying and underlying the oil shale, there is neither heat nor mass transfer across the boundaries of the model. At model initialization, the stock-tank oil initially in-place (STOIIP) is 9.1 x 10 5 m 3 .
The well models and configurations follow a standard SAGD process. The model comprises two horizontal injector-producer pairs, placed some 100m apart. As displayed on Fig. 6, the injectors are Inj1 and Inj2, while the corresponding producers are Prod1 and Prod2. Each well has a completion length of 500 m, radius of 0.16 m and perforated throughout their length. Both injectors and producers are characterized by a total skin factor of 5, assumed constant throughout the simulation times. The producers are placed 2 m above the base of the reservoir while the injectors are 5 m above the producers. Saturated steam is injected at 70% quality, and all the wells have 100%  Table 1.
Liquid constraints of 250 m 3 /D are imposed on the production wells, while the injection wells have a coldwater equivalent (CWE) steaming rate of 280 m 3 /D. By design, the steady-state voidage-replacement ratio is 1.1. Wellbore hydraulics are ignored. The producers and injectors are constrained to bottom-hole pressure (BHP) limits of 1000 and 4000 Pa, respectively. Limiting production water cut is 90%. In all simulation cases, geomechanical effects are ignored, and there is no provision for pre-heating of the reservoir. Initial reservoir temperature is fixed at 316 K.
The optimization study is set up to achieve two main objectives. These are (i) maximize oil recovery; and (iii) minimize cumulative steam-oil ratio (CSOR), over a twenty-year injection and production period. To account for time effects occasioned by differences in productivity and injectivity due to varied injection temperatures, this study includes discounted oil recovery and CSOR. Annual discount rate is 10%. In this work, discounted CSOR is simply the ratio of discounted oil recovery to the discounted cumulative volume of steam injected (in CWE). Therefore, as framed in this work, we seek the optimum operating temperature as that which best accomplishes all these objectives.
Oil recovery and CSOR as performance indicators: Considering absolute oil recovery and CSOR, performance plots for the production and injection wells are displayed in Figs. 7 and 8, respectively for various cases of constant injection temperatures in the range 350 -650 K. In all cases, the STOIIP is 9.0 x 10 5 sm 3 . As would be expected, oil recovery increases over time, until a limiting state is reached. Compared to other injection conditions, no appreciable oil recovery is achieved in the case of steaming at 350 K. This poor performance is attributed to the high viscosity (~3400 mPa.s) of in-situ oil at 350 K. To put this in context, the corresponding in-situ oil viscosities at 400, 450 and 550 K are approximately 280, 40 and 2 mPa.s, respectively. On the other hand, there is no material difference in the recovery rates and ultimate recoveries within 2.0 7 injection temperature (T inj ) cases are considered i.e. 350, 400, 450, 500, 550, 600 and 650 K. For clarity, we maintained each injection temperature throughout the simulation case, as against the optimization practice of gradually reducing injection temperature during the wind-down phase or substituting steam with cheaper thermal agents such as CO 2 and flue gas [33,34]. Use of these input data in Eq. 12 suggests 500.6 K as the optimum temperature for this example. This T opt <515 K, falls within the optimum envelope for total steam enthalpy ( Fig. 1) α (K −1 ) 0.002455 β −0.05325 T opt (K) 500.6 Fig. 7 Temporal variation of cumulative oil production for various steaming temperatures Fig. 8 Temporal variation of CSOR for various steaming temperatures the band 500-650 K. Clearly, viscosity effect is the main discriminator of the performances of the temperature band 350-450 K versus the range of higher injection temperatures in this example. Similar to the deductions reached from Fig. 7 in terms of production rate and ultimate recovery (UR), Fig. 8 suggests that the same two broad categeories of performance can be observed in this example problem. Considering CSOR as a performance indicator, system dynamics exhibit clear sensitivity to temperature in the range 350-450 K, but this sensitivity becomes negligible in the band 500-650 K. Given that CSOR is a reaosnable indicator of energy and environmental performances, it follows that the optimum steaming temperature is in the interval 500-650 K, within which energy and envonmental effciencies are maximized per unit volume of oil produced.
Effects of discounting: For better insights into the effects of time on the results, we present Figs. 9 and 10 as the absolute and discounted values of the performance measures, respectively. Fig. 8 shows that the upper limit of UR is essentially reached between 500 and 550 K, while the same interval also yields the lower limit of CSOR. Beyond this temperature range, no significant increase in oil UR and reduction in CSOR are obtained in this example problem. It is interesting that the use of discounted UR and CSOR, as displayed in Fig. 10, substantiates the optimality of the temperature band 500-550 K in this example.
It is worthy of note that the current thermal simulations account for limited heat losses to the surrounding formations vis-à-vis the overburden and underburden. As a result, the adverse effects of higher injection temperatures such as 550, 600 and 650 K on CSOR are not fully quantified in these cases. As demonstrated by other workers, such thermal losses to the surrounding rocks in the reservoir as well as the enviornment in the course of steam transport along the injection network can be significant and detrimental to the performance of a steam flood. These losses constitute another disincentive to high-temperature steaming operations [11,24,25,33]. This argument underscores the 500-550 K band as optimum for this example problem.
In comparison, using the same set of input data for the example problem in our analytical model yields an optimum injection temperature of 500.6 K ( Table 1). This estimate agrees with the optimum band 500-550 K obtained from detailed thermal simulations. Furthermore, this analytical model estimate is consistent with the theoretical maximum of 515 K obtained strictly from the standpoint of total enthalpy. Therefore, based on the results of the analytical model, thermal simulations and the upper limit of 515 K from total enthalpy, it is concluded that the optimum steam-injection temperature for the example case investigated is ~501-515 K.
Optimum operating temperature for productivity and thermal efficiency: The problem of determining the optimum operating temperature for a steam flood can be viewed from two broad perspectives. Based on the insights gained from this study, the two perspectives are related to the maximization of oil productivity and thermal efficiency (heat-transport capacity of steam). As shown in this study, the influence on thermal efficiency is generic because it is linked to total steam enthalpy, which is a universal property. In other words, providing the steaming temperature does not exceed 515 K, the operating condition of a steam flood remains optimum from the standpoint of thermal efficiency. Because of the positive correlation between energy consumption and environmental impacts, it follows that the same optimum temperature (< 515 K) is  favourable to achieving favourable environmental efficiency. In this context, environmental performance should be viewed in relation to the amounts of water and fuel consumed to generate steam, CO 2 emissions associated with steam generation as well as the amount of produced water handled per unit volume of oil produced.
Conversely, the optimum injection temperature required to maximize oil productivity is system specific. This optimum is controlled by the temperature characteristics of the rock and fluid properties, as described by the applicable relative-permeability and oil-viscosity models, respectively. Therefore, while the consideration of thermal (and environmental) efficiency yields 515 K as the upper limit of optimum operating temperature for a steam flood, the optimum range is further narrowed by the specific relative-permeability and oil-viscosity vs. temperature characteristics of the system of interest. While 515 K is the upper bound of optimum operating temperature for steam floods, the specific optimum operating condition would have to be evaluated on a case-by-case basis. This logical conclusion is an improvement over the positions advanced by Nguyen et al. [13] as well as Edmunds and Chhina [12], who recommended general optimum temperatures of 537 and 407-448 K for steam floods, respectively. We submit that such general recommendations by these previous workers do not appear to reflect the diverse and complex characteristics of the various heavy-oil and bitumen-bearing reservoirs, as well as the 515 K shown to be the universal upper limit to preserve the total enthalpy of saturated steam.

Conclusion
To aid the design and management of steam-flooding projects, a new semi-analytical model has been developed for conducting first-order estimation of the optimum temperature (pressure) at which oil productivity is maximized as a function of in-situ viscosity and relative permeability. Sensitivity tests suggest that it is sub-optimal to operate saturated-steam floods above 515 K (3.5 MPa).
As a robustness check, 3D multiphase thermal simulations were performed. Simulation results show that production rate, ultimate recovery and CSOR are sensitive to injection temperature in the range 350-450 K. However, for the cases examined in this work, these performance indicators do not exhibit strong dependency on injection temperature in the range 500-650 K.
A close examination of discounted oil recovery and CSOR from the simulation runs reveals an optimum injection temperature band of 500-550 K, which compares reasonably well with the ~501 K (< 515 K) from the new semianalytical model. It is noteworthy that these numerical and analytical results differ from the submissions of Nguyen et al. [13] as well as Edmunds and Chhina [12], who recommended general steam-flood optimum temperatures of 537 and 407-448 K, respectively. For practical applications, we recommend that the optimum operating conditions of a steam flood should be treated on a case-by-case basis, while the proposed model is used for guidance.
From the standpoint of thermal efficiency, the operating temperature (pressure) of a saturated-steam flood should not exceed 515 K (3.5 MPa). Otherwise, there is a non-zero risk of sub-optimal steam flooding, owing to lower injected total enthalpy and higher thermal losses to the surroundings. Operating steamfloods at an optimum temperature (< 515 K) is also favourable to the environment, as it offers maximum oil recovery while minimizing fuel costs, CO 2 emissions as well as the net amounts of water consumed and produced per unit oil recovery.
One limitation of this study is that the new analytical model does not account for potential geomechanical effects. In addition, it is worthwhile to examine the optimum operating temperatures and pressures for non-saturated steam floods.

Compliance with ethical standards
Conflicts of interest The authors declare that they have no conflict of interest.
Availability of data and material All the data pertaining to this study are available upon request.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.