Statics calculation that based on the refractions methods requires first breaks as mandatory input data set. Therefore, we need to do a picking of the first arrivals (guide-wave) on the input seismic traces. For this task g-Platform has Refraction FB picking - azimuthal solver guide / phase / aperture module. Next, we need to calculate statics corrections by using regression method solving or tomography algorithm (Tomo statics 2D module).

Create a new workflow 0070_Refraction_statics:

In this chapter we will pick first breaks, create velocity model via tomography algorithm and calculate (refraction) static corrections.

Notice, that tomography procedure requires high quality picks, so it is better to have accurate manual picking on the 2D line (manual picking if you have enough time). The seismic is vibro data and it is another obstacle of using auto picking method as well as trying to use tomography procedure, the result may be far from being sufficiently addressed in comparison with the conventional reciprocal method static solution.

Let’s look at the workflow with all necessary modules for static calculations:

1. Read seismic traces - loading seismic data for picking

2. Refraction FB picking – azimuthal solver guide / phase / aperture – manual or auto picking of the first breaks

3. Tomo statics 2D – tomography statics calculations

4. Static QC – quality control for statics from Tomo statics 2D

5. Sort traces - sort traces by CMP for stacking

6. Stack imaging - create a stack, make QC of refraction statics

1) Read seismic traces. For the first break picking step we should use seismic without any processing in order to have original first response from guide-waves. Load seismic data set 0001_Geometry.

2) Refraction FB picking - azimuthal solver guide / phase / aperture. This module computes the static correction based on the first breaks of the data. It allows the user to interactively pick or correct the automatically picked first breaks. Subsequently, the first breaks can be used for refraction static solution, generating static corrections for source and receiver positions. The module calculates only the refraction correction, that is, replacing the near surface low velocity layers from the surface to the bottom of the calculated near surface. The module does not include an elevation correction. Elevation static is a shift which moves the data from topography to final datum. This correction to final datum is performed with a separate module Shift Data. Since most processes in g-Platform that would require the statics applied (for example NMO) can be run from topography, being able to easily separate the refraction correction from the final datum correction is very useful. How to use floating and constant datum planes will be discussed in Datum planes chapter.

Refraction FB picking uses first breaks and  parameters for the solver to calculate the source and receiver static correction at each surface location. The solver includes options to calculate a smoothed solution, as well as to solve into multiple azimuths. Additionally, a residual correction can be calculated and applied at each source and receiver giving a short-wave correction. The replacement velocity is calculated internally by the module. Upcoming changes will allow specify a weathering velocity as well as a replacement velocity. With these parameters the module will build a near surface weathering model and then the user defined values will be used for static corrections calculation.

Fill the Refraction FB picking – azimuthal solver guide / phase / aperture parameters before performing calculation and displaying the data set. For the training seismic data set we will use the following example, but you should do some extra tests by changing parameters for better understanding:

Input data:

Parameters:

Let’s have a look at the main parameters we need for picking:

Load RCV gather

By default, it is FALSE. If checked, it will display the receiver gathers in the current receiver gather window.

Build offset step

Provide the offset step size. Depending on the offset step size, it will build the offset classes that can be chosen from the Offsets option parameter.

Offsets

Select the desired offset from the available drop down menu.

Number Azimuthal Direction

By default 1. In case of 3D, the user can define the azimuths.

Pick magnet phase

Provide the phase to pick the first breaks.

Min offset

Provide the minimum offset value for picking the first breaks.

Max offset

The user can limit the minimum and/or maximum offset values. Provide the maximum offset value for picking the first breaks.

Step offset

Offset step size. Create an offset class as a group of offsets for velocity estimation.

Sliding step

Provide the sliding value to jump to next source/receiver gather.

Super gather picking params:

Improve picks  – Set picks more accurate on phase if it possible.

Aperture

Super Gather Aperture, the aperture of locations that will be picked on either side of the current pick.

Min velocity

Minimum velocity range for super gather picking.

Max velocity

Maximum velocity range for super gather picking.

Step velocity

Step velocity between minimum and maximum to be used for super gather picking.

Run local solver to improve picking

Leave checked to use the results of the local refraction statics solver to improve the results of your super gather picking.

Improve picks

Set picks to be more accurate on the phase if possible.

Manual picking params:

Manual picking window

