Post-processing

Generic post-processing

GEOS Extract And Merge Blocks

PVGeosBlockExtractAndMerge is a Paraview plugin processing the input mesh at the current time in two steps:

  1. Extraction of domains (volume, fault and well) from a GEOS output multiBlockDataSet mesh

  2. Actions on each region of a GEOS output domain (volume, fault, wells) to:

    • Merge Ranks

    • Identify “Fluids” and “Rock” phases

    • Rename “Rock” attributes depending on the phase they refer to for more clarity

    • Convert volume meshes to surface if needed

    • Copy “geomechanics” attributes from the initial timestep to the current one if they exist

This filter results in 3 output pipelines with the vtkMultiBlockDataSet:

  • “Volume” contains the volume domain

  • “Fault” contains the fault domain if it exist

  • “Well” contains the well domain if it exist

Input and output meshes are vtkMultiBlockDataSet.

Important

  • The input mesh must be an output of a GEOS simulation or contain at least three blocks labeled with the same domain names:

    • “CellElementRegion” for volume domain

    • “SurfaceElementRegion” for fault domain if the input contains fault

    • “WellElementRegion” for well domain if the input contains well

    See more about Element region in GEOS documentation.

  • The filter detects automatically the three domains. If one of them is not in the input mesh, the associated output pipeline will be empty.

To use it:

  • Load the plugin in Paraview: Tools > Manage Plugins … > Load New … > …/geosPythonPackages/geos-pv/src/geos/pv/plugins/post_processing/PVGeosBlockExtractAndMerge

  • Select the Geos output .pvd file loaded in Paraview to process

  • Select the filter: Filters > 2- GEOS Post-Processing > GEOS Extract and Merge Blocks

  • Apply

GEOS Log Reader

PVGeosLogReader is a Paraview plugin that allows to read GEOS output log.

Input is a file and output is a vtkTable containing log information.

Warning

The reader is compliant with GEOS log before commit version #9365098. For more recent version, use the csv or hdf5 export options from GEOS.

To use it:

  • Load the plugin in Paraview: Tools > Manage Plugins … > Load New … > …/geosPythonPackages/geos-pv/src/geos/pv/plugins/post_processing/PVGeosLogReader

  • Open (File>Open…) and Select GEOS output log .out/.txt file

  • In the “Open data with…” window, Select PVGeosLogReader reader.

Geomechanics workflows

GEOS Geomechanics Calculator

PVGeomechanicsCalculator is a paraview plugin that allows to compute additional geomechanics properties from existing ones in the mesh.

To compute the geomechanics outputs, the mesh must have the following properties:

  • The Young modulus and the Poisson’s ratio named “youngModulus” and “poissonRatio” or bulk and shear moduli named “bulkModulus” and “shearModulus”

  • The initial Young modulus and Poisson’s ratio named “youngModulusInitial” and “poissonRatioInitial” or the initial bulk modulus named “bulkModulusInitial”

  • The porosity named “porosity”

  • The initial porosity named “porosityInitial”

  • The delta of pressure named “deltaPressure”

  • The density named “density”

  • The effective stress named “stressEffective”

  • The initial effective stress named “stressEffectiveInitial”

  • The pressure named “pressure”

The basic geomechanics properties computed on the mesh are:

  • The elastic moduli not present on the mesh

  • Biot coefficient

  • Compressibility, oedometric compressibility and real compressibility coefficient

  • Specific gravity

  • Real effective stress ratio

  • Total initial stress, total current stress and total stress ratio

  • Elastic stain

  • Real reservoir stress path and reservoir stress path in oedometric condition

The advanced geomechanics properties computed on the mesh are:

  • Fracture index and threshold

  • Critical pore pressure and pressure index

PVGeomechanicsCalculator paraview plugin input mesh can be a vtkUnstructuredGrid or a vtkMultiBlockDataSet of vtkUnstructuredGrid. If the input mesh is a vtkMultiBlockDataSet, the geomechanics properties will be computed on each vtkUnstructuredGrid. The output mesh has the same type than the input one.

To use it:

  • Load the plugin in Paraview: Tools > Manage Plugins … > Load New … > …/geosPythonPackages/geos-pv/src/geos/pv/plugins/post_processing/PVGeomechanicsCalculator

  • Select the mesh you want to compute geomechanics properties on

  • Select the filter: Filters > 2- GEOS Post-Processing > GEOS Geomechanics Calculator

  • Change the physical constants if needed

  • Select computeAdvancedProperties to compute the advanced properties

  • Apply

Plot Mohr’s circle

PVMohrCirclePlot is a ParaView plugin that allows to compute and plot Mohr’s circles of selected cells and times from effective stress attribute.

Input is a vtkMultiBlockDataSet or vtkUnstructuredGrid.

This filter results in opening a new Python View window and displaying Mohr’s circle plot.

