HowTo Analyse Angles

OIL & GAS / HOW TO
ANALYSE INCIDENCE ANGLES WITH ISATIS
AVO (Amplitude Versus Offset) and AVA (Amplitude Versus Angle) studies require angle stacks where the Common Middle Point (CMP) gathers are stacks according to the incidence angle at the main reflector . The incidence angle analysis provides maps to help the geophysicist assisted by the geologist to perform this stacking process. Areas of sub-stacks (typically near , mid, far) are defined from the angles and transcribed in offsets for effective implementation. Our current workflow is based on the Walden approximation to transform offset to incidence/reflection angle. Maps of extreme incidence angle related to mute function and fold map for sub-stacks are built. Each extreme angle is determined by a mute function which eliminates noise and/or possible refractions overlapping reflections.
Walden formula: = arcsin {VDIX / VRMS * sin [arctan (Offset / (T * VRMS)]}

Where: is the incidence angle, offset is the trace offset, T is the time, VRMS is the stacking velocity provided and VDIX is the Dix velocity computed as follows:
Dix formula: vDix = [(t2vrms22 t1vrms12) / (t2 t1)]1/2

where t1 and t2 is the seismic time from stratum between two reflectors.
1. Input Data and Parameters

RMS Velocities 2D or 3D Horizons in time Data of mute: To minimize air waves, ground roll and other early-arriving noise by removing the
corresponding part of the signal. Low-frequency traces and long-offset traces are typical targets for muting. Within the gather two types of mute laws are used: o o Internal mute law, used to determine the min offset to use External mute law, used to determine the max offset to use.
offset: distance between seismic source and receiver .
2. Workflow
In the case of tabular environments, the determination of the angles can be performed at any point in space by applying the Walden formula from a velocity field (VRMS), from a DIX velocity field and from the definition of the internal and external mute laws (offsetmin and offsetmax): min (X,Y ,T) = Walden [T , VRMS (X,Y ,T), VDIX (X, Y , T), offsetmin (X,Y ,T)] max (X,Y ,T) = Walden [T , VRMS (X,Y ,T), VDIX (X, Y , T), offsetmax (X,Y ,T)] min max are assigned at target horizons: min horizon = Walden (T horizon, VRMS horizon, VDIX horizon, offsetmin horizon) max horizon = Walden (T horizon, VRMS horizon, VDIX horizon, offsetmax horizon) Flowchart: Data Quality control: horizons (2D) and VRMS (3D) Dix Velocity calculation Velocity extraction at target horizons: VRMS and VDIX. Specification of offsetmin horizon and offsetmax horizon Angle determination with Walden Formula Fold map extraction Output maps analysis
3. Data Quality Control

The first step consists in checking the concordance between velocity cube VRMS and horizons. This step is essential when dealing with seismic data. Data control analysis allows detecting potential artifacts, problems of support, acquisition and interpolation during different processing steps. The first task is to carry out a proper data investigation:
Statistics Exploration Data Analysis
Display of a data basemap, in order to validate the data location; check the concordance between horizons
and velocity cube.
Compute a histogram to check the homogeneity of the statistical distribution of VRMS and time horizons and
between them.
Control the statistical characteristics: the mean (or average tendency) and the variance that measures the
dispersion of the values around the mean for time and velocity data.
Investigate correlation between several variables through the display of the scatter diagram.
The second task is to estimate large or medium scale structures in the velocity cube VRMS.
Calculation of the velocity trend:

Interpolate Quick Interpolation Least Square Polynomial fit More explanation about trend is given in the file How to compute quickly a reliable Trend that you can download from the Geovariances website.
Calculation of residual velocities:

File Calculator Residual velocities = Raw Data Trend velocities The variogram analysis is important to check the presence of potential anisotropies or if we need filtering some artifacts. For more information about filtering, please refer to white paper about filtering or the file How to remove quickly structured noise from a seismic that you can download from the Geovariances website.
4. Dix Velocity
Calculation of the Dix velocities:
File Calculator
VDix(x,y,[tn-1 ;tn] = sqrt {[(VRMS(x,y,tn)*tn)-(VRMS(x,y,tn-1)*tn-1)]/(t1-tn-1)}

where tn1 and tn (n=1max reflector) are reflection arrival times, in our case the sample time.
5. Extraction of RMS and Dix velocities at horizons:

Tools Time Sampling From Volume A simple verification of Dix velocity can be done with the ratio VDix/VRMS calculated with calculator . This ratio has to be between 1 and 2.
6. Determination of extremes offset

The mute law is provided by the geophysicist and defines offsetmax and offsetmax, therefore the offset range for each X,Y ,T position. The mute law is defined in order to discard the traces of poor quality , knowing that the signal degrades near limit angle (angle beyond which waves are refracted):
External Mute or offsetmax :

o offsetmax = constant for deep and ultra-deep environment (depth at which we have critical angle). o offsetmax = f (Time) for shallow environment (the trace at a given offset is only reliable below a given time). o offsetmax = f (x,y ,Time) for shallow environment with highly sloped seafloor (all offset are available only from a certain depth that varies with the seafloor).
Internal Mute or offsetmin :

o offsetmin = f (Time) for shallow environment (the trace at a given offset is only reliable above a given time). o offsetmin = f (x,y ,Time) for shallow environment with highly sloped seafloor (all offset are available only until a certain depth that varies with seafloor). The selected offset values are used in angles determination formula.
7. Angles determination
This step allows defining the two angles used to characterize the fold map. Use the Walden formula to calculate the minimum and maximum angles depending on offsetmax and offsetmin.
Angle calculation:
File Calculator min = asin (VDIX / VRMS * sin (atan (Offsetmin / (T * VRMS))) max = asin (VDIX / VRMS * sin (atan (Offsetmax / (T * VRMS)))
This can be done in 2D for each target horizon (where T is the horizon) or in 3D where T is the seismic time.
8. Output maps analysis

In 3D, we can do an extraction at the target horizon:
The angle min is extracted at shallow horizon and the angle max is extracted at deep horizon. If deep horizon
is not available we extract a horizon calculated as follow H = Top + thickness of the target (ms). Tools Time Sampling Time Sampling From Volume In theory, the minimum of max (min(max)) and maximum min (max(min)) angles define the interval from which to create sub-stacks (corresponding incidence angle ranges). Generally in practice, these statistics are not chosen as they often lead to a too low angle range, which limits the operational definition of sub-stacks. It is the geophysicist, in collaboration with the geologist in charge of the interpretation of the data, who decides the angle min(max) and max(min) values on the basis of the analysis of min horizon and max horizon maps. To help the decision making, cumulative curves are used to quantify the percentage of points discarded for the selected angles.
Cross plot angles / time, to position the envelope curves min and max.
Statistics Exploration Data Analysis: crossplot and / or cumulated histogram
1000 rho=-0.800
100.00
2000
75.00
3000 time
Cumulated Frequencies
50.00
4000
25.00
5000
6000
10
15
20
25
30 angle H1
35
40
45
50
55
60
0.00
6 angle H1
10
Minimum angle: 10% of H0 and 50% of H1 values are above 6. For more information, let us help you assessing the quality of your data and implementing this workflow. Contact our consultants at [email protected].
GEOVARIANCES offices France - Avon-Fontainebleau - +33 (0)1 6074 9090 Australia - Wynnum-Brisbane, QLD - +61 (0)7 3348 5333 [email protected] - www.geovariances.com

HowTo Analyse Angles

Uploaded by

Copyright:

Available Formats

HowTo Analyse Angles

Uploaded by

Document Information

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

HowTo Analyse Angles

Uploaded by

Copyright:

Available Formats

OIL & GAS / HOW TO

ANALYSE INCIDENCE ANGLES WITH ISATIS

Walden formula: = arcsin {VDIX / VRMS * sin [arctan (Offset / (T * VRMS)]}

Dix formula: vDix = [(t2vrms22 t1vrms12) / (t2 t1)]1/2

1. Input Data and Parameters

offset: distance between seismic source and receiver .

3. Data Quality Control

Statistics Exploration Data Analysis

Calculation of the velocity trend:

Calculation of residual velocities:

VDix(x,y,[tn-1 ;tn] = sqrt {[(VRMS(x,y,tn)tn)-(VRMS(x,y,tn-1)tn-1)]/(t1-tn-1)}

5. Extraction of RMS and Dix velocities at horizons:

6. Determination of extremes offset

External Mute or offsetmax :

Internal Mute or offsetmin :

8. Output maps analysis

You might also like