Skip to end of metadata
Go to start of metadata

Introduction

Example of a '.task' file used to conduct a magnetic geophysical inversion of a region in Queensland. This was performed as part of a project with the Geological Survey of Queensland. 

WARNING: Please note many of the parameters contained within this task file are specific to this region (eg. the IGRF, rock property values, project boundaries, stratigraphic column). All parameters should be reviewed if this file is used as a template for future processing.

Comments and annotations

Note annotations are written in italics, these should not be included in your task file, though they have commented-out just in case. If you would like to keep them for your own reference, or indeed if you want to add your own (highly recommeded) then precede each line with a '#'. A commented example is provided in the first line of the task. Commented lines will not be read by the geomodeller inversion engine, so therefore will not cause any errors. Be careful that you do not 'comment-out' lines that are necessary.

References to gravity parameters

While this .task has been specifically edited to execute a magnetic inversion, references to gravity-specific parameters have been retained for explanatory purposes and so this file can be used for gravity inversion. Please note the gravity parameters have not been edited, so have been left at their default values. Gravity parameters will need to be edited similar to the magnetic.

Editing task files

Task files can be edited using standard Windows text editors such as Wordpad and Notepad. Some more useful freeware text editors are Notepad++ , Programmer's Notepad and UltraEdit.

This is an example used for conducting magnetic inversion over the Lawn Hill region in Queensland. The setting for geometrical or property inversion is found in this file as well. It has been highlighted in the annotations in bold.

IMPORTANT: Throughout the file the output location for grids and voxets is defined (in the case of this example it is 'OutputGeology'). If this directory exists containing files from a previous run, the settings from the previous job will be incorporated into the new run. Any changes you make to the task file (rock property values, IGRF field etc.) will not be honoured. It is therefore important check whether the output directory is empty (either rename the output directory before re-running the job, or delete the contents) to ensure that the newer parameters will be incorporated and used in the new inversion job.

Task file

File begins below this line:

Task File
name: "new project from voxet"

description: "replicate the job file of the same name"

auditfile: "audit.rpt"

INTREPID_NODES=16

InversionTask {

    NewProjectFromVoxet {

    inputvoxet: "Model/Lawn_Hill_Voxet_subsampledxy.vo"; #_Defines the input voxet for the inversion. Ensure the voxet is the correct one and the path exists and is correct._

    new_project: "OutputGeology/imported.xml"; #_Defines the output directory and .xml file. This location will be automatically created once the process has been executed._

_\#The following lines define the stratigraphic column. It is important that the 'formations' or stratigraphic units defined here correspond correctly to other references throughout the file. Please note that these formations are specific to this particular region. They will need to be renamed for future work. _


        formation{

            index: 1;

            field: "a_BASEMENT";

            } #Note the name 'BASEMENT' is preceded by 'a_'. This sets this particular formation as the base layer. The next layer name is preceded with 'b_', the next with 'c_' and so in alphabetical order as the pile ascends.

        formation{

            index: 2;

            field: "b_BASE_LEICHHARDT_SUPERBASIN";

        }

        formation{

            index: 3;

            field: "c_CALVERT_SUPERBASIC";

        }

        formation{

            index: 4;

            field: "d_BASE_ISA_SUPERBASIN";

        }

        formation{

            index: 5;

            field: "e_BASE_RIVERSLEIGH";

        }

        formation{

            index: 6;

            field: "f_BASE_TERMITE";

        }

        formation{

            index: 7;

            field: "g_SOUTH_NICHOLSON_GROUP";

        }

        formation{

            index: 8;

            field: "h_CAMBRIAN";

        }

        formation{

            index: 9;

            field: "i_YELDHAM_GRANITE";

        }

        formation{

            index: 10;

            field: "j_WEBBERA_GRANITE";

        }

    }

}
_\#The following lines define some important parameters of the inversion.  _


