Interpretation of baseline surface seismic data at the Violet Grove CO 2 injection site, Alberta

Similar documents
Penn West Pembina Cardium CO 2 EOR seismic monitoring program

Baseline VSP processing for the Violet Grove CO 2 Injection Site

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

Passive seismic monitoring at a CO2 injection site, Violet Grove, Alberta, Canada

Delineating a sandstone reservoir at Pikes Peak, Saskatchewan using 3C seismic data and well logs

Time lapse view of the Blackfoot AVO anomaly

Post-stack inversion of the Hussar low frequency seismic data

Reservoir Characterization of Plover Lake Heavy-Oil Field

Inversion and interpretation of multicomponent seismic data: Willesden Green, Alberta

Pre-stack (AVO) and post-stack inversion of the Hussar low frequency seismic data

Experimental comparison of repeatability metrics

Seismic applications in coalbed methane exploration and development

Time-lapse surface seismic monitoring of injected CO 2 at the Penn West CO 2 -EOR site, Violet Grove, Alberta

Shallow P and S velocity structure, Red Deer, Alberta

Multicomponent seismic surveys at Sibbald Flats, Alberta

Seismic methods in heavy-oil reservoir monitoring

An empirical method for estimation of anisotropic parameters in clastic rocks

Synthetic Seismogram A Tool to Calibrate PP & PS Seismic Data

Analysis of multicomponent walkaway vertical seismic profile data

Four-D seismic monitoring: Blackfoot reservoir feasibility

Seismic modelling and monitoring of carbon storage in a shallow sandstone formation

A Petroleum Geologist's Guide to Seismic Reflection

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

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

INVERSION OF SEISMIC DATA FOR ASSESSING FLUID REPLACEMENT IN THE NISKU FORMATION PATRICIA GAVOTTI DON C. LAWTON

A 3C-4D surface seismic and VSP program for a coalbed methane and CO 2 sequestration pilot, Red Deer, Alberta

Porosity prediction using attributes from 3C 3D seismic data

Time-lapse seismic modelling for Pikes Peak field

Pluto 1.5 2D ELASTIC MODEL FOR WAVEFIELD INVESTIGATIONS OF SUBSALT OBJECTIVES, DEEP WATER GULF OF MEXICO*

AFI (AVO Fluid Inversion)

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

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

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

Rock physics and AVO applications in gas hydrate exploration

Detecting fractures using time-lapse 3C-3D seismic data

Multicomponent seismic survey at Spring Coulee: a data repeatability study

Interpretation of PP and PS seismic data from the Mackenzie Delta, N.W.T.

UNIVERSITY OF CALGARY. Reservoir property prediction from well-logs, VSP and. Natalia Soubotcheva A THESIS

An overview of AVO and inversion

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

Feasibility and design study of a multicomponent seismic survey: Upper Assam Basin

Reservoir Characterization using AVO and Seismic Inversion Techniques

The reason why acoustic and shear impedances inverted

Constrained inversion of P-S seismic data

Practical aspects of AVO modeling

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

New Frontier Advanced Multiclient Data Offshore Uruguay. Advanced data interpretation to empower your decision making in the upcoming bid round

Using multicomponent seismic for reservoir characterization in Venezuela

Use of Seismic and EM Data for Exploration, Appraisal and Reservoir Characterization

Demystifying Tight-gas Reservoirs using Multi-scale Seismic Data

Tim Carr - West Virginia University

Observation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story

Case study: AVO analysis in a high-impedance Atoka Sandstone (Pennsylvanian), North Arkoma Basin, McIntosh County, Oklahoma

BPM37 Linking Basin Modeling with Seismic Attributes through Rock Physics

AVAZ and VVAZ practical analysis to estimate anisotropic properties

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

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

Part I. Near surface characterization using a high-resolution 3C seismic survey

Cross-well seismic modelling for coal seam delineation

Application of nonlinear time-lapse AVO to the Pouce Coupe data set. Presented by Shahin Jabbari Supervisor Kris Innanen

URTeC: Abstract

AVO Attributes of a Deep Coal Seam

3D Converted Wave Data Processing A case history

