Source code for openquake.hazardlib.gsim.bozorgnia_campbell_2016_vh

# -*- coding: utf-8 -*-
# vim: tabstop=4 shiftwidth=4 softtabstop=4
#
# Copyright (C) 2014-2022 GEM Foundation
#
# OpenQuake 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.
#
# OpenQuake 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 OpenQuake. If not, see <http://www.gnu.org/licenses/>.

"""
Module exports :class:`BozorgniaCampbell2016VH`
               :class:`BozorgniaCampbell2016HighQVH`
               :class:`BozorgniaCampbell2016LowQVH`
               :class:`BozorgniaCampbell2016AveQJapanSiteVH`
               :class:`BozorgniaCampbell2016HighQJapanSiteVH`
               :class:`BozorgniaCampbell2016LowQJapanSiteVH`
"""
import numpy as np
from openquake.hazardlib.gsim.base import GMPE, CoeffsTable
from openquake.hazardlib import const, contexts
from openquake.hazardlib.imt import PGA, PGV, SA
from openquake.hazardlib.gsim import bozorgnia_campbell_2016
from openquake.hazardlib.gsim import campbell_bozorgnia_2014


def _get_tau_vh(C, mag, tau_v, tau_h):
    """
    Returns the inter-event random effects coefficient (tau) defined in
    Equation 10.
    """
    rhob1 = C['rhob1']
    rhob2 = C['rhob2']

    if mag <= 4.5:
        rhob = rhob1
    elif 4.5 < mag < 5.5:
        rhob = rhob2 + (rhob1 - rhob2)*(5.5 - mag)
    else:
        rhob = rhob2

    return np.sqrt(tau_v ** 2 + tau_h ** 2 - 2 * rhob * tau_v * tau_h)


def _get_phi_vh(C, mag, phi_v, phi_h):
    """
    Returns the intra-event random effects coefficient (phi) defined in
    Equation 11.
    """
    rhow1 = C['rhow1']
    rhow2 = C['rhow2']

    if mag <= 4.5:
        rhow = rhow1
    elif 4.5 < mag < 5.5:
        rhow = rhow2 + (rhow1 - rhow2) * (5.5 - mag)
    else:
        rhow = rhow2

    return np.sqrt(phi_v ** 2 + phi_h ** 2 - 2 * rhow * phi_v * phi_h)


