Flow an open source research tool for reservoir simulation

Transcription

1 Flow an open source research tool for reservoir simulation Robert Klöfkorn IRIS & The National IOR Centre of Norway April 27, 2016 Task 6 Reservoir Simulation Tools Pål (PostDoc UiS) Trine (PostDoc IRIS) Anna (PhD UiS/IRIS) Svein (Prof. UiS) Task Leader Steinar (Prof. UiS) Ove (IRIS) Robert (IRIS) Task Leader

2 Outline 1. Current status of OPM-Flow 2. Research Activities in Task 6 Higher order schemes for reactive advection-diffusion Modeling reactive advection-diffusion 3. Summary 2 / 23 Robert Klöfkorn

3 Outline 1. Current status of OPM-Flow 2. Research Activities in Task 6 Higher order schemes for reactive advection-diffusion Modeling reactive advection-diffusion 3. Summary 3 / 23 Robert Klöfkorn

4 Open Porous Media Initiative (OPM) OPM encourages open innovation and reproducible research for modeling and simulation of porous media processes. All OPM software is open-source. 4 / 23 Robert Klöfkorn

5 Open Porous Media Initiative (OPM) OPM encourages open innovation and reproducible research for modeling and simulation of porous media processes. All OPM software is open-source. Main contributors Statoil ASA (Alf B. Rustad, Joakim Hove, and many more) SINTEF ICT (group of K.A. Lie, Atgeirr Rasmussen and many more) IRIS Energy (Reservoir group) Ceetron Solutions AS A. Lauser (Poware) M. Blatt (Dr. Blatt HPC Simulation and Service) 4 / 23 Robert Klöfkorn

6 Open Porous Media Initiative (OPM) OPM encourages open innovation and reproducible research for modeling and simulation of porous media processes. All OPM software is open-source. Main contributors Statoil ASA (Alf B. Rustad, Joakim Hove, and many more) SINTEF ICT (group of K.A. Lie, Atgeirr Rasmussen and many more) IRIS Energy (Reservoir group) Ceetron Solutions AS A. Lauser (Poware) M. Blatt (Dr. Blatt HPC Simulation and Service) Other (current and former) partners: Uni Research CIPR University of Bergen Total University of Stuttgart University of Heidelberg 4 / 23 Robert Klöfkorn

7 Open Porous Media Initiative (OPM) 5 / 23 Robert Klöfkorn

8 Open Porous Media Initiative (OPM) Installation (releases twice a year): Ubuntu/Red Hat packing system Run OPM-Flow using a virtual machine Install from source on Linux and Mac OS X 5 / 23 Robert Klöfkorn

9 Open Porous Media Initiative (OPM) Installation (releases twice a year): Open-source (GPL v3+) Ubuntu/Red Hat packing system Run OPM-Flow using a virtual machine Install from source on Linux and Mac OS X 5 / 23 Robert Klöfkorn

10 OPM-Flow Features Fully implicit model formulation based on Automatic Differentiation: black-oil with dissolved gas and vaporized oil rock-dependent capillary and relative-permeability curves end-point scaling and hysteresis oil vaporization control 6 / 23 Robert Klöfkorn

11 OPM-Flow Features Fully implicit model formulation based on Automatic Differentiation: black-oil with dissolved gas and vaporized oil rock-dependent capillary and relative-permeability curves end-point scaling and hysteresis oil vaporization control EOR options: Todd-Longstaff type polymer model with adsorption, dead-pore space, permeability reduction, and shear effects (Flow-polymer) extra component equation(s), such as a solvent model (Flow-solvent) 6 / 23 Robert Klöfkorn

12 OPM-Flow Features Fully implicit model formulation based on Automatic Differentiation: black-oil with dissolved gas and vaporized oil rock-dependent capillary and relative-permeability curves end-point scaling and hysteresis oil vaporization control EOR options: Todd-Longstaff type polymer model with adsorption, dead-pore space, permeability reduction, and shear effects (Flow-polymer) extra component equation(s), such as a solvent model (Flow-solvent) Performance: 2014: 60 slower than Eclipse100 on SPE cases 6 / 23 Robert Klöfkorn