Search window used for manual first break picks.

Automatic picking window

Search window used for in automatic super gather picking mode.

Threshold

Threshold for manual picking.

Magnet type

Phase of first breaks to be picked.

Find magnet on click

Type of energy for first breaks.

Check lines on manual picking

By default TRUE.

Solve guide line picking params:

Aperture for picking around existed picks

Provide the aperture value for picking first breaks using the solve guide line with existing picks.

Magnet to phase

Choose the phase for picking the first breaks. By default, FALSE. If checked, the next two parameters Window and Threshold will be activated.

Window

Provide the picking window.

Threshold

Define the threshold. By default 3.

Magnet type in window

By default, Nearest.

Automatic picking params

Window

Search window for automatic first break picking.

Threshold

Provide the picking threshold value for automatic picking.

Min offset

Provide the minimum offset value for first break picking.

Max offset

Provide the maximum offset value for first break picking.

Magnet type in window

There are multiple magnet types are available for first breaking picking. By default, Nearest.

Regression correction

By default, FALSE. If checked, next two parameters will be activated.

Iterations to improve

Provide the number of iterations to improve the regression correction.

Radius of searching picks for regression

Provide the search radius.

LMO params

LMO type

Type of LMO to be applied to shot and receiver gathers for picking:

None - indicates no LMO will be applied;

V const - is a constant velocity LMO as indicated by the V LMO parameter;

Picking - will apply an LMO based on the bin gather picks.

Further step is setting up user work area, it means making configuration of all necessary Vista Windows for QC and interactive engagement. Click RMB on the module and open all Vista Groups:

On the first Vista Window, there is a Linear Move Out (LMO) picking on source gather. The LMO is one of the main parameters for FB picking, in view of the fact that it is used as a guide function for auto FB picking. Consequently, we should pick LMO guide function on a few source gathers for different locations through the acquisition area. Picking may be located exactly on the first arrivals or a little bit above them, it depends on what option (algorithm) and parameters you use for auto picking. For the training data it is enough to make LMO picking only for one source due to the fact the there is a flat relief.

All the visual setting were configured and we can propagate and save them as default:  

Next, open the source gather FB picking window and execute one source FB picking by pressing Auto picking current SRC<Alt+,S> button:

Then we are going to launch auto picking for the entire data set by clicking on Auto picking action button:

Check the result, it is obvious that picks became more accurate because the algorithm uses neighbor gathers for statistics:

Eventually we have all FB picks and need to check it, open the Map sources and receiver location and activate Selected source option:

Now we can choose any source from the location map and perform FB picks QC and editing. Other windows are connected to each other automatically, so it allows to perform interactive QC on fly:

Or choose another option for picks removing:

There are other QC windows available from the Refraction FB picking - azimuthal solver guide / phase / aperture module.

The final step is export FB picking to binary file that we will use for tomography process. Find Export picking action and choose Save in Binary, because tomo statics module requires binary format:

Define a name of the output binary file FBpick, we will use this file in the next step - tomography.

Now we have all FB picks and there two options for static calculation, by using the following modules:

• Refraction FB picking – azimuthal solver guide / phase / aperture;

•Tomo statics 2D.

Execute Refraction FB picking – azimuthal solver guide / phase / aperture, now we have refraction static correction and it is easy to do a QC. Go to vista window Current source gather and do a quick statics QC:

It was a fast QC of static corrections calculated by Refraction FB picking – azimuthal solver guide / phase / aperture. Of course we are going to calculate stacks, but before we do it there is another option for static calculation by using tomography - module Tomo statics 2D.

3) Tomo statics 2D. Tomography plays a key role in seismic imaging. Be it is a building depth velocity model or calculating the refraction statics by generating the near surface model. Refraction Tomography statics are used in the event of regular refraction statics are not giving optimum results or to improve the refraction statics solution. In this method, the weathering zone (WZ) is parametrized as a number of cells that used for filling up a depth velocity model. Travel times are goes through the cell-model, and residuals times are converted into velocity variations in 2D/3D cells. Here we are able to calculate vertical velocity gradient due to the fact that tomography is a nonlinear move out modeling of the first arrivals. 

The tomographic is sensitive to the initial model due to the fact that the nonlinear solution performs by iterating the ray trace calculation and velocity is updated by local linearization, so the solution of tomography is depending on quality of first breaks picking. Therefore, if we have 2D data set it is more reasonable to make FB picking more accurate by doing it manually.