Azimuthal Velocity Analysis of 3D Seismic for Fractures: Altoment-Bluebell Field

Synthetic seismic modelling and imaging of an impact structure

Sensitivity of interval Vp/Vs analysis of seismic data

Multiattribute seismic analysis on AVO-derived parameters A case study

Microseismicity applications in hydraulic fracturing monitoring

Robert R. Stewart, Henry C. Bland, Qi Zhang*, and Felix Guthoff**

Seismic Inversion on 3D Data of Bassein Field, India

A 3D seismic survey for mapping shallow targets

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

Elements of 3D Seismology Second Edition

Imaging complex structure with crosswell seismic in Jianghan oil field

Grand Rapids Oil Sands 3D Seismic Incorporating and Comparing Multiple Data Types for Reservoir Characterization

Oil and Natural Gas Corporation Ltd., VRC(Panvel), WOB, ONGC, Mumbai. 1

Crosswell Seismic Imaging for Horizontal Drilling of the High Porosity Sand, Crystal Viking Pool, Alberta, Canada.

Density estimations using density-velocity relations and seismic inversion

Fr Reservoir Monitoring in Oil Sands Using a Permanent Cross-well System: Status and Results after 18 Months of Production

A new S-wave seismic source

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

Synthetic seismograms, Synthetic sonic logs, Synthetic Core. Larry Lines and Mahbub Alam

Reprinted From. Volume 03 Issue 03-JulY 2010

Rock-Physics and Seismic-Inversion Based Reservoir Characterization of AKOS FIELD, Coastal Swamp Depobelt, Niger Delta, Nigeria

Rock physics and AVO analysis for lithofacies and pore fluid prediction in a North Sea oil field

Dynamic GeoScience Martyn Millwood Hargrave Chief Executive OPTIMISE SUCCESS THROUGH SCIENCE

Heriot-Watt University

Bulletin of Earth Sciences of Thailand

UNIVERSITY OF CALGARY. Processing and Interpretation of Time-lapse Vertical Seismic. Profile Data from the Penn West CO 2 Monitoring Project

INTEGRATED GEOPHYSICAL INTERPRETATION METHODS FOR HYDROCARBON EXPLORATION

Envelope of Fracture Density

Integrated well log and 3-D seismic data interpretation for the Kakinada area of KG PG offshore basin

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

Porosity prediction using cokriging with multiple secondary datasets

FUNDAMENTALS OF SEISMIC EXPLORATION FOR HYDROCARBON

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

SEG/New Orleans 2006 Annual Meeting

Prediction of Shale Plugs between Wells in Heavy Oil Sands using Seismic Attributes

The elastic properties such as velocity, density, impedance,

Chałupki Dębniańskie Field: Improving Drilling Success in Shallow Gas Reservoirs with VectorSeis

Transcription:

