Introduction

Thermal enhanced oil recovery (TEOR) is a family of tertiary processes defined as “any process in which heat is introduced intentionally into a subsurface accumulation of organic compounds for the purpose of recovering fuels through well.”

Thermal methods of EOR entail the application of heat to the oil well. Thermal methods have been tested since 1950s, and they are the most advanced among EOR methods, as far as field experience and technology are concerned. They are best suited for heavy oils (10–20° API) and tar sands (≤ 10° API). They have the ability to lower the viscosity of the oil and thus increase the mobility ratio. These methods are typically employed in relatively shallow oil wells with higher viscosity such as tar sands and heavy oil. Thermal methods of EOR have been highly successful in the USA, Canada, Venezuela and Indonesia and have also been used in China and Brazil. Within the USA, thermal methods account for around 40% of EOR production (Prats 1982).

Types of thermal EOR

  • Steam flooding.

  • Cyclic steam stimulation (CSS).

  • Steam-assisted gravity drainage (SAGD).

  • In situ combustion (ISC).

  • Solar EOR.

  • Electro-EOR (EEOR).

  • Heated annulus steam drive (HASD).

  • Steam and solvent processes.

FNE field is located about 9 km away from the North East of Fula field in Fula Sub-basin in Sudan; it is structurally subdivided into three main blocks, which are FNE-1, FNE-3 and FNE-N. There are three proved hydrocarbon-bearing formations (Aradaiba, Bentiu and AG), among which Bentiu is the main pay zone. Figure 1 shows the location of FNE.

Fig. 1
figure 1

FNE location map (Elbaloula et al. 2016)

FNE oil field is considered as heavy oil field and it has shallow reservoir, and Bentiu reservoir is affected by strong bottom water. The reservoir of FNE oil field is characterized by normal temperature and pressure system at 529 m depth. The average pressure is 576 psi and the average temperature is 43.9 °C. FNE has conventional heavy oil in both Aradaiba and Bentiu, and Bentiu reservoir has a relatively good quality with porosity ranging from 17 to 40%, permeability ranging from 100 to 11000 mD and viscosity of about 840 cP (Elbaloula et al. 2016).

The oil field put into production in 2010 with cold production and cyclic steam stimulation. By March 2018, a total of 117 wells had been drilled, including one horizontal well; 116 wells have been put into operation, of which 21 wells are producing as cold, 82 wells for steam stimulation and 13 wells for steam flooding. One hundred and eleven wells were opened with a daily oil production of 489 6bbl, a daily fluid production of 10,145 bbl, water cut of 54%, the total original oil in place of 332MM STB, and the up-to-date recovery factor of reservoirs of 0.8%.

Literature review and screening

Years of production shows that the daily oil production for steam stimulation is almost three times of the cold wells, indicating that injection steam into well bottom can heat the reservoir and reduce viscosity of the heavy oil to increase production capacity. Currently, the oil fields used for steam flooding include China’s Liaohe oil field, US Kern River oil field, Indonesian Dori oil field and so on (Jones et al. 1995).

From screening, the Qi 40 of Liaohe oil field which put into steam flooding has been selected early as a reference example to illustrate the feasibility of steam flooding in FNE oil field. Qi 40 reservoir was put into operation in 1983 by steam stimulation with oil take-off rate of 15.3%. During the period of 1990–2000, the steam flooding test was carried out. In 2001, intermittent steam flooding was officially started. The well pattern was anti-9 spot with well space of 100 meters. As of December 2003, the oil take-off rate increased to 41.4%, highlighting the superiority of steam flooding for heavy oil reservoirs (Jiu-quan et al. 2006).

By comparing the reservoir properties for each oil field and the applicability of steam flooding and comparing SPE table for EOR screening (Taber et al. 1997), as given in Table 1, it is considered that the Sudan FNE oil field is suitable for steam flooding, and this paper well covers the steam flooding screening, pilot area selection, pilot simulation and designing as well as optimization.

Table 1 Reservoir properties of each oil field using steam flooding (Elbaloula et al. 2016; Pearce and Megginson 1991; Jones et al. 1995; Taber et al. 1997)

Methodology

The geological data, reservoir data and production data for FNE oil field have been collected and used for analysis to identify the situation of the field, and it is suitable for conducting steam injection and full-field implementation.

Advanced thermal EOR simulator software has been used to build the dynamic model and to design and optimize steam flooding injection and production parameters for FNE oil field.

The steps in Fig. 2) are used in this paper.

Fig. 2
figure 2

The steps and workflow for this study

Results and discussion

The thermal recovery is essential for this field, and the first pilot has been conducted in FNE-16 well and the results showed that the CSS has the ability to double the production, and then, additional wells have been added at each phase. Up to 2016, the total CSS wells are 67 including 37 wells are under the first and second cycles, 24 wells under the third and fourth cycles and 6 wells under the fifth cycle.

