Heriot-Watt University

Similar documents
Calibration of the petro-elastic model (PEM) for 4D seismic studies in multi-mineral rocks Amini, Hamed; Alvarez, Erick Raciel

Th LHR2 08 Towards an Effective Petroelastic Model for Simulator to Seismic Studies

We G Quantification of Residual Oil Saturation Using 4D Seismic Data

Th LHR2 04 Quantification of Reservoir Pressure-sensitivity Using Multiple Monitor 4D Seismic Data

Integrating reservoir flow simulation with time-lapse seismic inversion in a heavy oil case study

Uncertainties in rock pore compressibility and effects on time lapse seismic modeling An application to Norne field

Time-lapse seismic modelling for Pikes Peak field

The elastic properties such as velocity, density, impedance,

Downloaded 09/16/16 to Redistribution subject to SEG license or copyright; see Terms of Use at

A E. SEG/San Antonio 2007 Annual Meeting. exp. a V. a V. Summary

Integration of Rock Physics Models in a Geostatistical Seismic Inversion for Reservoir Rock Properties

ROCK PHYSICS DIAGNOSTICS OF NORTH SEA SANDS: LINK BETWEEN MICROSTRUCTURE AND SEISMIC PROPERTIES ABSTRACT

HampsonRussell. A comprehensive suite of reservoir characterization tools. cgg.com/geosoftware

4D stress sensitivity of dry rock frame moduli: constraints from geomechanical integration

AFI (AVO Fluid Inversion)

Kondal Reddy*, Kausik Saikia, Susanta Mishra, Challapalli Rao, Vivek Shankar and Arvind Kumar

Interpretation and Reservoir Properties Estimation Using Dual-Sensor Streamer Seismic Without the Use of Well

Quantifying remaining oil saturation using seismic time-lapse amplitude changes at fluid contacts Alvarez, Erick; MacBeth, Colin; Brain, Jon

Feasibility study of time-lapse seismic monitoring of CO 2 sequestration

PRM on Johan Sverdrup - an unique Opportunity. Force seminar 2017 Stavanger, Maximilian Schuberth

P- and S-Wave Velocity Measurements and Pressure Sensitivity Analysis of AVA Response

PRELIMINARY REPORT. Evaluations of to what extent CO 2 accumulations in the Utsira formations are possible to quantify by seismic by August 1999.

SPE These in turn can be used to estimate mechanical properties.

Downloaded 11/02/16 to Redistribution subject to SEG license or copyright; see Terms of Use at Summary.

Integrating rock physics and full elastic modeling for reservoir characterization Mosab Nasser and John B. Sinton*, Maersk Oil Houston Inc.

Bertrand Six, Olivier Colnard, Jean-Philippe Coulon and Yasmine Aziez CGGVeritas Frédéric Cailly, Total

Summary. Seismic Field Example

Velocity Measurements of Pore Fluids at Pressure and Temperature: Application to bitumen

ScienceDirect. Model-based assessment of seismic monitoring of CO 2 in a CCS project in Alberta, Canada, including a poroelastic approach

The Weyburn Field in southeastern Saskatchewan,

We LHR3 06 Detecting Production Effects and By-passed Pay from 3D Seismic Data Using a Facies Based Bayesian Seismic Inversion

A021 Petrophysical Seismic Inversion for Porosity and 4D Calibration on the Troll Field

Downloaded 09/09/15 to Redistribution subject to SEG license or copyright; see Terms of Use at

The use of seismic attenuation to indicate saturation in hydrocarbon reservoirs: Theoretical study and modelling approach

Time-lapse (4D) seismic effects: Reservoir sensitivity to stress and water saturation

Time-Lapse Seismic: A Geophysicist s Perspective. Rob Staples Shell E & P Europe. SPE / EAGE workshop Copenhagen, March 2004

Simultaneous Inversion of Clastic Zubair Reservoir: Case Study from Sabiriyah Field, North Kuwait

Reservoir properties inversion from AVO attributes

A New AVO Attribute for Hydrocarbon Prediction and Application to the Marmousi II Dataset*