Linear (left) and nonlinear (right) solutions comparison:

The tomo-statics modules in the g-Platform system are based on the Fresnel seismic travel time tomography. A Fresnel volume approach is applied to represent wave propagation for seismic travel time tomography instead of rays. A Fresnel volume is defined as a set of many waves delayed after the shortest travel time by less than half a period. It is derived by calculating travel times both from a source and from a receiver. Tracing rays from sources to receivers is completely avoided.This considerably reduces computational time. We solved the eikonal equation by using a finite-difference method to calculate travel times. The advantage of this approach is as follows: first, the frequency of wave can be introduced into analysis. Therefore, we can evaluate the resolution of seismic tomography. Next, the smoothing feature can be naturally introduced. Finally, Fresnel volumes with finite bandwidth considerably reduces the sparseness of ray distribution.

The more physically-realistic representation of wave propagation is to treat a ray path as a beam with finite width. Using a Fresnel volume is a natural and an effective approach. A Fresnel volume is a set of many rays delayed after the shortest travel time by less than half the period of wave. The rays in a Fresnel volume are added constructively to form the first-arrival of wave. There have been several studies on the application of Fresnel volumes to seismic tomography since Harlan(1990). In this study, first, we discussed the characteristics of Fresnel volumes. Next, we formulated the inversion procedure. Then, we investigated the resolution of tomography with respect to the frequency (Toshiki W. Seismic travel time tomography using Fresnel volume approach).

Travel times scheme:

A module Tomo statics 2D is used for tomography and it requires the first break picks. In this module, it start with an initial model and performs updating the model during each iteration until it converges to the best solution.

Scheme of tomography steps:

Demo data also consists of the FB picking file on your disk in the same folder with the input seismic data set, but you just need to convert ASCII FB pick file into Binary one via Refraction FB picking - azimuthal solver guide / phase / aperture module (import->export). So, we can use the module Tomo statics 2D which requires first breaks as mandatory input data set.

Tomography procedure requires high quality picks, consequently preferable to have manual picking on the 2D line due to relatively small amount of source gathers. The seismic is vibro data and it is another obstacle of using auto picking method as well as trying to use tomography procedure, the result may be far from being sufficiently addressed in comparison with the conventional reciprocal method static solution.

For tomography procedure we will use Tomo statics 2D module, add it to the workflow RMB->Modules Next->Find module->Tomo statics 2D and set all the parameters.

Connect all the input data vectors from the previous module and import the FBpick binary file:

Usually, testing parameters starts with wider cells size of the model because initially the main goal is to find the most appropriate values for other parameters. When we have optimal result, it is reasonable to reduce cell size and increase the number of iterations. Eventually, a dense gradient depth velocity model is calculated as well as static corrections (ms) for sources and receivers. A depth velocity model may be used for the depth migration step as the upper part of the PSDM model.

Tomography parameters:

Tomography parameter explanation:

Velocity model params:

Definition cell-velocity model.

Deepest refractor velocity - the maximum values of guide-wave.

Datum – start value for tomography calculation, constant datum.

Max depth – maximum depth in meters of the model.

Step Z depth out – vertical cell size in meters.

Step X depth out – horizontal cell size in meters.

Step depth time table – vertical step for travel time solution, smaller value – more detailed result, but time consuming.

Horizontal smooth – horizontal velocity smoothing in meters.

Frequency:

Parameter for making solution more detailed in terms of travel time decomposition, bigger value – more detailed result (high spatial frequency), but time consuming. Pay attention on figure 3 Travel times scheme, green constrain depend on this parameter, less frequency - smoothie result.

Number of velocity updates

Tomography iterations.

Enable velocity inversion

Allows to have velocity inversion on the model.

Picking interpolation:

Options for interpolation missing first breaks picks:

Nearest – get picks from the nearest traces, use it in case of small gaps in FB picks.

Regression – model picks by using regression method.

Velocity min

Value for making minimum velocity constrain (m/s).

Velocity max

Value for making maximum velocity constrain (m/s).

Offset min

Value for making minimum offset constrain (m).

Offset max

Value for making maximum offset constrain (m).

Offset step

Step for splitting offset by classes (m).

