Sleipner Benchmark Dataset and Model Comparison

Transcription

1 Sleipner Benchmark Dataset and Model Comparison IEAGHG Modeling and Risk Assessment Network Meeting June 2013 Sarah Gasda, Center for Integrated Petroleum Research, Uni Research, Bergen, Norway

2 Sleipner CO 2 Injection Injection began in Mt/yr injection rate special topic Over 14 Mtons stored through C0 2 Sequestration Arts et al Figure 2 Scheme of the CO 2 injection at Sleipner. Arts et al. 2008

3 Injection with Data Collection Accumulated mass (million tons) CO 2 Injection rate ~0.9 Mt/yr Injected CO2 Cumulative volume stored >14 Mt Continuous monitoring programme seismic surveys Date seafloor Seismic mapping surveys highres seismic CSEM survey gravity surveys CSEM survey Gravity surveys Seafloor mapping

4 Depth Thickne Utsira Formation Shallow marine sandstones ~200-m thick saline aquifer Overlain by Nordland shale sequence Average porosity 34% Permeability 1-5 Darcy 98% sand content Depth Intra-formational shales Upper sand wedge above 5-m shale Extension of the Utsira Sand 100 km Thickness Rob Arts, Andy Chadwick & Ola Eiken E e reservoir consists of two reservoir units: a main lower 'Utsira Sand' unit and an upper 'Sand Wedge' unit, separated by a 6.5 m le unit (Figure 1). A number of thin intra-formational shale layers (approximately 0.5 to 2 m thick) are present within the Utsira nd unit. These appear as spikes on the gamma-ray log, and affect the flow behavior of the CO 2 in the reservoir. Chadwick et al Singh et al km Total area of the Utsira Sand is approximately km 2

5 Utsira Modeling and Monitoring Pre-injection simulation study Eclipse 100, 250 m gridblocks 20 yr injection plume < 3 km in diameter 30 bar wellhead pressure increase 18% dissolution Since then Arts et al Baklid and Korbøl, 1996

6 Layer 9 Benchmark Cumulative Mass (Mt) Singh et al. 2010

7 CO 2 Storage Benchmarks Pruess, et al. An Intercomparison Study of Simulation Models for Geologic Sequestration of CO 2, Energy Numerous problems covering a multitude of storage issues. Class et al. A benchmark-study on problems related to CO 2 storage in geologic formations, Comp. Geosci., Flow simple leaky well, and more realistic formations Nordbotten et al. Uncertainties in practical simulation of CO 2 storage. IJGCC, submitted. CO 2 flow and trapping in a simple brine aquifer

8 Layer 9 Simulations Conformance: History match observations and simulations observation flow simulations depth contours Core permeabilities 2-3 Darcy Well permeabilities 1-8 Darcy 3 Darcy 3 Darcy E-W 10 Darcy N-S 3/10 Darcy (higher temp) 1. Single feeder in model 2. Baseline depth map assumes laterally uniform seismic velocity in overburden. Introducing velocity changes within Chadwick known range of velocity et variation al could modify depth map to allow observed lateral migration. TOUGH2 Migration simulator Eclipse 100/300 Others SPE Singh et al. 2010

9 Sleipner Benchmark Downloads 38 downloads as of early countries Europe, N. America, and Asia Number of downloads US Norway Canada Australia Germany UK France Russia Denmark Spain Korea

10 Layer 9 Benchmark Results Model Institution Processes Assumptions Sensitivities Comments MUFITS (Multiphase Filtration Transport Simulator) Moscow State University CO 2 dissolution in water viscous and gravity forces nonisothermal (conduction and heat of dissolution) no NaCl solution in water, no capillary pressure, no mineralization initial formation temperature ( K) assumes CO 2 is heated upon injection VESA (Vertical Equilibrium w/subscale Analytical Model) Uni CIPR/UiB Vertical Equil Isothermal Sharp Interface Upscaled Convective Dissolution gravity segregated no capillary pressure porosity permeability topography CO 2 density MRST (Matlab Reservoir Simulator Toolbox) SINTEF Vertical Equil Black Oil Compositional depending on modeling choice temperature temporal available as opensource Matlab simulator (MRST- VE) download from SINTEF MPATH/ ECLIPSE Permedia/Statoil invasion percolation Eclipse 100/300 viscous flow is neglected Darcy-flow dominates pressure MPATH used to model full-field percolation More to come

11 Seismic Data Migration Simulator MUFITS Temperature=36.8 C, Mean CO 2 Density = kg/m 3 VESA Porosity = 36%, CO 2 Density = 700 kg/m Incomp seq exp 5 Fully implicit VE 0 Fully implicit comp5 Fully implicit VE 0 Fully implicit comp5 Fully implicit VE 0 Fully implicit comp Fully implicit VE MRST

12 Black Oil with pressure decay 2006 calibration 2010 prediction 2008 calibration 2012 prediction Cavanagh, 2012

13 Sensitivity to CO 2 density and porosity (VESA) Top of Utsira Base Case Low Porosity Low CO 2 density Low porosity and 36%, 700 kg/m 3 27% 400 kg/m 3 CO 2 density

14 Uncertainty in Temperature (MUFITS)

15 Sensitivity to Temperature and Simulation Method (MRST) T = rhow = 1200 T = rhow = T = rhow = 1200 T = rhow = Disolved Struct. residual Residual Free residual Struct. movable Free movable Leaked Volume (M m 3 ) Years since simulation start

16 Sensitivity to top surface topography (VESA) (i) (ii) (iii) (iv) Minimize traps and spillpoits Minimize traps Minimize spillpoints Smoothed surface

17 Fully implicit comp Fully implicit VE Fully implicit comp Fully implicit VE Fully implicit comp Fully implicit VE Fully implicit comp Fully implicit VE Fully implicit comp Fully implicit VE Volume (M m 3 ) Fully implicit comp Fully implicit VE Disolved Struct. residual Residual Free residual Struct. movable Free movable Leaked 30 years of injection 220 years of post-injection trapping increases, while a good portion of CO2 flows out of the domain Years since simulation start

18 Summary and Conclusions See a good scientific focus from the participants Gravity-dominated flow well-suited for percolation and VE models Understanding uncertainty is important for calibrating the simulators Progress with the Sleipner reference model is good, despite relatively few results to date. Feedback to Statoil can help improve the utility of the dataset in the future.

19 Acknowledgements The Sleipner Benchmark was provided by Statoil and is managed by IEAGHG Geomechanics benchmark, please contact if interested.