11.07.2014 Views

edinburgh anisotropy project applied seismic anisotropy 7 phase ...

edinburgh anisotropy project applied seismic anisotropy 7 phase ...

edinburgh anisotropy project applied seismic anisotropy 7 phase ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

EDINBURGH ANISOTROPY PROJECT<br />

APPLIED SEISMIC ANISOTROPY<br />

TH<br />

7 PHASE RESEARCH PROGRAMME<br />

2006-2009


THE EDINBURGH ANISOTROPY PROJECT<br />

APPLIED SEISMIC ANISOTROPY<br />

7 TH PHASE RESEARCH PROGRAMME 2006-2009<br />

EXCUTIVE SUMMARY<br />

The Edinburgh Anisotropy Project (EAP) conducts research into <strong>applied</strong> <strong>seismic</strong> <strong>anisotropy</strong>, and how<br />

to utilize such information for the understanding of hydrocarbon reservoirs. EAP is supported by a<br />

consortium of major oil, service and software companies. The EAP research programme is guided by<br />

industry requirements and is primarily focused on fracture characterization and <strong>anisotropy</strong> imaging.<br />

The programme covers a wide variety of topics including data processing, interpretation, modelling<br />

and acquisition design, as well as fundamental wave propagation issues.<br />

EAP was established in July 1988 with seventeen years of experience in <strong>seismic</strong> <strong>anisotropy</strong>. EAP has<br />

organised its research programme in three-year <strong>phase</strong>s. The <strong>project</strong> is now in its 6 th Phase, which ends<br />

in June 2006.<br />

The 7 th Phase Research Programme (2006-09)<br />

The programme proposed for the 7 th Phase is a continuation of the 6 th Phase, but with a particular<br />

emphasis on frequency-dependent <strong>anisotropy</strong>, arising from discussions with sponsors. The<br />

programme has three components:<br />

o AVD technology for fracture characterization<br />

Extension of traditional AVD (attributes versus direction) methods to exploit frequencydependent<br />

<strong>anisotropy</strong> for robust fracture density information together with fluid saturation and<br />

dominant scale length<br />

o Converted-wave analysis for improved imaging of reservoir structure<br />

Extension of current converted-wave imaging tools (e.g. CXtools) to account for diodic<br />

converted-wave effects and further development of depth-imaging tools for building a common<br />

anisotropic earth model<br />

o Anisotropic rock physics for analysing frequency-dependent <strong>anisotropy</strong><br />

Progress towards gaining more insights into the information content of frequency-dependent<br />

<strong>anisotropy</strong> through experimental studies, calibrated theoretical developments and analysis of<br />

time-lapse VSP data<br />

An important aspect of EAP research is that the development of theory, methodology and algorithms<br />

is combined with application to real data, supplied by the sponsors, to demonstrate the value of the<br />

technology produced.


Deliverables and Technology Transfer<br />

o Research results in the form of annual reports, technical papers and software releases<br />

o Courses on <strong>applied</strong> <strong>seismic</strong> <strong>anisotropy</strong>, tailored to individual sponsor’s requirements<br />

o Special consultancy and joint studies of company-specific <strong>anisotropy</strong>-related issues<br />

o PhD research studentships<br />

EAP believes in maintaining close links with each member of the consortium so that the technology<br />

developed is transferred efficiently, and the research programme remains dynamic and evolutionary.<br />

Fees<br />

For the 7 th Phase the annual fees will be set at a flat rate of £31k.<br />

For new sponsors there is, in addition to the annual subscription, a joining fee of £30k, or £10k plus<br />

some negotiable equivalent in the form of software or data.<br />

Meetings<br />

The main annual EAP Technical and Business Meeting is held in Edinburgh, usually immediately<br />

following the annual EAGE conference. A further meeting is held at the annual SEG conference.<br />

Consortium members<br />

EAP is currently sponsored by eighteen service and oil companies and one software company: Agip,<br />

BG, BGP, BP, Chevron, CNPC, ConocoPhillips, ExxonMobil, GX Technology, KerrMcGee,<br />

Marathon, Norsk Hydro, PDVSA, Schlumberger, Shell, SinoPec, Total, Veritas DGC and Landmark.<br />

Project Resources<br />

Staff: A team of 5 scientific staff and 6 PhD students (see Appendix A).<br />

Computing: a 24-node PC cluster, Sun Ultrasparc/Sparc Workstations, Intel-Based destop PCs,<br />

Terabyte Raid System For Data Storage.<br />

Datasets: EAP benefits from access to datasets supplied by the sponsors. A list of the datasets<br />

