Thursday, 20 October 2016

Investigating Spectral Enhancement using the Bedform Indicator Attribute

Spectral Enhancement is crucial process when analysing thin beds. The most common method to QC Spectral Enhancement is to use the slider in the Spectral Expression Tool. However, sometimes data quality can limit the true extent of enhancement that can be visualised, often characterised by washed out reflectors.

In order to get around this potential problem, we can use the Bedform Indicator attribute to help us. This attribute skeletonises the reflectors in the seismic data. It can be found in Workflows>Processes>Attributes>Trace Attributes>Bedform Indicator.

The image below shows the Bedform Indicator with the Peaks and Troughs option selected, which aids in visualising reflectors which appear more washed out.

This can sometimes lead to an over-cluttered image, so in the example below only Peaks will be outputted. The Bedform Indicator should be run on both the noise cancelled volume and the spectrally enhanced volume, and a slice can be created to compare the two results.

The red corresponds to positive reflectors while the green corresponds to positive doublets (often indicative of thin beds). As can be seen, some positive doublets now have been resolved fully into positive reflectors, helping to show the true extent of the spectral enhancement and allows the interpreter to explore the relationships between reflectors in more detail.

Tuesday, 4 October 2016

GeoTeric 2016.2

GeoTeric 2016.2 continues to improve the user experience by focusing on strengthening the core interpretation functionality. This includes:  

Faster performance in the Expression tools allowing you to work interactively with even larger datasets:
  • Fault Expression is up to 6 times faster
  • Noise Expression is up to 2-3 times faster.
  • Improvement seen will be dependent on data size, hardware and setup
  • The Expression tools also have Cursor location reporting feature added to help you navigate your volume.

Extended well functionality which includes: 

  • Ability to display logs as ribbons and up to 3 logs simultaneously for each well.
  • A new ‘Well Edit’ table allowing the user to QC the well data, and also to edit wells to remove spurious points, for example.
  • We have improved the handling of datum and depth measurements.
  • The 3rd party Links have been extended to allow control over which logs are imported.
  • The 3rd party Links also now clearly indicate when the well transfer is finished preventing the user from accidentally interrupting the transfer.
  • We have resolved several issues whereby the well would cause visualisation artefacts.
  • We have also improved the performance so visualising several wells with or without logs should not affect the overall interactivity of GeoTeric.

Left & Middle: GeoTeric’s new multi log display and improved multi well visualisation and performance. Right: Geobody, fault & horizon interpretation : in GeoTeric.

 Improved metrics & Hydrocarbon Calculator:

  • We have introduced further metrics to both Adaptive Geobodies and horizons. These can be found in the properties panel for each item, with the option to display it in the scene for easy snapshots.
  • To complete the Prospect Generation workflow in GeoTeric we have introduced a STOIIP/GIIP calculator to allow you to get a rough idea of the size of your reservoir. The hydrocarbon calculator can be accessed from the Tools menu, or alternatively by right-clicking on an Adaptive Geobody or a Horizon in either the project tree or 3D viewer.
  • The ruler in the 3D scene now has a show/hide option for better control when taking screenshots.

GeoTeric is now offered with a modular license option:

  • The modules are: Condition, Reveal, Interpret (which is the core license module) and Classify. For further information on what each module offers please contact Support.
  • Extended license functionality includes a “use offline” (borrow) feature.
  • We have improved the stability of GeoTeric when temporary interruptions to the license server occur.

Extended editing & transfer of fault sticks:
  • By selecting one or more points in a fault stick, the user can edit their location (either by mouse or keyboard).
  • Users can now delete points when picking Fault Sticks, either by stepping back when an erroneous pick is placed or by coming back later and deleting a misplaced pick.
  • We have also extended the 3rd party links to allow the transfer of fault sticks to and from GeoTeric.

Bug fixes and smaller enhancements include:
  • GeoTeric now fully supports projects with a negative start time or depth, including well data.
  • In the scatterplot tool we have added a save option, allowing you to compare data trends and determine the most optimal attribute to define the objective.
  • When exporting colour blends, there is an option to export a colour bar which is compatible with Jason*, allowing you to bring in your blends into Jason. Note: an option to import GeoTeric volumes is available in Jason 9.6.
  • Sometimes when interacting with the project tree in GeoTeric or moving through a volume, GeoTeric would slow down. This is resolved.
  • In GeoTeric 2016.1, the blends sometimes appeared with visual artefacts.  This is resolved.
  • Several issues have been resolved around Fault Stick picking including random crashes, faults disappearing, colours being wrong, ability to rename Fault Sets.
  • Fault Trends produced West-aligning faults for volumes orientated in certain way. This has been resolved.
  • Fault Expression sometimes failed to generate a 3D result.  This has been resolved.
  • Subset dialog now shows correct volume size
  • Sometimes the blend wasn’t showing the correct opacity applied. This has been fixed.
  • Now able to rename and delete items which are viewed in the 3D scene.
  • Expression tools crashed if the Z increment was not divisible with extent.  This has been resolved.
  • In some circumstances, volumes would be listed multiple times in selection boxes.  This has been resolved.
  • When a CMY blend was created using Fault Expression, it appeared in the tree with a red exclamation mark meaning it appeared invalid.  This has been fixed.