Pre-Stack Seismic Inversion and Amplitude Versus Angle Modeling Reduces the Risk in Hydrocarbon Prospect Evaluation

P235 Modelling Anisotropy for Improved Velocities, Synthetics and Well Ties

Edinburgh Anisotropy Project, British Geological Survey, Murchison House, West Mains

RESERVOIR SEISMIC CHARACTERISATION OF THIN SANDS IN WEST SYBERIA

Quantitative Interpretation

11282 Rock Physics Analysis and Time-lapse Rock Imaging of Geochemical Effects Due to CO2 Injection into Reservoir Rocks

The role of seismic modeling in Reservoir characterization: A case study from Crestal part of South Mumbai High field

THE USE OF SEISMIC ATTRIBUTES AND SPECTRAL DECOMPOSITION TO SUPPORT THE DRILLING PLAN OF THE URACOA-BOMBAL FIELDS

Downloaded 09/16/16 to Redistribution subject to SEG license or copyright; see Terms of Use at

We Simultaneous Joint Inversion of Electromagnetic and Seismic Full-waveform Data - A Sensitivity Analysis to Biot Parameter

We apply a rock physics analysis to well log data from the North-East Gulf of Mexico

23855 Rock Physics Constraints on Seismic Inversion

Four-D seismic monitoring: Blackfoot reservoir feasibility

Time-lapse seismic monitoring and inversion in a heavy oilfield. By: Naimeh Riazi PhD Student, Geophysics

Rock Physics Modeling in Montney Tight Gas Play

F003 Geomodel Update Using 4-D Petrophysical Seismic Inversion on the Troll West Field

Geological Classification of Seismic-Inversion Data in the Doba Basin of Chad*

C002 Petrophysical Seismic Inversion over an Offshore Carbonate Field

Shaly Sand Rock Physics Analysis and Seismic Inversion Implication

Toward an Integrated and Realistic Interpretation of Continuous 4D Seismic Data for a CO 2 EOR and Sequestration Project

A look into Gassmann s Equation

RP 2.6. SEG/Houston 2005 Annual Meeting 1521

Some consideration about fluid substitution without shear wave velocity Fuyong Yan*, De-Hua Han, Rock Physics Lab, University of Houston

Probabilistic seismic inversion using pseudo-wells

RC 2.7. Main Menu. SEG/Houston 2005 Annual Meeting 1355

Cold production footprints of heavy oil on time-lapse seismology: Lloydminster field, Alberta

IDENTIFYING PATCHY SATURATION FROM WELL LOGS Short Note. / K s. + K f., G Dry. = G / ρ, (2)

Time-lapse seismic modeling of CO 2 sequestration at Quest CCS project

Lithology prediction and fluid discrimination in Block A6 offshore Myanmar

Rock Physics & Formation Evaluation. Special Topic

Temperature Dependence of Acoustic Velocities in Gas-Saturated Sandstones

The Hangingstone steam-assisted gravity drainage

Linearized AVO and Poroelasticity for HRS9. Brian Russell, Dan Hampson and David Gray 2011

Towards Interactive QI Workflows Laurie Weston Bellman*

Rock Physics Perturbational Modeling: Carbonate case study, an intracratonic basin Northwest/Saharan Africa

Evaluation of Rock Properties from Logs Affected by Deep Invasion A Case Study

The Interaction of Reservoir Engineering and Geomechanics (a story)

SEISMIC PROFILE CGG SEISMIC INVERSIONS. by Lucia Levato, CGG. three offshore case studies show how one size does not fi t all. 18 seismic profile

QUANTITATIVE INTERPRETATION

Net-to-gross from Seismic P and S Impedances: Estimation and Uncertainty Analysis using Bayesian Statistics

Time lapse view of the Blackfoot AVO anomaly

Elements of 3D Seismology Second Edition

Hydrogeophysics - Seismics

Reservoir simulation and feasibility study for seismic monitoring at CaMI.FRS, Newell County, Alberta

AVO attribute inversion for finely layered reservoirs

Dynamic GeoScience Martyn Millwood Hargrave Chief Executive OPTIMISE SUCCESS THROUGH SCIENCE

