Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Paradigm Webinar Full Waveform Sonic Processing in Geolog - Cleaned

Download as pdf or txt
Download as pdf or txt
You are on page 1of 32

Optimize Full Waveform

Sonic Processing
Diego Vasquez – Technical Sales Advisor.
Paradigm Technical Session. May 18th, 2016.

©
© 2016,
2016, PARADIGM.
PARADIGM. ALL
ALL RIGHTS
RIGHTS RESERVED.
RESERVED.
AGENDA
• Introduction to Geolog.
• Introduction to Full Waveform Sonic Processing
in Geolog.
− Data Preparation.
− Pre-Processing
− Processing.
− Post-Processing.
• Best Practices – Custom menu.
• Live demo.
− DSI Processing.

• Wrap up.

© 2016, PARADIGM. ALL RIGHTS RESERVED. 2


Introduction to Geolog

Common platform
Single / multi well environment
Integrated data store (EPOS)

© 2016, PARADIGM. ALL RIGHTS RESERVED.


Acoustic measurement applications

• Formation compressional (P) and shear (S) wave slowness (inverse of velocity).
• Formation porosity estimation.
• Seismic applications.
− Time-depth curve
− Synthetic seismograms
− AVO

• Formation mechanical properties (Bulk, Shear and Young’s moduli,Poisson’s ratio).


− Borehole stability analysis
− Formation stress and fracture analysis
− Sanding potential

• Modelling (fluid substitution).


• Formation Shear anisotropy.

© 2016, PARADIGM. ALL RIGHTS RESERVED.


Sonic Processing in Geolog

• Vendor independent Full Waveform Sonic Log Processing,


supporting every tool on the market.
• Full set of tools for pre-processing and processing,
allowing to obtain the best semblance and picking
possible.
• Complete set of post-processing routines, including the
calculation of Mechanical Properties and Sonic Anisotropy.

© 2016, PARADIGM. ALL RIGHTS RESERVED. 5


Full Waveform Sonic Processing in Geolog

Anisotropy Anisotropy from Xdipole data, using Alford rotation.

© 2016, PARADIGM. ALL RIGHTS RESERVED.


Data preparation

©
© 2016,
2016, PARADIGM.
PARADIGM. ALL
ALL RIGHTS
RIGHTS RESERVED.
RESERVED.
Load Specification File
• Loads tool specific parameters required for sonic processing into the log attributes of the
sonic waveforms.
• The specification file contains all tool information required for sonic waveform processing.
• If the tool to be used is not handled by default, a new specification file can be easily
created.

© 2016, PARADIGM. ALL RIGHTS RESERVED.


Data Preparation - General
• If waveform recording did not start at time 0, set the correct start time for
waveforms for later processing.
• Exclude a particular receiver if data is bad or, replace bad data by missing values.
• All processing steps in Geolog require Waveform in packed form. Depending on
the data from the field, waveforms can be packed or unpacked as needed.
• Correct Receiver order can be easily check using the Array Data viewer in a
layout, and corrected using a dedicated module.

© 2016, PARADIGM. ALL RIGHTS RESERVED.


Pre-processing

©
© 2016,
2016, PARADIGM.
PARADIGM. ALL
ALL RIGHTS
RIGHTS RESERVED.
RESERVED.
Pre-Processing – Waveform Filtering
• Time-Average Filter:
− Remove regular event (e.g: low
background frequency), by calculating
an average in time direction and
extract it from the original value.

Original WF with Time averaged


low background WF with cleaner
noise background

© 2016, PARADIGM. ALL RIGHTS RESERVED.


Pre-Processing – Waveform Filtering
• Depth-Average Filter:
− Remove regular event (e.g: casing
signal, tool generated noise) by
calculating a moving average in depth
direction.
− User can then remove this average
from the original data (remove
coherent casing signal) or output the
mean (removes “chevron” patterns)

Original image with Casing noise removed


casing noise after moving depth
filtering

