NTNU THESIS Simulation Study of Enhanced Oil Recovery by ASP Flooding For Norne Field - Albadi - Farid
NTNU THESIS Simulation Study of Enhanced Oil Recovery by ASP Flooding For Norne Field - Albadi - Farid
Farid Abadli
Petroleum Engineering Submission date: July 2012 Supervisor: Odd Steve Hustad, IPT Co-supervisor: Jon Kleppe, IPT
Norwegian University of Science and Technology Department of Petroleum Engineering and Applied Geophysics
Simulation Study of Enhanced Oil Recovery by ASP (Alkaline, Surfactant and Polymer) Flooding for Norne Field C-Segment MASTER THESIS Farid Abadli
Supervisor: Odd Steve Hustad
NORWEGIAN UNIVERSITY OF SCIENCE AND TECHNOLOGY DEPARTMENT OF PETROLEUM ENGINEERING AND APPLIED GEOPHYSICS
2012
Abstract
This research is a simulation study to improve total oil production using ASP flooding method based on simulation model of Norne field C-segment. The black oil model was used for simulations. Remaining oil in the reservoir can be divided into two classes, firstly residual oil to the water flood and secondly oil bypassed by the water flood. Residual oil mainly contains capillary trapped oil. Water flooding only is not able to produce capillary trapped oil so that there is a need for additional technique and force to produce as much as residual oil. One way of recovering this capillary trapped oil is by adding chemicals such as surfactant and alkaline to the injected water. Surfactants are considered for enhanced oil recovery by reduction of oil water interfacial tension (IFT). The crucial role of alkali in an alkaline surfactant process is to reduce adsorption of surfactant during displacement through the formation. Also alkali is beneficial for reduction of oil-water IFT by in situ generation of soap, which is an anionic surfactant. Generally alkali is injected with surfactant together. On the other hand, polymer is very effective addition by increasing water viscosity which controls water mobility thus improving the sweep efficiency. In the first place, ASP flooding was simulated and studied for one dimensional, two dimensional and three dimensional synthetic models. All these models were built based on Csegment rock properties and reservoir parameters. Based on test runs, well C-3H was selected and used as a main injector in order to execute chemical injection schemes in the C-segment. Five studies such as polymer flooding, surfactant flooding, surfactant-polymer flooding, alkaline-surfactant and alkaline-surfactantpolymer flooding were considered in the injection process and important results from simulator were analyzed and interpreted. Sensitivity analyses were done especially focusing on chemical solution concentration, injection rate and duration of injection time. The polymer flooding project in this study has shown a better outcome compared to water flooding project. Economically best ASP solution flooding case is the flooding with concentration of alkaline at1.5kg/m3, surfactant at 15kg/m3 and polymer at 0.35 kg/m3 injecting for 5 years. AS flooding case for 4 years with alkali concentration at 0.5kg/m3 and surfactant concentration at 25 kg/m3 gave highest NPV value. It was found that surfactant flooding has a promising effect and it is more profitable than polymer flooding for the C segment in terms of NPV. Economic sensitivity analysis (Spider diagram) for low case, base case and high case at different oil prices, chemicals prices, and discount rate were also presented. It was found that change in oil price has significant effect on NPV compared to other parameters while polymer price has the least effect on NPV for high and low cases.
This work is dedicated to my respected parents and brother. May Allah always bless you!
ii
Acknowledgements
I would like to express my gratitude to my supervisor Odd Steve Hustad for his help, advices and academic guidance during my thesis work. I feel to thank to Professor Jon Kleppe at NTNU for his contribution and a lot of thanks go to Jan-Ivar Jensen for his technical help. I wish to give thanks to Statoil for provision of the Norne Field data for the purpose of research. My sincere thanks are for Jan ge Stensen (SINTEF), Dag Chun Standnes (Statoil), Richard Rwechungura (NTNU) for their suggestions and advices. I am grateful and say my deep thanks to my family for their support and encouragement. Above all, I thank to Almighty ALLAH on successful completion of this research work and thesis.
iii
Nomenclature
Mobility Ratio Resistance Factor Permeability Reduction Factor Residual Resistance Factor Water Mobility Oil Mobility Water Relative Permeability Oil Relative Permeability Oil Viscosity Water Viscosity Water Flow Rate Oil Flow Rate Porosity Permeability Water Saturation Oil Saturation Polymer Solution Concentration Enhanced Oil Recovery Alkaline, Surfactant and Polymer Alkaline and Surfactant Surfactant and Polymer Interfacial Tension Capillary Number Net Present Value Capillary Pressure Residual Oil saturation after water flooding Surfactant Viscosity Alkaline Concentration Transmissibility Discount Factor Pore Velocity
iv
Table of Contents
ABSTRACT ............................................................................................................................................ I DEDICATION ...................................................................................................................................... II ACKNOWLEDGEMENTS ................................................................................................................. III NOMENCLATURE ............................................................................................................................. IV LIST OF FIGURES .......................................................................................................................... VIII LIST OF TABLES ................................................................................................................................. X 1 INTRODUCTION ........................................................................................................................... 1 1.1 ENHANCED OIL RECOVERY METHODS ............................................................................. 2 1.2 ALKALINE FLOODING ............................................................................................................. 4 1.2.1 MECHANISMS OF ALKALINE ............................................................................................. 5 1.2.2 IN-SITU SOAP GENERATION ............................................................................................. 5 1.2.3 AQUEOUS REACTIONS ........................................................................................................ 6 1.2.4 ION EXCHANGE REACTIONS WITH CLAY ...................................................................... 6 1.2.5 DISSOLUTION AND PRECIPITATION REACTIONS ..................................................... 7 1.3 SURFACTANT FLOODING ...................................................................................................... 8 1.3.1 INTERFACIAL TENSION ...................................................................................................... 8 1.3.2 STRUCTURE OF SURFACTANT ......................................................................................... 8 1.3.3 SURFACTANTS TYPES AND SOME MATERIALS .......................................................... 9 1.3.4 CLASSIFICATION OF SURFACTANTS ........................................................................... 10 1.3.5 SURFACTANT FLOODING POTENTIAL IN THE NORTH SEA ................................ 11 v
1.4 POLYMER FLOODING ........................................................................................................... 11 1.4.1 POLYMER TYPES................................................................................................................ 12 1.4.2 POLYMER PARAMETERS................................................................................................. 13 1.4.3 POLYMER FLOODING EFFECTS ..................................................................................... 15 1.4.4 MOBILITY CONTROL .......................................................................................................... 17 2 NORNE FIELD ............................................................................................................................. 19 3 SIMULATION RESULTS........................................................................................................... 22 3.2 TWO DIMENSIONAL MODEL ............................................................................................. 24 3.3 THREE DIMENSIONAL MODEL.......................................................................................... 25 3.4 SIMULATION FOR C-SEGMENT FIELD MODEL ............................................................ 27 3.4.1 INJECTION WELL SELECTION ........................................................................................ 27 3.4.2 PLAN FOR CHEMICAL SIMULATION CASES ............................................................... 31 3.4.3 SURFACTANT STUDY ....................................................................................................... 32 3.4.3.1 CONCENTRATION SENSITIVITY ................................................................................ 32 3.4.3.2 INJECTION TIME PERIOD SENSITIVITY .................................................................. 35 3.4.4 POLYMER STUDY ............................................................................................................... 37 3.4.4.1 CONCENTRATION SENSITIVITY ................................................................................ 37 3.4.4.1 RATE SENSITIVITY ........................................................................................................ 41 3.4.5 ALKALINE-SURFACTANT STUDY.................................................................................. 43 3.4.6 SURFACTANT-POLYMER (SP) STUDY ......................................................................... 45 3.4.6 ALKALINE-SURFACTANT- POLYMER (ASP) STUDY ............................................... 47 vi
4 ECONOMIC EVALUATION ....................................................................................................... 50 4.1 ECONOMIC SENSITIVITY ANALYSES ............................................................................... 53 5 DISCUSSION ................................................................................................................................ 56 6 CONCLUSION .............................................................................................................................. 57 7 RECOMMENDATION ................................................................................................................ 58 REFERENCES .................................................................................................................................... 59 APPENDIX ......................................................................................................................................... 61 A. ALKALINE PROPERTIES ......................................................................................................... 61 B. SURFACTANT PROPERTIES .................................................................................................. 62 C. POLYMER PROPERTIES .......................................................................................................... 63 D. EXCEL TABLES FROM ECONOMIC EVALUATION............................................................ 64 E. ECLIPSE DATA FILE FOR SPECIFIC ASP CASE .................................................................. 73
vii
List of Figures Figures 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 Oil recovery mechanisms (1) Categorization of EOR methods (3) IOR methods implemented by lithology Schematic of alkaline recovery process (7) Structure of molecular surfactant (10) Critical Micelle Concentration (10) Comparison of oil saturation (13) Viscosity versus concentration for different salinities Viscosity versus concentration Residual oil in dead ends flushed by water, glycerin(13) Fingering effect with water flooding Decreased effect of fingering with polymer flooding Location of the Norne Field and Field Segments (16) Stratigraphical sub-division of the Norne reservoir NE-SW cross-section of fluid contacts for the Norne Field (20) The drainage strategy for the Norne Fielf(20) Overview of one dimensional Eclipse model Field Oil efficiency for different chemicals Reservoir pressure for different chemicals Overview of two dimensional Eclipse model Field Oil Efficiency for different chemicals Overview of three dimensional Eclipse model Field Oil Efficiency for different chemicals Waster-Cut for different chemicals C-Segment production and injection wells Recovery factor for different surfactant injection wells Surfactant Adsorption Total for different surfactant injection wells Surfactant Production Total for different surfactant injection wells C-segment total oil production with different polymer injection well Polymer production with different polymer injection wells C-segment 3D view with production and injection wells C-Segment oil production for different surfactant con. C-Segment water production for different surfactant con. Surfactant production for different concentration
viii
Page 1 2 3 6 9 10 11 13 14 17 18 18 19 20 21 21 22 23 23 24 24 25 26 26 27 28 28 29 29 30 30 33 33 34
35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57
Surfactant adsorption in Block (16, 23, and 18) C-segment oil production for different time period C-segment surfactant production for different continuous injection time period Well C-3H bottom-hole pressure for different continuous injection time period C-Segment oil production for different polymer concentration Polymer production total for different concentrations The overview of Polymer flooding to Buckley Leverett method Field reservoir pressure for different polymer concentration Well C-3H bottom-hole pressure for different polymer concentration C-Segment oil production for different polymer injection rates Reservoir pressure for different polymer injection rates Injector bottom-hole pressure for different polymer injection rates C-Segment oil production for different AS cases Alkali production for different AS cases C Segment oil production for different SP cases Field reservoir pressure for different SP cases Well C-3H bottom-hole pressure for different SP cases C-segment oil production for different ASP cases Reservoir pressure for different ASP cases Injector bottom-hole pressure for different ASP cases Net Present Value for Chemical Cases Different variable changes and their impact on NPV Spider Diagram for ASP flooding case with [email protected]/m3,SURF@15kg/m3 &[email protected] for 5yrs
34 35 35 36 38 38 39 40 40 41 42 42 44 44 45 46 46 48 49 49 52 54 55
ix
NPV for PLY flooding at concentration of 0.3kg/m3 for 4 yrs Values for variables at low, base and high case Incremental NPV for Surfactant flooding at concentration of 15kg/m3 for 4 yrs Incremental NPV for Surfactant flooding at concentration of 15kg/m3 for 6 yrs Incremental NPV for Surfactant flooding at concentration of 40kg/m3 for 4 yrs Incremental NPV for Polymer flooding at concentration of 0.3kg/m3 for 4 yrs Incremental NPV for Polymer flooding at concentration of 0.15kg/m3 for 4 yrs Incremental NPV for Polymer flooding at concentration of 11000kg/D for 4 years Incremental NPV for SP flooding with SURF@15kg/m3 & PLY @0.5kg/m3 Incremental NPV for SP flooding with SURF@40kg/m3 & PLY @0.2 kg/m3 SP flooding with SURF@15kg/m3 & [email protected]/m3 followed WI and PLYI Incremental NPV for AS flooding with [email protected]/m3 SURF @ 25 kg/m3 Incremental NPV for AS flooding with [email protected]/m3 & SURF@25kg/m3 Incremental NPV for AS flooding with [email protected]/m3 & SURF@25kg/m3 Incremental NPV for ASP with [email protected]/m3,SURF@15kg/m3 &[email protected] Incremental NPV for ASP with [email protected]/m3,SURF@15kg/m3 Incremental NPV for ASP with [email protected]/m3,SURF@15kg/m3 &[email protected] Incremental NPV for ASP flooding for 4 years slug AlK1.5_SURF15_PLY0.35 Incremental NPV for AS for 2 yrs following by 1yr WI and 2yrs PLY Incremental NPV for AS for 4 yrs following by 1yr WI and 2yrs PLY
x
1 Introduction
Oil production period mechanisms can be classified as primary, secondary and tertiary mechanisms. By the development of production time reservoir pressure is dropping, so different methods are used to control pressure and increase production. Most large oil fields are produced with some type of secondary pressure maintenance scheme, such as water flooding, gas flooding etc. Oil recovery mechanisms and their classifications are shown in Figure 1.
Right alkali is chosen based on some factors such as price and availability at the flooding area, the PH level, the temperature and mineralogy of the reservoir and composition of the mixed water. (6)
HAw - the concentration of acid in water, KD - the partition coefficient. By the time the acid in water will separate in the aqueous phase to generate soluble anionic surfactant (A-) referred to as soap according to the expression below:
The reaction above is one of the sources of alkaline consumption since the alkali uses hydrogen to generate soap by the following process:
As an example kaolinite, Al2Si2O5(OH)4, is found in most sandstone formations. The dissolution of kaolinite at high pH can result in generation of aqueous types such as
The most important preconditions for polymer flooding are reservoir temperature and the chemical properties of reservoir water. At high temperature or with high salinity in reservoir water, the polymer cannot be kept stabile, and polymer concentration will lose most of its viscosity. (12)
Figure 7: Comparison of oil saturation after polymer flooding and water flooding (13)
11
12
Figure 8: Viscosity versus concentration for different salinities Two factors should be considered when choosing polymer molecular weight. First, the polymer with highest Mw is practical to minimize the polymer volume. Second Mw must be small enough so that polymer can enter and propagate effectively through reservoir rock. (12)
13
Polymer Molecular Weight The effectiveness of a polymer flood is affected significantly by the polymer molecular weight (Mw). As illustrated in Figure 9, polymers with higher Mw provide greater viscosity. For many circumstances, larger polymer Mw also leads to improved oil recovery. Core flood simulation verifies this expectation for cases of constant polymer slug volume and concentration. (12)
Polymer Solution Concentration The polymer solution concentration dominates every index that changes during the course of polymer flooding. 1) Higher injection concentrations cause greater reductions in water cut and can shorten the time required for polymer flooding. For a certain range, they can also lead to an earlier response time in the production wells, a faster decrease in water cut, a greater decrease in water cut, less required pore volumes of polymer, and less required volume of water injected during the overall period of polymer flooding. As polymer concentration increases, enhanced oil recovery increases and the minimum in water cut during polymer flooding decreases. 2) Above a certain value, the injected polymer concentration has little effect on the efficiency of polymer flooding. For a pilot project, the economics of injecting higher polymer concentrations should be considered. The polymer solution concentration has a large effect on the change in water cut. However, consideration should also be given to the fact that higher concentrations will cause higher injection pressures and lower injectivity. For individual wells, the concentration can be adjusted to meet particular conditions. (12) 14
The amount of polymer adsorbed depends on the nature of the polymer and the rock surface. Both physical adsorption of polymer on solid surfaces and polymer retention by mechanical entrapment appear to play a big role in total polymer retention in a reservoir. Generally, three phenomena have been observed regarding polymer adsorption: (1) laboratory tests often indicate higher adsorption than field performance; (2) adsorption is significantly less in consolidated cores than in sand packs; and (3) adsorption increases with increasing water salinity. Typical field adsorption values range from 20 to more than 500Ib/acre-ft. Laboratory adsorption values range from 30 to several hundred jug m/gm.13 note that laboratory results often cannot be extrapolated to predict polymer adsorption in oil reservoirs, polymer retention is also important. (12, 15) Polymer Retention Retention of polymer in a reservoir can result from adsorption, entrapment or with improper application, physical plugging. Polymer retention tests are usually performed a polymer flood oil recovery test. If polymer retention tests are conducted with only water initially present in the core, a higher level of retention will result from the increased surface area available to the polymer solution in the absence of oil. Effluent samples from the core are collected both the polymer injection and a subsequent water flush. These samples are analyzed for polymer content. From a material balance, the amount of polymer retained in the core is calculated. Excessive retention will increase the amount of polymer that must be added to achieve the desired mobility control. The level of polymer retained in a reservoir depends on a number of variables: permeability of the rock, surface area, nature of the reservoir rock (sandstone, carbonate, minerals, or clays), nature of the solvent for the polymer (salinity and hardness), molecular weight of the polymer, ionic charge on the polymer, and the volume of porosity that is not accessible to the flow of polymer solution. (12, 15) Inaccessible Pore Volume Polymer solutions propagate through porous media at a velocity different from that of water because of adsorption and inaccessible pore volume. Adsorption tends to move the front edge of a polymer slug at a slower velocity than the water bank, and inaccessible pore volume tends to move the polymer slug at a higher velocity than the water bank. The combination of the two effects results in a smaller slug that is shifted forward. The phenomenon of inaccessible pore volume was first reported by Dawson and Lantz. They showed that all the pore space may not be accessible to polymer molecules and that this allow be accessible to polymer molecules and that this allows polymer solutions to advance and displace oil at a faster rate than predicted on the basis of total porosity. They also concluded from laboratory 15
results that about 30% of the total pore volume in the rock samples used was not contacted by the polymer solution. The inaccessible pore volume can have beneficial effects on field performance. The rock surface in contact with the polymer solution will be less than the total pore a volume, thus decreasing the amount of polymer adsorbed. More importantly, if connate water is present in the smaller pores inaccessible to the polymer, the bank of connate water and polymerdepleted injection water that precedes the polymer bank is reduced by the amount of inaccessible pore volume. However, movable oil located in the smaller pores will not be contacted by the polymer in some cases, and therefore it may not be displaced. (12) Permeability reduction Polymer reduces both the effective permeability of porous media and the mobility of displacing fluid the permeability reduction described by a reduction factor Rk:
Rk =
Kw Kp
Where Kw and Kp are the effective permeability of water and polymer. The mobility, change due to combined effect of increased viscosity and reduced permeability is the Resistance Factor Rr:
Rr =
R k * p w
R rr =
p p
The effect of permeability reduction persists even after the polymer solution has gone through the porous media. This effect is described by the residual resistance factor Rrr. (12) Where p and `p are the motilities before and after polymer solution, respectively. Dependence of Polymer behavior on shear rate Both types of polymer solutions are known to exhibition non-Newtonian, shear-thinning fluid behavior. Shear thinning properties usually can be determined in the laboratory using Brookfield viscometer.
ShearRate = C / K *
However, determination of the viscometric behavior of polymer solutions in reservoirs is more complex because shear rates are not well defined in the rock matrix.
16
At a given flow velocity or shear rate, the higher the polymer molecular weight is , the greater the mobility reduction will be, The degree of permeability reduction or residual resistance also increases with increasing polymer molecular weight. (12) Displacing residual Oil in Dead Ends Effect of polymer were studied in a laboratory using a glass etched core model with pore diameter of 250Urn, the oil in the core was first flushed by water until the water cut was 100%, then glycerin with viscosity of 30cP until 100% water cut, and then finally by polyacrylamide fluid (viscosity 30 cp) until 100 % water cut. The results showed that viscosity alone cannot mobilize the residual oil as shown by the glycerin flooding. However, the polymer fluid mobilized 4 times the amount of residual oil out of the dead ends than the glycerin. The polymer dragged the residual oil from the dead ends because of its elastic properties, where the fluids in front can pull the fluid behind and beside it. The elastic properties are lacking in the water and glycerin.
Figure 10: Residual oil in dead ends flushed by water, glycerin, and Polyacrylamide (13)
17
18
2 Norne Field
The Norne Field is situated approximately 200 km from the Norwegian coastline. It is located on the blocks 6608/10 and 6508/1 in the southern part of the Nordland II area. The water depth in the field zone is nearly 380 m. The main operator of the field is Statoil ASA and Eni Norge AS and Petoro are license partners.
Figure 13: Location of the Norne Field and field segments (16)
19
20
Figure 15: NE-SW cross-section of fluid contacts for the Norne Field (20)
Figure 16: The drainage strategy for the Norne Field from pre-start to 2014 (20)
21
3 Simulation Results
ASP simulation was done for both synthetic models and for the field model. Three different synthetic models were built based on the field data. Before doing ASP study analyses for field model, simulation was applied to one dimensional, two dimensional and three dimensional models.
0.4
0.3
FOE
0.2
0.1
06/11/1997
14/02/1998
25/05/1998
02/09/1998
11/12/1998
21/03/19
300
RESERVOIR PRESSURE
ALK_SUR ALK_SUR_POLY BASE CASE POLYMER SURFACTANT
250
Pressure (BARSA)
200
150
100
50 06/11/1997
14/02/1998
25/05/1998
02/09/1998
11/12/1998
21/03/19
Time, years
23
FOE
13/09/1998
03/10/1998
23/10/1998
12/11/1998
02/12/1998
Time, Years
24
25
FOE
1998
1998
1998
1999
Water cut
1997
1998
1998 TIME,year
1998
1998
1999
26
Figure 25: C-Segment production and injection wells As it can be seen from Figure 25 we have 4 injectors in C-segment, in the first place, surfactant properties and relevant keywords were included to eclipse data file and continuous surfactant was injected for 4 years through each injection well to compare wells effect based on filed simulation models production performance. The wells were added in 2005 injecting surfactant with concentration of 15kg/m3 until the end of 2008, and then only water was injected for rest of simulation life. Figure 26 indicates that surfactant mixture with water has positive effect with all injectors on oil recovery compared to base case performance. Obviously, injection through C-4H has smallest effect among other wells but production is higher and almost the same when C-2H and C-3H were used as a chemical injection well. 27
0.5
BASE CASE Surfactant injection Surfactant injection Surfactant injection Surfactant injection in in in in Well Well Well Well C1H C2H C4H C3H
FOE
0.5
0.5
0.5
0.5 2021
2021
2021
2021
2021
2022
TIM E, Ye ar
Figure 26: Recovery factor for different surfactant injection wells When we look at surfactant adsorption and production plots, injection with C-3H leads to smallest adsorption in the formation and highest chemical production as a result of higher oil production. This can be because of area around Well C-3H is high permeable zone and certain amount of chemical reaches producers in a short time with low adsorption so we lose less chemical in the reservoir but recovery is high. Based on surfactant implementing on wells C-3H can be use as a main chemical injector.
200M M
150M M
FTADSUR (KG)
100M M
50M M
0 2005
2008
2011
2014 TIM E, Ye ar
2016
2019
2022
Figure 27: Surfactant Adsorption Total for different surfactant injection wells
28
FTPTSUR (KG)
100M M 80M M 60M M 40M M 20M M 0 2005 2008 2011 2014 TIM E, Ye ar 2016 2019 2022
Figure 28: Surfactant Production Total for different surfactant injection wells For optimal injection well selection polymer was used as well. Flooding was implemented starting in 2013 lasting 4 years with polymer solution concentration at 3kg/m3 through all 4 injection wells (C-1H, C-2H, C-3H, C-4H). From Figure 29 and Figure 30 we can see that polymer injection in well C-3H gives higher oil production compared polymer injection in other injection wells and in this case production of polymer is relatively higher.
42808000
41808000
40808000
PLY injection with Well PLY injection with Well PLY injection with Well PLY injection with Well Group Base Case
39808000
38808000
37808000 2008
2011
2014
2016
2019
2022
Time,Year
Figure 29: C-segment total oil production with different polymer injection wells
29
200000
150000
100000
50000
0 1997
2000
2003
2006
2008
2011
2014
2017
2019
2022
TIME,Year
Figure 30: Polymer production with different polymer injection wells Based on the proper well selection process in terms of surfactant and polymer flooding, well C-3H gave better results among all injectors and it was selected an optimal well in order to make all predictions, chemical sensitivity cases and ASP combination scenarios.
30
Surfactant Study
Case 1: Case 2: Case 3: Case 4: Case 5: Case 6: Surfactant flooding at concentration of 15kg/m3 Surfactant flooding at concentration of 40kg/m3 Continuous Surfactant flooding for 2 years Continuous Surfactant flooding for 3 years Continuous Surfactant flooding for 4 years Continuous Surfactant flooding for 6 years
Polymer Study
Case 1: Case 2: Case 3: Case 4: Case 5: Case 6: Case 7: Polymer flooding at concentration of 0.15kg/m3 Polymer flooding at concentration of 0.3kg/m3 Polymer flooding at concentration of 0.5kg/m3 Polymer flooding at concentration of 0.8kg/m3 Polymer flooding at injection rate of 5000kg/D Polymer flooding at injection rate of 8000kg/D Polymer flooding at injection rate of 11000kg/D
31
32
44733000
44233000
43733000
43233000
42733000
2017
2018
2019
TIME, Year
2020
2021
GROUP WATER PRODUCTIN TOTAL 6.0x107 Group Base Case SURF Concentration 15 Kg/m3 SURF Concentration 25 Kg/m3 SURF Concentration 40 Kg/m3
5.0x107
Water production total (SM3)
4.0x107
3.0x107
2.0x107
1.0x107
1997
2000
2003
2006
2008
2011
2014
2017
2019
2022
Time, Year
33
SURFACTANT PRODUCTIONTOTAL
60000000.0
50000000.0
40000000.0
GTPTSUR (KG)
30000000.0
20000000.0
10000000.0
1997
2000
2003
2006
2008
2011
2014
2017
2019
2022
TIME,Year
BLOCK (16, 23, 18) SURF. ADSORBTION TOTAL 0.00020 SURF Concentration 15 Kg/m3 SURF Concentration 25 Kg/m3 SURF Concentration 40 Kg/m3 0.00015
BTADSUR (KG/KG)
0.00010
0.00005
1997
2000
2003
2006
2008
2011
2014
2017
2019
2022
TIME,Year
Figure 35: Surfactant adsorption in Block (16, 23, and 18) for different concentration
34
44868000
44818000
44768000
44718000
44668000
44618000 2020
2020
2021
2021
TIME, Years
2021
2021
Figure 36: C-segment oil production for different continuous injection time period
20000000.0
2 3 4 6
GTPTSUR (KG)
15000000.0
10000000.0
5000000.0
1997
2000
2003
2006
2008
2011
2014
2017
2019
2022
TIME,Years
Figure 37: C-segment surfactant production for different continuous injection time period 35
Figure 38: Well C-3H bottom-hole pressure for different continuous injection time period
From Figure 36 and Figure 37, it can be seen that longer injection period in early life of simulation leads to higher oil production and obviously higher chemical mixture in the flooding water. At different stages, injection for the period of 2, 3 and 4 years gave very close recovery values. Based on the analyses, to maximize production and at the same time to optimize the economic recovery long injection period is not profitable where injection for 2 or 3 years can be optimal choice. Even it can be wasteful injecting surfactant for 6 years when taking economic evaluation into consideration which will be discussed in details later. That is why better approach would be short term injection time period. The injection well bottom-hole pressure is important parameter for EOR projects. There is a need for detailed analyzing of pressure around the injection well to keep the formation stability in the reservoir and avoid formation damage. Pressure graph demonstrates that with the same injection rate and the same solution chemical concentration injector bottom-hole pressure is going up a certain constant pressure limit for all injection time period sensitivity cases. To sum up, surfactant flooding is recommended for the Norne C-Segment especially when the concentration is low and injection occurs in the early years. Injection of surfactant at a later time might not be profitable.
36
to increase which causes improvement in mobility ratio that leads to flow spread in the reservoir by diverting injected water from zones that have been swept and widening the flow pattern zone by decreasing effect of fingering. As a result, more adsorption takes place because of wide sweeping area. Taking these concepts, production of polymer is lower for higher concentration in our case.
42839000
42339000
2016
2017
2020
2021
300000
200000
100000
1997
2000
2003
2006
2008
2011
2014
2017
2019
2022
TIME,years
It is believed that polymer flooding cannot reduce the Sor, but it is still an efficient way to reach the Sor more quickly or/and more economically which is found in Buckley Leverett solution as described in figure BL.
The graphs in Figure 42 and Figure 43 shows how pressure near the injection well change due to polymer flooding at different concentrations. In the period of 4 years continuous injection the incremental pressure is observed around injector but reservoir pressure starts to drop. Study of concentration effect shows the rate of reservoir pressure is opposite dependent with polymer concentration where concentration at 0.8kg/m3 leads to lower reservoir pressure. However, consideration should also be given to the fact that higher concentrations will cause higher injection pressures and lower injectivity. For the concentration at 0.8kg/m3 recovery is highest but highest bottom-hole pressure where incremental pressure is more than 50 bars makes it unfavorable case. Based on the analyses, the case at 0.8 kg/m3 is unfavorable though the case at 0.5kg/m3 can be optimal proposal.
39
RESERVOIR PRESSURE
330 320 310 300 290 280 270 260 250 2011 2012 2012 2013 2014 2015 2015 2016 2017 2017 2018 2019 2019 2020 2021 2021 TIME,years PLY PLY PLY PLY concentration 0.3 kg/m3 concentration 0.5 kg/m3 concentration 0.8 kg/m3 concentration 0.15 kg/m3
Figure 43: Well C-3H bottom-hole pressure for different polymer concentration
Pressure (BARSA)
40
42945000
41945000
40945000 2015
2016
2017
2017
2018
2019
2019
2020
2021
2021
TIM E,years
Figure 44: C-Segment oil production for different polymer injection rates
41
RESERVOIR PRESSURE
450
400
Pressure (BARSA)
350
300
250
200 1997 1999 2000 2001 2003 2004 2006 2007 2008 2010 2011 2012 2014 2015 2017 2018 2019 2021 2022 TIME,years
400
300
PLY Rate 11000 sm3/D PLY Rate 8000 sm3/D PLY Rate 5000 sm3/D
200
100
0 1997
2000
2003
2006
2008
2011
2014
2017
2019
2022
TIME,years
Figure 46: Injector bottom-hole pressure for different polymer injection rates
42
There are some benefits of alkali which is applicable most of the flooding process. The main concept of alkali surfactant flooding is the reduction of oil water Interfacial tension (IFT) and alkali plays a big role in an alkaline surfactant process to reduce adsorption of the surfactant during displacement. Important point is that AS reduces the amount of surfactant required and causes higher production in an economic way. As it can be seen from Figure 47 simulator has produced lower oil production for the AS flooding with alkali at higher concentration. Obviously, production development is almost the same with alkali concentration at 1.5kg/m3 and 2.5kg/m3. Simulation with lower concentration at 0.5kg/m3 gives relatively better result. The reason can be that if alkali concentration is too much it can decrease viscoelasticity of AS solution which can reduce sweep efficiency. And we need to consider limit for alkali concentration due to the concept of minimum achievable value of IFT which is optimal value for most positive effect on movable oil. In our case, according to simulation result alkali is necessary for AS solution but high alkali concentration is not so effective. Figure 48 tells us that more alkali concentration more chemical production. As a result, AS solution is an effective flooding process making more contribution with attendance of alkali but key observation from simulation results is that high alkali concentration low oil production. So Eclipse results for particularly alkali effect are matched with lab results discussed in literature part.
43
2021
2021
2021
2021 TIME,years
2021
2021
2022
GTPTALK (KG)
2000000
1000000
1997
2000
2003
2006
2008
2011
2014
2017
2019
2022
TIME, years
44
44860000
44660000
44460000
44260000
44060000 2017
2018
2019
2019
2020
2021
2021
TIME,years
RESERVOIR PRESSURE
450
400
SURF ( 15 kg/m3 ) and PLY (0.5 kg/m3) SURF & PLY followed by water and PLY SURF ( 40kg/m3) and PLY (0.2 kg/m3)
Pressure (BARSA)
350
300
250
200 1997
2000
2003
2006
2008
2011
2014
2017
2019
2022
TIME,years
500
Bottom hole pressure (BARSA)
400
300
SURF ( 15 kg/m3 ) and PLY (0.5 kg/m3) SURF & PLY followed by water and PLY SURF ( 40kg/m3) and PLY (0.2 kg/m3)
200
100
0 1997
2000
2003
2006
2008
2011
2014
2017
2019
2022
TIME,years
46
47
From Figure 52 we can see that longer duration of ASP flooding gives better recovery. During this time residual oil saturation is decreased while long and continuities flooding leads to much trapped oil to move. As illustrated on the graph, periodical slug injection efficiency is very close to the efficiency of longer continuous injection where slug injection demands less chemical compared to continuous flooding. So the advantage of slug injection is to reduce chemicals costs by means of lower amounts of chemicals. Obviously, based on the last two cases injecting AS and polymer separately gives lower recovery. Considering this point, attendance of all chemicals in the solution is more beneficial and would be better choice. We get higher incremental reservoir pressure and injector well pressure when followed polymer slug injection is used. But polymer in the solution with other chemicals does not have so big effect on pressure increase. Due to the constant concentration of each chemical, maximum injection pressure level while chemical flooding is stable for all cases since rate does not change as well. From what have been discussed above, ASP as an effective chemical flooding process and it can be the potential EOR method for the C-segment.
44805000
Oil production total (SM3)
44605000
44405000
44205000
44005000
43805000
2017
2017
2018
2019
2019
TIME,years
2020
2021
2021
48
RESERVOIR PRESSURE
312 302 292 282 272 262 252 242 2012
Pressure (BARSA)
ASP flooding for 2 yrs ASP flooding for 4 yrs ASP flooding for 5 yrs Periodical ASP slug flooding for 5 yrs 2 yrs AS inj. followed by WI & PLY inj. 4 yrs AS inj. followed by WI & PLY inj.
2013
2014
2014
2015
2016
2016
2017
2018
2019
2019
2020
2021
2021
2022
TIME, years
500
400
300
200
ASP flooding for 2 yrs ASP flooding for 4 yrs ASP flooding for 5 yrs Periodical ASP s lug flooding for 5 yrs 2 yrs AS inj. followed by WI & PLY inj. 4 yrs AS inj. followed by WI & PLY inj.
100
0 1997
2000
2003
2006
2008
2011
2014
2017
2019
2022
Time, years
49
4 Economic Evaluation
Profitability analysis of the project is very important evaluation after the technical analysis is completed. This analysis is the basis for investment decision of the project. For this EOR project economic analysis and proper decisions have been done based on Net Present Value (NPV). NPV takes into consideration the value of cash earned in the future and converts the money to the present value. The equation below is used to calculate NPV.
NPV =
Re venue Expense (1 + r )t
Where r represents discount rate and t represents projects year. In this simulation study implementation of ASP method started from 2013 and simulation lasted until 2022. The year 2013 is considered as a start year of the project therefore it is zero as a project year. So for economic evaluation years started from zero till nine which is 2021. DF is discount factor (for determining time value of money) that is taken to account in NPV formula. A project is categorized as profitable when NPV is positive. And we are always encouraged to choose higher NPV when we have many alternatives with positive NPV. A high NPV indicates a high profitability from the project. In our case, only revenue is got from selling of crude oil. It is assumed that expenses only include surfactant costs, polymer costs and alkali costs. For a specific case to calculate NPV and make decision, incremental revenue between chemical case and base case (no chemical, only water injection) has been used as a basis cash flow. For the estimation of income, annual total production data were taken from the Eclipse and included to excel. The following prices for oil and chemicals were used. And conversion of units was done to get NPV with NOK. The crude oil price for current year is referred to http://www.oil-price.net/
50
Table 1: NPV for Polymer flooding at concentration 0.3kg/m3 for 4 yrs Table 1 represents the procedure how to calculate NPV for one of the polymer cases. Likewise, the procedure as described in the table was done for all chemical cases. Inflation affected only oil price while assuming chemical prices are constant for predicted years. We get present value by discounting net cash flow. Incremental oil production has been found by subtracting annual base case production from production with polymer flooding. NPV for each year was calculated cumulating of annual present value. We can see that cumulative incremental NPV is 1.77E+12 NOK in the end of simulation life which is key element for decision making of investment. Obviously, it is positive in this case and can make project being accepted. Economic evaluation process demonstrated in the Table 1 was done for all typical chemical flooding cases and NPV values versus years were plotted in figure 55. By considering the following plots, it is evident that the highest NPV is from AS flooding for 4 years with alkali and surfactant at concentration of 0.5kg/m3 and 25 kg/m3 respectively. We have almost the same expected net present value at 1.09E+13 NOK for SP case with concentration of 15kg/m3 and 0.5kg/m3 respectively. The incremental NPV for polymer cases are lower than all other cases. Interesting observation is that in terms of NPV both AS and SP cases are better than ASP solution cases. Even adding only surfactant to water makes higher expected incremental NPV. Similarly, attendance of surfactant is visible in all high profitable cases. As long as the expected net present value is positive, the project should be accepted but considering many alternatives, the cases with surfactant and its solution leads to much profit and based on the profitability evaluation and from investment point of view injecting only surfactant should be better than injecting combination of chemicals for C segment. Surfactant is very favorable chemical for the Norne C segment while polymer flooding is not encouraged to implement because that single polymer cases cause lower income. Purposefully, more analysis and studies are required to investigate surfactant effect for this field segment.
51
52
The graphs in Figure 57 were plotted based on change values (%) for each variable shown in Figure 56. From the results of the sensitivity analysis we can see that (+11.76) % and (-11.76) % change in oil price causes (+11.77) % and (+11.77) % change respectively in NPV. This means that change of oil price has a big effect on NPV value. Compared to chemical price effect, percentage change of discount rate has great effect on NPV where change of (-37.5) % and (+37.5) % in discount rate corresponds nearly (-9.44) % and (+11.31) % change in NPV respectively. Chemical prices have much lower effect on NPV compared to oil price. Surfactant has highest impact of (+0.2/-0.2) % on NPV when it is deviated by (+39/-39) % from base case surfactant price. From analyses and observations, it can be concluded that oil price has highest impact on NPV in terms of percentage change. On the contrary NPV is least sensitive to the polymer price.
53
OIL PRICE Low Case 2810.601853 -11.76% i 6.37E+12 -11.77% ALKALINE PRICE Low Case 5.958 -33.33% i 7.22E+12 0.000779 %
High Case 11.916 33.33% L 7.22E+12 -0.000779 % High Case 27.4068 39.39% L 7.22E+12 -0.020264 %
SURFACTANT PRICE Low Case Base Case Surfactant price(NOK/kg) 11.916 19.6614 Change (%) -39.39% 0.00% i k NPV 7.22E+12 7.22E+12 Change (%) 0.020264 % 0.000000 % POLYMER PRICE Low Case 14.895 -37.50% i 7.22E+12 0.000546 % DISCOUNT RATE Low Case 0.11 37.50% i 6.54E+12 -9.443865 %
54
Figure 57: Spider diagram for ASP flooding case with [email protected]/m3,SURF@15kg/m3 &[email protected] injecting for 5yrs
55
5 Discussion
In this project main purpose was to investigate ASP flooding as a tertiary recovering method in Norne field C segment through simulation study. The work was carried out using Eclipse simulation. Norne C segment has 13 active wells including 9 producers and 4 injectors. We needed to choose one injection well to implement injection scheme and sensitivity analyses. Purposefully, the same chemical injection strategy was applied through each of injector and the impact on field production performance, chemical adsorption and production was studied. Finally, the well C-3H was selected as an optimal injector based on its location and numerical simulation results. This well was used as a chemical injector in the further ASP cases. Firstly, only surfactant was added to water with different concentration and injection time period cases. It was discovered that, flooding is profitable when the concentration is low and injection occurs in the early years. Injection of surfactant at a later time might not be profitable. On the other hand, long injection period might not be profitable. So injection for 2 or 3 years can be optimal choice. From economic evaluation results, surfactant is very favorable chemical for the Norne C segment that all cases with surfactant gave high net present value. The next step was polymer flooding study. Analyses showed that polymer injection as EOR method is more efficient than water injection because of improved mobility ratio, increasing the cumulative oil production and decreasing the total water production. Oil recovery factor increases with an increase in polymer solution concentration because of increase in displacing fluids viscosity, but decreases with increase of rate so lower injection rate causes much better volumetric sweep efficiency. The incremental NPV for polymer cases are lower than all other cases which make polymer flooding as unfavorable method for C Segment. Moreover, analyses were done for combination of 2 and 3 chemicals. From interpretations, AS solution was an effective flooding process with attendance of alkali. Main observation was that higher alkali concentration, lower oil production. SP solution caused higher oil production compared to single chemical flooding with positive sides that surfactant decreases IFT at the same time polymer improves mobility by increasing fluid viscosity. It was evident that the highest NPV value belongs to AS flooding case. SP case with concentration for surfactant at 15kg/m3 and polymer at 0.5kg/m3 was second best economical case. Simulation study showed that ASP as an effective chemical flooding process can be the potential EOR method for Norne C-segment by considering right chemical combinations with proper concentrations and duration of injection time. Spider diagram clearly proved that oil price and discount rate has highest impact on NPV in terms of percentage change in price.
56
6 Conclusion
The injector C-3H is optimal well for ASP flooding in Norne field C-segment. The ASP flooding is applicable EOR method for C segment. Economically, surfactant is very favorable chemical for C-segment. It is more profitable to inject surfactant for 2 and 3 years and it is wasteful injecting for 6 years.
Polymer did not perform efficiently for both synthetic and field models. Change of chemical prices has lower effect on NPV compared to oil price. NPV was least sensitive to polymer price while it was highest sensitive to surfactant price.
Production of polymer is lower for higher concentration in our case. Simulation with lower alkali concentration gives relatively better result. AS flooding for 4 years with alkali and surfactant at concentration of 0.5kg/m3 and 25 kg/m3 respectively gives highest NPV value.
Polymer flooding at rate of 11000kg/D for 4 years gives lowest NPV which is most unfavorable case.
Economically best ASP solution flooding case is the flooding with concentration of alkaline at 1.5kg/m3, surfactant at 15kg/m3 and polymer at 0.35 kg/m3 injecting for 5 years.
57
7 Recommendation
For more accurate and real results it is recommended detailed laboratory work focused on Csegment formation properties should be done considering ASP flooding. There is a need for up-scaling strategy to be development for our field model. Moreover, selection of proper time and duration of injection for chemical is important. Eclipse assumes that shear rate is proportional to the flow viscosity and does not show shear effect in the reservoir. Taking the possibility of share effect in the reservoir wide studies needs to be carried out.
58
References
1. Schmidt, R.L. Thermal Enhanced Oil RecoveryCurrent Status and Future Needs. January 1990. 2. Modern Chemical Enhanced Oil Recovery-Theory and Practice (www.knovel.com) 3. S.M Farouq Ali & S. Thomas, University of Alberta, The Promise and Problems of Enhanced Oil Recovery Methods JCPT September 1996, volume 35, No. 7. 4. Dr. Leonid Surguchev,Vice President,Rogaland Research, Dr. Eduardo Manrique, Questa Engineering Corp., Golden, CO, USA, Prof. Vladimir Alvarado, Catholic Pontifical University of Rio de Janeiro, Rio de Janeiro, Brazil. Improved Oil Recovery: Status and Opportunities. 5. Hourshad Mohammadi, SPE, Mojdeh Delshad and Gary A. Pope. Mechanistic Modeling of Alkaline/Surfactant/Polymer Floods, the University of Texas at Austin, SPE 110212 6. R.C. Nelson, J.B. Lawson, D.R. Thigpen and G.L. Stegemeier, Shell Development Co. Co-surfactant Enhanced Alkaline Flooding, SPE 12672 7. Sheng, James J. Modern Chemical Elsevier Inc, USA. 2011. Enhanced Oil Recovery (Theory and Practice),
8. Y.P. Zhang, S.G. Sayegh, S. Huang. The Role of Effective Interfacial Tension in Alkaline/Surfactant/Polymer Flooding, Petroleum Society, Canadian Institute of Mining, Metallurgy and Petroleum, Paper 2005-022 9. Chinenye Clara Emegwalu, Enhanced Oil Recovery For Norne Fields E-Segment Using Surfactant Flooding, Master Thesis, NTNU, Spring 2010 10. A.K. Flaaten, Q.P. Nguyen, J.Zhang, H. Mohammadi, G.A. Pope. ASP Chemical Flooding Without the Need for Soft Water, SPE, The University of Texas at Austin 11. SVEIN M. SKJVELAND & JON KLEPPE. Recent Advances in Improved Oil Recovery Methods for North Sea Sandstone Reservoirs SPOR Monograph 12. Farid Abadli. Investigation of Polymer Flooding for Improved oil Recovery, Semester Project, NTNU, December 2011. 13. Wang Demin, Li Qun, SPE, Gong Xiaohong, Wang Yan / Daqing Oil Company, The Engineering and Technical Aspects of Polymer Flooding in Daqing Oil Field, SPE 64722 14. Y. Du, L. Guan. Field-Scale Polymer Flooding: Lessons Learnt and Experiences Gained During Past 40 Years, SPE, New Mexico Institute of Mining and Technology, Texas A&M University, SPE 91787 15. Eldar Sadikhzadeh. Evaluation of polymer flooding for improved recovery in the Statfjord Formation of the Gullfaks Field, Semester Project, NTNU, December 2006 59
16. Statoil (2004) Annual Reservoir Development Plan, Norne Field 17. Statoil (1994) PDO-Reservoir Geology, Support Documentation 18. Statoil (2006) Annual Reservoir Development Plan, Norne and Urd Field. 19. N. Arihara, Lu XiangGuo, T. Yoneyama, Y. Akita. Oil Recovery Mechanisms of Alkali-Surfactant-Polymer Flooding, SPE, Waseda University, JNOC and Daqing Petroleum Institute, SPE 54330 20. Study of Enhanced-Oil-Recovery Mechanisms of Alkali/Surfactant/Polymer flooding in porous media from experiments 21. Yugal Kishore Maheshwari. A Comparative Simulation Study of Chemical EOR Methodologies (Alkaline, Surfactant and/or Polymer) Applied to Norne Field E-Segment, Master Thesis, NTNU, August 2011 22. Schlumberger. Eclipse Technical Description Manual. 2009.2 23. Schlumberger. Eclipse Reference Manual. 2009.2
60
APPENDIX
A. Alkaline Properties
ALSURFST --Alkaline --conc --Kg/m3 0.0 6.0 15.0 20.0 30.0 ALPOLADS --Alkaline --conc --Kg/m3 0.0 3.0 6.0 9.0 0.0 3.0 6.0 9.0 ALSURFAD --Alkaline --conc --Kg/m3 0.0 3.0 6.0 9.0 0.0 3.0 6.0 9.0 ALKADS --Alkaline --conc --Kg/m3 0.0 3.0 6.0 9.0 10.0 0.0 3.0 6.0 9.0 10.0 Water/oil Surface Tension Multiplier 1.0 0.5 0.3 0.1 0.0 /
Polymer adsorption Multiplier 1.0 0.7 0.5 0.3 / 1.0 0.7 0.5 0.3 /
Adsorption Multiplier 1.0 0.7 0.5 0.0 / 1.0 0.7 0.5 0.0 /
Alkaline Adsorbed on rock (kg/kg) 0.000000 0.000005 0.000007 0.000008 0.000009 / 0.000000 0.000005 0.000007 0.000008 0.000009 /
61
B. Surfactant Properties
SURFVISC --Surf conc --Kg/m3 0.0 0.318 5.0 0.449 10.0 0.503 15.0 0.540 20.0 0.630 / Water Viscosity Centipoises ---0.42
SURFST --Surf conc Water/oil Turface Tension --Kg/m3 (N/m) 0 30.0E-03 0.1 10.0E-03 0.25 1.60E-03 0.5 0.40E-03 1.0 0.07E-03 2.0 0.01E-03 3.0 0.006E-03 5.0 0.004E-03 10.0 0.006E-03 15.0 0.008E-03 20.0 0.01E-03 / SURFADS --Surf conc Adsorbed mass --Kg/m3 (kg/kg) = kg surf /kg rock 0.0 0.00000 1.0 0.00017 5.0 0.00017 10.0 0.00017 / 0.0 0.00000 1.0 0.00017 5.0 0.00017 10.0 0.00017 / SURFCAPD --Log10 (capillary --number) -8 -7 -6 -5.0 -2.5 0 5 10 SURFROCK --Adsorption --index -2 2
Miscibility function 0 = immisible, 1= miscible 0.0 0.0 0.0 0.0 1.0 1.0 1.0 1.0 /
62
C. Polymer Properties
--Polymer shear thinning data -- Wat. Velocity Visc reduction -- M/day CP --0.0 1.0 --2.0 1.0 / -- Polymer solution Viscosity Function PLYVISC -- Ply conc. Wat. Visc. mult. -- Kg/m3 0.0 1.0 0.1 1.55 0.3 2.55 0.5 5.125 0.7 8.125 1.0 21.2 / -- Polymer Adsorption Function PLYADS -- Ply conc. Ply conc. -Adsorbed by rock -- Kg/m3 kg/kg 0.0 0.0 0.5 0.0000017 1.0 0.0000017 / 0.0 0.5 1.0 0.0 0.0000017 0.0000017 /
-- Todd-Longstaff Mixing Parameters TLMIXPAR 1 1* / -- Polymer-Salt concentration for mixing -- Maximum polymer and salt concentration PLYMAX -- Ply conc. Salt conc. -- Kg/m3 kg/m3 1.0 0.0 / --Polymer-Rock Properties PLYROCK --dead residual --pore resistance --space factor 0.16 0.16 1.0 1.0
mass density
Ads. Index
2650.0 2650.0
2 2
63
Table 3: Incremental NPV for Surfactant flooding at concentration of 15kg/m3 for 4 yrs
Table 4: Incremental NPV for Surfactant flooding at concentration of 15kg/m3 for 6 yrs
64
Table 5: Incremental NPV for Surfactant flooding at concentration of 40kg/m3 for 4 yrs
Table 6: Incremental NPV for Polymer flooding at concentration of 0.3kg/m3 for 4 yrs
65
Table 7: Incremental NPV for Polymer flooding at concentration of 0.15kg/m3 for 4 yrs
Table 8: Incremental NPV for Polymer flooding at concentration of 11000kg/D for 4 years
66
Table 10: Incremental NPV for SP flooding with SURF@40kg/m3 & [email protected]/m3
67
Table 11: SP flooding with SURF@15kg/m3 & [email protected]/m3 followed WI and PLYI
Table 12: Incremental NPV for AS flooding with [email protected]/m3 & SURF@25kg/m3
68
Table 13: Incremental NPV for AS flooding with [email protected]/m3 & SURF@25kg/m3
Table 14: Incremental NPV for AS flooding with [email protected]/m3 & SURF@25kg/m3
69
70
Table 18: Incremental NPV for ASP flooding for 4 years slug AlK1.5_SURF15_PLY0.35
71
Table 19: Incremental NPV for AS for 2 yrs following by 1yr WI and 2yrs PLY
Table 20: Incremental NPV for AS for 4 yrs following by 1yr WI and 2yrs PLY
72
73
130 36 20 84 / LGR --MAXLGR MAXCLS MCOARS 0 0 848 / --WSEGDIMS -- 3 30 3 / --mlLGR -- maxlgr maxcls mcoars mamalg mxlalg lstack interp -4 2000 0 1 4 20 'INTERP' / TABDIMS --ntsfun ntpvt nssfun nppvt ntfip nrpvt ntendp 107 2 33 60 16 60 / -- WI_VFP_TABLES_080905.INC = 10-20 VFPIDIMS 30 20 20 / -- Table no. -- DevNew.VFP =1 -- E1h.VFP =2 -- AlmostVertNew.VFP = 3 -- GasProd.VFP =4 -- NEW_D2_GAS_0.00003.VFP = 5 -- GAS_PD2.VFP = 6 -- pd2.VFP = 8 (flowline south) -- pe2.VFP = 9 (flowline north) -- PB1.PIPE.Ecl = 31 -- PB2.PIPE.Ecl = 32 -- PD1.PIPE.Ecl = 33 -- PD2.PIPE.Ecl = 34 -- PE1.PIPE.Ecl = 35 -- PE2.PIPE.Ecl = 36 -- B1BH.Ecl = 37 -- B2H.Ecl = 38 -- B3H.Ecl = 39 -- B4DH. Ecl= 40 -- D1CH.Ecl = 41 -- D2H.Ecl = 42 -- D3BH.Ecl = 43 -- E1H.Ecl = 45 -- E3CH.Ecl = 47 -- K3H.Ecl = 48
74
----------------------------------------------------------Input of grid geometry ---------------------------------------------------------GRID NEWTRAN GRIDFILE 2 / -- optional for postprocessing of GRID MAPAXES 0. 100. 0. 0. 100. 0. / GRIDUNIT METRES / -- do not output GRID geometry file --NOGGF -- requests output of INIT file INIT MESSAGES 8*10000 20000 10000 1000 1* / PINCH 0.001 GAP 1* TOPBOT TOP/ NOECHO COARSEN -- I1 I2 J1 J2 K1 K2 NX NY NZ -6 9 43 92 1 3 1 1 1 / 10 12 45 96 1 3 1 1 1 / 13 18 48 100 1 3 1 1 1 / 19 19 58 90 1 3 1 1 1 / 20 25 49 70 1 3 1 1 1 / 26 30 50 90 1 3 1 1 1 / 31 41 65 91 1 3 1 1 1 / 31 38 61 64 1 3 1 1 1 / 31 34 56 60 1 3 1 1 1 / 38 41 92 100 2 3 1 1 1 / 31 33 54 55 1 3 1 1 1 / 31 32 53 53 1 3 1 1 1 / 35 36 59 60 1 3 1 1 1 / 39 39 63 64 1 3 1 1 1 / 35 37 93 98 2 3 1 1 1 / 36 37 99 99 2 3 1 1 1 / 33 34 94 96 2 3 1 1 1 / 30 34 93 93 1 3 1 1 1 / 30 37 92 92 1 3 1 1 1 /
75
39 41 102 102 2 3 1 1 1 / 19 19 48 57 1 3 1 1 1 /
6 9 43 92 19 22 1 1 1 / 10 12 45 96 19 22 1 1 1 / 13 18 48 100 19 22 1 1 1 / 19 19 58 90 19 22 1 1 1 / 20 25 49 70 19 22 1 1 1 / 26 30 50 90 21 22 1 1 1 / 31 41 65 91 22 22 1 1 1 / 31 38 61 64 21 22 1 1 1 / 31 34 56 60 21 22 1 1 1 / 38 41 92 100 22 22 1 1 1 / 31 33 54 55 21 22 1 1 1 /
76
31 32 35 36 39 39 35 37 36 37 33 34 30 34 30 37 39 41 19 19 31 35 37 40 31 35 37 40
53 53 21 22 1 1 1 / 59 60 21 22 1 1 1 / 63 64 21 22 1 1 1 / 93 98 22 22 1 1 1 / 99 99 22 22 1 1 1 / 94 96 22 22 1 1 1 / 93 93 22 22 1 1 1 / 92 92 22 22 1 1 1 / 102 102 22 22 1 1 1 / 48 57 21 22 1 1 1 / 52 58 60 64 52 58 60 64 1 1 1 1 22 22 22 22 1 1 1 1 1 1 1 1 1 1 1 1 / / / /
/ /
---------------------------------------------------------Grid and faults ----------------------------------------------------------- Simulation grid, with slooping faults: --- file in UTM coordinate system, for importing to -------------DecisionSpace INCLUDE './INCLUDE/GRID/IRAP_1005.GRDECL' / -- '/project/norne6/res/INCLUDE/GRID/IRAP_0704.GRDECL' / -INCLUDE './INCLUDE/GRID/ACTNUM_0704.prop' / --- Faults --INCLUDE './INCLUDE/FAULT/FAULT_JUN_05.INC' / -- Alteration of transmiscibility by use of the 'MULTFLT' -------keyword -INCLUDE './INCLUDE/FAULT/FAULTMULT_AUG-2006.INC' / -- '/project/norne6/res/INCLUDE/FAULT/FAULTMULT_JUN_05.INC' / -- Additional faults --Nord for C-3 (forlengelse av C_10) EQUALS MULTY 0.01 6 6 22 22 1 22 / / -- B-3 water EQUALS 'MULTX' 0.001 9 11 39 39 1 22 / 'MULTY' 0.001 9 11 39 39 1 22 /
77
26 29 39 39 1 22 /
---------------------------------------------------------Input of grid parametres ---------------------------------------------------------INCLUDE './INCLUDE/PETRO/PORO_0704.prop' / -INCLUDE './INCLUDE/PETRO/NTG_0704.prop' / -INCLUDE './INCLUDE/PETRO/PERM_0704.prop' / -- G segment north EQUALS PERMX 220 32 PERMX 220 33 PERMX 220 34 PERMX 220 35 PERMX 220 36 PERMX 220 37 PERMX 220 38 PERMX 220 39 PERMX 220 40 PERMX 220 41 /
32 33 34 35 36 37 38 39 40 41
2 2 2 2 2 2 2 2 2 2
2/ 2/ 2/ 2/ 2/ 2/ 2/ 2/ 2/ 2/
-- C-1H MULTIPLY PERMX 4 21 29 39 49 16 18 / PERMX 100 21 29 39 49 19 20 / / COPY PERMX PERMY / PERMX PERMZ / / -- Permz reduction is based on input from PSK -- based on same kv/kh factor -- ****************************************** -- CHECK! (esp. Ile & Tofte) -- ****************************************** MULTIPLY 'PERMZ' 0.2 1 46 1 112 1 1 / Garn 3 'PERMZ' 0.04 1 46 1 112 2 2 / Garn 2 'PERMZ' 0.25 1 46 1 112 3 3 / Garn 1 'PERMZ' 0.0 1 46 1 112 4 4 / Not (inactive anyway) 'PERMZ' 0.13 1 46 1 112 5 5 / Ile 2.2
78
'PERMZ' 0.13 1 46 1 112 6 6 / Ile 2.1.3 'PERMZ' 0.13 1 46 1 112 7 7 / Ile 2.1.2 'PERMZ' 0.13 1 46 1 112 8 8 / Ile 2.1.1 'PERMZ' 0.09 1 46 1 112 9 9 / Ile 1.3 'PERMZ' 0.07 1 46 1 112 10 10 / Ile 1.2 'PERMZ' 0.19 1 46 1 112 11 11 / Ile 1.1 'PERMZ' 0.13 1 46 1 112 12 12 / Tofte 2.2 'PERMZ' 0.64 1 46 1 112 13 13 / Tofte 2.1.3 'PERMZ' 0.64 1 46 1 112 14 14 / Tofte 2.1.2 'PERMZ' 0.64 1 46 1 112 15 15 / Tofte 2.1.1 'PERMZ' 0.64 1 46 1 112 16 16 / Tofte 1.2.2 'PERMZ' 0.64 1 46 1 112 17 17 / Tofte 1.2.1 'PERMZ' 0.016 1 46 1 112 18 18 / Tofte 1.1 'PERMZ' 0.004 1 46 1 112 19 19 / Tilje 4 'PERMZ' 0.004 1 46 1 112 20 20 / Tilje 3 'PERMZ' 1.0 1 46 1 112 21 21 / Tilje 2 'PERMZ' 1.0 1 46 1 112 22 22 / Tilje 1 /
---------------------------------------------------------Barriers ---------------------------------------------------------- MULTZ multiplies the transmissibility between blocks -- (I, J, K) and (I, J, K+1), thus the barriers are at the -- bottom of the given layer. -- Region barriers --INCLUDE './INCLUDE/PETRO/MULTZ_HM_1.INC' / --- Field-wide barriers -EQUALS 'MULTZ' 1.0 1 46 1 112 1 1 / Garn3 - Garn 2 'MULTZ' 0.05 1 46 1 112 15 15 / Tofte 2.1.1 - Tofte 1.2.2 'MULTZ' 0.001 1 46 1 112 18 18 / Tofte 1.1 - Tilje 4 'MULTZ' 0.00001 1 46 1 112 20 20 / Tilje 3 - Tilje 2 -- The Top Tilje 2 barrier is included as MULTREGT = 0.0 / -- Local barriers -INCLUDE './INCLUDE/PETRO/MULTZ_JUN_05_MOD.INC' /
79
NOECHO MINPV 500 / EQUALS 'MULTZ' 0.00125 26 29 30 37 10 10 / better WCT match for B-2H 'MULTZ' 0.015 19 29 11 30 8 8 / better WCT match for D-1CH 'MULTZ' 1 'MULTZ' .1 / 6 12 16 22 8 11 / for better WCT match for K-3H 6 12 16 22 15 15 / for better WCT match for K-3H
EDIT --------------------------------------------------------------------------------- modification related to HM of G-segment aug-2006 MULTIPLY 'TRANX' 0.1 30 46 72 112 2 2 / 'TRANX' 0.1 30 46 72 112 3 3 / 'TRANY' 5 30 46 72 112 2 2 / 'TRANY' 10 30 46 72 112 3 3 / -'TRANX' 10 29 29 67 70 1 3 / 'TRANY' 10 30 41 67 67 1 3 / -'TRANX' 0.05 34 34 76 95 1 3 / 'TRANY' 0.001 30 41 67 67 1 3 / Open against the main field -'TRANY' 0.5 30 30 90 93 1 3 / Increase TRANY against the well 'TRANY' 0.5 31 32 94 94 1 3 / Increase TRANY against the well --'TRANY' 0.5 31 31 87 93 1 3 / --'TRANY' 0.5 30 30 85 89 1 1 / 'TRANY' 2 30 30 72 82 1 3 / 'TRANY' 0.8 30 30 82 93 1 3 / --'TRANX' 10 34 34 92 95 1 3 / Increase TRANX trough the fault against the well 'TRANX' 0 34 34 90 91 1 3 / 'TRANX' 2 34 38 88 89 1 3/ --'TRANX' 2 35 36 93 95 1 3 / 'TRANX' 0.1 35 36 90 91 1 3 / 'TRANX' 10 35 38 95 98 1 3 / 'TRANX' 5 31 31 91 92 1 3 / Increase TRANX against the well --'TRANX' 2 31 33 92 95 1 3 / -'TRANY' 2 30 31 79 86 3 3 / 'TRANY' 3 30 30 86 86 2 2 / --'TRANY' 0.7 34 41 72 80 1 3 /
80
'TRANX' 2 31 31 87 94 1 3 / -'TRANY' 0.0004 37 41 71 71 1 3 / 'TRANY' 2 30 31 87 93 2 3 / 'TRANX' 5 34 34 88 90 1 3 / -'TRANY' 1.5 33 35 94 96 2 3 / -'TRANX' 2 30 41 68 70 1 3 / Increase trans around F-4H -/
EQUALS 'TRANY' 20 31 31 85 85 1 3 / SET TRANY ulik 0 trougth the fault 'TRANY' 30 30 30 93 93 2 2 / 'TRANY' 30 32 32 84 84 1 3 / 'TRANY' 30 30 30 93 93 3 3 / --'TRANY' 30 31 32 95 95 2 3 / 'TRANY' 30 31 32 94 94 1 1 / 'TRANY' 20 33 33 96 96 2 3 / 'TRANY' 20 34 34 97 97 2 3 / --'TRANX' 0 33 33 71 81 1 3 / set the fault tight 'TRANX' 0 34 34 76 85 1 3 / -'TRANY' 0 33 33 71 81 1 3 / Set the fault tigt 'TRANY' 0 34 34 76 85 1 3 / -'TRANY' 0 33 36 71 71 1 3 / 'TRANX' 0 34 41 71 71 1 3 / -'TRANY' 0 33 33 71 72 1 3 / Decrease TRANY trougth the fault -'TRANX' 0 34 34 73 75 1 3 / Set the fault tight 'TRANY' 0 34 34 71 75 1 3 / -/ -------------------------------------------------------------------------------PROPS ---------------------------------------------------------------------------------- Input of fluid properties and relative permeability ---------------------------------------------------------NOECHO
-- Input of PVT data for the model -- Total 2 PVT regions (region 1 C,D,E segment, region 2 Gsegment) -INCLUDE './INCLUDE/PVT/PVT-WET-GAS.DATA' / INCLUDE
81
'./INCLUDE/ASP.inc' / TRACER 'SEA' 'WAT' / 'HTO' 'WAT' / 'S36' 'WAT' / '2FB' 'WAT' / '4FB' 'WAT' / 'DFB' 'WAT' / 'TFB' 'WAT' / / ------------------------------------------------------------ initialization and relperm curves: see report blabla ------------------------------------------------------------ rel. perm and cap. pressure tables --INCLUDE './INCLUDE/RELPERM/HYST/swof_mod4Gseg_aug-2006.inc' / -- '/project/norne6/res/INCLUDE/RELPERM/HYST/swof.inc' / --Sgc=10 0.000000or g-segment -INCLUDE './INCLUDE/RELPERM/HYST/sgof_sgc10_mod4Gseg_aug-2006.inc' / -- '/project/norne6/res/INCLUDE/RELPERM/HYST/sgof_sgc10.inc' / ---INCLUDE -- './INCLUDE/RELPERM/HYST/waghystr_mod4Gseg_aug-2006.inc' / -- '/project/norne6/res/INCLUDE/RELPERM/HYST/waghystr.inc' / --RPTPROPS -- 1 1 1 5*0 0 /
--------------------------------------------------------------------------------
REGIONS -INCLUDE './INCLUDE/PETRO/FIPNUM_0704.prop' / -INCLUDE './INCLUDE/PETRO/SATNUM_0704.prop' / EQUALS 'SATNUM' 102 30 41 76 112 1 1 / 'SATNUM' 103 30 41 76 112 2 2 / 'SATNUM' 104 30 41 76 112 3 3 / / MISCNUM 113344*1/
82
-INCLUDE './INCLUDE/PETRO/IMBNUM_0704.prop' / EQUALS 'SATNUM' 102 30 41 76 112 1 1 / 'SATNUM' 103 30 41 76 112 2 2 / 'SATNUM' 104 30 41 76 112 3 3 / / -INCLUDE './INCLUDE/PETRO/PVTNUM_0704.prop' / EQUALS 'PVTNUM' 1 1 46 1 112 1 22 / /
INCLUDE './INCLUDE/PETRO/EQLNUM_0704.prop' / -- extra regions for geological formations and numerical layers INCLUDE './INCLUDE/PETRO/EXTRA_REG.inc' / --------------------------------------------------------------------------------SOLUTION RPTRST BASIC=6 / RPTSOL FIP=3 / ---------------------------------------------------------------------------------- equilibrium data: do not include this file in case of RESTART --INCLUDE './INCLUDE/PETRO/E3.prop' / -- restart date: only used in case of a RESTART, remember to use SKIPREST --RESTART -- 'BASE_30-NOV-2005' 360 / AT TIME 3282.0 DAYS ( 1-NOV-2006) THPRES 1 2 0.588031 / 2 1 0.588031 / 1 3 0.787619 / 3 1 0.787619 / 1 4 7.00083 / 4 1 7.00083 / / -- initialise injected tracers to zero TVDPFSEA 1000 0.0 5000 0.0 / TVDPFHTO 1000 0.0 5000 0.0 /
83
TVDPFS36 1000 0.0 5000 0.0 / TVDPF2FB 1000 0.0 5000 0.0 / TVDPF4FB 1000 0.0 5000 0.0 / TVDPFDFB 1000 0.0 5000 0.0 / TVDPFTFB 1000 0.0 5000 0.0 / ------------------------------------------------------------------------------SUMMARY RUNSUM SEPARATE EXCEL -INCLUDE './INCLUDE/SUMMARY/summary.data' / -INCLUDE './INCLUDE/SUMMARY/extra.inc' / -INCLUDE './INCLUDE/SUMMARY/tracer.data' / -INCLUDE './INCLUDE/SUMMARY/gas.inc' / -INCLUDE './INCLUDE/SUMMARY/wpave.inc' / --------------------------------------------------------------------------------
SCHEDULE -- use SKIPREST in case of RESTART --SKIPREST -- No increase in the solution gas-oil ratio?! DRSDT 0 / -- Use of WRFT in order to report well perssure data after first -- opening of the well. The wells are perforated in the entire reservoir -- produce with a small rate and are squeesed after 1 day. This pressure -- data can sen be copmared with the MDT pressure points collected in the -- well.
84
NOECHO ---------------------------------------------=======Production Wells========---------------------------------------------INCLUDE './INCLUDE/VFP/DevNew.VFP' / -INCLUDE './INCLUDE/VFP/E1h.VFP' / -INCLUDE './INCLUDE/VFP/NEW_D2_GAS_0.00003.VFP' / -INCLUDE './INCLUDE/VFP/GAS_PD2.VFP' / -INCLUDE './INCLUDE/VFP/AlmostVertNew.VFP' / -INCLUDE './INCLUDE/VFP/GasProd.VFP' /
-- 01.01.07 new VFP curves for producing wells, matched with the latest well tests in Prosper. lmarr -INCLUDE './INCLUDE/VFP/B1BH.Ecl' / -INCLUDE './INCLUDE/VFP/B2H.Ecl' / -INCLUDE './INCLUDE/VFP/B3H.Ecl' / -INCLUDE './INCLUDE/VFP/B4DH.Ecl' / -INCLUDE './INCLUDE/VFP/D1CH.Ecl' / -INCLUDE './INCLUDE/VFP/D2H.Ecl' / -INCLUDE './INCLUDE/VFP/D3BH.Ecl' / -INCLUDE './INCLUDE/VFP/E1H.Ecl' / -INCLUDE './INCLUDE/VFP/E3CH.Ecl' / -INCLUDE './INCLUDE/VFP/K3H.Ecl' /
--------------------------------------------
85
--=======Production Flowlines========----------------------------------------------- 16.5.02 new VFP curves for southgoing PD1,PD2,PB1,PB2 flowlines -> pd2.VFP -INCLUDE './INCLUDE/VFP/pd2.VFP' / --- 16.5.02 new VFP curves for northgoing PE1,PE2 flowlines -> pe2.VFP -INCLUDE './INCLUDE/VFP/pe2.VFP' /
-- 24.11.06 new matched VLP curves for PB1 valid from 01.07.06 -INCLUDE './INCLUDE/VFP/PB1.PIPE.Ecl' / --24.11.06 new matched VLP curves for PB2 valid from 01.07.06 -INCLUDE './INCLUDE/VFP/PB2.PIPE.Ecl' / --24.11.06 new matched VLP curves for PD1 valid from 01.07.06 -INCLUDE './INCLUDE/VFP/PD1.PIPE.Ecl' / --24.11.06 new matched VLP curves for PD2 valid from 01.07.06 -INCLUDE './INCLUDE/VFP/PD2.PIPE.Ecl' / --24.11.06 new matched VLP curves for PE1 valid from 01.07.06 -INCLUDE './INCLUDE/VFP/PE1.PIPE.Ecl' / --24.11.06 new matched VLP curves for PE2 valid from 01.07.06 -INCLUDE './INCLUDE/VFP/PE2.PIPE.Ecl' /
---------------------------------------------=======INJECTION FLOWLINES 08.09.2005 --------------------------------------------- VFPINJ nr. 10 Water injection flowline WIC -INCLUDE './INCLUDE/VFP/WIC.PIPE.Ecl' / -- VFPINJ nr. 11 Water injection flowline WIF -INCLUDE './INCLUDE/VFP/WIF.PIPE.Ecl' / ---------------------------------------------======= INJECTION Wells 08.09.2005 --------------------------------------------
========--
========--
86
-- VFPINJ nr. 12 Water injection wellbore Norne C-1H -INCLUDE './INCLUDE/VFP/C1H.Ecl' / -- VFPINJ nr. 13 Water injection wellbore Norne C-2H -INCLUDE './INCLUDE/VFP/C2H.Ecl' / -- VFPINJ nr. 14 Water injection wellbore Norne C-3H -INCLUDE './INCLUDE/VFP/C3H.Ecl' / -- VFPINJ nr. 15 Water injection wellbore Norne C-4H -INCLUDE './INCLUDE/VFP/C4H.Ecl' / -- VFPINJ nr. 16 Water injection wellbore Norne C-4AH -INCLUDE './INCLUDE/VFP/C4AH.Ecl' / -- VFPINJ nr. 17 Water injection wellbore Norne F-1H -INCLUDE './INCLUDE/VFP/F1H.Ecl' / -- VFPINJ nr. 18 Water injection wellbore Norne F-2H -INCLUDE './INCLUDE/VFP/F2H.Ecl' / -- VFPINJ nr. 19 Water injection wellbore Norne F-3 H -INCLUDE './INCLUDE/VFP/F3H.Ecl' / -- VFPINJ nr. 20 Water injection wellbore Norne F-4H -INCLUDE './INCLUDE/VFP/F4H.Ecl' /
TUNING 1 10 0.1 0.15 3 0.3 0.3 1.20 / 5* 0.1 0.0001 0.02 0.02 / --2* 40 1* 15 / / -- only possible for ECL 2006.2+ version ZIPPY2 'SIM=4.2' 'MINSTEP=1E-6' / /
--WSEGITER --/
87
-- PI reduction in case of water cut -INCLUDE './INCLUDE/PI/pimultab_low-high_aug-2006.inc' / -- History and prediction --INCLUDE './INCLUDE/ASP4YRS.SCHEDULE' /
END
88