InversionTask {

    NewCase {

    filename: "OutputGeology/imported.xml"; #this is the path where the inversion output will be sent to.

    case: "Case_Magnetic"; #this is a subdirectory of the above referenced location describing the nature of the modelling. In the case of joint inversion both 'Case_Magnetic' and 'Case_Gravity' will be created.

    externalvoxet: "Model/Lawn_Hill_Voxet_subsampledxy.vo"; #this is where your original voxet is located.

#the following lines describe the input geophysical grid that is being used for modelling. This example is a magnetic grid, as the fields suggest.
        
ObservedGridList {

            ObservedGrid{

            type: Magnetism;

            mean_elevation: 500; #mean elevation height of the survey

            precision: 1.0;

            Match_Trend: true;

            Match_Trend_Degree: 1;

            Match_Trend_Rate: 0;

            grid: "../InputGrids/LH_project_tmi.out.ers"; #location of the input grid.

            }

        }

    }

}
#this section assigns the following parameters to the case defined in the previous section


InversionTask {

    CaseControl {

        filename: "OutputGeology/imported.xml"; #keep this parameter consistent with that defined in the previous section

        case: "Case_Magnetic"; #keep this parameter consistent with that defined in the previous section

# this section describes the IGRF (International Geomagnetic Reference Field). These parameters will need to be checked and adjusted - see link in the above introduction. These settings are ignored in the case of gravity forward modelling and inversion.
        SetReferenceMagneticField{

            Magnitude:51093.0;

            Inclination: \-51.385;

            Declination: 6.609; 

            }


        SetReferenceDensity:2.67; #reference density for gravity inversion - ie what you think the average density of the rocks are in your study area.

        SetAboveTopoDensity:0.0; #set density for the volume above the top rock surface. On-shore assumption is that the volume is air, therefore density = '0.0'. Can be set to accommodate water or other substances if necessary.


#the following lines determine the rock properties and data distributions for modelling. Note the PDF (probability density function) is defined by the type ('Normal' - typically used for densities, 'LogNormal' - typically used for mag. susecptibility and 'Weibull' - a distrubtion typically used in extreme value analysis and change detection)
#the parameters for the PDF are set by specifying the mean and standard deviation. A bi-modal distribution can be designated by using a third variable that determines the proportion of the distribution that the PDF corresponds to.
#an example of a bi-modal distribution is shown in the magnetic parameters for the Base of the Leichhardt Superbasin (b_BASE_LEICHHARDT_SUPERBASIN) - extracted and shown below
#LithologyProperty {lithology:"b_BASE_LEICHHARDT_SUPERBASIN";property:Susceptibility;ProbabilityDistributionFunction {type:LogNormal;mean: 0.0500   ;stddev: 0.00001; percentage: 66};ProbabilityDistributionFunction2 {type:LogNormal;mean: 0.0800   ;stddev: 0.00001; percentage: 34}}

        SetLaw {

            LithologyProperty {lithology:"h_CAMBRIAN";property:Density;ProbabilityDistributionFunction {type:Normal;mean: 2.500  ;stddev: 0.005;}}

            LithologyProperty {lithology:"g_SOUTH_NICHOLSON_GROUP";property:Density;ProbabilityDistributionFunction {type:Normal;mean: 2.650  ;stddev: 0.04;}}

            LithologyProperty {lithology:"f_BASE_TERMITE";property:Density;ProbabilityDistributionFunction {type:Normal;mean: 2.600  ;stddev: 0.04;}}

            LithologyProperty {lithology:"e_BASE_RIVERSLEIGH";property:Density;ProbabilityDistributionFunction {type:Normal;mean: 2.8118 ;stddev: 0.005;}}

            LithologyProperty {lithology:"d_BASE_ISA_SUPERBASIN";property:Density;ProbabilityDistributionFunction {type:Normal;mean: 2.6935 ;stddev: 0.08;}}

            LithologyProperty {lithology:"c_CALVERT_SUPERBASIC";property:Density;ProbabilityDistributionFunction {type:Normal;mean: 2.597  ;stddev: 0.03;}}

            LithologyProperty {lithology:"b_BASE_LEICHHARDT_SUPERBASIN";property:Density;ProbabilityDistributionFunction {type:Normal;mean: 2.6647 ;stddev: 0.12;}}

            LithologyProperty {lithology:"a_BASEMENT";property:Density;ProbabilityDistributionFunction {type:Normal;mean: 2.8    ;stddev: 0.16;}}

            LithologyProperty {lithology:"i_YELDHAM_GRANITE";property:Density;ProbabilityDistributionFunction {type:Normal;mean: 2.650  ;stddev: 0.02;}}

            LithologyProperty {lithology:"j_WEBBERA_GRANITE";property:Density;ProbabilityDistributionFunction {type:Normal;mean: 2.650  ;stddev: 0.01;}}

            LithologyProperty {lithology:"h_CAMBRIAN";property:Susceptibility;ProbabilityDistributionFunction {type:LogNormal;mean: 0.00019  ;stddev: 0.00004;}}

            LithologyProperty {lithology:"g_SOUTH_NICHOLSON_GROUP";property:Susceptibility;ProbabilityDistributionFunction {type:LogNormal;mean: 0.0019   ;stddev: 0.00004;}}

            LithologyProperty {lithology:"f_BASE_TERMITE";property:Susceptibility;ProbabilityDistributionFunction {type:LogNormal;mean: 0.004    ;stddev: 0.00202;}}

            LithologyProperty {lithology:"e_BASE_RIVERSLEIGH";property:Susceptibility;ProbabilityDistributionFunction {type:LogNormal;mean: 0.00089   ;stddev: 0.00006;}}

            LithologyProperty {lithology:"d_BASE_ISA_SUPERBASIN";property:Susceptibility;ProbabilityDistributionFunction {type:LogNormal;mean: 0.005    ;stddev: 0.00013;}}

            LithologyProperty {lithology:"c_CALVERT_SUPERBASIC";property:Susceptibility;ProbabilityDistributionFunction {type:LogNormal;mean: 0.0002     ;stddev: 0.00001;}}

            LithologyProperty {lithology:"b_BASE_LEICHHARDT_SUPERBASIN";property:Susceptibility;ProbabilityDistributionFunction {type:LogNormal;mean: 0.0500   ;stddev: 0.00001; percentage: 66};ProbabilityDistributionFunction2 {type:LogNormal;mean: 0.0800   ;stddev: 0.00001; percentage: 34}}

            LithologyProperty {lithology:"a_BASEMENT";property:Susceptibility;ProbabilityDistributionFunction {type:LogNormal;mean: 0.00001  ;stddev: 0.000005;}}

            LithologyProperty {lithology:"i_YELDHAM_GRANITE";property:Susceptibility;ProbabilityDistributionFunction {type:LogNormal;mean: 0.01     ;stddev: 0.00001;}}

            LithologyProperty {lithology:"j_WEBBERA_GRANITE";property:Susceptibility;ProbabilityDistributionFunction {type:LogNormal;mean: 0.00001  ;stddev: 0.00001;}}

            LithologyProperty {lithology:"h_CAMBRIAN";fixed:false;}

            LithologyProperty {lithology:"g_SOUTH_NICHOLSON_GROUP";fixed:false;}

            LithologyProperty {lithology:"f_BASE_TERMITE";fixed:false;}

            LithologyProperty {lithology:"e_BASE_RIVERSLEIGH";fixed:false;}

            LithologyProperty {lithology:"d_BASE_ISA_SUPERBASIN";fixed:false;}

            LithologyProperty {lithology:"c_CALVERT_SUPERBASIC";fixed:false;}

            LithologyProperty {lithology:"b_BASE_LEICHHARDT_SUPERBASIN";fixed:false;}

            LithologyProperty {lithology:"a_BASEMENT";fixed:false;}

            LithologyProperty {lithology:"i_YELDHAM_GRANITE";fixed:false;}

            LithologyProperty {lithology:"j_WEBBERA_GRANITE";fixed:false;}

            LithologyProperty {lithology:"h_CAMBRIAN";preserve: {node:false node:false node:false node:true node:true node:true node:true node:true node:true node:true;}}

            LithologyProperty {lithology:"g_SOUTH_NICHOLSON_GROUP";preserve: {node:false node:false node:false node:false node:true node:true node:true node:true node:true node:true;}}

            LithologyProperty {lithology:"f_BASE_TERMITE";preserve: {node:false node:false node:false node:false node:false node:true node:true node:true node:true node:true;}}

            LithologyProperty {lithology:"e_BASE_RIVERSLEIGH";preserve: {node:false node:false node:false node:false node:false node:false node:true node:true node:true node:true;}}

            LithologyProperty {lithology:"d_BASE_ISA_SUPERBASIN";preserve: {node:false node:false node:false node:false node:false node:false node:false node:true node:true node:true;}}

            LithologyProperty {lithology:"c_CALVERT_SUPERBASIC";preserve: {node:false node:false node:false node:false node:false node:false node:false node:false node:true node:true;}}

            LithologyProperty {lithology:"b_BASE_LEICHHARDT_SUPERBASIN";preserve: {node:false node:false node:false node:false node:false node:false node:false node:false node:false node:true;}}

            LithologyProperty {lithology:"a_BASEMENT";preserve: {node:false node:false node:false node:false node:false node:false node:false node:false node:false node:false;}}

            LithologyProperty {lithology:"i_YELDHAM_GRANITE";preserve: {node:false node:false node:true node:true node:true node:true node:true node:true node:true node:true;}}

            LithologyProperty {lithology:"j_WEBBERA_GRANITE";preserve: {node:false node:true node:true node:true node:true node:true node:true node:true node:true node:true;}}

            LithologyProperty {lithology:"h_CAMBRIAN";property: CommonalityVolume;ProbabilityDistributionFunction {type:Normal;mean: 0.0;stddev: 1.;}}

            LithologyProperty {lithology:"g_SOUTH_NICHOLSON_GROUP";property: CommonalityVolume;ProbabilityDistributionFunction {type:Normal;mean: 0.0;stddev: 1.;}}

            LithologyProperty {lithology:"f_BASE_TERMITE";property: CommonalityVolume;ProbabilityDistributionFunction {type:Normal;mean: 0.0;stddev: 1.;}}

            LithologyProperty {lithology:"e_BASE_RIVERSLEIGH";property: CommonalityVolume;ProbabilityDistributionFunction {type:Normal;mean: 0.0;stddev: 1.;}}

            LithologyProperty {lithology:"d_BASE_ISA_SUPERBASIN";property: CommonalityVolume;ProbabilityDistributionFunction {type:Normal;mean: 0.0;stddev: 1.;}}

            LithologyProperty {lithology:"c_CALVERT_SUPERBASIC";property: CommonalityVolume;ProbabilityDistributionFunction {type:Normal;mean: 0.0;stddev: 1.;}}

            LithologyProperty {lithology:"b_BASE_LEICHHARDT_SUPERBASIN";property: CommonalityVolume;ProbabilityDistributionFunction {type:Normal;mean: 0.0;stddev: 1.;}}

            LithologyProperty {lithology:"a_BASEMENT";property: CommonalityVolume;ProbabilityDistributionFunction {type:Normal;mean: 0.0;stddev: 1.;}}

            LithologyProperty {lithology:"i_YELDHAM_GRANITE";property: CommonalityVolume;ProbabilityDistributionFunction {type:Normal;mean: 0.0;stddev: 1.;}}

            LithologyProperty {lithology:"j_WEBBERA_GRANITE";property: CommonalityVolume;ProbabilityDistributionFunction {type:Normal;mean: 0.0;stddev: 1.;}}

            LithologyProperty {lithology:"h_CAMBRIAN";property: Commonality;ProbabilityDistributionFunction {type:Weibull;mean: 0.001;stddev: 1.;}}

            LithologyProperty {lithology:"g_SOUTH_NICHOLSON_GROUP";property: Commonality;ProbabilityDistributionFunction {type:Weibull;mean: 0.001;stddev: 1.;}}

            LithologyProperty {lithology:"f_BASE_TERMITE";property: Commonality;ProbabilityDistributionFunction {type:Weibull;mean: 0.001;stddev: 1.;}}

            LithologyProperty {lithology:"e_BASE_RIVERSLEIGH";property: Commonality;ProbabilityDistributionFunction {type:Weibull;mean: 0.001;stddev: 1.;}}

            LithologyProperty {lithology:"d_BASE_ISA_SUPERBASIN";property: Commonality;ProbabilityDistributionFunction {type:Weibull;mean: 0.001;stddev: 1.;}}

            LithologyProperty {lithology:"c_CALVERT_SUPERBASIC";property: Commonality;ProbabilityDistributionFunction {type:Weibull;mean: 0.001;stddev: 1.;}}

            LithologyProperty {lithology:"b_BASE_LEICHHARDT_SUPERBASIN";property: Commonality;ProbabilityDistributionFunction {type:Weibull;mean: 0.001;stddev: 1.;}}

            LithologyProperty {lithology:"a_BASEMENT";property: Commonality;ProbabilityDistributionFunction {type:Weibull;mean: 0.001;stddev: 1.;}}

            LithologyProperty {lithology:"i_YELDHAM_GRANITE";property: Commonality;ProbabilityDistributionFunction {type:Weibull;mean: 0.001;stddev: 1.;}}

            LithologyProperty {lithology:"j_WEBBERA_GRANITE";property: Commonality;ProbabilityDistributionFunction {type:Weibull;mean: 0.001;stddev: 1.;}}

            LithologyProperty {lithology:"h_CAMBRIAN";property: VolumeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.01;}}

            LithologyProperty {lithology:"g_SOUTH_NICHOLSON_GROUP";property: VolumeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.01;}}

            LithologyProperty {lithology:"f_BASE_TERMITE";property: VolumeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.01;}}

            LithologyProperty {lithology:"e_BASE_RIVERSLEIGH";property: VolumeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.01;}}

            LithologyProperty {lithology:"d_BASE_ISA_SUPERBASIN";property: VolumeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.01;}}

            LithologyProperty {lithology:"c_CALVERT_SUPERBASIC";property: VolumeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.01;}}

            LithologyProperty {lithology:"b_BASE_LEICHHARDT_SUPERBASIN";property: VolumeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.01;}}

            LithologyProperty {lithology:"a_BASEMENT";property: VolumeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.01;}}

            LithologyProperty {lithology:"i_YELDHAM_GRANITE";property: VolumeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.01;}}

            LithologyProperty {lithology:"j_WEBBERA_GRANITE";property: VolumeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.01;}}

            LithologyProperty {lithology:"h_CAMBRIAN";property: ShapeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.05;}}
    
            LithologyProperty {lithology:"g_SOUTH_NICHOLSON_GROUP";property: ShapeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.05;}}
    
            LithologyProperty {lithology:"f_BASE_TERMITE";property: ShapeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.05;}}

            LithologyProperty {lithology:"e_BASE_RIVERSLEIGH";property: ShapeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.05;}}

            LithologyProperty {lithology:"d_BASE_ISA_SUPERBASIN";property: ShapeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.05;}}

            LithologyProperty {lithology:"c_CALVERT_SUPERBASIC";property: ShapeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.05;}}

            LithologyProperty {lithology:"b_BASE_LEICHHARDT_SUPERBASIN";property: ShapeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.05;}}

            LithologyProperty {lithology:"a_BASEMENT";property: ShapeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.05;}}

            LithologyProperty {lithology:"i_YELDHAM_GRANITE";property: ShapeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.05;}}

            LithologyProperty {lithology:"j_WEBBERA_GRANITE";property: ShapeRatio;ProbabilityDistributionFunction {type:LogNormal;mean: 0;stddev: 0.05;}}

        }

    }

}