Shear Wave Velocity Estimation Utilizing Wireline Logs for a Carbonate Reservoir, South-West Iran

SEG Houston 2009 International Exposition and Annual Meeting. that the project results can correctly interpreted.

Seismic reservoir characterization in offshore Nile Delta.

Adding Value with Broadband Seismic and Inversion in the Central North Sea Seagull Area

ROCK PHYSICS MODELING FOR LITHOLOGY PREDICTION USING HERTZ- MINDLIN THEORY

Fluid-property discrimination with AVO: A Biot-Gassmann perspective

Rock Physics and Quantitative Wavelet Estimation. for Seismic Interpretation: Tertiary North Sea. R.W.Simm 1, S.Xu 2 and R.E.

Bandlimited impedance inversion: using well logs to fill low frequency information in a non-homogenous model

Sensitivity Analysis of Pre stack Seismic Inversion on Facies Classification using Statistical Rock Physics

BPM37 Linking Basin Modeling with Seismic Attributes through Rock Physics

Reducing Uncertainty through Multi-Measurement Integration: from Regional to Reservoir scale

An Overview of the Tapia Canyon Field Static Geocellular Model and Simulation Study

Summary. Introduction

An empirical method for estimation of anisotropic parameters in clastic rocks

We Density/Porosity Versus Velocity of Overconsolidated Sands Derived from Experimental Compaction SUMMARY

Transcription:

Heriot-Watt University Heriot-Watt University Research Gateway 4D seismic feasibility study for enhanced oil recovery (EOR) with CO2 injection in a mature North Sea field Amini, Hamed; Alvarez, Erick Raciel; Wilkinson, David; Lorsong, J.; Slater, J. ; Holman, G. ; Timofeeva, O. Publication date: 2014 Document Version Publisher's PDF, also known as Version of record Link to publication in Heriot-Watt University Research Gateway Citation for published version (APA): Amini, H., Alvarez, E. R., Wilkinson, D., Lorsong, J., Slater, J., Holman, G., & Timofeeva, O. (2014). 4D seismic feasibility study for enhanced oil recovery (EOR) with CO2 injection in a mature North Sea field. Poster session presented at 4th EAGE CO2 Geological Storage Workshop, Stavanger, Norway.

P33 4D Seismic Feasibility Study for Enhanced Oil Recovery (EOR) with CO2 Injection in a Mature North Sea Field H. Amini (Senergy Ltd), E. Alvarez (Senergy Ltd), D. Wilkinson (Senergy Ltd), J. Lorsong (2Co Energy Limited), J. Slater (Talisman Sinopec Energy UK Limited), G. Holman (Talisman Sinopec Energy UK Limited) & O. Timofeeva* (Senergy Ltd) SUMMARY A mature North Sea field with a long history of water flooding has been modelled to predict the potential for enhanced oil recovery (EOR) with CO2 injection. This technique has the potential to further extend the life of the field by extracting more oil and to store the injected CO2 in the reservoir. A feasibility study was conducted to assess whether 4D (time-lapse) seismic could be used as a monitoring tool during CO2 injection. Different approaches were employed in this project; the initial examination using a 4D risk assessment table (Lumley et al., 1997) suggests that 4D seismic can be applied. For a more quantitative analysis, a petro-elastic model was designed to capture both the saturation and pressure signatures from the injected CO2. 1D modelling using fluid substitution at well locations, and simulator to seismic modelling in conjunction with the compositional simulation model was performed to predict the threedimensional 4D signatures within the simulation model grid. Based on this analysis a seismic impedance change of up to 9% is predicted which is above the limit that is typically considered for 4D applications. The discrimination between pressure and saturation signals is also discussed and the possibility to use 4D seismic to detect pressure changes was evaluated. This study makes suggestions regarding the seismic survey characteristics and appropriate timing to acquire the monitor surveys.

