Navigation lines, shot point numbers and stacking velocities

The Navigation plug-in Geocap handles seismic navigation lines, stacking velocities and depth conversion. The plug-in is part of the Oil & Gas module.

The following describes the command syntax for the functionality collected in this plug-in.
Together with general mapping and modeling functionality in Geocap it governs a flexible platform for velocity analysis and horizon based depth conversion.

The functionality can be operated through scripts, commands or a user dedicated menu.



On this page:

Navigation lines are also called stacking velocities or shot point lines and will be used interchangeably.

The navigation lines go in one direction. Each navigation line contains many different shot point numbers.
Shot points of same number go in the other direction. Underneath each shot point sits the stacking velocities in a gather.


Shot points and navigation lines


Detail of navigation lines and id and some shot point numbers

The following ASCII shot point file represents an example:

        Line name     Shot point  x coord   y coord    time   velocity
        3359I          1280     530196.6  6757992.7   548    1517
        3359I          1280     530196.6  6757992.7   608    1566
        3359I          1280     530196.6  6757992.7   664    1585
        .  .  .  .  .  .  .  .  .
        .  .  .  .  .  .  .  .  .
        3359I          1280     530196.6  6757992.7  3774    3318
        3359I          1280     530196.6  6757992.7  4014    3464
        3359I          1280     530196.6  6757992.7  4986    3961
        3359I          1320     530064.0  6758475.0     0    1470
        3359I          1320     530064.0  6758475.0   475    1470
        3359I          1320     530064.0  6758475.0   562    1533
        3359I          1320     530064.0  6758475.0   615    1573
        .  .  .  .  .  .  .  .  .
        .  .  .  .  .  .  .  .  .
        3439I          1280     531160.9  6758258.8     0    1470
        3439I          1280     531160.9  6758258.8   472    1470
        3439I          1280     531160.9  6758258.8   519    1476
        3439I          1280     531160.9  6758258.8   587    1525
        .  .  .  .  .  .  .  .  .
        .  .  .  .  .  .  .  .  .

This special format is called a navigation velocity stacking format in Geocap and the whole file is converted to a Geocap vtk file and placed in the project.

This file is read into a project folder by Import->Stacking Velocities. Shot point are in field data SP and Line Name are in field data 2Dline.

Check the import afterwards in Table View to see that the field data columns are properly assigned.

The navigation file is converted to a cell structure making a new cell whenever the shot point number changes.

The command syntax for manipulation this dataset is placed under the command nav for navigation, because this kind of file traditionally belongs to a navigation file.

All functionality connected to the nav command is placed in the plug-in called Navigation. This plug-in is included in the Oil & Gas module and must be loaded into Geocap in order to operate the functionality.

The following commands require that the stacking velocity dataset is the active data.

nav mak cells ;# make cells of stacking velocity. Each vertical stacking column is a cell. Should be performed if the stacking velocity file has no cell structure; i.e. contains only one cell.

nav mak shotpoints ;# make one line set of all the first top shot points.

nav mak timevel selshot cur | _shotpoint_number ;# select time and velocity data for shot point closest to the cursor position (cur) or for a certain shot point number shotpoint_number. Produces also an additional dataset navcross in workspace which can be triangulated by nav gen trg.

nav mak timevel cur ;# make time and velocity column at the cursor position.

The situation is that stacking velocities are present as a SEGY file which is read into a seismic folder as usual. The task is to convert the stacking velocities in the seismic format to the standard stacking velocity format.

The following command require that the outline is in workspace outline and the amplitude grid in the SEGY file is in active.

nav lin [resolution_in_xresolution_in_z] [xstart x zstart z] [name SEGY_file_name ]

resolution_in_x

The sample resolution in x direction is resolution_in_x value. Similar for depth or z direction.

xstart x

Specifies where to start the sampling in x direction. Similar in z direction.

name SEGY_file_name

Specifies the name in the field data 2Dline. This should normally be the name of the SEGY dataset.

Example of converting SEGY stacking velocity format to standard navigation stacking velocity format

nav lin ;# makes vertical lines for all nodes (grid in active, outline in outline)
nav lin 1 4 ;# makes vertical lines of sample 1 in length direction (x) and 4 in depth (z).
# The result is placed in workspace seislines.

The following commands require that the stacking velocity dataset is the active data.

nav dis shotpoint | shotline | shotlinepoint cur [ref:arguments];# display shot line and/or point values together at cursor position.

arguments:
col color_name – specify color.
sel selection – select every selection.
size size_value – size of numbers.

The following commands require that the stacking velocity dataset is the active data.

nav cut z_level [nofd] [ver] [ext] ;# cut the navigation file at a certain level

The navigation file will be cutted a level z_level. The new created dataset will be present in active after operation.

Default is that present field data SP and 2Dline are generated in the new dataset and line structure is kept.

arguments:

z_level - value of cut level. Must be present directly after nav cut.
nofd – no field data is generated, but line structure is kept if field data SP and 2Dline are present in the navigation file.
ver – interpolation is done only if the line to be cut is completely vertical. This will prevent cutting nonvertical lines if the cell structure is incorrect. A correct cell structure has one vertical line per cell.
ext – extrapolation is performed when necessary at start and stop of the line if the cut level is outside the line. In that case a curve called extrapolationcurve shall be present in workspace. The extrapolation curve has entry on X and output on Y coordinates.

The following command require that the selected velocity stack along a line is active data.