available to the <strong>project</strong> during the 6 th Phase is given in Appendix B.<br />

Further Information<br />

Enquiries about EAP activities and membership details should be addressed to:<br />

Professor Xiang-Yang Li, Technical Director<br />

Edinburgh Anisotropy Project, British Geological Survey,<br />

Murchison House, West Mains Road, Edinburgh EH9 3LA, UK<br />

Tel: +44 (0)131 667 1000/Fax: +44 (0)131 667 1877/E-mail: xyl@bgs.ac.uk<br />

Further information can also be found at http://www.eap.bgs.ac.uk


7 th PHASE RESEARCH PROGRAMME<br />

Introduction<br />

During the 6 th Phase (2003-2006), we have focused on consolidating our ideas and findings for<br />

technology transfer. This led to the release of three software packages, i.e. CXtools for convertedwave<br />

imaging, RU for rock physics and fracture modelling, and GeoLab for depth model building and<br />

imaging. Major progress has been made in the following areas in the past five years:<br />

o Converted-waves analysis<br />

- New theories and algorithms for converted-wave imaging are consolidated into CXtools,<br />

which makes it possible to perform anisotropic imaging with high accuracy but low cost<br />

- A number of successful case studies (e.g. Lomond 2D 4C data, Afflect 3D 4C data)<br />

- Ideas for decoupled migration to perform multiwave type prestack depth imaging<br />

o Fracture characterization<br />

- Frequency-dependent P- and S-wave attributes to characterize fluid saturation and dominant<br />

scale length<br />

- A number of successful case studies (e.g. Clair and Emilio 3D 4C data )<br />

o Time-lapse <strong>seismic</strong> data analysis<br />

- Improved dynamic rock physical model with multi-scale fluid flow<br />

- Ideas of frequency-dependent time-lapse <strong>anisotropy</strong> for discriminating the effects of pore<br />

pressure and saturation<br />

These developments have been endorsed and well received by all sponsors. During this Phase, EAP<br />

has also attracted a record number of sponsors, which encourages us to continue developing our work<br />

through real data applications during the next <strong>phase</strong> from 2006 to 2009, with a primary focus on<br />

integrated studies and the delivery of user-friendly tools.<br />

EAP activities seek to address the needs of industry (fractures, converted waves, imaging in<br />

complex areas), while retaining a component of basic strategic science for future developments<br />

(frequency-dependent <strong>anisotropy</strong>, fluid flow). It is EAP’s intention to achieve a healthy balance<br />

between <strong>applied</strong> and basic research as well as ‘blue-sky’ research.<br />

In addition to the industrial sponsorship, every year EAP also receives matching funds from the<br />

Natural Environment Research Council through its science budget. This provides great value for<br />

money for EAP’s members. From time to time EAP also receives government grants through<br />

competitive tender and review exercises. These matching funds and research grants enable us to carry<br />

out more basic and curiosity-driven research. As a whole, EAP research builds on a solid foundation<br />

of theoretical expertise and field-data analysis in <strong>seismic</strong> <strong>anisotropy</strong> developed over the past<br />

seventeen years. The following sections describe the major components of the 7 th Phase programme.


Converted wave analysis<br />

Interest in converted-waves was rekindled by the emergence of 4C OBC technology in the mid<br />

1990’s. Since then, more than a decade has elapsed. Despite an intensive effort in research and<br />

development of the converted-wave technology, the anticipated breakthrough into a mainstream<br />

technology has still to be realized. A general consensus is that the processing and interpretation of<br />

converted-wave data appears to lag behind the corresponding acquisition technology that has evolved<br />

from analogue phones to digital sensors, from line receivers to point receivers, from dragged cables to<br />

permanently installed trench-buried systems. The difficulties in processing the converted waves lie in<br />

two main areas. In onshore applications, the near-surface is still the main obstacle for obtaining highquality<br />

and high resolution converted-wave data; whilst offshore, converted-wave imaging suffers<br />

from the wide occurrence of <strong>anisotropy</strong> in marine sediments together with the asymmetric raypaths<br />

and the diodic velocity effects due to shallow inhomogeneities, which has even made the data<br />

acquired with permanently installed systems un-interpretable. Due to these problems, there is a<br />

general lack of user-friendly and reliable tools for processing and interpreting the converted-waves.<br />

Since the onset of Phase 5 (2000-03), we have made significant progress in bridging the above gap for<br />

converted-wave imaging. During Phase 5 (2003-03), we have developed a common theoretical<br />

framework to describe the converted-wave behaviour in anisotropic media in terms of convertedwavefield<br />