*Trademark of CGG

Wednesday, 21 September 2016

Optimisation of Amplitudes and Scaling for Visualisation

One common question for GeoTeric users is how best to scale amplitude data for optimal visualisation and analysis. In this article, we present some useful definitions and tips for ensuring the best experience when interpreting seismic data in GeoTeric.

1.    GeoTeric Data Types

GeoTeric only supports integer values for volume data, meaning that each amplitude within a seismic cube being displayed must be a whole number. To facilitate this, GeoTeric volume formats support 8, 16 and 32 bits, such that the following dynamic ranges are possible:
Bit depth
Maximum Dynamic range
Unsigned 8-bit
Signed 16-bit
-32768 to +32767
Signed 32-bit
-2147483648 to 2147483647
This means that for example, a 16-bit seismic cube can contain approximately 65,000 different discrete amplitude values. In this notation, n-bit means 2n values, so 16-bit is equivalent to 216 degrees of freedom.

You can see what amplitudes are present in the interval of interest by looking at the cursor reporting in GeoTeric scenes, the Opacity Editor, and the colour bar compression metrics.

2.   Colour Bars and Visualisation

While most modern GPUs are capable of displaying 24-bit colour depth (8-bits per channel in red, green and blue), it can be extremely resource intensive and impractical to render seismic data or attributes in this manner. LCD monitors are also not always good at displaying 24-bit true colour. For this reason, different visual scaling and interpolation methods are used to balance quality of display with computer resource consumption.

The choice of interpolation method can also have an influence, as selected in Options > Volume Interpolation.  Standard interpolation utilizes full 24-bit rendering, at the expense of memory, however Enhanced and Smooth interpolation render at 12-bit.

In GeoTeric, there are two aspects in play: ensuring that the seismic amplitudes are scaled appropriately for the task in hand, and that the colour bar is interacting with those amplitudes suitably. The GeoTeric colour bars do not have fixed values: when displaying data, the colour map will be normalized and extrapolated across the range of values present, and the visual fidelity depends on the interpolation method used.

In general, GeoTeric will autofit colour maps across a range of 2n that is the smallest bracket around the detected amplitude range. So if a seismic cube has an amplitude range of approximately +/- 8000, GeoTeric will map the colour bar across a 14-bit range (16384 total values, or +/-8192).

A commonly encountered issue is in the visualization of very low amplitudes, relative to those elsewhere in the seismic cube. In these scenarios, the lowest amplitudes might look faint or indistinct as they represent a very small segment within the entire range, as in the below figure. Some solutions are presented below.

These seismic amplitudes appear blocky and washed out where they are low, because they are occupying a relatively small portion of the colour map space.
3.   Scaling on Import

The most effective step is to ensure the data has been scaled appropriately when loaded into GeoTeric, especially if it originally comprised floating point values, e.g. a range of -1.476 to 1.928.

Search in the GeoTeric software help file for “Data Loading Tips” and “Load Floating Point Data” for tutorials on these topics.

4.   Scaling within GeoTeric

Once data is loaded in GeoTeric however, the amplitudes of interest might still be very low and difficult to see. This can be because there is a strong contrast between high and low amplitudes, peculiarities of gain correction, or low impedance variation. In this case rescaling can be applied within GeoTeric. Note however that if amplitudes have been compressed during the loading process, and information has been lost, the data might need to be reloaded with attention paid to the preservation of those amplitudes. 

Please search the GeoTeric help file for “Data Scaling”, “Convert”, and “Parser” for more information on carrying this out.

5.   Specialist Scaling

In some cases, tricks using the GeoTeric colour map compression can be used to help with low amplitudes. If a 16-bit volume has high amplitudes in redundant parts of the data, for example on the seafloor reflection, but low amplitudes in the reservoir section, it can be clipped to take advantage of the colour map compression. The Parser tool is used to achieve this.

