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


def _compute_mean(C, mag, repi, hypo_depth):
    """
    Compute mean value for MSK-64.
    """
    return C['a1'] * mag + C['a2'] + _get_term01(C, repi, hypo_depth)


def _get_term01(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


[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'} #: Required distance repi REQUIRES_DISTANCES = {'repi'} fixedh = None
[docs] def compute(self, ctx: np.recarray, imts, mean, sig, tau, phi): """ See :meth:`superclass method <.base.GroundShakingIntensityModel.compute>` for spec of input and result values. """ for m, imt in enumerate(imts): C = self.COEFFS[imt] if self.fixedh: ctx = ctx.copy() ctx.hypo_depth = np.full_like(ctx.hypo_depth, self.fixedh) mean[m] = _compute_mean(C, ctx.mag, ctx.repi, ctx.hypo_depth) sig[m] = C['sigma']
#: 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'} fixedh = 15. #: 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 """)