Introduction After 30 years of production in a mature North Sea field (Field-X), where water flooding has been the major drive mechanism, CO 2 injection has been considered as an enhanced oil recovery and storage strategy. Senergy has performed a feasibility study, covering reservoir engineering and geophysical aspects of this project. The reservoir engineering side of the project covers the optimum scenario regarding the number of required wells (producers and injectors), well placement and completion intervals to maximise utilisation of injected gas and minimise inefficient recycling. The use of the existing wells was also investigated. The proposed development plan consists of a continuation of the current water flooding for 5 years, followed by CO 2 enhanced oil recovery. As part of the reservoir engineering studies, a compositional simulation model (Eclipse 300) was developed. This allows the reservoir performance under gas+co 2 EOR for the next 30 years for different scenarios of well placements to be evaluated. The geophysical studies define the extent of the storage complex, identify the potential spill points from structural closures, characterise the subsequent migration of CO 2 and potential secondary/tertiary seals. A 4D seismic feasibility study was also completed. This 4D seismic feasibility study is the subject of this paper and included three stages: 1) preliminary 4D screening and petro-elastic modelling, 2) well-based modelling, and 3) 3D simulator to seismic modelling. The details of these studies are outlined in the following sections. 4D risk assessment (Lumley et al., 1997) The 4D technical risk spreadsheet (Lumley et al., 1997) was examined for the field. Comparatively, based on the rock and fluid characteristics, Field-X falls (Score 16) into the range of a typical North Sea reservoir (score 15), with a good possibility for monitoring gas+co2 injection by the use of 4D seismic. scores Property Field-X 5 4 3 2 1 0 Dry rock bulk modulus ~13 GPa <3 3-5 5-10 10-20 20-30 30+ Fluid compressibility contrast Up to 95% change 250+ 150-200 100-150 50-100 25-50 0-25 Fluid saturation change Water replaced by gas 50+ 40-50 30-40 20-30 10-20 0-10 Porosity ~23 % 35+ 25-35 15-25 10-15 5-10 0-5 Impedance change Up to 9% change 12+ 8-12 4-8 2-4 1-2 0 Preliminary 4D screening Petro-elastic model (PEM) parameters were calibrated using an optimisation algorithm (Amini & Alvarez, 2014) based on well-log data from one well. Gassmann s (1951) fluid substitution model and MacBeth s (2004) model for stress-sensitivity was used to take pore-pressure effects into account. The magnitude of the impedance change due to gas+co 2 injection is predicted to be up to 9% (Figure 1). Well-based modelling The pressure and saturation information were extracted from the compositional simulation model at the well location at different time steps. Using the calibrated PEM parameters, the elastic parameters were calculated at each time step and the 4D AVO seismic response was modelled. The 4D AVO analysis shows that to separate pressure and saturation responses based on angle stacks is challenging (Figure 2). 3D simulator to seismic modelling Even though well-based 4D seismic modelling is useful in establishing the feasibility and evaluation of 4D AVO responses, they might be inadequate in understanding the lateral-sweep efficiency and the timing of the repeat surveys (Johnston, 2013). Additionally, the well data may not represent the variability of reservoir properties (wells are targeted to high net-to-gross areas) and hence bias the modelling results. Simulator to seismic modelling is a valuable tool in overcoming these challenges. Static and dynamic properties from the simulation model at the desired time-steps are extracted, and using the calibrated PEM parameters the elastic parameters were calculated at each time-step