attributes. This paved the way for the use of standard procedures such as semblance<br />

analysis and Kirchhoff summation to implement converted-wave anisotropic processing with high<br />

accuracy but low cost. During Phase 6 (2003-06), we extended the previous developments into 3D-<br />

4C data for both anisotropic imaging and parameter estimation, and consolidated our achievement<br />

into an user-friendly software, referred to as CXtools. Using this tool, we have reprocessed a number<br />

of 2D and 3D-4C dataset from the North Sea , and in each case, a substantially improved subsurface<br />

image was obtained that has made significant impact on the development of the relevant oilfield.<br />

During Phase 7, our emphasis is to extend CXtools to account for the diodic converted-wave effects<br />

and further development of our depth-imaging tools for building a common anisotropic earth model,<br />

and our aim is to fill the wide acknowledged gap in converted-wave processing with CXtools.<br />

Areas of interest and deliverables for the 7 th Phase:<br />

o diodic decomposition for 3D 4C data;<br />

o the corresponding depth imaging work flows using PP-independent PS imaging approach;<br />

o decoupled migration and associated work flows for anisotropic inversion;<br />

o conversion point signature from 3D 4C data;<br />

o anisotropic estimation in TTI media;<br />

o case studies of 3D-4C data from permanently installed systems in the North Sea;<br />

o user-friendly fully-documented algorithms/software for PC cluster.<br />

Datasets available for the 7 th Phase:<br />

o Valhall LoFS (Life of Field Seismic) 3D-4C data, North Sea;<br />

o BG 4C data<br />

o Shenli 3D-3C data, Yellow River Delta


TH<br />

7 PHASE RESEARCH PROGRAMME<br />

CONVERTED-WAVE ANALYSIS<br />

2006-09<br />

Research developments<br />

-Diodic PS decomposition for building a<br />

common earth model from 3D-4C data<br />

- Anisotropic estimation in orthorhombic<br />

or TTI media<br />

- Optimizing anisotropic depth model<br />

through a prestack time migration<br />

- Case studies of 3D 4C <strong>seismic</strong> data<br />

- User-friendly algorithms and software<br />

CXtools for anisotropic model building<br />

Affleck: EAP Line 60, through Tor penetration (PS data in PP time)<br />

(a) EAP processing<br />

Tor<br />

Top Tor Tie<br />

(b) isotropic/conventional<br />

The Afflect 3D-4C experience (courtesy of KerrMcGee): (a) Reprocessed by<br />

EAP using the Cxtools; and (b) conventional isotropic processing.


AVD technology for fracture characterization<br />

EAP has pioneered the use of <strong>seismic</strong> <strong>anisotropy</strong> for characterization of natural fractured reservoirs.<br />

Currently, we are extending our expertise to go beyond what the conventional <strong>anisotropy</strong> concept for<br />

fracture imaging can offer by developing innovative methods to evaluate the dominant scale length<br />

from 3D VSP, 3D surface and 3D 4C-OBC <strong>seismic</strong> data. Our new theories and methods allow the<br />

<strong>seismic</strong> response across a broad band of <strong>seismic</strong> frequencies to be computed, making it possible to<br />

observe and analyse the directional spectral/dispersion/attenuation variation and frequency-dependent<br />

<strong>anisotropy</strong> in <strong>seismic</strong> data and then to invert these variations for the properties of the fracture<br />

distribution. This will require advanced techniques for extracting spectral/dispersion and attenuation<br />

characteristics from <strong>seismic</strong> data using, for example, multi-instantaneous attribute analysis and<br />

modern spectral decomposition techniques.<br />

EAP is also developing advanced numerical methods to simulate wave propagation in media with<br />

discrete fracture systems using the discrete dynamic rock physics model, which allows us to develop<br />

inversion techniques to obtain information about discrete formation or meso-scale fractures. The work<br />

so far focused on 2D will be extended to 3D with an aim to develop a discrete fracture network (DFN)<br />

model by using <strong>seismic</strong> <strong>anisotropy</strong> attributes as inputs to combine with other geological/ geophysical<br />

data.<br />

Areas of interest and deliverables for the 7 th Phase:<br />

o Characteristics of directional dispersion/spectral/attenuation variations in P- and converted-wave<br />

data and their applications for determining fracture properties such as the dominant scale length,<br />

and fluid saturation;<br />

o Generation of realistic fracture networks and patterns at the meso-scale for <strong>seismic</strong> modelling, and<br />

computation of 3D synthetic VSP and reflection dataset using realistic reservoir models with<br />

discrete fracture distributions;<br />