In this example, the amplitudes above and below 4000 are clipped such that the data fits within a 12-bit (+/- 4096) dynamic range, and will then make better use of the full range of the colour map. The amplitudes in the interval of interest remain unchanged.

Volume type
Highest amplitudes
+/- ~30,000
Amplitudes of interest
+/- ~3000
Parser expression
((im1>=4000)*4000) + ((im1<=-4000)*-4000) + (((im1<4000)&(im1>-4000))*im1)

The amplitudes of interest in this data now occupy the full useful extent of the colour bar, and are readily visible without further colour bar compression.

Non-linear scaling can also be applied where there are very extreme amplitude contrasts, although this should not be applied in cases where amplitude preservation is sensitive. In this example a square root function has been applied to the amplitudes, such that the low values are boosted, and the high values suppressed. Note that the power law function in the GeoTeric Parser retains the correct sign, to preserve polarity in volume calculations.

Volume type
Parser expression

In this example, the highest amplitude reflector is preserved, while the background signal which was previously too weak to see has been restored.

For further Parser examples and explanations of the terms, expressions and syntax applied, search for topics containing “Parser” in the help.

For help and more information on any of these topics, feel free to contact your local GeoTeric office, or email

Thursday, 8 September 2016

Effect of Post-Stack Noise Attenuation in Acoustic Inversion

It's well known that there are many methods of data conditioning that improve the signal to noise ratio. In this case study the post stack noise attenuation process applied to the data was carefully performed to preserve the signal associated with geological features within the Pre-Salt context of the Santos Basin. The objective of the noise attenuation is to improve the input volume to the inversion process. 
Figure 1 - common noise in the pre-salt context

Mothodology - Noise Attenuation
The workflow to reduce the noise is out lined below and illustrated in Figure 2:

  1. Determine the grid oriented filter size, a large filter in the X&Y directions and small in Z, will reduce the continuity of the vertical noise associated with the salt/migration smiles. In order to create an input for the Steering Volumes, the Smoothing Attribute is calculated. Recommended filter size for this case: 9x9x5.
  2. The Steering Volumes, of Dip and Azimuth are calculated from the Smoothing Attribute, this becomes the reference volume.
  3. Finally the SO Noise filter which is an structurally oriented mean noise filter, is applied to the original seismic volume, which is guided by the reference volume.
Figure 2 – Noise Attenuation Workflow for Pre-Salt

Methodology – Acoustic Inversion
In geosciences, through reflection seismic data, the seismic inversion process estimates the quantitative rock-petrophysical properties such as porosity, lithology and fluid saturation of a reservoir. This process, in a simply way, is based on a straightforward convolution of a reflectivity model and an estimated wavelet to produce synthetic traces, which are compared to the seismic input. The reflectivity model is updated to ensure there is a high correlation between the synthetic and seismic data (figure 3).
Figure 3 – Deterministic Inversion workflow
Noise Attenuation

As can be seen in the figure below (figure 4), the SO Noise filter has removed a lot of the migration smiles associated with pre-salt seismic data. The input data is on the left and the noise attenuated data is on the right. The reflectors are more continuous and the faults are still visible.
Figure 4 – Original Seismic (Left) x Noise attenuated seismic (right)

Acoustic Inversion
In this case study, both the input data and the noise attenuated data were inverted.  On the left is the inversion results using the original seismic data and on the right the inversion results using the noise attenuated data.  The inversion results using the noise attenuated as an input contains more continuous layers, which is more consistent with what well data would contain.  This result can be correlated to petrophysical properties and the reservoir mapper more easily as compared to the result on the left.
Figure 5 – Acoustic inversion over raw seismic (left) x Acoustic inversion over noise attenuated seismic (Right). the colors nearer to yellow have low impedance and nearer to blue are the high ones.
By applying a noise attenuation filter to noisy pre-salt data, it allows for better interpretation of events on the output data, but also provides a better input into other attributes, such as an inversion. 

Wednesday, 31 August 2016

GeoTeric Technology Forums

This autumn GeoTeric is busy hosting Technology Forums in two different locations. 


The 27th of September we’ll be hosting our first ever Technology Forum in Norway. Join us for a day packed with technical talks, the latest showcasing of the new release of GeoTeric & our collaborative workflows as well as special topics presented by our guest speakers.

@Cyviz AS, Vestre Svanholmen 6, 4313 Sandnes

For full details please click here


On the 13th of October, it's Houston's turn.
We will be demonstrating how the industry uses Cognitive Interpretation to reveal geology from seismic data, via informative case studies, live demos and discussion.  
@Omni Houston Hotel at Westside, 13210 Katy Freeway, Houston, TX 77079
For full details on the Houston forum- please see here