© 2016, PARADIGM. ALL RIGHTS RESERVED.


Pre-Processing – Waveform Filtering
• Frequency Filter:
− To remove unwanted frequency
component.

Original image Filtered waveform

© 2016, PARADIGM. ALL RIGHTS RESERVED.


Pre-Processing – Waveform Filtering
• F-K Filter:
− To remove or isolate reflections (i.e.
chevron patterns), which might be
caused by fractures, dipping bed
boundaries or abrupt changes in hole
diameter

~ 22
m

© 2016, PARADIGM. ALL RIGHTS RESERVED.


Processing

©
© 2016,
2016, PARADIGM.
PARADIGM. ALL
ALL RIGHTS
RIGHTS RESERVED.
RESERVED.
Semblance processing

… …
Moving time step

Slowness end Semblance map at


Moving time each depth
window length

Slowness step

Slowness start
… …
Processing start time Processing end time

© 2016, PARADIGM. ALL RIGHTS RESERVED.


Semblance processing

Semblance map at each depth Projection log

© 2016, PARADIGM. ALL RIGHTS RESERVED.


Semblance Processing
• Semblance processing allows controlling both the processing time and
slowness limits, as well as the windows length and step.

© 2016, PARADIGM. ALL RIGHTS RESERVED.


Picking
• Auto picking with seed or existing curve (manual curve insert in layout)
• Data could be adjusted manually in a layout, and/or smoothed as needed.

© 2016, PARADIGM. ALL RIGHTS RESERVED.


Post-Processing

©
© 2016,
2016, PARADIGM.
PARADIGM. ALL
ALL RIGHTS
RIGHTS RESERVED.
RESERVED.
Post-processing and Utilities
• Calculate Mean:
− Calculate mean of 2 slowness logs.
• Dispersion Correction:
− Calculates the dispersion curve of the flexural wave and applies the correction at each
depth position.
• Traveltime:
− Calculates the travel time of a selected arrival at a given receiver.
• Create other attributes:
− For in-depth analysis of waveform, instantaneous phase, amplitude and frequency
component can be extracted.
• First Arrival Detection:
− This module is to detect the first arrival and pick the onset time automatically.
• Mechanical Properties:
− Calculate VP/VS, Poisson ratio and elastic parameters.

© 2016, PARADIGM. ALL RIGHTS RESERVED.


Mechanical properties
• Calculate VP/VS, Poisson ratio and elastic parameters.

© 2016, PARADIGM. ALL RIGHTS RESERVED.


Anisotropy Analysis

©
© 2016,
2016, PARADIGM.
PARADIGM. ALL
ALL RIGHTS
RIGHTS RESERVED.
RESERVED.
Anisotropy analysis – Alford Rotation
• Cross Dipole tools have 2 dipole transmitters and 2
sets of multiple receivers placed 90 degrees apart.
• This configuration allows recording 4 sets of
waveforms, 2 inline (XX, YY) and 2 cross-line (XY and
YX).
Y

X Receivers Receivers
(X) (Y)

Transmitter Transmitter
(X) (Y)

© 2016, PARADIGM. ALL RIGHTS RESERVED. 24


Anisotropy analysis – Alford Rotation
• When the tool is aligned with the anisotropy axis, the
shear wave splits on fast and slow directions:
− Energy of the inline waveforms would be maximum.
− Energy of the Cross-line waveforms would be 0.
• As the orientation of the tool is recorded, the orientation Y
of the fast and slow shear wave propagation can be
obtained. X
• The Alford method: mathematical technique that allows
reconstructing measurements done at any orientation.
Fast
Shear

Slow
Shear

© 2016, PARADIGM. ALL RIGHTS RESERVED. 25