nav gen trg ;# generate triangle model from selected velocity stack data along a line.

An alternative is to run a prestep that makes an equal number of points for adjacent shot point gathers. This option is safer because it may work with a high resolution in all cases.

mak bandsca fix ;# run prestep and generate triangle model from selected velocity stack data along a line.

After this command one gets a solid vertical section that can be mapped and contoured.

The stacking velocities are active data.

nav gen cube stack | cel nx ny nz ;# generate 3d cube model nx ny nz from stacking velocities using time slice method and gridding.

stack – use original stacking velocities.
cel – use stacking velocities transferred to cells.

nav set hicoreid – sets an id that can be used in other nav commands.

It can work like this:

set job <something>

nav set hicoreid $job

and then, later on

nav cal dix stack will need the following data sets to be present in workspace:

  • topTimeVelocity$job
  • bottonTimeVelocity$job
  • topTimeSurface$job
  • bottonTimeSurface$job

This is a common feature in all nav commands which use hicore (workspace).

Preset hicore names that can be adjusted by a hicoreid

            topTimeSurface

            topTimeVelocity

            topTimeVelocityStack

            bottomTimeSurface

            bottomTimeVelocity

            bottomTimeVelocityStack

            stackingVelocityData

            stackingVelocityCube

            intervalVelocityDix

            intervalVelocityStackDix

            highSurface

The velocity cube is active data.

nav acu wells_in_workspace ;# adjust cube towards check shot wells in wells_in_workspace. The format of the check shot well is: x y depth time. The scalar part is in fact time_in_milliseconds.

The active data is a time grid in milliseconds to be depth converted. The general formula for depth conversion is Depth = Time x Velocity . The time grid is probed with values from the velocity cube to get the velocity at the time grid positions.

nav depcon velocity_cube ;# perform depth conversion of time grid using velocity_cube in workspace. The result depth grid will be the new active data and is also saved as workspace depthgrid.


Velocity cube used to depth convert time grid to depth grid

The stacking velocities are active data.

nav gen cube stack | cel nx ny nz ;# generate 3d cube model nx ny nz from stacking velocities using time slice method and gridding.

stack – use original stacking velocities.
cel – use stacking velocities transferred to cells.

The time surface is active data.

nav cal velocity stack ;# calculating velocity data for time surface in workspace. The navigation stacking file is in workspace: stackingVelocity

nav cal dix stack ;# calculate interval velocity after Dix formula.

The following datasets must be present in workspace under exact these names:

topTimeVelocity – polydata derived by nav cal velocity stack.
bottomTimeVelocity - polydata derived by nav cal velocity stack.
topTimeSurface – top time surface grid.
bottomTimeSurface – bottom time surface grid.

The stacking velocities are active data.

nav cal int cel value ;# calculate interpolated velocity at value

nav cal int surface ;# calculate interpolated velocity for surface in workspace highSurface

nav cal ave [mid | bot] ;# calculate average depth velocities.

Default will use interval value above calculation point. mid: use middle interval. bot: use interval below. The result is saved in depthAverageValues[Top | Mid | Bot] .



The stacking velocities are active data.

nav select cel select [ref:start start_number*] ;# select every select vertical stacking cell and optionally start at _start_number.

nav sto dix filename [ref:topbottomvelocity] ;# store interval velocity data after having it calculated so that everything is in proper place in workspace. 6 columns are written out by default. If topbottomvelocity is specified, two extra columns (7 and 8) are written out.

Standard column output:

          Shotline  Shotpoint   x           y         time     int.Velocity
          3359I     2360    526618.38   6771010.5   2319.050   2389.451
          3359I     2284    526751.00   6770528.5   2284.787   2363.030
          3359I     2279    526883.50   6770046.0   2279.813   2399.876
          .  .  .  .  .  .  .  .  .
          .  .  .  .  .  .  .  .  .

The stacking velocities are active data.

nav set diagonalnumber diagonalnumber ;# set number for division of diagonals. To be used for very dense stacking velocity datasets. Default number is 700. This may be increased and represents a temporary grid resolution when calculating Dix interval velocity.

nav set hicoreid hicoreid ;# all fixed workspace names can have an additional id tag attached to the name making it possible for several projects with different datasets to operate simultaneously.

These scripted examples can be used as command objects with the project manager or as script commands activated in the shell.

Stacking velocity dataset is the active data.

          nav dis shotpoints sel 1 size .7  ; # display shot points
        


Detail of shot point numbers

Stacking velocity dataset is the active data.

          nav dis shotlines size 1  ; # display shot line numbers
          col bla  ; # black color
          poll ; # draw cell lines


Detail of shotlines and id numbers

Stacking velocity dataset is the active data.

          nav mak timevel selline cur  ; # make time and velocity section at cursor position
          mhi timevel ; # save data
          mlo navcross ; # get saved navcross data
          nav gen trg ; # generate triangulate model
          map ; # map time and velocity

The velocity cube itself is displayed with an opacity: dis opa .4


Selected stacking velocities along a line and triangulated into a vertical section


Composite display of navigation lines, stacking velocity sections and time grids


Detail of navigation lines and stacking velocities


Stacking velocities mapped and navigation lines

These menus are using scripted syntax described in this reference. They can be activated directly on a stacking velocity dataset in the project.

Stacking velocity dataset is the active data.

          nav dis shotpoints sel 1 size .7  ; # display shot points
        


Menu for mapping stacking velocities


Menu for displaying shot lines and shot points