InversionTask {

    NewRun {

        filename: "OutputGeology/imported.xml";

        case: "Case_Magnetic";

        run: "Run_Magnetic";

        NumberOfIterations:100000;

    }

}

InversionTask {

    RunControl {

        filename: "OutputGeology/imported.xml";

        case: "Case_Magnetic";

        run: "Run_Magnetic";

        NumberOfIterations:4000000;

        AllowNeighbourPropertyDifferentCheck: false;

        PropertyChangeOnFixedCell: true;

        DensityInitialiseStyle: FromLawMean;

        SusceptibilityInitialiseStyle: FromLawMean;

        RemenantMagnetizationInitialiseStyle: FromLawMean;

        LithologyInitialiseStyle: FromVoxet;

        property_voxet: "Model/Lawn_Hill_Voxet_subsampledxy.vo";

        PreserveTopology: false;

        ProbabilityOfPropertyChangeOnly: 100;

        SetLargeMisfitRMSThresholdFactor:1.0;

        SetLargeMisfitTemperature:1.0;

        SetDualTemperatureMode: false;

        Mode1TemperatureFactor: 1.0;  # if dual temperature mode is true

        Mode1Duration: 100.0;

        Mode2TemperatureFactor: 500.0;

        Mode2Duration: 4500.0;

        AllowGeophysicalTests: true;

            \#CommonalityVolumeTest {

            \#SetWeights: none;

        \#}

            \#VolumeRatioTest {

            \#SetWeights: LocalAbundance;

        \#}

            \#CommonalityTest {

            \#SetWeights: none;

        \#}

            \#ShapeRatioTest {

            \#SetWeights:LocalAbundance;

            \#IncludeAboveTopoFaces: true;

            \#IncludeOutsideFaces: true;

        \#}

        SetOutputMisfits: true;

        PreserveVerticalRelationshipTest: false;

        \#Seed:1135314798;

    }

}