Anisotropy analysis in Geolog
Geolog uses the Alford Rotation for Sonic Anisotropy. Three
indicators of formation anisotropy are calculated:
• Energy anisotropy:
− Principal measure of anisotropy. It is less affected by processing, thus
more reliable than the other anisotropy indicators.
− Energy anisotropy is the percentage of the cross-components (XY, YX)
relative to all four components (XX, XY, YY, YX). The minimum and
maximum cross-line energy for each depth frame are computed.
− The difference between the maximum and minimum cross-energy is a
measure of the strength of anisotropy.
• Slowness anisotropy:
− Difference between the fast and slow slownesses calculated by
semblance processing from the rotated waveforms. It is calculated by
dividing the slowness difference by the average of the fast and slow shear
slowness to get a percentage difference.
• Travel time anisotropy:
− Difference in arrival time between the fast and slow shear waves at the
receivers. The travel time is obtained by first arrival picking from the fast
and slow waveforms.

© 2016, PARADIGM. ALL RIGHTS RESERVED. 27


Anisotropy analysis Workflow
• Calculate DTSM and Travel time for first receiver.
• Amplitude Recovery:
− Converts waveform amplitudes back to its original form from gain control and/or
normalization factor applied to data during data acquisition. Waveforms may be
unpacked after recovery.
• Cross Dipole Processing:
− Calculates Fast and slow Waveforms, Fast Shear Azimuth, Energy Map and
Anisotropy Map.
• TT and DTS Anisotropy:
− Calculates Traveltime and Shear Anisotropy, using the shear wave picked from the
fast and slow waveforms.
• Analysis:
− Automatic analysis to select the anisotropic zones, by analysing the three indicators.

© 2016, PARADIGM. ALL RIGHTS RESERVED.


Best Practice – Custom Menu
• Creating a custom menu per tool type can be
an invaluable tool for processing sonic data in
Geolog.

• The same module can be included different


times, pointing to a different spec file relevant
to such step.

© 2016, PARADIGM. ALL RIGHTS RESERVED. 29


Full Waveform Sonic – Live Demo

©
© 2016,
2016, PARADIGM.
PARADIGM. ALL
ALL RIGHTS
RIGHTS RESERVED.
RESERVED.
Wrap Up

• Vendor independent Sonic Log Processing, supporting


every tool on the market.
• Full set of tools for pre-processing and processing,
allowing to obtain the best semblance and picking
possible.
• Complete set of post-processing routines, including the
calculation of Mechanical Properties and Sonic Anisotropy.

© 2016, PARADIGM. ALL RIGHTS RESERVED. 31


Thank You!
Diego Vasquez
Technical Sales Advisor – Formation Evaluation
Diego.Vasquez@pdgm.com
+61-8-9327-1800

©
© 2016,
2016, PARADIGM.
PARADIGM. ALL
ALL RIGHTS
RIGHTS RESERVED.
RESERVED.
References
• The Well Logging Handbook
O. Serra - Editions Technip, 2008 (ISBN 978-2-7108-0912-8)

• Quantitative Borehole Acoustic Methods


X.-M. Tang and A. Cheng – Handbook of Geophysical Exploration, Elsevier 2004 (ISBN 0-08-044051-7)

• The Rock Physics Handbook (2nd ed.)


G. Mavko, T. Mukerji, J. Dvorkin – Cambridge University Press, 2009 (ISBN 978-0-521-86136-6)

• Slowness-frequency projection logs: A new QC method for accurate sonic slowness evaluation
T. Plona, M. Kane, J. Alford, T. Endo, J. Walsh and D. Murray (Schlumberger) - SPWLA 46th Annual Logging Symposium, June 26-
29, 2005

• Isotropic AVO Methods to Detect Fracture Prone Zones in Tight Gas Resource Plays
Bill Goodway, John Varsek and Christian Abaco (EnCana Corp., Calgary, AB, Canada) – 2007 CSPG CSEG convention

• Shear data in the presence of azimuthal anisotropy


Alford, R. M. - 1986, 56th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 476-479.

© 2016, PARADIGM. ALL RIGHTS RESERVED.

You might also like