Thursday, 25 August 2016

Testing on Small Volumes, Processing on Large

With the acquisition of ever larger and denser surveys, it can be challenging to efficiently optimise and run processes on your data. GeoTeric solves this problem by allowing you two different methods to rapidly test your processes on small volumes, before running the processing on the full dataset. An example will be shown below using the Noise Expression tool. The steps assume you have the full dataset loaded into GeoTeric.

  1. Visualise the full dataset in the main viewer.
  2. Click the ‘Extents’ button in the bottom right.
  3. Trim the volume to focus on a target area.
  4. Enter a name in the ‘Subset’ box and click ‘Save Subset’.

  1. Open Noise Expression, selecting the subset from the dropdown.
  2. Optimise your process as normal.
  3. Once you are satisfied with the results, hit ‘Generate 3D Volumes’.
  4. Firstly, note the location and name of the batch job. This can be applied to any dataset at a later point. Instructions for this will be given later.
  5. If you wish to apply the current process to the full volume, simply select it from the dropdown at the top.
  6. If you choose the full volume from the dropdown and click ‘Process’, the processing will be applied to that dataset.
  7. This method can be used in Noise Expression, Spectral Expression, HDFD and Fault Expression.

Using the Batch Job

This will allow you to load up any batch job that has been saved previously, and apply it to any datasets you have available in your project. Care should be taken to only apply processes to appropriate datasets.

  1. From the menus in the top left of the main viewing window, select ‘Workflows’ > ‘Processes and Workflows…'
  2. From the top left of the Processes and Workflows window, select ‘File’ > ‘Open…’ (The window is actually titled ‘Batch Processing Framework’)
  3. Browse to the location of your batch job and open it. This will generally be in the ‘batchjobs’ directory within your project.
  4. On the left hand tab, select the input dataset from the drop down.
  5. The output dataset names can also be changed at the bottom of each tab. This is not essential however.

Thursday, 11 August 2016

Combining different Noise Cancellation/Spectral Enhancement volumes using Time Variant

The aim of this blog post is to describe how to combine different Noise Cancellation/Spectral Enhancement volumes using a time variant method.
This workflows can be used for:
  1. Creating a combined Noise Cancellation volume in which different noise filter were applied to target different noise contents that vary vertically and can be bounded by a horizon.
  2. Creating a Spectral Enhancement volume with different enhancement parameters that can be divided by a horizon.
The following workflows is for two Noise Cancellation volumes where there is a need to apply a more aggressive noise filter below the horizon and only a gentle filter above the horizon, but the same workflow applies for the Spectral Enhancement volumes too.

Input needed for this workflows are horizons which separate the areas with different noise content vertically, Noise Cancellation volumes, which target the different intervals and a time volume, that is used for the combination. Example in this post will show the workflows to combine 2 areas which is separated by a horizon.
Figure 1: Original Seismic Volume

The first step is to create two Noise cancellation volumes, a gentle noise filter above the horizon and an aggressive filter below the horizon. The next step is to combine these two volumes into one. To combine these two volumes, a time volume is needed, which is unflatten using the horizon and formatted for the combination purpose.

The following is the workflow to create the a formatted Unflatten Time Volume.
  1. Go to Workflows>Processes and Workflows>Processes>Utilities>Time Volume.Input your seismic volume to create the time volume.
    Figure 2: Time volume displayed using the spectrum colour bar, where green is the higher value.  
  2. Use the Horizon Tools to unflatten the Time Volume. Go to Tools> Horizon Tools> Flatten/Unflatten. Toggle Seismic as volume type. Then specify the input horizon, the input volume will be the newly created Time Volume. The output volume will need to be named as Unflatten Time Volume. Change the mode to Un-flatten Horizon and Click Apply. Once you input the horizon, a Set Flattened value, which is the mean of the input horizon, will be displayed at the bottom of the Horizon Tool menu. Write down the value as it will be used to format the Unflatten Time Volume later.
    Figure 3: Unflatten Time Volume   
  3. The next step is to format the Unflatten Time Volume so that one of the Noise Cancellation volumes can be assigned to the data above the horizon and the other Noise Cancellation can be assigned to below the horizon. Go to Workflows>Processes and Workflows>Processes>Volume Math>Parser. Input the Unflatten Time Volume, and use the following parser expression: ((im1>0)*(im1-A)) where A is the “set flattened value to” number that was given when unflattening the time volume. Now you have the formatted Unflatten Time Volume.

    Figure 4: Formatted Unflatten Time Volume displayed as Redwhiteblue colourbar. Red are the negatives values, whites are 0 values and blues are the positive values. It will be observed that along the horizon, values are all 0, this will be used as a transition zone between the 2 volumes.