The average daily oil production for this field has been increased from 5300 bbl./day as of December 2014 to 8300 bbl./day as of September, 2016, and the peak oil production has been recorded in 2016 as 9000 bbl./day. The CSS is very successful, and the average oil rate is almost 1.6 times compared to cold production. The CSS only can increase the recovery percent from 32.5 to 34.2%, which makes it more attractive as development scenario for FNE oil field (Elbaloula 2017).

The steam flooding pilot has been implemented in September 2015, and the results encourage to go for the second steam flooding pilot in this field in order to get full picture of full-field implementation as given in Table 2.

Table 2 Screening criteria for thermal recovery and development (Elbaloula et al. 2016)

Pilot area selection analysis

In order to select the best location for pilot, screening of FNE oil field with steam flooding has been conducted; the following factors are considered:-

  1. 1.

    The pilot area must have abundant reservoir resources, and then, the properties can represent the general level of the whole oil field.

  2. 2.

    The sand body between the producer and the injector must have good connectivity and a uniform oil–water contact, and the interlayer should be fully considered to make sure that the injected steam plays a more desirable role.

  3. 3.

    At present, the daily oil production of some wells is in a decreasing stage, and the water cut is relatively high. Dynamic analysis needs to be conducted to assist optimization.

  4. 4.

    The producers have experienced at least one cycle of high-temperature steam injection and soaking, and even four–six cycles. The cementing facility of some wells becomes damaged. Therefore, it is necessary to check the well cementing and propose workover before steam flooding.

Block 1 is located at a high structure with 16 wells, including 2 cold wells and 14 steam simulation wells. Two cold wells and three CSS wells were put into production in 2010 and 11 steam simulation wells in 2015. The main perforation layers are B1a and B1b. Block 2 is located in a relatively low structure with 11 wells, including three cold wells commissioned in 2010 and eight steam wells commissioned in 2015.

According to above screening principle and considering the block structure, the effective thickness of the reservoir (Fig. 3), the cementing quality and the current production status, the selected block 1 is selected as the steam flooding pilot area to carry out the next numerical simulation and development plan prediction.

Fig. 3
figure 3

Reservoir properties and pilot area

Building the numerical model

Advanced thermal process reservoir simulator has been used which includes options such as chemical/polymer flooding, thermal applications, steam injection, horizontal wells, dual porosity/permeability, directional permeability, flexible grids, fireflood and many more.

The pilot area has been cut as sector model from static model and then used for initialization, history matching and prediction using advanced thermal EOR simulator. The grid number of model is 49 × 48 × 119 = 279,888, cell size is dx = dy = 20 m, dz = 1 m, and vertical formation includes Ara-D, AE, AF, B1a, B1b, B1c, B1d and B1e. Figure 4 shows 3D view of porosity distribution and wells location in FNE oil field.

Fig. 4
figure 4

3D view of porosity distribution and wells location in FNE oil field

After finding the difference between calculated production and observed production, some parameters of the dynamic model are modified. The process is continuous in order to make dynamic model match closely to the actual reservoir performance and ensure the reliability of the remaining oil distribution.

History matching

During history matching of oil field, the matched parameters are pressure, liquid rate, oil rate and water cut. The reservoir pressure and production of the field were firstly matched. Then, the flowing bottom-hole pressure and production of single well were matched (Fig. 5).

Fig. 5
figure 5

History matching result of block and single well

For the history matching result, 10 wells taking 50 percent of SF2 pilot area are well matched, 5 wells are fairly matched and 1 well is poorly matched. The history matching fitting rate is 94%.

Scenarios design for pilot test well group

Essentially four types of well arrangements are used in fluid injection projects:

  1. 1.

    Irregular injection patterns.

  2. 2.

    Peripheral injection patterns.

  3. 3.

    Regular injection patterns.

  4. 4.

    Cristal and basal injection patterns.

A wide variety of injection–production well arrangements have been used in injection projects. The most common patterns, as shown in Fig. 6, are the following (Ahme 2018):

Fig. 6
figure 6

Flood patterns (Ahme 2018)

  • Direct line drive.

  • Staggered line drive.

  • Five spots.

  • Seven spots.

  • Nine spots.

In terms of cumulative oil production, Eq. 1 can be written as:

$${\text{NP}} = {\text{NS}}\;{\text{ED}}\;{\text{EA}}\;{\text{EV}}$$
(1)

Applying the above methods for this field, it has been found that totally four different well groups including 5 spots, 7 spots and 9 spots can be optimized (Fig. 7). In these cases, liquid rate is the last point of time; injection rate is 100 m3/day. For the forecast result, taking the new wells needed and economic condition into account, option 3 (7 spots) including 6 new producers is recommended in this case (Table 3).