Violet Grove seismic interpretation Interpretation of baseline surface seismic data at the Violet Grove CO 2 injection site, Alberta Fuju Chen and Don Lawton ABSTRACT Time-lapse seismic technology has been implemented in the Penn West pilot project, in west-central Alberta, Canada, to monitor the CO 2 storage in the Cardium Formation. A multi-component 3D surface seismic baseline survey was acquired in March 2003. An interpretation using both PP and PS seismic data has been made. Four horizons of Ardley, Cardium, Blackstone, and Viking Formations have been picked. At well 102/7-11-48-9, the top of Cardium Sand correlates to a weak peak at approximately 1043 ms and 1690 ms in the PP and PS survey, respectively; the top of the Blackstone Formation correlates to a strong trough at 1060 ms and 1714 ms in the PP and PS survey, respectively; There exists a small amplitude anomaly in the Cardium event around the observation well and the injection pad in the PP data. V p /V s values were calculated from the compression and shear sonic logs: The good reservoir of the Cardium Formation has a low V p /V s value (1.6-1.8); the shale above and below the Cardium Formation has a relatively high V p /V s value (1.8-2.0). The average V p /V s between Ardley and Viking horizons is approximately 2.0. INTRODUCTION The time-lapse seismic technique has been applied successfully to monitor the movement of fluid, pressure fronts, and water-oil contact during hydrocarbon production (e.g. Meyer, 2001; Gouveia, 2004). Some case studies on monitoring CO 2 storage in geological formations by time-lapse seismic also have been described in recent years (Davis, 2003; Raef, 2004). In Canada, an innovative CO 2 monitoring strategy has been implemented in the Penn West pilot in Alberta, which involved a multi-component 3D surface seismic program integrated with active and passive monitoring using geophones permanently cemented into an observation well (Lawton, 2005). The pilot site is located at the Pembina Oil Field, operated by Penn West Petroleum, in west-central Alberta, Canada, with CO 2 injection into the Upper Cretaceous Cardium Formation at a depth of approximately 1650 m below surface. The Cardium Formation is bounded by black shales of the Blackstone Formation below and the Wapiabi Formation above. It consists of two lithostratigraphic units: the lower Pembina River Member (an upwardly coarsening sequence, from shale through sandstone to conglomerate), and the upper Cardium Zone Member (a black shale and siltstone interval with pebbly stringers). The total thickness of the Cardium Formation at the site is approximately 20 m. The formation pressure and temperature are 19 MPa and 50 C, respectively. The produced oil in this area is approximately 40 API (Nielsen, 1984). There are 15 wells in the 3D survey area in total. Well 100/7-11-48-9 was chosen as the observation well. The baseline surface seismic data was shot in March 2005 using a 2 kg dynamite charge. It consists of two parallel, multi-component 2D lines, 400 m apart and oriented east west, and one orthogonal multi-component 2D north-south line of 3 km length, CREWES Research Report Volume 17 (2005) 1

Chen and Lawton respectively. Two additional short north-south receiver lines were also included (Figure 1). During acquisition, all lines were live with a receiver interval of 20 m, and a source interval of 40 m. Veritas DGC of Calgary processed the 2D lines as individual lines and also as a sparse 3D volume, which covers the injection pad. Pre-stack Kirchhoff migration datasets of vertical, radial and transverse components were delivered. Only the filtered/scaled noise attenuated PP and PS 3D data were interpreted here. FIG. 1. Map showing the Penn West CO 2 injection site. Multi-component seismic lines are shown in red, receivers-only lines are shown in green. The 3D survey is shown by a semi-transparent red rectangle. The observation (VSP) well location is shown by a yellow circle, and the CO 2 injection well locations are shown by blue circles. INTERPRETATION Synthetic seismogram and PP-PS correlation A dipole sonic log was available from well 102/7-11-48-9, which penetrated the top of the Blackstone Formation. The PP and PS wavelets extracted from the seismic data at the well location were used to generate the PP and PS synthetic seismograms. The top of the Cardium sand correlates to a weak peak at approximately 1043 ms and 1690 ms in the PP and PS survey, respectively; the top of Blackstone Formation correlates to a strong trough at 1060 ms and 1714 ms in the PP and PS survey, respectively (Figure 2). In order to correlate the PP and PS data well, the strong peak at 360 ms and 685 ms in PP and PS data, respectively, which correlates to a depth of approximately 440m at the well, was picked. This reflection is from the Ardley Coal Zone. Figure 3 shows the un-stretched, un-squeezed PP seismic and synthetic seismogram and the squeezed PS seismic and synthetic seismogram used to correlate the PP and PS datasets. 2 CREWES Research Report Volume 17 (2005)