[docs]class BozorgniaCampbell2016VH(GMPE): """ Implements the GMPE by Bozorgnia & Campbell (2016) vertical-to-horizontal ratio for ground motions from the PEER NGA-West2 Project This V/H model is combined from VGMPE by Bozorgnia and Campbell (2016) as the vertical model, and HGMPE by Campbell and Bozorgnia (2014) as the horizontal model. **Reference:** Bozorgnia, Y. & Campbell, K. (2016). Ground Motion Model for the Vertical-to-Horizontal (V/H) Ratios of PGA, PGV, and Response Spectra *Earthquake Spectra*, 32(2), 951-978. Implements the global model that uses datasets from California, Taiwan, the Middle East, and other similar active tectonic regions to represent a typical or average Q region. Applies the average attenuation case (Dc20=0) """ VGMPE = bozorgnia_campbell_2016.BozorgniaCampbell2016() HGMPE = campbell_bozorgnia_2014.CampbellBozorgnia2014() #: Supported tectonic region type is active shallow crust DEFINED_FOR_TECTONIC_REGION_TYPE = const.TRT.ACTIVE_SHALLOW_CRUST #: Supported intensity measure types are spectral acceleration, peak #: ground velocity and peak ground acceleration DEFINED_FOR_INTENSITY_MEASURE_TYPES = {PGA, PGV, SA} #: Supported intensity measure component is the #: :attr:`~openquake.hazardlib.const.IMC.VERTICAL_TO_HORIZONTAL_RATIO` DEFINED_FOR_INTENSITY_MEASURE_COMPONENT = ( const.IMC.VERTICAL_TO_HORIZONTAL_RATIO) #: Supported standard deviation types are inter-event, intra-event #: and total; see the section for "Aleatory Variability Model". DEFINED_FOR_STANDARD_DEVIATION_TYPES = { const.StdDev.TOTAL, const.StdDev.INTER_EVENT, const.StdDev.INTRA_EVENT} #: Required site parameters are taken from the V and H models REQUIRES_SITES_PARAMETERS = ( VGMPE.REQUIRES_SITES_PARAMETERS | HGMPE.REQUIRES_SITES_PARAMETERS) #: Required rupture parameters are taken from the V and H models REQUIRES_RUPTURE_PARAMETERS = ( VGMPE.REQUIRES_RUPTURE_PARAMETERS | HGMPE.REQUIRES_RUPTURE_PARAMETERS) #: Required distance measures are taken from the V and H models REQUIRES_DISTANCES = ( VGMPE.REQUIRES_DISTANCES | HGMPE.REQUIRES_DISTANCES)
[docs] def compute(self, ctx, imts, mean, sig, tau, phi): """ See :meth:`superclass method <.base.GroundShakingIntensityModel.compute>` for spec of input and result values. """ mean_, sig_, tau_, phi_ = contexts.get_mean_stds( [self.VGMPE, self.HGMPE], ctx, imts) for m, imt in enumerate(imts): # V/H model, Equation 1 and 12 (in natural log units) mean[m] = mean_[0, m] - mean_[1, m] # Get standard deviations C = self.COEFFS[imt] t = _get_tau_vh(C, ctx.mag, tau_[0, m], tau_[1, m]) p = _get_phi_vh(C, ctx.mag, phi_[0, m], phi_[1, m]) sig[m] = np.sqrt(t ** 2 + p ** 2) tau[m] = t phi[m] = p
#: Table of regression coefficients obtained from supplementary material #: published together with the EQS paper COEFFS = CoeffsTable(sa_damping=5, table="""\ IMT rhow1 rhow2 rhob1 rhob2 0.010 0.783 0.718 0.916 0.893 0.020 0.785 0.718 0.917 0.891 0.030 0.784 0.703 0.919 0.884 0.050 0.782 0.678 0.931 0.877 0.075 0.781 0.681 0.933 0.909 0.100 0.778 0.657 0.944 0.900 0.150 0.775 0.653 0.946 0.875 0.200 0.774 0.630 0.946 0.826 0.250 0.770 0.642 0.945 0.784 0.300 0.757 0.658 0.949 0.819 0.400 0.750 0.661 0.953 0.719 0.500 0.742 0.666 0.955 0.716 0.750 0.730 0.688 0.965 0.718 1.000 0.721 0.684 0.967 0.765 1.500 0.707 0.663 0.967 0.796 2.000 0.691 0.664 0.966 0.799 3.000 0.652 0.629 0.956 0.822 4.000 0.669 0.613 0.945 0.839 5.000 0.665 0.586 0.921 0.860 7.500 0.586 0.573 0.898 0.685 10.00 0.639 0.547 0.854 0.720 PGA 0.782 0.720 0.915 0.893 PGV 0.754 0.680 0.882 0.699 """)
[docs]class BozorgniaCampbell2016HighQVH(BozorgniaCampbell2016VH): """ Implements the GMPE by Bozorgnia & Campbell (2016) vertical-to-horizontal ratio for ground motions from the PEER NGA-West2 Project Applies regional corrections in path scaling term for regions with low attenuation (high quality factor, Q) (e.g. eastern China) """ VGMPE = bozorgnia_campbell_2016.BozorgniaCampbell2016(sgn=+1) HGMPE = campbell_bozorgnia_2014.CampbellBozorgnia2014HighQ()
[docs]class BozorgniaCampbell2016LowQVH(BozorgniaCampbell2016VH): """ Implements the GMPE by Bozorgnia & Campbell (2016) vertical-to-horizontal ratio for ground motions from the PEER NGA-West2 Project Applies regional corrections in path scaling term for regions with high attenuation (low quality factor, Q) (e.g. Japan and Italy) """ VGMPE = bozorgnia_campbell_2016.BozorgniaCampbell2016(sgn=-1) HGMPE = campbell_bozorgnia_2014.CampbellBozorgnia2014LowQ()
[docs]class BozorgniaCampbell2016AveQJapanSiteVH(BozorgniaCampbell2016VH): """ Implements the GMPE by Bozorgnia & Campbell (2016) vertical-to-horizontal ratio for ground motions from the PEER NGA-West2 Project Incorporates the difference in linear Vs30 scaling for sites in Japan by activating the flag variable in shallow site reponse scaling Applies the average attenuation case (Dc20=0) """ VGMPE = bozorgnia_campbell_2016.BozorgniaCampbell2016(SJ=1) HGMPE = campbell_bozorgnia_2014.CampbellBozorgnia2014JapanSite()
[docs]class BozorgniaCampbell2016HighQJapanSiteVH( BozorgniaCampbell2016AveQJapanSiteVH): """ Implements the GMPE by Bozorgnia & Campbell (2016) vertical-to-horizontal ratio for ground motions from the PEER NGA-West2 Project Incorporates the difference in linear Vs30 scaling for sites in Japan by activating the flag variable in shallow site reponse scaling Applies regional corrections in path scaling term for regions with low attenuation (high quality factor, Q) """ VGMPE = bozorgnia_campbell_2016.BozorgniaCampbell2016(SJ=1, sgn=+1) HGMPE = campbell_bozorgnia_2014.CampbellBozorgnia2014HighQJapanSite()
[docs]class BozorgniaCampbell2016LowQJapanSiteVH( BozorgniaCampbell2016AveQJapanSiteVH): """ Implements the GMPE by Bozorgnia & Campbell (2016) vertical-to-horizontal ratio for ground motions from the PEER NGA-West2 Project Incorporates the difference in linear Vs30 scaling for sites in Japan by activating the flag variable in shallow site reponse scaling Applies regional corrections in path scaling term for regions with high attenuation (low quality factor, Q) """ VGMPE = bozorgnia_campbell_2016.BozorgniaCampbell2016(SJ=1, sgn=-1) HGMPE = campbell_bozorgnia_2014.CampbellBozorgnia2014LowQJapanSite()