13 OPM-Flow Features Fully implicit model formulation based on Automatic Differentiation: black-oil with dissolved gas and vaporized oil rock-dependent capillary and relative-permeability curves end-point scaling and hysteresis oil vaporization control EOR options: Todd-Longstaff type polymer model with adsorption, dead-pore space, permeability reduction, and shear effects (Flow-polymer) extra component equation(s), such as a solvent model (Flow-solvent) Performance: 2014: 60 slower than Eclipse100 on SPE cases 2015: 8 slower than Eclipse100 on Norne field model 6 / 23 Robert Klöfkorn

14 OPM-Flow Features Fully implicit model formulation based on Automatic Differentiation: black-oil with dissolved gas and vaporized oil rock-dependent capillary and relative-permeability curves end-point scaling and hysteresis oil vaporization control EOR options: Todd-Longstaff type polymer model with adsorption, dead-pore space, permeability reduction, and shear effects (Flow-polymer) extra component equation(s), such as a solvent model (Flow-solvent) Performance: 2014: 60 slower than Eclipse100 on SPE cases 2015: 8 slower than Eclipse100 on Norne field model 2016: 2.2 slower than Eclipse100 on Norne field model 6 / 23 Robert Klöfkorn

15 OPM-Flow Features Fully implicit model formulation based on Automatic Differentiation: black-oil with dissolved gas and vaporized oil rock-dependent capillary and relative-permeability curves end-point scaling and hysteresis oil vaporization control EOR options: Todd-Longstaff type polymer model with adsorption, dead-pore space, permeability reduction, and shear effects (Flow-polymer) extra component equation(s), such as a solvent model (Flow-solvent) Performance: 2014: 60 slower than Eclipse100 on SPE cases 2015: 8 slower than Eclipse100 on Norne field model 2016: 2.2 slower than Eclipse100 on Norne field model Other features: Code review via github merge requests Performance monitoring on linuxbenchmarking.org Parallel version of flow working for limited number of cores 6 / 23 Robert Klöfkorn

16 NORNE field model (SOIL) Visualization with ResInsight 7 / 23 Robert Klöfkorn

17 NORNE field model (SWAT) Visualization with ResInsight 8 / 23 Robert Klöfkorn

18 NORNE field model (D-2H) WBHP ECL OPM-Flow WOPR ECL OPM-Flow days days ECL OPM-Flow 3x x10 6 ECL OPM-Flow WWPR WGPR 2x x10 6 1x days days 9 / 23 Robert Klöfkorn

19 NORNE field model (E-3H) ECL OPM-Flow ECL OPM-Flow WBHP WOPR days days ECL OPM-Flow ECL OPM-Flow WWPR WGPR days days 10 / 23 Robert Klöfkorn

20 OPM-Flow future directions Ongoing work: More field cases (Model 2,...) Enhanced solvent model (CO 2,...) Use OPM-Flow in history matching activities in Task 7 Coupling with IORSim Improved overall performance Parallelization to many cores 11 / 23 Robert Klöfkorn

21 OPM-Flow future directions Ongoing work: More field cases (Model 2,...) Enhanced solvent model (CO 2,...) Use OPM-Flow in history matching activities in Task 7 Coupling with IORSim Improved overall performance Parallelization to many cores Long term: New non-linear solvers and splitting methods Higher order methods discretization methods Seamless integration with history matching and UQ 11 / 23 Robert Klöfkorn

22 OPM-Flow future directions Ongoing work: More field cases (Model 2,...) Enhanced solvent model (CO 2,...) Use OPM-Flow in history matching activities in Task 7 Coupling with IORSim Improved overall performance Parallelization to many cores Long term: New non-linear solvers and splitting methods Higher order methods discretization methods Seamless integration with history matching and UQ Except a bumpy ride. 11 / 23 Robert Klöfkorn