Violet Grove seismic interpretation Compression sonic PP Seismic data Shear sonic PS Seismic data Ardley Ardley Cardium Blackstone Cardium Blackstone FIG. 2. PP (left) and PS (right) seismic correlation at well 102/7-11-48-9. The blue trace is synthetic seismogram, the red trace is the seismic trace extracted from the seismic data at the well location. In order to correlate the deep formations with the seismic events, the digital compression-sonic log of well 102/08-14-48-9W5/0, which penetrated the top of the Banff Formation, was projected into the 3D survey. The synthetic seismogram ties the seismic data very well, and the horizons of Cardium and Blackstone picks defined by the synthetic seismogram of this well match very well with the horizons identified from the observation well. The top of the Viking Formation correlates to a strong peak at 1232 ms in PP data. The top of the Viking Formation also correlates to a strong peak at approximately 1950 ms for PS data, which was defined by comparing the reflection character of the PP and PS data (Figure 4). PP and PS data interpretation Four horizons, which correspond to the Ardley Coal Zone, the tops of Cardium, Blackstone, and Viking Formations, were picked, and their time structural maps are shown (Figures 5, 6, 7, and 8). All the maps are very flat in PP and PS data. This means the formations are very flat in this area. On the Cardium and Blackstone structural maps, the central part of the 3D survey is slightly lower than areas near the edges of the survey. A small amplitude anomaly was shown around the observation well and the injection pad on the PP amplitude maps of the Cardium and Blackstone horizons (Figures 9 and 10). CREWES Research Report Volume 17 (2005) 3

Chen and Lawton Compression sonic PP Seismic data PP synthetic seismogram PS Seismic data PS synthetic seismogram Ardley Ardley Cardium Blackstone Cardium Blackstone FIG. 3. PP and PS synthetic seismogram at well 102/7-11-48-9. Gamma Ray Density P-wave sonic PP seismic Cardium Blackstone Viking FIG. 4. PP seismic correlation at well 102/08-14-48-9W5/0 and horizon picking. The blue trace is the synthetic seismogram, the red trace is the seismic trace extracted from the seismic data at the well location. 4 CREWES Research Report Volume 17 (2005)

Violet Grove seismic interpretation FIG. 5a. PP time structural map of the Ardley horizon. The observation well is shown by a solid circle, the injection wells are shown by open circles; the gaps are due to the absence of data. FIG. 5b. PS time structural map of the Ardley horizon. The observation well is shown by a solid circle, the injection wells are shown by open circles. CREWES Research Report Volume 17 (2005) 5

Chen and Lawton FIG. 6a. PP time structural map of the Cardium horizon. The observation well is shown by a solid circle, the injection wells are shown by open circles. FIG. 6b. PS time structural map of the Cardium horizon. The observation well is shown by a solid circle, the injection wells are shown by open circles. 6 CREWES Research Report Volume 17 (2005)

Violet Grove seismic interpretation FIG. 7a. PP time structural map of the Blackstone horizon. The observation well is shown by a solid circle, the injection wells are shown by open circles. FIG. 7b. PS time structural map of the Blackstone horizon. The observation well is shown by a solid circle, the injection wells are shown by open circles. CREWES Research Report Volume 17 (2005) 7

Chen and Lawton FIG. 8a. PP time structural map of the Viking horizon. The observation well is shown by a solid circle, the injection wells are shown by open circles. FIG. 8b. PS time structural map of the Viking horizon. The observation well is shown by a solid circle, the injection wells are shown by open circles. 8 CREWES Research Report Volume 17 (2005)

Violet Grove seismic interpretation FIG. 9a. Maximum amplitude map of the Cardium horizon in PP data. The observation well is shown by a solid circle, the injection wells are shown by open circles (A 10ms time gate centered at the Cardium horizon was used). FIG. 9b. Maximum amplitude map of the Cardium horizon in PS data. The observation well is shown by a solid circle, the injection wells are shown by open circles (A 10ms time gate centered at the Cardium horizon was used). CREWES Research Report Volume 17 (2005) 9

Chen and Lawton FIG. 10a. Maximum negative amplitude map of the Blackstone horizon in PP data. The observation well is shown by a solid circle, the injection wells are shown by open circles. (A 10ms time gate centered at the Blackstone horizon was used). FIG. 10b. Maximum negative amplitude map of the Blackstone horizon in PS data. The observation well is shown by a solid circle, the injection wells are shown by open circles. (A 10ms time gate centered at the Blackstone horizon was used). 10 CREWES Research Report Volume 17 (2005)