(Figures 3-4). Simulator to seismic modelling was performed using a Senergy in-house software programme. It is observed that the saturation changes due to injection results in impedance changes above the threshold of what is considered a detectable 4D signal. This is due to the fact that most of the reservoir is water flooded after a long period of water injection, and the contrast between water and a gas+co 2 mixture is large enough to generate an evident 4D signal. 3D modelling indicates that both pore pressure (pressure build-up) and saturation changes (mainly water being replaced by gas+co 2 ) result in a decrease in impedance (softening). However, due to the choice of a relatively low stress sensitivity curve, the pressure signal due to CO 2 injection will be too small to be detected. Conclusions The 4D risk spreadsheet shows that Field-X is comparable to any typical North Sea field in terms of the rock and fluid characteristics needed to have a measurable 4D response. The results of the well-based modelling and simulator to seismic modelling show that 4D seismic can be used effectively to monitor the gas+co 2 injection, even at low seismic resolution (dominant frequency 22 Hz). The expected magnitude of the 4D seismic response is estimated to be about 9%, which is above the threshold of what is considered a detectable 4D signal. Our modelling also indicates that the pressure signal will be too small to be detected, and it will be coupled with the effect of an increase in gas saturation. Our analysis can be used to plan future 4D seismic monitor surveys. The largest uncertainty in our study is related to the dry rock characterisation which could not be calibrated due to a lack of shear sonic data and laboratory measurements on core data; this will mainly affect our estimation of the pressure signal. Acknowledgments We thank 2Co Energy Limited, Talisman Sinopec Energy UK Limited and Senergy (GB) Limited for permission to publish this study. We would also like to thank John McCurry (Senergy (GB) Limited) for technical support and useful discussions. References Amini, H., Alvarez, E., 2014, Calibration of the petro-elastic model (PEM) for 4D seismic studies in multi-mineral rocks, submitted to EAGE/FESM Joint Regional Conference Petrophysics Meets Geoscience, Kuala Lumpur, Malaysia. Gassmann, F., 1951, Uber die elastizitat poroser medien: Vierteljahrsschrift der Natur. Gesselschaft, 96, 1-23. (English translation from http://sepwww.stanford.edu/sep/berryman/ps/gassmann.pdf). Johnston, D. H., 2013, practical applications of time-lapse seismic data, 2013 DISC, SEG. Lumley, D. E., Behrens, R. A., and Wang, Z., 1997, Assessing the technical risk of a 4-D seismic project: The Leading Edge, 16, 1287 1291. MacBeth, C., 2004, A classification for the pressure-sensitivity properties of a sandstone rock frame: Geophysics, 69. a) b) Paleocene (West of Shetland (WoS)) Forties (Nelson) Rotliegend (Southern North Sea) P-Impedance (Km.cc)/(gr.s) Oil-Water Change in Kdry (%) Applied to Field-X Gas-Water Pressure build-up 100 - S W Change in Pore Pressure (Mpa) Figure 1. a) The curves represent a theoretical trend showing only the saturation effect (no pressure effect) for a rock with ϕ=23%, b) published rock stress sensitivity curves (MacBeth, 2004) for different samples of sandstones with a similar porosity as Field-X. In the absence of rock-stress sensitivity measurements, the curve with the lowest stress-sensitivity (West of Shetland) was selected. Since the reservoir in Field-X is deeper than the West of Shetland sample, the rock stress sensitivity is expected to be lower.

DEPTH (FT) a) Start of injection After 5 years After 10 years After 15 years After 20 years After 25 years After 30 years 4D gather 4D gather 4D gather 4D gather 4D gather 4D gather shale sand oil water gas b) c) 4D change (monitor-baseline) 4D Amplitude Angle of incidence 5 10 15 20 25 30 Figure 2. Pseudo-logs at well were extracted from the simulation model at 7 time steps over 30 years of gas+co 2 injection. Calibrated petro-elastic model parameters were used to calculate V P and density. A Greenberg-Castagna relationship was used to estimate V S. Ricker wavelet (22 Hz) was used to calculate synthetic seismograms. a) the AVO gathers at 7 time steps; 4D signal is a function of both gas thickness and gas saturation. b) the AVO response at different time steps, c) extracted amplitude variations at the well location. Change in pore pressure Change in gas saturation Figure 3. The vertical sections show the change in gas saturation and pore pressure (from the compositional simulation model) and the corresponding expected change in impedance after 15 years of gas+co 2 injection (from simulator to seismic modelling). Change in P-impedance

a) Pore pressure change after 2 years Gas saturation change after 2 years P-Impedance change after 2 years P change (psi) S g change (v/v) IP change (%) b) Pore pressure change after 15 years Gas saturation change after 15 years P-Impedance change after 15 years P change (psi) S g change (v/v) IP change (%) c) Pore pressure change after 30 years Gas saturation change after 30 years P-Impedance change after 30 years P change (psi) S g change (v/v) IP change (%) Figure 4. The maps show the change in gas saturation and pore pressure (from the compositional simulation model) and the corresponding expected change in impedance at different time steps (from simulator to seismic modelling).