o Analysis and inversion procedure to quantify fracture size and fracture porosity (storability) from<br />

3D VSP, 3D surface and 3D-4C OBC data including both synthetic and real data;<br />

o Multi-Attribute Analysis, including coherency analysis, multi-component instantaneous<br />

polarization analysis, and spectral decomposition of <strong>seismic</strong> data to reveal anomalies that are<br />

related subsurface fracture distributions;<br />

o Comprehensive software packages for fracture modelling and <strong>seismic</strong> imaging including<br />

enhanced 3D modelling capability to simulate <strong>seismic</strong> wave propagation in fractured media and<br />

software to extract information on fracture properties (size/storability, density, orientations, and<br />

fluid properties).<br />

Datasets available for the 7 th Phase:<br />

o Emilio 3D 4C <strong>seismic</strong> dataset<br />

o The Clair field OBC data and associated VSPs<br />

o Physical modeling 3D datasets (CNPC Key Labs)<br />

o Wind River 3D surface data.


TH<br />

7 PHASE RESEARCH PROGRAMME<br />

Fracture characterization<br />

2006-09<br />

Research developments<br />

- Analysis of P-wave attenuation <strong>anisotropy</strong>/shear-wave splitting<br />

- Dominant scale length and fluid saturation from 3D-4C <strong>seismic</strong><br />

- Techniques for modeling the <strong>seismic</strong> response of large fractures<br />

- multi-scale decomposition and attributes analysis<br />

- analysis procedures/user-friendly algorithms & software<br />

Azimuthal frequency anomalies in 3D <strong>seismic</strong> data<br />

East-west azimuth<br />

North-south azimuth<br />

Wind River Basin (Queen and Lynn 2003)<br />

0.44<br />

39.039.239.439.639.8<br />

0.46<br />

Ti me<br />

0.48<br />

0.50<br />

Frequency-independent Frequency-dependent Perpendicular Parallel<br />

0.52


Anisotropic rock physics for frequency-dependent <strong>anisotropy</strong><br />

Frequency-dependent geophysical properties have attracted more and more interests from the research<br />

communities and the industry, which potentially links the subsurface <strong>seismic</strong> response to permeability<br />

and the dominant scale length of the heterogeneities in the subsurface. The situation regarding this<br />

subject is comparable to that of research on <strong>seismic</strong> <strong>anisotropy</strong> in the 80’s: the potential is high yet<br />

more work needs to be done towards both the understanding of the phenomenon and its implications.<br />

During the 6 th Phase, we have developed a dynamic rock physical model accounting for changes in<br />

both pore pressure and saturation. A range of observable effects was predicted, including the<br />

frequency dependence of shear-wave time delays, differential attenuation between fast and slow<br />

shear-waves and P-wave attenuation <strong>anisotropy</strong>. We have analysed a series of VSP datasets in order<br />

to evaluate these effects, including multicomponent VSPs, walkway VSPs and time-lapse VSPs.<br />

These studies have confirmed that these effects are detectable in field data and can be related to<br />

fracture properties such as dominant scale length and fluid saturation, using robust modelling and<br />

inversion procedures.<br />

Phase 7 will develop these ideas across a spectrum of activities, encompassing further theoretical<br />

development, laboratory calibration, development of processing techniques and application to field<br />

data. The goal of the research will be to use the modern ideas in spectral decomposition and<br />

poroelastic theory to establish a robust link between <strong>seismic</strong> <strong>anisotropy</strong> and the petrophysical<br />

properties of rock, particularly the anisotropic variation of permeability. We will continue to analyse<br />

VSP surveys including time-lapse VSP data, to consolidate our understanding, but the new direction<br />

of Phase 7 will be the emphasis on applying our ideas to surface <strong>seismic</strong> data.<br />

Areas of interest and deliverables for the 6 th <strong>phase</strong><br />

o improved rock physics models, calibrated against field and laboratory data, which relate the<br />

frequency dependent anisotropic response to fracturing, fluid saturation and stress state;<br />

o finite difference modeling for wave propagation through a discrete fracture network, to establish<br />

the effect of scattering from fractures of different sizes on the anisotropic response;<br />

o use of attenuation measurements for reservoir characterisation, in particular a comprehensive<br />

theory for interpretation of azimuthal AVO in the presence of frequency dependent <strong>anisotropy</strong>;<br />

o Analysis of time-lapse multi-azimuth VSP data to establish the relationship between the<br />

anisotropic <strong>seismic</strong> response and pressure and saturation variations.<br />

o a fully tested comprehensive rock physics subroutine library<br />

Datasets available for the 6 th Phase:<br />