Violet Grove seismic interpretation Rock properties The P-wave and S-wave velocities of different lithologies were analyzed by using the sonic logs of well 102/7-11-48-9. The gamma ray values of clean sands range from 33 to 65 API, for P-wave velocities from 4400-4800 m/s and S-wave velocities from 2700-2800 m/s, respectively; the gamma ray values of shales range from approximately 100 to 130 API, with P-wave velocities of 3500-4200 m/s and S-wave velocities of 1700-2500 m/s, respectively. Figure 11 shows a plot of V p and V s versus the gamma ray value over the depth interval from 1500 to 1660 m. Clean sand Shale FIG. 11. Cross plot map of Gamma ray value and V p, V s over the interval from 1500 to 1660m for well 102/07-11-048-09W5/0. The V p /V s values were also calculated and the cross plot map of measured depth, gamma, and V p /V s were shown in Figure 12, and 13. Over the interval from 1500 m to 1665 m, the V p /V s value ranges from 1.61 to 2.11. The good reservoir of the Cardium Formation has a low V p /V s value (1.6-1.8); the shale above and below the Cardium Formation has a relatively high V p /V s value (1.8-2.0). The average V p /V s ratio map was computed from the PP and PS travel times of the interval between Ardley and Viking by using the following equation (Xu, 2004): V p 2* Tps / Vs = 1 T pp Figure 14 shows an almost constant V p /V s ratio of 1.95-2.0, that means the lithology is constant laterally in the 3D survey. This value is similar to that derived from the dipole sonic log (Figure 12). CREWES Research Report Volume 17 (2005) 11

Chen and Lawton Gamma Ray P-Wave Sonic S-Wave Sonic V p /V s Cardium Blackstone FIG. 12. Logs from well 102/7-11-48-9. Shale Clean sand FIG. 13. Cross plot map of gamma ray versus V p /V s over the interval from 1500 to 1660 m of well 102/7-11-48-9. 12 CREWES Research Report Volume 17 (2005)

Violet Grove seismic interpretation FIG. 14. Average V p /V s map between Ardley and Viking horizons. The gaps in the map are due to the absence of data. CONCLUSIONS Four horizons (Ardley, Cardium, Blackstone, and Viking) were picked. At the location of well 102/07-11-48-9W5/0, the top of Cardium sand correlates to a weak peak at approximately 1043 ms and 1690 ms in the PP and PS survey, respectively; the top of the Blackstone Formation correlates to a strong trough at 1060 ms and 1714 ms in the PP and PS survey, respectively; There exists a small amplitude anomaly in the Cardium reservoir around the observation well and the injection pad in the PP data. V p /V s values were calculated from the compression and shear sonic logs: The good reservoir of the Cardium Formation has a low V p /V s value (1.6-1.8); the pure shale above and below the Cardium Formation has a relatively high V p /V s value (1.8-2.0). The average V p /V s between Ardley and Viking horizons is approximately 2.0. FUTURE WORK In the near future, some CO 2 detection methods such as acoustic impedance inversion, AVO etc. will be tested on the baseline data in preparation for analysis of monitor survey data. In early 2006, the first monitor survey will be shot and processed. After that, the monitor survey will be interpreted and compared with the baseline survey. CREWES Research Report Volume 17 (2005) 13

Chen and Lawton ACKNOWLEDGEMENTS This program was funded through grants from the Alberta Energy Research Institute (AERI), Western Economic Diversification (WED) and the Consortium for Research in Elastic Wave Exploration Seismology (CREWES) at University of Calgary. Hampson- Russell Software Services and Landmark Graphics Corporation donated the soft-wares used in the interpretation. We thank Penn West Petroleum for providing the observation well and logistic support. We also thank Kevin Hall, Han-Xing Lu, Henry Bland, and Marcia Coueslan for their technical support and data processing. REFERENCES Lawton, Don, et al., 2005, Seismic survey design for monitoring CO 2 storage: integrated multi-component surface and borehole seismic surveys: Penn West Pilot, Alberta, Canada, Abstract, GHGT-8 Nielsen, A and Porter, J., 1984, Pembina oil field in retrospect: CSPG Memoir 9, P. 1-13 Xu, Chuandong and Stewart, R, 2004, Using 3D multi-component seismic data, logs and VSP to interpret a sand-shale oil reservoir: EAGE 66 th Conference & Exhibition Paris, France, 7-10 June 2004 14 CREWES Research Report Volume 17 (2005)