Source code for openquake.hazardlib.gsim.bindi_2011_ipe

# coding: utf-8
# The Hazard Library
# Copyright (C) 2012 GEM Foundation
#
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU Affero General Public License as
# published by the Free Software Foundation, either version 3 of the
# License, or (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU Affero General Public License for more details.
#
# You should have received a copy of the GNU Affero General Public License
# along with this program.  If not, see <http://www.gnu.org/licenses/>.
"""
Module exports :
class:`BindiEtAl2011Repi`,
class:`BindiEtAl2011RepiFixedH`,

"""
from __future__ import division
import numpy as np
from openquake.hazardlib.gsim.base import CoeffsTable, GMPE
from openquake.hazardlib import const
from openquake.hazardlib.imt import MMI


[docs]class BindiEtAl2011Repi(GMPE): """ Implements IPE developed by Dino Bindi et al. 2011 and published as "Intensity prediction equations for Central Asia" (Geo-physical journal international, 2011, 187,327-337). Model implemented by laurentiu.danciu@gmail.com """ DEFINED_FOR_TECTONIC_REGION_TYPE = const.TRT.ACTIVE_SHALLOW_CRUST DEFINED_FOR_INTENSITY_MEASURE_TYPES = {MMI} DEFINED_FOR_INTENSITY_MEASURE_COMPONENT = const.IMC.HORIZONTAL DEFINED_FOR_STANDARD_DEVIATION_TYPES = {const.StdDev.TOTAL} REQUIRES_SITES_PARAMETERS = set() REQUIRES_RUPTURE_PARAMETERS = {'mag', 'hypo_depth'} REQUIRES_DISTANCES = {'repi'}
[docs] def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types): """ See :meth:`superclass method <.base.GroundShakingIntensityModel.get_mean_and_stddevs>` for spec of input and result values. """ # extract dictionaries of coefficients specific to required # intensity measure type C = self.COEFFS[imt] mean = self._compute_mean(C, rup.mag, dists.repi, rup.hypo_depth) stddevs = self._get_stddevs(C, stddev_types, num_sites=dists.repi.shape) return mean, stddevs
def _compute_mean(self, C, mag, repi, hypo_depth): """ Compute mean value for MSK-64. """ return (C['a1'] * mag + C['a2'] + self._get_term01(C, repi, hypo_depth)) def _get_term01(self, C, repi, hypo_depth): h = hypo_depth term_repi = np.sqrt((repi**2+h**2)/h**2) term_h = np.sqrt(repi**2+h**2)-h return -C['a3']*np.log10(term_repi)-(C['a4']*term_h) def _get_stddevs(self, C, stddev_types, num_sites): """ Return total standard deviation. """ stddevs = [] for stddev_type in stddev_types: assert stddev_type in self.DEFINED_FOR_STANDARD_DEVIATION_TYPES stddevs.append(C['sigma'] + np.zeros(num_sites)) return stddevs #: Coefficient table constructed from the electronic suplements of the #: original paper.Table 1 .page 331 COEFFS = CoeffsTable(table="""\ IMT a1 a2 a3 a4 sigma MMI 0.898 1.215 1.809 0.003447 0.737 """)
[docs]class BindiEtAl2011RepiFixedH(BindiEtAl2011Repi): """ Implements IPE developed by Dino Bindi et al. 2011 and published as "Intensity prediction equations for Central Asia" (Geo-physical journal international, 2011, 187,327-337). for a fixed depth of 15 km and epicentral distance (equation 5 in the paper) Implements the Repi with fixed depth at 15km /coeff on Table 1 Model implmented by laurentiu.danciu@gmail.com """ REQUIRES_SITES_PARAMETERS = set() REQUIRES_DISTANCES = {'repi'} REQUIRES_RUPTURE_PARAMETERS = {'mag', 'hypo_depth'}
[docs] def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types): """ Extends the Bindietal2011 IPE for Repi without a fixed depth of 15Km """ C = self.COEFFS[imt] rup.hypo_depth = 15. mean, stddevs = super(BindiEtAl2011RepiFixedH, self).\ get_mean_and_stddevs(sites, rup, dists, imt, stddev_types) return mean, stddevs
#: Coefficient table constructed from the electronic suplements of the #: original paper. COEFFS = CoeffsTable(table="""\ IMT a1 a2 a3 a4 sigma MMI 1.049 0.686 2.706 0.0001811 0.689 """)