o Laboratory measurements of synthetic fractured rocks with known fracture geometries<br />

o Multi-azimuths VSPs from the Wyoming Basin<br />

o Time lapse 3D VSPs from Weyburn<br />

o Time lapse VSPs from Quarnalam/PDO.


APPENDIX A – PROJECT PERSONNEL<br />

The EAP team presently comprises five members of BGS staff, a senior consultant in the Unversity of<br />

Edinburgh, and six PhD Students. It also benefits from its close association with BGS staff working in<br />

other areas of seismology, and associated academic staff in the University of Edinburgh.<br />

EAP Staff<br />

Technical director:<br />

Professor Xiang-Yang Li (Multicomponent seismology)<br />

Principal Geophysicist:<br />

Dr. Enru Liu (Associate director and Fracture imaging)<br />

Senior Geophysicists:<br />

Dr. Hengchang Dai (Anisotropic imaging)<br />

Dr. Mark Chapman (Rock physics model and fluid-rock interaction)<br />

Senior Consultant:<br />

Professor Stuart Crampin (Fundamentals of <strong>anisotropy</strong>)<br />

(University of Edinburgh)<br />

Associated Staff<br />

Senior Seismologist:<br />

Technical Engineer:<br />

Dr David Booth (Project support)<br />

Mr Brian Bainbridge (Technical and computing support)<br />

Associated Academics<br />

University of Edinburgh:<br />

Professor Ian Main (Rock Physics)<br />

Dr Roger Scrutton (Log data intepretation)<br />

Dr Andrew Curtis (Exploration seismology)<br />

EAP Students<br />

Jinghua Zhang (4C <strong>seismic</strong> and converted-wave)<br />

Lifeng Wang (4C <strong>seismic</strong> and fractures)<br />

Isabel Varela (Fracture modelling)<br />

Zhongping Qian (4C <strong>seismic</strong> and fractures)<br />

Adam Wilson (Frequency-dependent <strong>anisotropy</strong>)<br />

Joao Marcos Ferreira (Micro <strong>seismic</strong> monitoring)


APPENDIX B - EAP ACTIVITIES AND DATASETS<br />

Converted-waves COMPLETED ONGOING PROMISED<br />

Marine 4C <strong>seismic</strong>s<br />

BP – Valhall(2D)<br />

Agip – Johnd Field 2D<br />

BP - Lomond 2D<br />

Agip – Emilio 3D 4C<br />

BG – 4C<br />

BP-Vauhall Life<br />

of Field Seismic<br />

Mahagony 2D<br />

Chevron Alba 2D<br />

Afflect 3D 4C<br />

Shell – Guilemot 2D<br />

Chevron – Alba 3D<br />

Land 3C data Changqing 2D SinoPec 3D 3C<br />

Petro-China 2D 3C<br />

Offset VSPs/<br />

3D VSPs<br />

PetroChina – Ordos<br />

VSPs<br />

Fractures COMPLETED ONGOING PROMISED<br />

Agip – Emilio 3D<br />

Clair – 3D 4C<br />

Marine 4C <strong>seismic</strong>s<br />

2D/3D <strong>seismic</strong>s<br />

Saga 2D (Geophone<br />

coupling)<br />

Shell – Guilemot 2D<br />

Shenli 3D wide azimuth<br />

Fina – 2D +3D lines<br />

Norsk hydro<br />

Grange 3D<br />

Physical modelling 3D wide DOE – Wyoming;<br />

azimuth data<br />

wind river basin<br />

Walkaway VSPs /<br />

3D VSPs<br />

Bluebell – VSPs<br />

San Juan Basin VSPs<br />

NW Colorado VSPs<br />

BP Clair VSP<br />

BP Egypt VSP<br />

TIME LAPSE<br />

SEISMIC<br />

Repeat VSPs/<br />

COMPLETED ONGOING PROMISED<br />

Mobil – Camelot<br />

Shell – Quarn Alam<br />

Weyburn VSPs<br />

Saga –Snorre, 15 nearoffset<br />

VSPs<br />

Micro-<strong>seismic</strong>ity<br />

USGS micro<strong>seismic</strong>ity<br />

data<br />

Marine 4C <strong>seismic</strong>s Teal South 3D Vauhall LoFS data<br />

COMPLEX AREAS COMPLETED ONGOING PROMISED<br />

Surface <strong>seismic</strong> Amera Hess – long offset Marathon 3D<br />

two-boat shoot (Flare)<br />

Veritas DGC - 12km<br />

cable<br />

Low frequency data Texaco-2D lines<br />

Veritas – 2D lines

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!