NHLIB Integration: Implement Event Based calculator

Bug #1021338 reported by Damiano Monelli
6
This bug affects 1 person
Affects Status Importance Assigned to Milestone
OpenQuake (deprecated)
Fix Released
High
Lars Butler

Bug Description

With the stochastic event set calculator (see pr: https://github.com/gem/nhlib/pull/48) and ground motion field calculator (see pr: https://github.com/gem/nhlib/pull/47), NHLIB has all the core calculators for supporting in OpenQuake the Event Based calculator.

The configuration file for a `event based' calculation should be something similar to the following:

[general]

calculation_mode = event based
random_seed = 23

[geometry]

region = 6.5 45.8, 6.5 46.5, 8.5 46.5, 8.5 45.8
region_grid_spacing = 20.0

or
sites = ......

[logic_tree]

number_of_logic_tree_samples = 2

[erf]

rupture_mesh_spacing = 1
width_of_mfd_bin = 0.3
area_source_discretization = 10

[site_params]

reference_vs30_type = measured
reference_vs30_value = 760.0
reference_depth_to_2pt5km_per_sec = 5.0
reference_depth_to_1pt0km_per_sec = 100.0

or

site_model= ...

[calculation]

source_model_logic_tree_file = source_model_logic_tree.xml
gsim_logic_tree_file = gmpe_logic_tree.xml
investigation_time = 50.0
truncation_level = 3
maximum_distance = 200.0
number_of_stochastic_event_sets_per_logic_tree_sample: Integer, representing number of stochastic event sets to be generated for each logic tree sample
ground_motion_correlation_model: name of ground motion correlation model to be used in ground motion field calculation. It can be empty -- meaning that no correlation is taken into account
ground_motion_correlation_parameters: dictionary containing ground motion model parameter names, and values, needed to properly set the ground motion correlation model. Example {'vs30_clustering': False}. It can be empty meaning that the correlation model does not require extra parameters.

[output]

complete_logic_tree_ses: True or False. If true a stochastic event set obtained by joining the stochastic event sets for each sampled tree path and associated to a time span equal to:

investigation_time * number_of_stochastic_event_sets_per_logic_tree_sample * number_of_logic_tree_samples

is computed and saved.

ground_motion_fields: True or False. If True ground motion fields are saved.

The workflow for the event based calculator is as follows:

1) Given source model and gsim logic tree, n (=number_of_logic_tree_samples) logic tree paths are sampled and corresponding source models and gsims dictionaries (gsims = {'tectonic region type 1': gsim1, 'tectonic region type 2': gsim2, ...}) are computed.

2) For each tuple (source_model, gsim), the following calculations are done:

      2a) generate m (=number_of_stochastic_event_sets_per_logic_tree_sample) stochastic event sets from source_model

      2b) for each rupture in each stochastic event set, compute a ground motion field using the appropriate gsim in gsims.

3) For each logic tree sample, save the m stochastic even sets, and the ground motion fields generated by each rupture in each stochastic event set.

4) Give option to serialize stochastic even sets and ground motion fields to xml file.

The calculator should report to the user the status of the calculations (calculation done vs. calculation requested)

description: updated
description: updated
description: updated
Changed in openquake:
milestone: none → 0.8.2
assignee: nobody → Lars Butler (lars-butler)
importance: Undecided → High
status: New → Confirmed
description: updated
Changed in openquake:
status: Confirmed → In Progress
Changed in openquake:
status: In Progress → Fix Committed
Changed in openquake:
status: Fix Committed → Fix Released
To post a comment you must log in.
This report contains Public information  
Everyone can see this information.

Other bug subscribers

Remote bug watches

Bug watches keep track of this bug in other bug trackers.