Fig. 7
figure 7

Steam flooding well group options

Table 3 SF designed and forecast result of each option

Optimization

The optimization is done for the injection parameters such as injection rate, injection volume, injection pressure and temperature, also for the timing of converting from CSS to steam flooding besides the production parameters as well.

Timing of converting to steam flooding

There are four cases designed in order to know the timing of converting to steam flooding in this pilot area, including converting in January 2020, January 2021, January 2022 and January 2023. For the forecast result of these four cases, the oil recovery factor decreases along with the converted timing delay (Table 4, Fig. 8). Based on the temperature distribution in 2020 and 2021, the heat connectivity is already formed and most of the wells are in the fourth cycle until January 2020. According to experience of SF#1 block, SF#2 pilot area could be converted for steam flooding in January 2020.

Table 4 Forecast index of converting to SF timing cases
Fig. 8
figure 8

Cumulative oil and injection steam along with injection timing delay

Liquid rate and injection rate

In this case, the producers’ product intensity is 5–40 bbl/day m according to the production history of each well and the injector’s injection index is 20–70 bbl/day m same as SF#1 (Tables 5, 6). Based on reasonable liquid–steam ratio which ranges from 0.8 to 2.0, ten cases are selected for optimization (Table 7).

Table 5 Summary of product intensity for producers
Table 6 Summary of injection index for injectors of SF#1
Table 7 Liquid rate and injection rate of each well in selected ten cases

For the forecast result of these ten cases, the oil recovery factor is increasing progressively along with the injection rate increasing (Table 8). Boiler condition is as follows: Steam outlet temperature is 300 °C, boiler efficiency is 88%, and if not considered the heat loss, the enthalpy of oil is 41900 kJ/kg. Enthalpy of saturated steam under 300 °C is 2750 kJ/kg, which means 0.5 barrel oil exchanges 1 m3 steam. For example, compared with case 6, case 10 needs to spend 0.65 MMB oil to get 1.3 MMm3 steam. Considering the economic problem (water cut 98%), case 10 needs to get net oil at least 0.276 MMbbl. So it needs to add oil 0.926 MMstb totally; however, it only needs to add oil 0.73 MMbbl (Table 9).

Table 8 Liquid rate and injection rate of each well in selected ten cases
Table 9 Liquid rate and injection rate of each well in selected ten cases

In conclusion, the liquid rate and injection rate of case 6 are chosen, producers’ product intensity is 20 bbl./day m per well, steam injection index is 7 m3/day m per well, and the average injection rate is 96 m3/day per well.

Steam quality and steam temperature

There are totally 18 cases in this case, and the steam quality of injectors ranges from 0.6 to 0.9 and temperature ranges from 200 to 340 °C (Table 10). According to the forecast result, recovery factor of each case is different, but the difference is very small (Table 11). Taking the rationality of steam temperature and quality into consideration, the current case (case 8) is chosen although the recovery factor of case 12 is the highest.

Table 10 Steam quality and seam temperature optimization cases
Table 11 Forecast result of steam quality and steam temperature optimization cases

SG intermittent open

There are totally 18 cases in this case, and the SG stop time ranges from 1 to 4 weeks and SG open time ranges from 1 to 6 weeks (Table 12). According to the forecast result, recovery factor of case 15 (open 6 weeks and stop 1 week) is the highest, but the value of current case (no stop) is higher (Table 13). Finally, the current case (no stop) could be kept.

Table 12 SG open and stop times optimization cases
Table 13 Forecast result of SG open and stop time optimization cases

Product layer and injection layer

There are totally 15 cases in layer optimization case, and oil zones include B1a, B1b and upper B1c (Table 14). There is no clear interbred between B1a and B1b in the sector model also (Figs. 9, 10), and the wells grouping has been selected as shown in Fig. 11. The shale barrier has three types including clear shale barrier, thin shale barrier and no clear shale barrier after checking the logging interpretation result (Figs. 12, 13, 14). So case 12 was designed to simulate this situation, in which wells’ perforation has three types (Table 15). Wells with clear shale barrier perforate B1a and B1b. Wells with thin shale barrier perforate B1a and upper B1b. Wells with no clear shale barrier perforate B1a. Some wells need to adjust perforation based on actual production situation such as FNE-5 and FNE-73.

Table 14 Product layer and injection layer optimization cases
Fig. 9
figure 9

Well type distribution

Fig. 10
figure 10

Net gross of (B1a and B1b) layers

Fig. 11
figure 11

Well group map

Fig. 12
figure 12

Sand section around FNE-66

Fig. 13
figure 13

Sand section around FNE-64 and FNE-66

Fig. 14
figure 14

Sand section around FNE-86 and FNE-87