To use it:

This plugin requires the presence of a stressEffective attribute in the mesh. Moreover, several timesteps should also be detected.

Warning

The whole ParaView pipeline will be executed for all timesteps present in the initial PVD file. Please be aware that the number of pipeline filters and timesteps should be as limited as possible. Otherwise, please consider going to get a cup of coffee.

  • Load the plugin in Paraview: Tools > Manage Plugins … > Load New … > …/geosPythonPackages/geos-pv/src/geos/pv/plugins/post_processing/PVMohrCirclePlot

If you start from a raw GEOS output, execute the following steps before moving on. - First, consider removing some unnecessary timesteps manually from the PVD file in order to reduce the calculation time and resources used in the following steps. - Load the data into ParaView, then apply the PVGeosBLockExtractAndMerge plugin on it. - Select the filter output that you want to consider for the Mohr’s circle plot.

  • Extract a few number of cells with the ExtractSelection ParaView Filter, then use the MergeBlocks ParaView Filter

  • Select the resulting mesh in the pipeline

  • Select the filter: Filters > 0- Generic-Processing > Plot Mohr’s Circle

  • Select the cell Ids and time steps you want

  • (Optional) Set rock cohesion and/or friction angle

  • Apply

Note

After a first application, select again cells and time steps to display, then

  • Apply again

  • Click on Refresh Data (you may have to click twice to refresh the Python view correctly).

GEOS Surface Geomechanics

PVSurfaceGeomechanics is a Paraview plugin that allows to compute additional geomechanical attributes from the input surfaces, such as shear capacity utilization (SCU).

Input and output are vtkMultiBlockDataSet.

Important

  • Please refer to the GEOS Extract and Merge Blocks plugin to provide the correct input.

  • This filter only works on triangles at the moment. Please apply a triangulation algorithm beforehand if required.

To use it:

  • Load the plugin in Paraview: Tools > Manage Plugins … > Load New … > …/geosPythonPackages/geos-pv/src/geos/pv/plugins/post_processing/PVSurfaceGeomechanics

  • Select any pipeline child “Fault” from “GEOS Extract and Merge Blocks” plugin

  • Select the filter: Filters > 0- Generic-Processing > GEOS Surface Geomechanics

  • (Optional) Set rock cohesion and/or friction angle

  • Apply

GEOS Geomechanics Workflow

PVGeomechanicsWorkflow is a Paraview plugin that executes multiple plugins:

  1. PVGeosBlockExtractAndMerge

  2. PVGeomechanicsCalculator

  3. PVSurfaceGeomechanics (if the input mesh contains faults)

PVGeosBlockExtractAndMerge is a Paraview plugin processing the input mesh at the current time in two steps:

  1. Extraction of domains (volume, fault and well) from a GEOS output multiBlockDataSet mesh

  2. Actions on each region of a GEOS output domain (volume, fault, wells) to:

    • Merge Ranks

    • Identify “Fluids” and “Rock” phases

    • Rename “Rock” attributes depending on the phase they refer to for more clarity

    • Convert volume meshes to surface if needed

    • Copy “geomechanics” attributes from the initial timestep to the current one if they exist

PVGeomechanicsCalculator is a paraview plugin that allows to compute basic and advanced geomechanics properties from existing ones in the mesh. This is donne on each block of the volume mesh.

The basic geomechanics properties computed on the mesh are:

  • The elastic moduli not present on the mesh

  • Biot coefficient

  • Compressibility, oedometric compressibility and real compressibility coefficient

  • Specific gravity

  • Real effective stress ratio

  • Total initial stress, total current stress and total stress ratio

  • Elastic stain

  • Real reservoir stress path and reservoir stress path in oedometric condition

The advanced geomechanics properties computed on the mesh are:

  • Fracture index and threshold

  • Critical pore pressure and pressure index

PVSurfaceGeomechanics is a Paraview plugin that allows to compute additional geomechanical attributes from the input surfaces, such as shear capacity utilization (SCU). This is donne on each block of the fault mesh.

This filter results in 3 output pipelines with the vtkMultiBlockDataSet:

  • “Volume” contains the volume domain

  • “Fault” contains the fault domain if it exist

  • “Well” contains the well domain if it exist

Input and output meshes are vtkMultiBlockDataSet.

To use it:

  • Load the plugin in Paraview: Tools > Manage Plugins … > Load New … > …/geosPythonPackages/geos-pv/src/geos/pv/plugins/post_processing/PVGeomechanicsWorkflow

  • Select the Geos output .pvd file loaded in Paraview to process

  • Select the filter: Filters > 2- GEOS Post-Processing > GEOS Geomechanics Workflow.

  • Change the physical constants if needed

  • Select computeAdvancedProperties to compute the advanced properties on volume mesh

  • Apply