23 Outline 1. Current status of OPM-Flow 2. Research Activities in Task 6 Higher order schemes for reactive advection-diffusion Modeling reactive advection-diffusion 3. Summary 12 / 23 Robert Klöfkorn

24 Reactive Advection-Diffusion Problems Motivation Unlike water fronts, polymer fronts described by linear waves are not self-sharpening. φ tc + (uc D(u) c ) = S(c) 13 / 23 Robert Klöfkorn

25 Reactive Advection-Diffusion Problems Motivation Unlike water fronts, polymer fronts described by linear waves are not self-sharpening. φ tc + (uc D(u) c ) = S(c) Low Order Finite Volume 13 / 23 Robert Klöfkorn

26 Reactive Advection-Diffusion Problems Motivation Unlike water fronts, polymer fronts described by linear waves are not self-sharpening. φ tc + (uc D(u) c ) = S(c) Low Order Finite Volume Higher Order Finite Volume 13 / 23 Robert Klöfkorn

27 Reactive Advection-Diffusion Problems Motivation Unlike water fronts, polymer fronts described by linear waves are not self-sharpening. φ tc + (uc D(u) c ) = S(c) Low Order Finite Volume Higher Order Finite Volume R. Klöfkorn, D. Kröner, and M. Ohlberger. Local adaptive methods for convection dominated problems. Int. J. Numer. Methods Fluids, 40(1-2):79 91, / 23 Robert Klöfkorn

28 Fully implicit higher order schemes for polymer flooding Motivation Unlike water fronts, polymer fronts described by linear waves are not self-sharpening. tρφs o + (ρ ou o(s w,o) ) = 0 tρφs w + (ρ wu w(s w,o) ) = 0 t R(c, s w) + (cρ wu wp(s o, s w) ) = 0 Challenges: Fully implicit fully coupled formulation Integration of slope limiter techniques into OPM s Automatic Differentiation framework 14 / 23 Robert Klöfkorn

29 Fully implicit higher order schemes for polymer flooding Grid = Grid = Challenges: Fully implicit fully coupled formulation Integration of slope limiter techniques into OPM s Automatic Differentiation framework Trine S. Mykkeltvedt (IRIS), Xavier Raynaud (SINTEF), Knut-Andreas Lie (SINTEF). Fully implicit higher-order scheme applied to polymer flooding. In preparation. 14 / 23 Robert Klöfkorn

30 Higher order schemes in 3D on regular meshes tc + (uc) = 0 Error Change w.r.t. CPU on Cartesian Grid Logarithm of the error L1-norm First Order Second Order EOC = 1 EOC = Logarithm of the CPU time Results on the Cartesian grid for the second order scheme (top-left) and the first order scheme (bottom-right). 15 / 23 Robert Klöfkorn

31 Higher order schemes in 3D on regular meshes tc + (uc) = 0 Error Change w.r.t. CPU on Cartesian Grid Logarithm of the error L1-norm First Order Second Order EOC = 1 EOC = Logarithm of the CPU time A. Dedner and R. Klöfkorn. A Generic Stabilization Approach for Higher Order Discontinuous Galerkin Methods for Convection Dominated Problems. J. Sci. Comput. 47(3): , Results on the Cartesian grid for the second order scheme (top-left) and the first order scheme (bottom-right). 15 / 23 Robert Klöfkorn

32 Higher order schemes on polyhedral meshes Results on a grid consisting of hexagonal prims for the second order scheme (top-left) and the first order scheme (bottom-right). Hexagonal prism grid 16 / 23 Robert Klöfkorn

33 Higher order schemes on polyhedral meshes Results on a grid consisting of hexagonal prims for the second order scheme (top-left) and the first order scheme (bottom-right). Hexagonal prism grid Anna Kvashchuk (UiS/IRIS), R. Klöfkorn (IRIS), and M. Nolte (Uni Freiburg). Higher Order Finite Volume Schemes on Polygonal and Polyhedral Grids. In preparation. Anna Kvashchuk (UiS/IRIS), IOR Norway 2016, Poster. 16 / 23 Robert Klöfkorn