Table 15 Perforation of each well in Case 12

According to the forecast result, the recovery factor of case 12 is the highest, 1% enhanced by the case with current perforation (Fig. 15). Case 12 for the shale barrier is recommended.

Fig. 15
figure 15

Recovery factor curve of each case

Converting time from CSS to SF

There are many reasons to convert CSS to SF such as changes in the temperature, pressure, and saturation of the well group, and from the previous study for this field and from numerical simulation results for turning to steam flooding and after comparing recovery percent of reservoirs, switching to steam flooding after the second period of steam stimulation is also the best (Elbaloula et al. 2016).

Considering the previous studies and experience of CSS implementation in FNE oil field, the assumption for the simulation was to start the SF as soon possible since most wells has more than two CSS cycles and many wells start production with CSS at different times.

Final results

Base on the above optimized result, FNE-03 block SF#2 pilot area can choose 7 spots with 6 new producers and 4 wells converting to injectors in January 2020 (18producers + 4 injectors). In the recommended case, the producers’ production intensity and injectors’ injection index are 20 bbl/day m and 7 m3/day m, respectively (Tables 16, 17). Steam quality and temperature need to maintain current situation of 0.7 and 270 °C, respectively. The SG should not stop to inject steam, and considering shale barrier when adding new perforation is necessary.

Table 16 Liquid rate of producers in recommended case
Table 17 Injection rate of injectors in recommended case

For the recommended case forecast result, by the end of 2035, cumulative oil is 9.69 MMbbl, oil recovery factor of Bentiu is 32.3% enhanced by 15% compared with do nothing case (Tables 18, 19 and Fig. 16). Compared the optimized oil rate and liquid rate in the recommended case with the actual current data (May. 2018), the production performance will be greatly improved. Also Aradaiba D formation can consider for future production and could get additional cumulative oil of 1.47 MMstb.

Table 18 Forecast result of DNC and recommended case
Table 19 Recommended case development indicator prediction
Fig. 16
figure 16

Oil rate, water cut and recovery factor curve of recommended case

Economic evaluation

Although this work is a scenario-based optimization, simple economic evaluation model has been done for this research in order to fully picture the feasibility study of steam flooding in this area, and the assumptions for the sensitivity analysis are as follows:

Evaluation period is 2016–2027, under the oil price of 50$/bbl, we do the IRR sensitivity analysis of Contractor, and the variation scopes are 10% and 20%, and the detailed sensitivity analysis results are given in Table 20 and Fig. 17.

Table 20 The net cash flow for the contractor (50$/bbl)
Fig. 17
figure 17

IRR sensitivity analysis (50$/bbl)

From sensitivity analysis, the oil price is the most sensitive factor. When the oil price decreases 20%, the IRR is 12.85% still greater than 10%, which means the contractor has the strong ability to resist risks, and under the present oil price, the contractor is profit.

Also when the oil price is 40$/bbl, the IRR of the four basic cases is 32.45%, 31.48%, 25.12% and 32.47%. The state NPV of the four basic cases is 153.85 MM$, 208.29 MM$, 206.43 MM$ and 212.46 MM$, and in all of the cases except DNC case, from the sensitivity analysis, the oil price got the highest sensitivity of IRR, in the break-even point, the oil price is 34.78$/bbl of recommended case, all the cases are profit, but we still should pay more attentions to the international oil price to avoid investment risks.

Conclusions and recommendation

  • The feasibility study to design the second steam flooding pilot test patterns in Sudanese oil field has been done, and the pilot area has been selected and designed using advanced thermal EOR simulator.

  • After the analogy with the development effect of similar international oil fields, combined with international standards for thermal recovery, the preliminary result from the screening shows that FNE oil field is suitable for steam flooding.

  • According to the geological and dynamic facts of FNE oil field, the steam flooding pilot test area and well groups have been selected in FNE-03 block SF#2 pilot area and the optimized result to choose 7 spots with 6 new producers and 4 wells to be converted to injectors in 2020.

  • The results show that the optimum steam flooding parameters for this area are steam injection temperature of 270 °C, with 5–7 MPa injection pressure, steam injection quality of more than 0.6 and steam injection rate of 1.6 m3/day/ha/m, and the final recovery ratio will reach 32.3%.

  • The optimum production rate for producers is 20 bbl/day m and for injectors’ injection index is 7 m3/day m).

  • From simple economic analysis, the IRR for the recommended case is 32.47%, NPV is 212.46 MM$, and in the break-even point, the oil price is 34.78$/bbl.

  • It is highly recommended for this area to convert from CSS to steam flooding producing in 2020 in order to produce high recovery factor, and it is necessary to consider shale barrier when adding new perforation.

  • It is strongly recommended to conduct detail economic evaluation for this area for field implementation.