Add script to calculate mean, stdev for continuous 3D field parameter…

…s and estimated facies probabilities for discrete 3D parameter

pyproject.toml

@@ -117,6 +117,8 @@ sw_model_utilities = "subscript.sw_model_utilities.sw_model_utilities:main"
 sunsch = "subscript.sunsch.sunsch:main"
 vfp2csv = "subscript.vfp2csv.vfp2csv:main"
 welltest_dpds = "subscript.welltest_dpds.welltest_dpds:main"
+field_statistics = "subscript.field_statistics.field_statistics:main"
+
 
 [project.entry-points.ert]
 subscript_jobs = "subscript.hook_implementations.jobs"

src/subscript/field_statistics/WF_FIELD_PARAM_STATISTICS

@@ -0,0 +1,21 @@
+-- Workflow job for ERT to calculate
+-- - mean and stdev of ensemble of continuous 3D parameters with name <name> saved for geogrid 
+--   under <ensemble_path>/realization-*/iter-*/share/results/grids/geogrid--<name>.roff
+-- - estimate facies probabilities of discrete 3D parameters with name <name> saved for geogrid
+--   under <ensemble_path>/realization-*/iter-*/share/results/grids/geogrid--<name>.roff
+INTERNAL   False
+EXECUTABLE  ../scripts/field_statistics.py
+
+MIN_ARG   6
+ARG_TYPE    0   STRING
+ARG_TYPE    1   STRING
+ARG_TYPE    2   STRING
+ARG_TYPE    3   STRING
+ARG_TYPE    4   STRING
+ARG_TYPE    5   STRING
+ARG_TYPE    6   STRING
+ARG_TYPE    7   STRING
+
+
+
+

src/subscript/field_statistics/field_param_stat.yml

@@ -0,0 +1,59 @@
+# Configuration file for script wf_field_param_statistics.py
+field_stat: 
+    # Number of realizations for specified ensemble
+    # Required.
+    nreal: 100
+
+    # Iteration numbers from ES-MDA in ERT (iteration = 0 is initial ensemble,
+    # usually iteration=3 is final updated ensemble)
+    # Required.
+    iterations: [0, 3]
+
+    # Selected set of zone names to use in calculations of statistics.
+    # Must be one or more of the defined zones.
+    # Require at least one zone to be selected.
+    use_zones: ["Valysar", "Therys", "Volon"]
+
+    # Zone numbers with zone name dictionary
+    zone_code_names: 
+        1: "Valysar"
+        2: "Therys"
+        3: "Volon"
+
+    # For each zone specify either Proportional, Top_conform or Base_conform
+    # as grid conformity.
+    # Conformity can be checked by opening the RMS job that has created
+    # the geogrid and check the grid settings for grid layers.
+    # Proportional means that number of layers is specified.
+    # Top or base conform means that grid cell thickness is specified.
+    # Required (but only for zones you want to use)
+    zone_conformity:
+        "Valysar": "Proportional"
+        "Therys": "Top_conform"
+        "Volon": "Proportional"
+
+    # For each zone specify which discrete parameter to use to calculate
+    # facies probability estimates.
+    # Possible names are those found in the
+    # share/results/grids/geogrid--<name>.roff
+    # files that are of discrete type.
+    # This key can be omitted or some of the lines specifying parameters
+    # for a zone if you don't want to use it.
+    discrete_property_param_per_zone: 
+        "Valysar": ["facies"]
+        "Therys": ["facies"]
+        "Volon": ["facies"]
+
+    # For each zone specify which continuous parameter to use to
+    # calculate estimate of mean and stdev over ensemble.
+    # Possible names are those found in the
+    #  share/results/grids/geogrid--<name>.roff
+    # files that are of continuous type
+    # This key can be omitted or some of the lines specifying
+    # parameters for a zone if you don't want to use it.
+    continuous_property_param_per_zone:
+        "Valysar": ["phit"]
+        "Therys": ["phit"]
+
+    # Size of ertbox grid for (nx, ny, nz)
+    ertbox_size: [92, 146, 66]