Introduction to IORSIM The National IOR Centre of Norway, University of Stavanger 2 The National IOR Centre of Norway, IRIS 3 The National IOR Centre of Norway, IFE 4 The National IOR Centre of Norway, Schlumberger Stavanger Research
Recovery Factor RF = E PS E S E D E C IOR A. Muggeridge et al., Recovery rates, enhanced oil recovery and technological limits, Phil. Trans. R. Soc. A 372: 2232, 23
MORE WELLS MORE OIL?
Why IORSim? We know fluid chemistry affects flow properties on core scale (~ cm). Compaction and wettability in chalk 2. Water diverging How to translate core scale knowledge/processes to field scale? Implementation of IOR processes Interpretation of field data
NaCl Effect of sulphate on oil production NaCl NaCl Na 2 SO 4 Na 2 SO 4 Na 2 SO 4 R. Ahsan, M. V. Madland, F. Bratteli, A. Hiorth A STUDY OF SULPHATE IONS - EFFECTS ON AGEING AND IMBIBITION CAPILLARY PRESSURE CURVE SCA, 34 (22)
Molality Ion (Total).28 - HCO 3.24 2- CO 3.27 Cl -.56 SO 4 2- % 8% 6% 4% 2% % Anion complexes In Seawater SO4-- HCO3- CO3-- Cl- Free Ion Ca-anion Mg-anion Na-anion K-anion
IORSim overview Originally developed as BUGSIM, in order to interpret microbial enhanced oil recovery (MEOR) experiments 3D grids (Cartesian, radial), based on a black oil formulation Includes several capabilities : Bacteria, substrates and nutrients (MEOR) Silicate Polymers Surfactants Geochemical components (e.g., ions, aqueous Tracers complexes, minerals, charged surfaces) Aims to provide physically sound models for capturing EOR mechanisms
D integration Global level Flow Local (Block) level non linear chemistry Block-by-block upstream integration Chemical reactions Mass in Dissolution, precipitation, ion exchange, surface complexation Mass out 2 N c in, (t + Δt) c in,2 (t + Δt) c in,n (t + Δt)
Can we extend to 2D? Yes, if there are no flow loops (tear up the loops) 3x3 Example 7 28 39 4 5 6 Renumbering of ECLIPSE grid implicit integration of the concentrations if no flow loops are detected If flow loops, explicit couplings are inserted to tear up the flow loops The algorithm is fast ( O(N) ) 7 82 93 Similar ideas: Natvig, J.R. and Lie, K.A. [28]. J. Comp. Phys. 227(24), 8 24.
The Approach IORSim Technology Eclipse Reservoir simulator Restart Files Sw, Po, Pw, qw IORSim advect components Modular interface Geochemistry Silicate Oil Rate Water Rate Mg SO4
IORSim Backward Coupling Stop Eclipse process Update SATNUM X = EOR response, e.g.: Low sal reduce Sor Silicate reduce Perm SATNUM -2
Sodium Silicate Chemistry Gelation time gelation 5wt% Sodium Silicate Na 2 O:(SiO 2 ) n + HCl Oligomeric silicic acid Phase separation Nano sized aggregation A. Stavland, H. Jonsbråten, O. Vikane, K. Skrettingland and H. Fischer, In-depth Water Diversion Using Sodium Silicate Preparation www.clker.com/ for Single Well Field Pilot on Snorre, 6th European Symposium on Improved Oil Recovery Cambridge, UK, 2-4 April 2 Sebastian Wilhelm and Matthias Kind, Polymers 25, 7, 254 252; doi:.339/polym72528
Sodium silicate model 4 aqueous species + mineral phase. Silica 2. HCl 3. Mobile Gel 4. Calcium 5. Immobile Gel Sodium Silicate Activator Nano size gel particles Divalent ions Rigid gel Perm Reduction
Example: Silicate Injection Permeability field Temperature profile md 4mD 2C 95C
Designing the silicate slug Producer no treatment Poor sweep Main parameters Slug volume Fracture size Silicate System Gelation kinetics Position of the Temperature front BHP
5wt% silica injection for 65 days Water saturation Water inj Oil prod..7 Water prod Precipitation of silicate Permeability modification
5wt% silica injection for 65 days Water saturation Water inj Oil prod..7 Water prod Precipitation of silicate Permeability modification
5wt% silica injection for 65 days Water saturation Water inj Oil prod..7 Water prod Precipitation of silicate Permeability modification
5wt% silica injection for 65 days Water saturation Water inj Oil prod..7 Water prod Precipitation of silicate Permeability modification
5wt% silica injection for 65 days Water saturation Water inj Oil prod..7 Water prod Precipitation of silicate Permeability modification
5wt% silica injection for 65 days Water saturation Water inj Oil prod..7 Water prod Precipitation of silicate Permeability modification
5wt% silica injection for 65 days Water saturation Water inj Oil prod..7 Water prod Precipitation of silicate Permeability modification
5wt% silica injection for 65 days Water saturation Water inj Oil prod..7 Water prod Precipitation of silicate Permeability modification
Tracer curve after treatment After treatment After Tracer Injection
Convergence Oil production water production ECLIPSE ΔT > IORSim ΔT Convergence dependent on both Eclipse report step and IORSim timestep
Does it converge to the right solution? IORCoreSim: In-house reservoir simulator
Heterogeneous case Permeability field Immobile Gel
Summary IORSim and ECLIPSE coupling works well Block sorting solve implicit & modular Can be extended to other simulators INTERSECT OPM Applied to realistic cases Speed Not critical, need adaptive time steps 2. order scheme Local grid refinement Positive technical potential for silicate Fracture case Heterogeneous case
Acknowledgement The authors acknowledge the Research Council of Norway and the industry partners; ConocoPhillips Skandinavia AS, Aker BP ASA, Eni Norge AS, Maersk Oil Norway AS, Statoil Petroleum AS, ENGIE E&P NORGE AS, Lundin Norway AS, Halliburton AS, Schlumberger Norge AS, Wintershall Norge AS, DEA Norge AS of The National IOR Centre of Norway for support.