Effects of Solvent Properties and Injection Strategies on Solvent-Enhanced Steam Flooding for Thin Heavy Oil Reservoirs with Semi-Analytical Approach

Compared with conventional steam flooding and Steam-Assisted Gravity Drainage (SAGD), Solvent-Enhanced 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 semi-analytical 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 multi-solvent 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.


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 steam-flooding strategy has been regarded as a significantly superior technology compared to Steam-Assisted Gravity Drainage (SAGD) operations in terms of energy efficiency [3] and has been applied widely in thin heavy oil reservoirs in China [4][5][6].However, significant heat losses to the understrata and overburden can render the steam-flooding process uneconomic because of the small thickness of the reservoir [7].Solvent-Enhanced 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 three-dimensional 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 steam-CO 2 injection strategies in a scaled model and revealed that steam-CO 2 continuous injection results in improved performance compared to steam-CO 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 well-designed 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 solvent-steam process.The results indicate that mediumalkane solvents provide the greatest increase in total oil production, whereas heavy-alkanes 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 steam-solvent 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 time-consuming 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 semi-analytical 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 co-injection chamber.Gupta and Gittins [13] developed a semi-analytical approach to estimate oil drainage in the co-injection 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 semi-analytical model is proposed based on mass and energy balance principles, diffusion theory, and the theory of Buckley-Leverett.First, the mathematical model is solved in a semi-analytical 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.

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: 1. 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 condensate-producing stage).The pressure inside the steam zone is constant during all stages [16,17]; 2. 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 steam-flooding 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; 3. 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; 4. the gas-phase 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 VAPor-EXtraction (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; 5. 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 co-injecting with medium solvents.Therefore, the solvents used in this model are assumed to be medium solvents and a combination of medium and light solvents; 6. 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 rock-fluid 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 condensate-producing stage).These are introduced in detail in the following.

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: where H o is the heat injection rate per unit well length; A is the area of the steam front; q r , q o , q s , and q 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 unsteady-state solvent concentration distribution along the steam front can be described by Equation (3) [24]: The corresponding boundary conditions are: where e 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 e; D is the solvent molecular diffusion coefficient; a 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: where d 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: 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: where A d c is the solvent diffusion area at time t and " C s 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.
see equation (10) at the bottom of this page

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: 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: According to assumption (2), the geometry of the water front is assumed to be a diameter-constant circular shape; thus, the change of the heated area can be transformed to the change of distance: Carslaw and Jaeger [27] built and solved the heat loss model of a semi-infinite plate with a constant temperature boundary.As the hot water mainly contacts the understrata and overburden, the heat loss rate is written as: where T s and T w are the temperature of the steam and hot water, respectively; q cap is the density of the understrata and overburden rock; k 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: Based on the mass conservation law, the generating rate of the cold steam condensate per unit well length can be given as: where v m c 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: Shutler and Boberg [11] proposed that piston-like flow through a thin layer is assumed; therefore, the isothermal two-phase 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: where f 0 w ðS wf Þ 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): where L e is the effective operating thickness of the solvent; " S w 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.
The solvent diffusion rate and oil production rate can also be calculated by Equations ( 6) and (7).

Steam Condensate-Producing Stage
The steam condensate-producing 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]: where l o , l s , and l 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: The oil production rate can be rewritten from Equation ( 7) by considering the water and solvent production: where S wp is the water saturation at the production well, which is determined by f 0 w ðS wo Þ and l m , and f w ðS wp Þ is the fractional flow of water at the production well, which can be acquired by considering l m and S wp .

Model Treatment
A semi-analytical 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): 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: where

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: 1. 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 dd c /dt, and solvent penetration depth d c ; 2. the production rate in the first stage is calculated using Equations ( 7)-( 10); 3. repeat steps 2 and 3 substituting d 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; 4. the steam front velocity dx/dt and the steam front displacement x in the second stage are calculated using Equations ( 15) and ( 27); 5. the cold steam condensate generating rate v m c is calculated using Equations ( 16)- (18).Then, Equations ( 5), (6), and ( 19)-( 22) are coupled and solved for the solvent dd c /dt, d c , and the cold steam condensate distribution beyond the edge of the steam front; 6. the production rate in the second stage is also calculated using Equations ( 7)-(10); 7. repeat steps 5-7 substituting d c calculated in step 6 into Equation ( 6) until the cold steam condensate reaches the production well, when the process enters the third stage; 8. v m c , dd c /dt, and d 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); 9. repeat steps 9 substituting d c calculated in step 9 into Equation ( 6) until f w (S wo ) = 1, when the process ends.

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.
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 Buckley-Leverett theory in the steam-flooding process in a thin layer.
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].RMSRE (Rooting-Mean-Square 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]: where NPTS is the number of data points; y cal i is the calculated value at point i by the proposed model; and y sim i 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 steam-flooding 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 one-dimensional assumption of the Buckley-Leverett 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.
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 9 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 9 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.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 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.

Effects of Solvent Properties on 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.Multi-solvent 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.RE distribution along oil production rate and time.

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 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.
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 semi-analytical 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: 1. 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; 2. 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; 3. 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 multi-solvent injection using different values of operating thickness and solubility; 4. 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.

Figure 1 Figure 2
Figure 1 Illustration of steam front rising stage.a) Three-dimensional diagram of the stage; b) two-dimensional diagram of the stage.

Figure 3
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.

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

Figure 5
Figure 5Schematic of solvent operating thickness.

Figure 12 Figure 6 Figure 9 RE distribution along cumulative oil production. 7 Comparison
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

Figure 11 Comparison
Figure 13Effects of solvent solubility on production performance of SESF.

Figure 14 EffectsFigure 15 Effects
Figure 14Effects of injection rate on production performance of SESF.

TABLE 1
Average property parameters of X oil field in China.