Regression aperture

Value of aperture for regression algorithm (m), bigger value – smoother solution.

Maximum deviation (picking to model)

Maximum variations of observed picks according into the model trend.

TopoFilling:

Visual parameters for displaying velocity model, it is not used for calculations:

Fill topo

Values above the relief elevations for filling.

Velocity filling – define exact value for filling.

Samples topo fill – how many samples to fill.

Solver params

Residual statics solver parameters:

Step X depth time for solver - Vertical step for residual travel time solution, smaller value – more detailed result, but time consuming.

Regression aperture residual - Value of aperture for regression algorithm (m), bigger value – smoother solution.

Max pick time deviation - Maximum picking variations (ms) for residual static solution.

Global iteration count - Number of global iterations for residual static solution.

Local iteration count - Number of local iterations for residual static solution.

By the way, for 3D type tomography Tomo statics 2D/3D there is an option for define the low velocity zone by using Window for calculation V0 parameter which determines an average velocity in the interval (datum - window).

Now we have defined all the parametrization and the next step is launching the tomography calculation by double LMB click on the module or pressing Solve action. Use GPU and distribution mode for time consuming procedures or big acquisition. Those functions must be configures on a cluster by your system administrator:

To visualize tomography velocity model and convergence, user should generate vista items. The module generates many QC windows: the location map, convergence refraction tomo, tomography velocity view and Current source/receiver/bin gather.

For tomography velocity view, the user can define how many number of iteration required in the parameters tab. Once the solution is computed we can do a QC of the statics result by going to the any one of the current source/receiver/bin gather and enable/disable the Current source gather and Current source gather applied options from the View properties of the Current source/receiver/bin gather. Interactive fast QC by applying statics correction on gathers with different sorting: source, receiver and bin.

Check the source gather before and after statics applying. Choose by LMB clicking on any source on the location map:

And look at the source gather ViewVista window and check gather before and after statics applying:

Source gather before (left) and after (right) tomography statics applying:

Source gather and different types of picks for QC: User picks (purple), tomography (green) and residual (red):

Another important QC tool is convergence graphs where you can understand how many tomography iterations should be enough. Every iteration must reduce an error/mistie between observed and modeled times. If the error is too high or the tomo convergence is not matching with the residual convergence then the user can adjust the solver parameters and/or other parameters and click on the Solve option from the action item menu. It will recalculate the solution and update the displays.

Graphs of the main-tomography iteration (purple) and residual (red) solution:
 

Depth velocity model after tomography (try to change default color palette):

4) Static QC. The next step is checking statics corrections by using Statics QC module.

       ----------------------------------------------------------------------------------------------------------------------------------------------------------

For static QC we can use two modules:

• Statics QC if input statics correction from Tomo statics 2D

• Statics azimuthal QC if input statics correction from Refraction FB picking - azimuthal solver guide / phase / aperture

----------------------------------------------------------------------------------------------------------------------------------------------------------

Make input data item connection:

Use default parameters:

Open all vista groups. There are many QC tools like source, receiver statics graphs and maps, source and receiver gathers before and after static applying, tables with static corrections etc. Statics QC vista windows: source statics graph (left top), source statics map (bottom), source gathers overlaying before and after static correction applying (right):

      -----------------------------------------------------------------------------------------------------------------------

An internal data base of the g-Platform includes seismic data sets and all types of libraries.

Library is all other types of data except the seismic, i.e. velocity, statics, muting, etc.

----------------------------------------------------------------------------------------------------------------------

Before we continue, let's save library of tomo statics into data base. Add new module to the workflow: Save item, it saves any libraries into DB. Define a name of statics library - Refraction_statics and execute the module. After saving a library we can remove Save item module or save it in the workflow:

Check that your static library was saved into DB. Open g-Platform main window and go to Project data tab and write a mask for filtering static library:

Now we can do a quick QC on a stack section. Apply statics corrections to seismic traces and build a stack section via Stack imaging.

Stack section before (left) and after (right) tomography statics applying:

We will apply refraction statics correction during residual static and noise attenuation steps.

Next step >>> Velocity analysis (iteration 1, manual picking)

If you have any questions, please send an e-mail to: support@geomage.com    

First Break Picking and Refraction Statics Tutorial - Geomage g-Platform - YouTube