InversionTask {

    RunControl {

        filename: "OutputGeology/imported.xml";

        case: "Case_Magnetic";

        run: "Run_Magnetic_Forward";

        LithologyInitialiseStyle: FromVoxet;

        property_voxet: "Model/Lawn_Hill_Voxet_subsampledxy.vo";

    }

}

InversionTask {

    Run {

        filename: "OutputGeology/imported.xml";

        case: "Case_Magnetic";

        run: "Run_Magnetic";

    }

}

InversionTask {

    MakeSuperSummaryStats {

        filename: "OutputGeology/imported.xml";

        case: "Case_Magnetic";

        RunInterval {

            run: "Run_Gravity";

            start_iterations: 400000;

            end_iterations: 4000000;

        }

        threshold_percentage: 90;

        statistical_voxet: "Results/StatsReport/Case_Magnetic_SuperSummary_90";

    }

}

InversionTask {

    MakeEvolutionMovie {

        filename: "OutputGeology/imported.xml";

        case: "Case_Magnetic";

            RunInterval {

            run: "Run_Magnetic";

            start_iterations: 400000;

            end_iterations: 4000000;

        }

        Section: "MYSECTION";

        MaxPixelCols: 900;

        MaxPixelRows: 900;

        VerticalExaggeration: 3;

        QualityPercentage: 100;

        FramesPerSecond: 1;

        IterationIncrement: 1000;

        lut: LithologyColours;

        Scaling {

            doLog: false;

        }

        show: Geology;

        movie: "Results/Movies/EvolutionMovie_Case_Magnetic_Lithology_MYSECTION";

    }

}
Labels:
None
Enter labels to add to this page:
Please wait 
Looking for a label? Just start typing.