34 UiS Postdoc activity: Paper 1 Issue: How would lab scale experiment behave on higher scale? Consider MgCl2 which reacts with chalk in similar way as seawater: Dissolution + precipitation of minerals Exchange of ions at the surface Obtain reaction kinetic parameters from core scale experiment

35 Scaled model for MgCl2 injection into fracture-matrix geometry Advection along fracture Diffusion between matrix and fracture Ion exchange and dissolution + precipitation in matrix

36 Important dimensionless numbers Alfa: time scale ratio of diffusion to advection Beta: volume ratio of matrix to fracture Gamma: time scale ratio of reaction to advection Time scale given in injected fracture volumes (FV)

37 Ex 3: Flow in different regions 3 regions with different fracture spacing and apertures, same volume Case with similar flux to each region (depends on parameters) Produced compositions from regions 1-3 and combined composition Even if the regions receive same amount of brine the interaction is very different in each region (alfa * beta and beta * gamma) P. Ø. Andersen and S. Evje. A Model for reactive flow in fractured porous media. Chem. Eng. Sci., 2016

38 Outline 1. Current status of OPM-Flow 2. Research Activities in Task 6 Higher order schemes for reactive advection-diffusion Modeling reactive advection-diffusion 3. Summary 21 / 23 Robert Klöfkorn

39 Summary SINTEF ICT Collaborations outside of the IOR Centre Cluster of Excellence in Simulation Technology (University of Stuttgart) Collaboration through DUNE, especially EXA-DUNE: Flexible PDE Solvers, Numerical Methods, and Applications 1. RWTH Aachen 2. Heidelberg University 3. University of Freiburg 4. University of Heidelberg 5. University of Münster 6. University of Stuttgart 7. University of Warwick and Imperial College London... Colorado School of Mines 22 / 23 Robert Klöfkorn

40 Summary SINTEF ICT Collaborations outside of the IOR Centre Cluster of Excellence in Simulation Technology (University of Stuttgart) Collaboration through DUNE, especially EXA-DUNE: Flexible PDE Solvers, Numerical Methods, and Applications (University of Heidelberg and others)... Colorado School of Mines Task 6 research contributes to improved reservoir simulation capabilities Research is supposed to be integrated into OPM 22 / 23 Robert Klöfkorn

41 Summary SINTEF ICT Collaborations outside of the IOR Centre Cluster of Excellence in Simulation Technology (University of Stuttgart) Collaboration through DUNE, especially EXA-DUNE: Flexible PDE Solvers, Numerical Methods, and Applications (University of Heidelberg and others)... Colorado School of Mines Task 6 research contributes to improved reservoir simulation capabilities Research is supposed to be integrated into OPM Flow (OPM-SIMULATORS) allows for simulation of field scale models Performance is currently addressed (right now within reach of the commercial simulators) Improve research transfer to industry relevant cases Upcoming release ( OPM Meeting at SINTEF in Oslo, June 1-2, / 23 Robert Klöfkorn

42 Summary SINTEF ICT Collaborations outside of the IOR Centre Cluster of Excellence in Simulation Technology (University of Stuttgart) Collaboration through DUNE, especially EXA-DUNE: Flexible PDE Solvers, Numerical Methods, and Applications (University of Heidelberg and others)... Colorado School of Mines Task 6 research contributes to improved reservoir simulation capabilities Research is supposed to be integrated into OPM Flow (OPM-SIMULATORS) allows for simulation of field scale models Performance is currently addressed (right now within reach of the commercial simulators) Improve research transfer to industry relevant cases Upcoming release ( OPM Meeting at SINTEF in Oslo, June 1-2, 2016 Thank you for your attention. 22 / 23 Robert Klöfkorn

43