PETREL Shell - Refresher - Day3
PETREL Shell - Refresher - Day3
Process Manager
Process Manager
Synthetics
Utility Systems
Any Core
Well Design
Data Analysis
OpenSpirit Plug-in
Geophysics
Seismic Volume Rendering and Extraction
Multitrace Attributes
Well Correlation
Combined Core
Geology
Facies Modeling
Petrophysical Modeling
Engineering
Reservoir
Engineering
Reservoir
FrontSim Locked
Core
Geoscience Core AND Seismic Interpretation AND/OR Seismic Vol. Rend. & Extract
3rd
1
3
Similar options for posting 3D and 2D seismic are available for 3D windows as well.
The arbitrary polyline intersection is composed of 2 parts, the polygon trace and the vertical
seismic intersection. The polygon trace is stored under the General Intersection object in the
Input tab. Remember these “General Intersections” can be renamed to something more
meaningful.
When entering the points for the polygon, the movement of the mouse will be restricted to either
the inline or crossline direction.
1 3 4
5
2
Ghost
The Ghost curve is used to create a small bitmap of reflector events on a seismic line. The area
can be moved to other parts of the same seismic line (or other seismic lines) to compare signal
pattern and identify the same horizons across a fault e.g.
When a ghost is created it will be located in the seismic interpretation window folder in the
Windows tab.
The seismic in the interpretation window and the ghost share the interpretation window display.
If you click on the ghost in the interpretation window you will notice that the ghost is activated in
the windows tab, if you click on the seismic displayed in the interpretation window you will notice
that the Interpretation 1 under the Interpretation window is activated. In practice this means that
you can have different settings for the ghost and the seismic in the interpretation window, e.g.
the ghost can be displayed with wiggles or another attribute while the seismic in the
interpretation window is displayed in a conventional mode.
From the windows tab the settings and administrations for the ghost is accessed. When the
ghost is bolded/activated the display options for the ghost can be changed.
In the windows tab you can toggle the ghost on/off, delete ghosts and access the settings for the
ghost. Note that you can have as many ghosts you like, but it is recommended to clean up the
ghosts once in a while.
1 4
3
2
If changing the increment or navigating to a line not falling on the used increment, the
Next/Previous interpretation will not display.
The Survey manager is a tool used to sort and manage the seismic
data in Petrel. Within this the user can control all parameters
regarding 2D lines, like CDP, SP, Trace number, navigation file
assignment, sort and move data into folders etc. The Survey
Manager can also be used to visualize data.
To use it, right click on the survey folder and select Survey Manager.
Mis-tie manager
The Mis-tie manager is an interactive tool for managing the mis-ties in Petrel. With this the user
can calculate or specify corrections, select reference lines, deselect lines not to be used for mis-
tie analysis and also toggle to display 2D lines in the active window.
Mis-tie Analysis allows the user to extract and apply mis-tie value corrections for 2D seismic
surveys. All lines contained inside a Survey Folder or a sub-folder will be used when a Mis-tie
Set is created.
2
1
3
4
5
6
7
There are three options when creating an attribute map. That is to create a new surface, add the
attribute to an already existing surface or overwrite an attribute.
Static attribute maps are calculated based on the input as it is at the time of creation. A dynamic
attribute map will be recalculated if any input data are changed.
Note: this process only works with surfaces, not horizon interpretations directly.
2
1
3
4
5
6
8
7
1
3
4
6
5
Process Manager
The Process manager is a very useful tool to do looping processes and operations, or redoing
operations on an already existing 3D Grid.
This example shows how to make surfaces from several interpreted events and generate
attribute maps from the resulting surfaces. The Make/Edit Surface and Make attribute map
processes need parameter inputs and must be opened (double click) to give in appropriate data.
1
4 3
Generated attribute maps are attached to surfaces as various attributes. To make the attribute
map as a ‘stand alone’ map, it has to be converted.
3
1
4
2
5
A Variogram map displays the separation distance in X-direction and Y-directions, not ordinary
XY coordinates. This means that the variance map must be displayed in a map window since the
center of the map is 0 (zero) X- and Y separation with increasing separation in X- and Y-
directions respectively.
4
3
The combined variogram information derived from the attribute map can be used in facies or
petrophysical modeling. If a relation excists between the attribute map and distributions of facies
or any petrophysical parameter, the sample variogram can be used for establishing a modeled
variogram that can be used as input.
This exercise chapter describes the new functionalities that are part of the Shell Global Deal. We
will open and use a Petrel project where necessary data is stored or import the missing data
when needed.
Exercise overview:
• Seismic data visualization in base map window
• Create random intersections
• Seismic Ghost
• Neighboring interpretation
• Survey Manager
• Mis-tie analysis
• Make Attribute Map process
• Seismic Workflows
• Variograms from seismic attribute maps
One Inline and Crossline will automatically be generated when importing a 3D seismic cube but
they are not visualized before they are toggled on.
Exercise Steps
1. Open the Petrel project Shell_Seismic_startup_project.pet. It contains realized 2D and
3D seismic data that will be used in the following exercises.
2. Expand and then double click on the 3D seismic folder .
3. A Settings window for the Survey folder pops up. Here you can define the Base Map
Annotation. Review the default settings and Cancel the window. Note the settings will apply
to both the 3D and the Base Map (2D) windows.
9. Toggle on the checkbox in front of the 2D Lines folder. The 2D lines are now displayed in
the 2D window.
Exercise Steps
1. Activate the previously used 3D Window. Right click on mig.sgy [Crop] [Realized] and
select Load into memory when needed. Right click again on the seismic volume and select
Insert Time Slice Intersection.
2. Expand the window to use the full screen. Deselect the other seismic lines if any.
3. Make sure the Seismic Interpretation process is active.
4. Select the “Create arbitrary polyline intersection” icon in the Function bar. Draw a
polygon using the left mouse button and double click to stop the polygon line. A seismic
random line is automatically generated and is stored under the mig.sgy [Crop] [Realized]
folder.
5. To edit the Random line, open the “General Intersection” folder, display and activate the
Arbitrary Polygon.
6. Activate the Make/Edit Polygons process under Utilities in the Process Diagram.
7. Select the “Select and Edit/Add Points” icon . Choose a point on the polygon and drag
it to a new position. Notify that the seismic random line will automatically be updated.
Exercise Steps
1. Open a new Interpretation window (preferably a crossline) and identify a discontinuity.
2. Activate Insert Seismic Ghost . Drag a rectangle over the zone of interest.
3. Place the cursor inside the seismic ghost area. Move it to the desired position.
4. Identify the corresponding events.
Neighboring interpretation
Another help when interpreting in an Interpretation window is the Show Previous/Next
interpretation functionality. This displays the interpretation from the incremental (controlled from
the Select plane step increment in the lower left corner of Petrel) neighboring seismic lines.
Exercise Steps
1. Activate the Seismic Interpretation process under Geophysics. In a 3D window, display an
inline and crossline from the realized seismic data and use Seeded 3D Autotracking to
interpret on a strong event.
2. Open an Interpretation window and display the crossline.
3. Open Settings for the newly interpreted horizon. Go to the Style tab. Under Neighboring
toggle on Show Previous and Show Next.
4. Select color and size for the display and observe the interpretation window.
5. If needed, change Select plane step increment from the Player
to fit an interpretation.
Exercise Steps
1. Right click on the 2D Lines folder and select Survey Manager from the appearing menu.
2. View the active 2D window while toggling any of the checkboxes in the eye column.
This column turns on/off the line display.
3. Click and drag the cursor over several rows in the numbered column to select lines.
10. If none of the options on the left is toggled, the search engine will search for the words
containing the specification in the field. The search starts from the top of the Survey
Manager.
Comments
The Survey Manager works for the survey folder it is started from. To easily manage all 2D lines
in a project, keep all surveys in the same folder. Sort the different surveys into sub-folders at
30 • New Seismic functionalities - Exercises Shell Refresher - Day 3
several levels. If a project with separate survey folders exists it is possible to drag and drop lines
between folders.
A survey Manger can also be opened from a 3D survey folder. Instead of separate lines 3D
volumes will be displayed.
Exercise Steps
1. Right click on the 2D Lines folder and select Insert Mis-tie Set from the appearing menu.
2. The Mis-tie Manager opens up displaying the 2D lines crossing each other.
3. Toggle on Gain and Vertical leave Phase toggled off. Click on Compute Mis-ties and
observe the values entered into the different columns.
4. Use Constant as method for calculating mis-tie corrections then click Compute Corrections.
5. The Gain correction and Vertical correction columns are filled in. These computed values
could be edited manually. In the Vertical correction column, type in –10.0 (10 ms downward
relative to the crossing lines) for all entries of slb2.sgy [Realized].
6. Different 2D lines can be used and locked in this process by toggling the appropriate boxes
in the and lock columns. Use all lines, but do not lock any of them. When satisfied
with the setup, click Realize to generate physical versions of the mis-tied lines.
7. Click Yes in the appearing message window. Use Create new folder and name the folder
mis-ties.
8. These lines can now be found in the mis-ties folder.
9. Open a 3D window and display the original lines by toggling boxes in the first eye column
to the left in the Mis-tie Manager. The realized mis-tied lines are displayed by toggling from
the second eye column. Click OK to save and close the Mis-tie manager.
Comments
A Mis-tie set is created for the survey folder or sub-folder it is inserted from. All lines in that
folder will be displayed in the mis-tie manager including all sub-folder lines.
Attribute Maps
Generating an attribute map is now treated as a proper process in Petrel. Functionalities related
to this process is essentially the same as in the previous releases of Petrel, but the resulting
attribute map is associated with a surface.
Exercise Steps
1. Expand the Surfaces folder and display the HOUSTON surface in a 3D window.
2. Open the Make Attribute Map process under Geophysics.
3. Use Create new surface (at seismic resolution). Use the default name or give in your own.
4. Use Static map type and drop in the mig.sgy [Crop] [Realized] as Seismic. Use
Instantaneous, select and drop in surface HOUSTON. Apply to generate the attribute map.
Exercise Steps
1. Make sure an active 3D window is open.
2. Click on Insert from the Menu bar and select New Workflow.
3. An empty workflow opens up. Change the workflow name and Author if needed. Type in a
description in the available field if preferred.
10. From Utilities, find and insert Move. Variable A is the Object to be moved. Insert a new
folder into the Input tab, give it a name and input it as Folder to move it into.
11. Insert Clear All Visualizations and View All into the workflow.
12. Also insert Set visible and use Variable A as input.
13. Insert Set Reference into the workflow, Variable B should equal Variable A in the
expression.
7. Use default settings for Transform and XY range. Under the Output tab, make sure As a
surface is toggled. Execute.
8. Open a new Map window and display the Variance map.
9. Observe the elongated nature of the (variance) contours.
10. Use the Measure Distance tool and follow the direction of lowest variance (ESE –
WNW). Read out the direction from the information area in the lower right corner of Petrel. It
will be used in the next exercise step.
11. Reopen Settings for the attribute point dataset. Go to the Variogram tab.
Schlumberger Confidential
• Correct sonic log and establish Time-Depth relation
• Make Acoustic Impedance log and Reflection Coefficient Series
• Generate Synthetic and Extracted Wavelets
• Produce Synthetic Seismogram
• Compare seismogram with real well seismic
• Manual adjustment and event picking
• Display of synthetic trace in 3D window
1 2
3. Expand the Global well 4. Color fill logs by using the Create/edit curve fill tool
logs folder and select the
available logs to use 4
When selecting logs from the Global well logs folder, all wells containing the selected logs will
have them displayed in the Well section panel.
The other option is to select the logs under each well that will be a more inefficient procedure to
follow.
3. Sonic correction
2
4. Time/Depth lookup
3
5. Acoustic Impedance
4
6. Reflection Coefficient 5
6
7. Generate Wavelet 7
8
8. Synthetic Seismogram
9
9. Correlate Seismic 10
10. Manual Adjustment
Synthetic Seismogram
Synthetics is a separate process in the Process Diagram tab, found under Stratigraphic
Modeling.
Input
Fundamental input in the Synthetics process is a sonic log, a density log and check shots. If
density logs are missing, estimated density logs can be derived from sonic logs.
In addition must a mathematical wavelet be defined based on known parameters from the
seismic data.
The sonic log is corrected by using check shots, giving a velocity log, which is then multiplied
with the density log. The resulting accoustic impedance log gives a reflection series which, when
convolved with a wavelet gives the synthetic seismogram. Displayed together with conventional
seismic data from around the well position, a well tie to seismic can be made.
Manual adjustments based on known well tops can be performed.
1
2 3
This log calculator procedure shown, is used to produce a copy of the original input DT log to
avoid overwriting it.
1
2
Despiking the sonic log is used to remove spike values. Since these error readings will
accumulate down through the well when integrating the sonic values, they need to be removed
before sonic correction is done.
5
1
6
Log estimation in Petrel can produce a range of different logs based on real data. The quality of
any estimated log is only as good as the input permits.
Here the density log is estimated by Gardner’s Approximation. Default parameters are used, but
will typically change with the sediments depth of burial.
1. Double click on the 2. Move the mouse over each of the different icons to see hints
Synthetics process under how to use them
Stratigraphic Modeling
1
2
3
Schlumberger Confidential
4. OK (one file) or OK For All (two or
more files if in same format) other time/depth information
6
2
3
B
4
Time
1
2
Depth
4
1
2
4 3
1 2
Phase
Ormsby
3
Wavelet shape
Amplitude Klauder
Butterworth
Ormsby – wavelet designed from four points forming a trapezoid bandpass filter.
Butterworth – defined by amplitude, frequency, high/low cut frequency and high/low order.
Bartlett
1. Select Extracted type of 1
seismogram
2. From Input tab, select 2 Phase
4
seismic data and Well. Drop
Cosine 3
both ino the dialog box
using the drop in arrows Amplitude
3. Set Neighbourhood extent
4. Select Taper window type
Hamming
Wavelet shape
5. Alternative Type of
seismogram is loaded from 5
file Hanning
6
6. Browse for the file and open
it in the dialog box
7
7. Set Time unit
Papoulis
Extracting wavelets in Petrel is a method of statistical describing the shape of the wavelet.
The extraction option gives statistical wavelets from the seismic data at a borehole based on a
portion of a 3D volume of seismic traces.
The taper window is essentially a function of time that is multiplied by a data segment. Its
purpose is to smooth or shape the resulting spectrum.
Signal:
S=R*Wavelet
1 2
=
3
*
Extracted Well Seismic in Petrel follows the well path of the different wells. For a deviated well,
this means that extracted traces for the deviated part is different from the traces higher up in the
well.
1
4
2
When Manual Adjustment is performed (toggled), the calculated TWT Auto (se previous slide) is
recalculated using the updated TWT Picked values. This is a reversible process since toggling
off the Manual Adjustment restores the TWT Auto to the original that is again based on Time
settings for the well.
5 7
6
9
Manual adjustment of the synthetic to the real seismic is based on well tops. If no well tops are
available to the project, users can create a set of artificial well tops at levels where easily
identified events can be observed.
1
2
3
EXERCISE
Synthetic Seismogram
The fundamental difference between well data (in depth) and seismic data (in time) is bridged in
the process of generating synthetic seismograms. When using velocity information from wells, a
time – depth relation can be established for each well. All logs are then attached to both depth
and time domain and by combining different log data, synthetic seismic traces following the well
paths can be made. Comparing the synthetic seismogram for a well with conventional seismic, a
starting point for seismic interpretation can be defined.
Exercise overview:
• Make Well Section
• Import check shots
• Correct sonic log and establish Time - Depth relation
• Make Acoustic Impedance log and Reflection Coefficient series
• Generate Wavelets and make Synthetic Seismogram
• Compare seismogram with real well seismic
• Manual adjustment and event picking
• Display the synthetic trace in a 3D and Interpretation Window
Exercise Steps
1. Go to Window in the main menu and select the New Well Section Window.
2. Toggle on wells Apatite-E13, Copper-6, Dolomite-B1 and Diamond-14.
3. Open the Global well logs folder and toggle on GR, DT, RHOB and CALI. Make a visual
examination of the different logs in the Well Section window.
4. Save your project.
Comments
Although only the sonic and density logs are used for generation of synthetic seismograms,
Gamma Ray and Caliper logs are helpful when preparing the data for this process. The GR log
Exercise Steps
1. The Synthetics process is found in the Process Diagram under Stratigraphic Modeling.
Activate the Synthetics process dialog by double clicking on it. A description of the process
is found by resting the cursor over the icon.
2. All wells or one by one well can be selected from the Selected Well list. By toggling on and
selecting a Well Section, changes and generated objects can by displayed from the process
diagram. Use the Well Section created in this chapters first Exercise (Apatite-E13, Copper-
6, Dolomite-B1 and Diamond-14).
3. Use the edited sonic log (Sonic_despiked) as Original Sonic input.
6. Check statistics from the settings icon and open the Spreadsheet by clicking on the Show
Spreadsheet icon . Change the depth column to TVD. Notice that the two first columns
can be edited if needed and that data for all loaded wells are found in the same
spreadsheet.
Comments
In the Checkshots spreadsheet box all grayed out columns are calculated. Avg. Vel and Int. Vel
are derived from the input time – depth pairs, while Sonic Time and Sonic Int. Vel come from the
input Sonic Log. Drift is the difference between the two data sources. Missing rows for sonic
derived values implies that the sonic log do not exist at these levels.
TVD for Checkshots in Petrel means TVDSS (i.e. referenced to MSL).
Exercise Steps
1. Click on the Create new object icon behind the greyed out field for Sonic Correction.
Open the settings dialog box for the corrected sonic and see that the input data can be
changed. Leave it as it is. Notice the two options for correction curve fitting. Cancel the
settings dialog box.
Exercise Steps
1. Create the Acoustic Impedance log. Open the settings dialog box, make sure that RHOB is
used for Density log and Sonic_despiked is used for Sonic log. Click OK.
2. Make the Reflection Coefficients series. Under settings, observe that Acoustic Impedance 1
is input to Reflection Coefficients 1.
Exercise Steps
1. Make a new wavelet by clicking on the icon. Wavelet 1 is created. Open the settings
dialog box and view its content. Under the Info tab, rename it to Ricker Wavelet.
2. Go to the Settings tab, leave the Type to Synthetic, Sample Rate and Length as set. Make
sure Synthetic is Ricker and use Constant Phase 0 (assumes zero-phase seismic). Set the
Peak Frequency to 30 Hz. Make sure the SEG Polarity norm is toggled ON. Click OK.
Comments
The Ricker, Ormsby, Klauder and Butterworth wavelets try to express in a mathematical way the
expected waveform as it passes through the earth. Parameters are set according to known
properties (or estimated) from the seismic data cube. Extracted wavelets on the other hand, are
stochastically defined wavelets that try to minimize the difference between the synthetic
seismogram and the well seismic. Independent of wavelet generation method, adjustments may
be necessary to give a better result when comparing with conventional seismic (next exercise).
Exercise Steps
1. Create a new object for Well Seismic and open Settings. From the Input tab in Petrel
Explorer, find the Seismic folder and the mig.sgy seismic. Dump it in as Seismic Volume in
the settings window, using the blue arrow . Leave Sample Rate at 0.5. Use Radial
extraction with default settings.
2. Click OK.
3. Toggle on the seismic data in the Well Section window. Compare the synthetic and real
seismic data.
4. Can any matching events be found in the first well? Is there any recognizable pattern in the
seismic that can be seen at a different level in the synthetic seismogram? What about the
other well?
Exercise Steps
1. Toggle on the Manual Adjustment Checkbox in front of Well Tops 1 in the Synthetics dialog
box. Observe that no adjustment is done. Toggle off Manual Adjustment again.
2. Right click on Well Tops 1 and open the Spreadsheet. Identify well Diamond-14 and
Surface Kobe. Type in the value for TWT Picked for this well top to 2150 ms. Click Ok, and
toggle on Well Top 1 in the Synthetics dialog box.
3. Toggle on and off the Manual Adjustment Checkbox. Observe the adjustment of the Well
Seismic.
4. Graphically the manual adjustment can be done in the Well Section Window in the following
way. Make sure Well Tops 1 is displayed in the Well Section Window and that the Manual
Adjustment Checkbox is toggled off.
5. Click on the Create/Edit Well Tops icon , then on the Show Well Top Time icon
and finally on the Edit Well Top Time icon . The Well Tops turns into stippled lines.
6. Click and drag the stippled Kobe Fm top in well Copper-6 to another position. If the Well
Tops 1 spreadsheet is open, observe that the TWT Picked for this well top changes
interactively. By toggling the Manual Adjustment on and off from the Synthetics dialog box,
this edit is taken into account.
7. Open a new 3D window and display all wells together with a Crossline. Position the
Crossline close to (but behind) well Copper-6 using the Manipulate Plane icon. Toggle
on Well Tops 1. From the Global Well Logs folder, display the GR log. From the tool bar,
change depth settings for all wells by opening the drop down menu and changing from Any
to TWT. Observe the changes to the position of well paths, log data and well tops relative to
the seismic line. Observe that only well paths for the four wells used for generating
Synthetic Seismograms is displayed. Reset the depth setting to TVD. Observe the effect.
Exercise Steps
1. Open a new 3D window and display all wells. Use the same procedure as described in the
previous exercise to position a Crossline close to any of the wells.
2. From the Global Well Logs list, toggle on Synthetic 1 and view it in the 3D window. Change
between the Ricker Wavelet and Extracted Wavelet in the drop down menu in the
Synthetics Process dialog box. Observe the difference. Also change between TVD and
TWT in the drop down menu over the graphic window.
3. In the Input tab, right click on the Crossline selected in step 1 (this exercise) and select
Create Interpretation Window. Display wells, well tops and Synthetic 1. Change between
TVD and TWT. Observe that in depth you cannot display conventional seismic data, but all
data related to wells (logs, well tops, synthetic traces and extracted seismic along the well
path) can be displayed in TVD.
4. Open Setting for Synthetic 1 and go to the 2D Style tab. Toggle on Show wiggles, set
Trace repeat count to 2 and if needed toggle on and set Manual gain [dB] to 4. Apply the
changes and view the result in the active Interpretation window. You might want to change
the well path thickness to have a better view. Go to Settings for the Wells folder, Style tab
and type in a smaller Pipe width number in the Path tab. Click Apply.
Comments
Display of data in different domains is dependent on whether the data has information of position
of it in that particular domain. Well paths and well tops without time information cannot be
displayed in TWT, just as seismic data without depth information cannot be displayed in TVD.
After Time – Depth relationship for wells are established and velocity models are made for which
seismic can be depth converted, display of data can be done across domains.
• Ant-tracking
• Neural Nets
• Bricked seismic
The new workflow is a top-down approach, where the interpreter interacts with automatically
extracted fault surfaces in a 3D canvas. Several properties are connected with the surfaces,
which the interpreter can use for organizing the data. For example, the surfaces can be split into
groups representing fault systems. The faults that make a system have common strike, meaning
that the same stress field has created them in the same time period.
In most of the cases the interpreter has to use some of the editing tools to edit the fault patches
before converting them to fault interpretation or fault pillars.
Structural Smoothing
Smoothing of the input signal guided by the local structure to increase the continuity of the
seismic reflectors (Randen, 2002). Principal component dip and azimuth computation are used
to determine the local structure. Gaussian smoothing is then applied parallel to the orientation of
this structure.
Chaos
The chaotic signal pattern contained within seismic data is a measure of the “lack of
organization” in the dip and azimuth estimation method. Chaos in the signal can be affected by
gas migration paths, salt body intrusions, and for seismic classification of chaotic texture. No
parameters are required.
Other attributes may be used (e.g. variance, coherency etc). Some testing may be required to
find the optimum attribute and parameters.
Ant-tracking
This unique algorithm is part of an innovative workflow that introduces a new paradigm in fault
interpretation. It emulates the behavior of ant colonies in nature and how they use pheromones
to mark their paths in order to optimize the search for food. Similarly, virtual ants are put as
seeds on a seismic discontinuity volume to look for fault zones. Virtual pheromones deployed by
the ants capture information related to the fault zones in the volume. The result is an attribute
volume that shows fault zones very sharp and detailed. A second optional output is the
automated extraction of “fault patches” out of the volume.
The Ant-tracking attribute cube and the Fault Patches will be put in the Input tab
On the settings for fault Patches, select the Style tab. You can define the method to use to color
the patches on the 3D Window. The default is “Different”, which will assign a different color for
each fault patch. You can display a different property:
1. Click on the Edit Extracted Faults process to step in Process diagram. The corresponding
tools will then become available on the function to the right of the Display window.
2. Activate the Fault Patch object inside the Fault Patch folder you want to work with in the Petrel
explorer.
3. Toggle on the object to display all the fault patches available in the folder.
4. Double-click on the Fault Patch object or use the right mouse-button to access the settings.
This will open the settings dialog for the Fault Patches and allow you to interact with the
displayed fault patches. Select the tab to perform the desired operation on the patches.
Seismic – Ant-tracking
Settings for Fault Patches
2. Go to the Histogram
2
tab and filter on the
Fault Patch Size.
Click on the histogram
representing the small
faults and click Apply
You can filter fault patches using the properties and the histogram tool. This is a powerful
method to discriminate fault patches using any combination of properties.
Filter settings work in an additive mode: so you can remove patches based on size and then
change the filter to another property and the second filter will work with the result of the first.
1
2
Stereonet Tab
The Stereonet tool can be used to filter a priori area where the Ant agents should not look for
discontinuities. This powerful tool can be very useful to filter correlated noise like acquisition
footprint or to isolate specific fault systems within the reservoir for detailed structural studies.
You can access the Stereonet tool by selecting the “Stereonet” Tab that is available when the
Ant-tracking Attribute is active. Otherwise this tab option is grayed out.
How does the Stereonet work? (After Stein-Inge Pedersen et. al, Stavanger Research)
The surface orientation is represented by its normal given in spherical coordinates. Surfaces will
be close to planar, so the normal is a good representation of its orientation. On the Stereonet the
normal is plotted as a pole. The azimuth [0 – 360] is read along the perimeter of the plot, and the
dip [0 –90] is read as the distance from the centre.
A horizontal surface will be plotted in the center of the plot. Two surfaces having the same
azimuth, but conjugate dips will be plotted symmetrically around the center of the stereonet.
From structural geology, it is known that faults appear in systems, that is, they have the same
strike (azimuth) and conjugate dips. The same stress field, indicating that faults have been
created in the same time period, has created a fault system. On the polar plot, fault systems can
be identified as clusters symmetrical around the center.
Seismic – Ant-tracking
Edit Extracted Faults process
1. Activate the Edit Extracted Faults process under Geophysics in the Process diagram
2. Double click on the Edit Extracted Faults process to get information about the Automatic Fault
Extraction Process
2
1
The Edit Extracted Faults process provides a set of interactive tools to display, analyze and edit
extracted fault-patches. Fault patches must have been created first by the Ant Tracking attribute.
Fault patches can be filtered based on properties. Initially, all the extracted surfaces are
displayed. This gives the interpreter an overview of the possible fault systems in the data. The
interpreter will then, using the surface orientations, split the surfaces into their corresponding
systems. Each of these systems can then be studied and validated separately. This represents
an interaction with the result at a fault system level. If the interpreter finds that a system is due to
e.g. acquisition footprints, or other coherent noise, the system is simply deleted. If the system is
real, the interpreter will then verify, and if necessary, edit the individual surfaces within the
system. This corresponds to working on the surface level.
Seismic – Ant-tracking
Edit Extracted Faults process
1. Activate the Edit Extracted 1
Faults process under
Geophysics in the Process
diagram
Selected faults will turn gray; multiple faults can be selected by holding down the Ctrl key while
selecting. Selected faults will also be annotated by red dots in the stereonet.
1 2 3
1
3
2
The settings for the conversion of the fault patches can be adjusted on the settings tab of the
Edit Extracted Faults process. The number of fault sticks, the sampling, the height above and
below given min/max points and the number of fault pillars can be defined.
3
1
2
1
4
3
6 5
Neural net evaluation is a good tool for automatically finding relationships between multiple
known parameters and a single potentially unknown parameter. Exemplified this can be used to
identifying the relationship between facies and other measured logs or between unknown
facies/attribute classes and generated seismic attributes. Once found the relationship can be
used in any position (wells or 3D grids) where the input logs or attributes are present to make an
estimate of the unknown parameter. This can be done in two ways:
Untrained analysis
The user provides the attribute cubes to use in the analysis and chooses how many classes to
create (X). The algorithm plots all the input data together and tries to split the data points into X
different seismic classes which each show similar relationship.
Trained analysis
The user provides an interpreted attribute at a location (wells) and the associated attribute cubes
for that data. The algorithm tries to recognize patterns that separate the different assigned
classes. Trained analysis can also be used to estimate continuous data.
Once the relationship has been modeled it can be applied to areas where only the input data is
present.
7 9
10
Unsupervised – Specify the desired number of classes. The algorithm will separate the data as
logically as possible into this number of classes (classification only).
Supervised – Specify the training data to be used in the classification or the target for the
estimation. This data must be present at all of the points used as input or on the points used as
Training Data if that option was selected (seismic and surfaces only). If your desired training
data is not available in this drop down list then check that your selected input data has the
training data.
Next specify the training options:
Max Number of Iterations – The algorithm will stop at this number even if an adequate result
has not been reached.
Error Limit (%) – When the number of points classified incorrectly is below this limit, model is
assumed to be trained and will stop.
Cross Validation (%) (supervised only) – This is the percentage of the input data which is used
to test the result and give the error. The remaining part is used to train the model.
Probability Threshold (%) (classification only) – The algorithm finds the probability of each
class being found at the point, then assigns the most probable class to the point. If the
probability for the most probable class is below this value then the point will be left undefined.
Exercise overview:
• Introduction to Ant-tracking
• Generate the Structural Smoothing attribute cube
• Generate Chaos attribute cube as input for Ant-tracking attribute cube
• Generate the Ant-tracking cube and fault patches
• Settings for fault patches
• Edit Extracted Faults process
• Convert fault patches to fault interpretation
• Train Estimation Model
• Bricked seismic
Introduction to Ant-tracking
This unique algorithm is part of an innovative workflow that introduces a new paradigm in fault
interpretation. It emulates the behaviour of ant colonies in nature and how they use pheromones
to mark their paths in order to optimize the search for food. Similarly, virtual ants are put as
seeds on a seismic discontinuity volume to look for fault zones. Virtual pheromones deployed by
the ants capture information related to the fault zones in the volume. The result is an attribute
volume that shows fault zones very sharp and detailed. A second optional output is the
automated extraction of “fault patches” out of the volume.
The input for the ant attribute is an edge enhanced volume, like variance or chaos.
Ant-tracking is a cpu-intensive process. It is therefore advisable to use a small-cropped volume
first to test the parameters before attempting to calculate larger areas. As a rule of thumb,
account for 20 min. of calculation for a 100 Mbyte volume using a machine with 2 Gbyte of RAM
and 1.7 GHz of cpu-clock.
The time it takes to generate a volume will heavily depend also on the amount of discontinuities
that are present in the input volume and on the mode you are running the Ant-tracking.
If the Generate Fault Patches option is active, calculation time will essentially double due to the
patch generation and merging operations.
Ant-tracking can only be run in Realized mode.
3. Click Apply.
4. Toggle on an Inline and Crossline and quality check the result. Change the parameters and
Apply to view the effect on the displayed lines.
5. Specify the other settings the same as in the figure above. If there is a name in the Output
seismic box, click on it and delete it using the delete key on the keyboard. Toggle on Set
output name and specify the name. Click OK.
The Structural Smoothing attribute cube has been generated and will be found in the Input tab in
the Explorer window. This attribute volume can be Realized to create a physical version on disk,
put is not needed.
Exercise Steps
1. Right click on the Structural smoothing attribute cube and click on Attribute Generation.
2. Choose Chaos under the Stratigraphic methods in the Attribute Library button. Toggle
on Virtual to render the attribute on-the-fly.
3. Click Apply.
4. Toggle on an Inline and Crossline and quality check the result. This volume can be Realized
if preferred, but is not necessary.
5. Set up the other settings the same as in the figure above. If there is a name in the Output
seismic box, click on it and delete it using the delete key on the keyboard. Toggle on Set
output name and specify the name. Click OK.
Exercise Steps
1. Right click on the attribute cube (Chaos) and click on Attribute Generation using the LMB.
2. Choose AntTracking under the Structural methods in the Attribute Library button. Ant-
tracking can only be run in realized mode.
4. Click on the Show Fault Patch Setting Options box . Leave the settings as
default.
5. Click on the Show Ant Tracking Settting Options box to come back to the Ant Tracking
settings
6. Click on the Passive Ants and leave the AntTracking settings as default.
7. Toggle on Set output name and specify the name (AntTracking).
8. Click OK to generate the AntTracking attribute cube and the Fault Patches.
Both the AntTracking attribute cube and the Fault patches will be found in the Input tab in the
Explorer window.
Exercise Steps
1. Right click on the Fault Patches 1 and open the Settings.
2. Go to the Histogram tab and filter on the Fault Patch Size. Click with the LMB on the first
histogram bar representing the smallest faults. The color will change from blue to pink to
indicate the selection. In this way the small faults will not be displayed.
3. Click Apply.
4. Go to the Stereonet tab and filter on the Fault Patch dip and azimuth.
The Stereonet tool can be used to filter a priori area where the Ant agents should not look for
discontinuities. This powerful tool can be very useful to filter correlated noise like acquisition
footprint or to isolate specific fault systems within the reservoir for detailed structural studies.
The circles represent the dip. The inner black circle represent faults with less than 15 degrees
dip and they are filtered away. The next circle represents the faults with less than 30 degrees
etc. if the dip increment is set to 15. It is possible to filter away the faults with certain dips and
azimuths by clicking on the sectors with using the left mouse button as shown in the figure
below.
Exercise Steps
1. Activate the Edit Extracted Faults process under Geophysics in the Process diagram
2. A set of tools will be available in the Function bar. Activate the Set select/Pick mode
icon and select a fault that needs to be smoothed.
3. Click on the Smooth Selected Patches icon and the selected fault is smoothed. If the
fault patch need to be even more smooth click on the Smooth Selected Patches icon
ones more.
4. To merge faults click on the Set select/Pick mode icon and select the two faults to be
merged. Click on the Merge Selected Patches icon and the selected faults will be
merged. The merged fault needs often to be smoothed.
5. Convert all the fault patches to Fault Interpretation by selecting all the Fault Patches in the
3D Window using the Ctrl-A key combination on the keyboard.
Comments
To convert the Fault Patches directly to Fault Pillars select one or several Fault Patches and
click on the Convert to Fault in Fault Model icon . The selected fault patches will be
converted to Fault Pillars and added to the active Fault model. It is displayed automatically. The
settings for the conversion of the fault patches can be adjusted on the settings tab of the Edit
Extracted Faults process. The number of fault sticks, the sampling, the height above and below
given min/max points and the number of fault pillars can be defined.
Exercise Steps
1. Open the Attribute Generation process found under Geophysics. Use mig.sgy [Crop]
[Realized] as Input Seismic. Use Graphic Equalizer to band-pass filter all frequencies
except 10, 20 and 30 Hz.
Comments
Only virtual attributes that are associated with the parent cube can be used in the estimation
model.
Bricked seismic
Any seismic volume, either loaded into or generated in Petrel can easily be exported to the hard
disk for further use. Any cropped, realized cube or attribute cube will be exported using the
SEGY format or the ZGY bricked format.
To save disk space and time, use your knowledge of cropping down the already cropped and
realized cube to less than 10 MB (use statistics to verify) before export.
Exercise Steps
1. Right click on the seismic volume of choice and go to Export in the appearing menu.
2. Specify which directory (i.e. Student) and give in an appropriate name. Select Save as type
to be SEG-Y seismic data. A progress bar in the lower right corner of Petrel tells how much
of the process is performed.
3. In the student folder, find your SEG-Y file, right click on it and open Properties. Observe the
size of the file.
4. In Petrel, import the newly generated SEG-Y file as described previously. Rename it to
mig_small.sgy and look at Statistics, the SEG-Y Settings tab and check the ASCII header
information.
5. Open the Operations tab. Toggle on ‘Optimize for 3d Visualization’, set Alignment and
Resolution at your preference. Realize the volume.
6. Give in an appropriate name, use *.zgy as format and save it to disk. Click Yes in the pop
up box.
7. In the Input tab, find the new ZGY object. Observe that by default is Inline Crossline and Z
available (it resides in RAM by default).
8. Open Settings for the new object, under 3D Style toggle on Volume render. Go to the Colors
tab and set the opacity curve to only display extreme amplitude values.
Comments
Note that 2D seismic lines can only be exported in SEG-Y format.
Exercise overview:
• Create a simple grid
Exercise Steps
1. Open the Make Simple Grid process found under Utilities.
2. Make sure Create new is selected and give in an appropriate name.
3. With Insert horizons toggled on in the Input data tab, find and drop in the surfaces found in
the Surfaces folder.
4. Open the Geometry tab, click on the HOUSTON surface in the Input tab, then on Get limits
from selected.
5. Change Grid increment to 250 X 250 and click OK.
6. Find and display horizons and edges form the newly created 3D Grid.
Comments
The simple grids created this way could be used in velocity modeling, sampling seismic data into
the grid to create a seismic property or to have regular grids for upscaling.
Velocity model
Time domain Depth domain
1. In the Process tab, open the Geophysics folder and double-click on the Make Velocity Model
2. Click the ”Append item in the table” icon as many times as available input objects
3. Change Bottom Interval to desired type.
4. Select the levels to model, from Input or Models tab, click on its name and drop it into the dialog box
using the drop in arrow .
5. Keep the Correction to None
2
1
3 5
4
1. Select type of velocity model for each interval. When using Linvel function for
modeling, constants, surfaces or Well TDR can be used as input
2. Set V0 to Well TDR – Surface/ Constant
3. Set k to Well TDR – Surface/ Constant
3
2
This exercise uses the established Time/Depth relationship in the different wells to derive V0
and k points for the specific intervals (each well with a T/D relation gives V0/k values for each
interval). These points are gridded or averaged to make surfaces or constant values for the
intervals used in creating the velocity model.
When creating velocity models, the results can be output in different formats. If a correction is
performed when modeling, both original and final velocity models as well as the correction can
be output.
3 1
2 2. Possible to set MD increment, depth
and time tolerance and clipping
parameters for V0 and k
Output settings
Gridding:- The chosen gridding algorithm will be used for creating all the surfaces made during
the velocity modeling process. For more control over this gridding e.g. using different methods
for different zones, create the output as points, grid them using the make surface process and
input surfaces instead of the TDR.
Resolution:- Output from the velocity model (SEGY, surfaces and logs) will be created at the
resolutions specified here.
1. Open Settings for Wells and go to the Time tab. Select Checkshots, Well Tops or
corrected sonic logs as source for time/depth (move the prefered source to the top)
2. Keep the parameters for the velocity model as already set up
3. Apply or OK to generate V0 and k parameters
As previously stated, an alternative to use Well Tops as surce for time/depth relation in all wells
is to substitute with Checkshots, velocity logs or sonic logs for each well.
This is done in Settings for each well, Time tab where Override global settings is toggled on and
the correct T/D source is selected and put at the top of the list using the up arrow.
1 3
4
5
2
The next example shows how to use velocity data in the form of a SEG-Y velocity cube in the
modeling process. A velocity cube must be of type Average velocity (i.e. each point contains the
velocity from MSL to that point). The input volume is converted to have the correct type of data
for velocity modeling. This approach is not covered in the exercises section.
3 4
This example is different from the previous one since well data is not used to make the velocity
model, but Well Tops are used for correcting the final model.
2
1
4
2
The next example shows how to use Stacking velocity data points to model velocities.
1
3
5
4
2
1 3 4
3
1
Using Relative (%) for Min/Max sets the output to be ralitive to the input in %.An alternative
exists where Absolute is used. The data range needs to be estimated by clicking the Estimate
button. All zones must be unlocked and output data range estimated for each.
5 3
1
3
5
1 3
6
7 8
5
9
EXERCISE
Domain Conversion
Commonly called depth conversion, or more accurately described, conversion between domains
is a two-step process in Petrel. Velocity modeling is a pre-requisite for and is separated from
Domain conversion. The velocity model is defined by surfaces or constants, well time-depth
relations, velocity cubes, 3D grids or other velocity information available. The same model can
be used to convert from time to depth or depth to time. It is also possible to convert within the
same domain (i.e. time to time or depth to depth). This is done for correction purposes.
After building a velocity model, objects can be converted independent of 3D grids. It is then
possible to interactively switch between time and depth domains for that object. If more than one
velocity model exists, several depth versions of the same object can co-exist. Domain
conversion is integrated with and can be performed within the process manager in Petrel.
In this section we will create three different velocity models that later can be used for domain
conversion.
Exercise overview:
• Define velocity model intervals
• Create a velocity model using well data
• Create a velocity model using stacking velocities
• Depth convert a 3D grid
• General depth conversion
Exercise Steps
1. From the Process tab, open the Make Velocity Model process found under Geophysics.
2. Click the Append item in the table icon once. Change Bottom interval type to Surface,
leave Correction to None and V=V0 for Model type.
The V0 and k output can be found under the Velocities [Final] folder for the model. These
objects can be used for quality control of the velocity model as well. The time and velocity logs
are found under the Wells folder under the Input tab.
Exercise Steps
1. Go to the Vstack_points found under Vstack folder. Expand it together with the Attributes
sub-folder. Display the data set in a 3D window. Open Settings for the Vstack_points, then
go to the Operations tab.
2. Expand the Eliminate Where folder and select the function Z < Constant. Type in –3600
(remember to use minus) and click on Execute. Observe the displayed data points in the
3D window.
5. Apply or OK will assign values to the grid cells containing any point data. After completion,
open the Properties folder and display Average Velocity [U] in a 3D Window.
6. The next step is to populate the upscaled values out into the total 3D grid. From the
Process tab, open the Petrophysical Modeling process found under Property Modeling.
Make sure the Use existing property is toggled on and use Average velocity [U].
7. For the first zone, unlock it with the lock icon .
8. Select Moving Average (Interpolation) as Method for zone/facies. For Output data range,
keep Min and Max to 0 (zero) and Relative(%).
9. Copy settings from the selected zone , and then Paste settings to all zones . In this
way, all zones will be populated using the same settings.
10. Apply the settings and view the result in a 3D window.
14. Make sure that all toggles are set for Well report under the Output tab.
15. Apply the new settings and view the result in a 3D window after completion.
Exercise Steps
1. Go the Models tab and activate Coarse 3D Grid for VelMod (bold text). In the Process
tab, open Depth Convert 3D Grid.
2. From the Velocity model drop-down menu, select the Velocity Model 1 to use.
3. Alternatively reduce pillar geometry complexity by toggling off geometries for both faulted
and non faulted pillars (this 3D grid has simple geometries, so this option will have limited
effect if any).
4. Apply the settings.
5. A new Coarse 3D Grid for VelMod (DC) will appear in the Models tab. View and check it
in a 3D window.
Exercise Steps
1. First we will depth convert the cropped and realized seismic volume used throughout this
course. In the Models tab, expand the Velocity Models folder and activate Velocity
Model 2.
2. Go to the Input tab, right click on mig.sgy [Crop] [Realized] and use Depth Convert by
active Velocity Model from the menu.
8. Click OK. Open a new 3D Window and display all objects used in the previous step.
9. Change between TWT and TVD with the Change time/depth settings for the window
drop-down menu in the Tool bar. Observe the differences. Also note that these object now
have an Attributes sub-folder available. Expand one of the depth-converted objects and
observe what information there is in the Attributes sub-folder.
Comments
Another way of domain converting (i.e. DC) a seismic volume is to attach it to a 3D grid depth
converted inside Petrel. The conversion from time to depth for the 3D grid is stored with it and
any attached seismic volume will be stretched to fit the grid. The result is a seismic volume that
is defined inside the 3D grid and undefined outside the grid interval.