The final steps now is to combine the two noise cancellation volume using the formatted Unflatten Time Volume. To combine the volumes, go to Workflows>Processes and Workflows>Processes>Volume Math>Parser. The first volume should be the formatted Unflatten Time Volume as im1, Gentle Noise Cancellation as im2 and Aggressive Noise Cancellation as im3. The Parser expression used will be:

((im1<0)*im2) + ((im1>100)*im3) + (((im1>=0) & (im1<=100))*(im1*im3/100)) + (((im1>=0) & (im1<=100))*((100-im1)*im2/100))
As a part of QC, you can compare the Combined Noise Cancellation result with the Aggressive and Gentle Noise Cancellation volume, or look at the Difference volume of each of the Noise Cancellation volumes (Original-Noise Cancellation Volume) to see the related noise in the Combined Noise Cancellation volume.
Combining these two volumes in this way, will avoid abrupt changes along the horizon. Instead there will be a transition zone of 100ms taken from both volumes. This is not limited to only two areas of Noise Cancellation, more volumes can be combined, in a sequence, using the same workflow.This workflows can also be used to combine different level of Spectral Enhancement as long as there are horizons that separate the areas with different Spectral Enhancement.

Monday, 1 August 2016

Multi-volume phase alignment QC

GeoTeric offers a number of workflows that allow a QC of the phase alignment between multiple volumes. These workflows can be used to analyse the similarity between angle stacks or 4D vintages. This enables the user to assess the phase alignment of the data to ensure an accurate analysis of the data.

Peaks and troughs replacement
The peaks of stack 1 can be replaced by the data from stack 2. If they are phase aligned, the result will be a peak. If they are misaligned, the result will be a trough. The troughs from stack one will be ignored and replaced by zeros.

This volume can be produced by selecting the two stacks in the Parser and using the expression (im1>0)*im2

The Parser is available in Processes & Workflows -> Processes -> Volume Maths -> Parser.

In a similar way, the troughs of stack 1 can be replaced by the data from stack 2. If they are aligned, the result will be a trough. If they are misaligned, the result will be a peak.

This volume can be produced by using the expression (im1<0)*im2 in the Parser.

Example from the Parihaka dataset, offshore New Zealand. On the left, peaks in the near have been replaced by the far, so troughs indicate misalignment. On the right, troughs in the near have been replaced by the far, so peaks indicate misalignment.

Bedform stack
A much more detailed similarity analysis can be achieved by using a Bedform stack.
The Bedform indicator attribute can be computed for each of the stacks, and then multiple Bedform attributes can be combined using a Parser expression to create the Bedform stack.
To combine three Bedform indicators, we can use the following Parser expressions depending on the data type.

For 32bit:

For 16bit:

For 8bit:

These Parser expressions are simply (im1+im2+im3)/3, but removing the doublet values in each of the Bedform attributes.

If the output is visualised using the Azimuth colour map, good alignment will be seen in black, areas with dispersion along a peak in red, and areas with dispersion along a trough in blue. The lighter the blue and red colours, the worse the alignment.

Example from the Parihaka dataset, offshore New Zealand. The bedform stack has been computed using the near, mid and far angle stacks.

Alignment quantification workflow
To quantify the alignment between two stacks (im1 and im2), the following Parser expression can be applied:


This will produce a volume with
+1 when there is good phase alignment between the stacks
-1 when there is bad phase alignment between the stacks

Note that any blank areas in any of the stacks will bias the results, so it is recommended to apply this workflow on a subset around the area of interest.

By using the Voxel Pick tool, available in the properties panel when highlighting the volume, a voxel pick can be placed for a +1 value and another one for a -1 value. Going to the Picked voxels tab, the two picked voxels will be displayed. Right clicking on each of them will bring up the Measure option, which allows a measurement of the number of voxels within the dataset that have the value of the picked voxel. Using the default options in the measure tool we can eaily get the number of voxels with a +1 value, and we can repeat the process for the voxels with a -1 value.

This way, as the number of voxels with a +1 value and the number of voxels with a -1 value is known, a simple ratio or alignment percentage can be manually calculated.

If the alignment is close to 50%, that is a bad alignment, as it means a peak in stack 1 has almost the same likelihood of being a peak or a trough in stack 2. Anything close to 90% or 100% means there is a good phase alignment.

Using GeoTeric we are able to check if the stacks are aligned properly or not. If there is a large amount of misalignment, the user should consider reprocessing or recreating the angle stacks before continuing with the AVO or 4D analysis.