def iter_ruptures(self, temporal_occurrence_model):
"""
See :meth:`nhlib.source.base.SeismicSource.iter_ruptures`.
Uses :func:`_float_ruptures` for finding possible rupture locations
on the whole fault surface.
"""
whole_fault_surface = ComplexFaultSurface.from_fault_data(
self.edges, self.rupture_mesh_spacing)
whole_fault_mesh = whole_fault_surface.get_mesh()
cell_center, cell_length, cell_width, cell_area = (
whole_fault_mesh.get_cell_dimensions())
for (mag, mag_occ_rate) in self.get_annual_occurrence_rates():
rupture_area = self.magnitude_scaling_relationship.get_median_area(
mag, self.rake)
rupture_length = numpy.sqrt(rupture_area *
self.rupture_aspect_ratio)
rupture_slices = _float_ruptures(rupture_area, rupture_length,
cell_area, cell_length)
occurrence_rate = mag_occ_rate / float(len(rupture_slices))
for rupture_slice in rupture_slices:
mesh = whole_fault_mesh[rupture_slice]
# XXX: use surface centroid as rupture's hypocenter
# XXX: instead of point with middle index
hypocenter = mesh.get_middle_point()
surface = ComplexFaultSurface(mesh)
yield ProbabilisticRupture(mag, self.rake,
self.tectonic_region_type,
hypocenter, surface, type(self),
occurrence_rate,
temporal_occurrence_model)
def test_mesh_spacing_more_than_two_lengths(self):
edge1 = Line([Point(0, 0, 0), Point(0, 0.1, 0)])
edge2 = Line([Point(0, 0, 10), Point(0, 0.1, 20)])
with self.assertRaises(ValueError) as ar:
ComplexFaultSurface.from_fault_data([edge1, edge2],
mesh_spacing=27)
self.assertEqual(
str(ar.exception),
'mesh spacing 27.0 km is too big for mean length 13.0 km'
)
def __init__(self, source_id, name, tectonic_region_type,
mfd, rupture_mesh_spacing,
magnitude_scaling_relationship, rupture_aspect_ratio,
# complex fault specific parameters
edges, rake):
super(ComplexFaultSource, self).__init__(
source_id, name, tectonic_region_type, mfd, rupture_mesh_spacing,
magnitude_scaling_relationship, rupture_aspect_ratio
)
NodalPlane.check_rake(rake)
ComplexFaultSurface.check_fault_data(edges, rupture_mesh_spacing)
self.edges = edges
self.rake = rake
def test_2(self):
edge1 = Line([Point(0, 0, 1), Point(0, 0.02, 1)])
edge2 = Line([Point(0.02, 0, 0.5), Point(0.02, 0.01, 0.5)])
edge3 = Line([Point(0, 0, 2), Point(0, 0.02, 2)])
surface = ComplexFaultSurface.from_fault_data([edge1, edge2, edge3],
mesh_spacing=1)
self.assert_mesh_is(surface=surface, expected_mesh=[
[(0.00000000e+00, 0.00000000e+00, 1.00000000e+00),
(0.00000000e+00, 1.00000000e-02, 1.00000000e+00),
(0.00000000e+00, 2.00000000e-02, 1.00000000e+00)],
[(8.70732572e-03, 5.33152318e-19, 7.82316857e-01),
(8.67044753e-03, 7.83238825e-03, 7.83238813e-01),
(8.57984833e-03, 1.57100762e-02, 7.85503798e-01)],
[(1.74146514e-02, 1.06630462e-18, 5.64633714e-01),
(1.73408950e-02, 5.66477632e-03, 5.66477626e-01),
(1.71596963e-02, 1.14201520e-02, 5.71007595e-01)],
[(1.47979150e-02, 3.18525318e-19, 8.90156376e-01),
(1.48515919e-02, 6.28710212e-03, 8.86130608e-01),
(1.49871315e-02, 1.25064345e-02, 8.75965152e-01)],
[(7.39895749e-03, 7.71565830e-19, 1.44507819e+00),
(7.42579599e-03, 8.14355113e-03, 1.44306530e+00),
(7.49356586e-03, 1.62532174e-02, 1.43798258e+00)],
[(0.00000000e+00, 0.00000000e+00, 2.00000000e+00),
(0.00000000e+00, 1.00000000e-02, 2.00000000e+00),
(0.00000000e+00, 2.00000000e-02, 2.00000000e+00)],
])
def iter_ruptures(self, temporal_occurrence_model):
"""
See :meth:`nhlib.source.base.SeismicSource.iter_ruptures`.
Uses :func:`_float_ruptures` for finding possible rupture locations
on the whole fault surface.
"""
whole_fault_surface = ComplexFaultSurface.from_fault_data(
self.edges, self.rupture_mesh_spacing
)
whole_fault_mesh = whole_fault_surface.get_mesh()
cell_center, cell_length, cell_width, cell_area = (
whole_fault_mesh.get_cell_dimensions()
)
for (mag, mag_occ_rate) in self.get_annual_occurrence_rates():
rupture_area = self.magnitude_scaling_relationship.get_median_area(
mag, self.rake
)
rupture_length = numpy.sqrt(rupture_area
* self.rupture_aspect_ratio)
rupture_slices = _float_ruptures(rupture_area, rupture_length,
cell_area, cell_length)
occurrence_rate = mag_occ_rate / float(len(rupture_slices))
for rupture_slice in rupture_slices:
mesh = whole_fault_mesh[rupture_slice]
# XXX: use surface centroid as rupture's hypocenter
# XXX: instead of point with middle index
hypocenter = mesh.get_middle_point()
surface = ComplexFaultSurface(mesh)
yield ProbabilisticRupture(
mag, self.rake, self.tectonic_region_type, hypocenter,
surface, occurrence_rate, temporal_occurrence_model
)
def test(self):
edges = [Line([Point(-3.4, 5.5, 0), Point(-3.9, 5.3, 0)]),
Line([Point(-2.4, 4.6, 10), Point(-3.6, 4.9, 20)])]
polygon = ComplexFaultSurface.surface_projection_from_fault_data(edges)
elons = [-2.4, -3.6, -3.9, -3.4]
elats = [4.6, 4.9, 5.3, 5.5]
numpy.testing.assert_allclose(polygon.lons, elons)
numpy.testing.assert_allclose(polygon.lats, elats)
def test_surface_with_variable_width(self):
edges = [Line([Point(0.1, 0.1, 0),
Point(0.459729190252, 0.0999980290582, 0.0),
Point(0.819458380482, 0.0999960581553, 0.0)]),
Line([Point(0.1, 0.1, 20.0),
Point(0.459729190252, 0.0999980290582, 20.0),
Point(0.819458380482, 0.0999960581553, 43.0940107676)])]
surface = ComplexFaultSurface.from_fault_data(edges, mesh_spacing=5.0)
self.assertAlmostEqual(surface.get_width(), 26.0, places=0)
def test_edges_with_nonuniform_length(self):
edges = [Line([Point(1, -20, 30), Point(1, -20.2, 30),
Point(2, -19.7, 30), Point(3, -19.5, 30)]),
Line([Point(1, -20, 50), Point(1, -20.2, 50),
Point(2, -19.7, 50)])]
polygon = ComplexFaultSurface.surface_projection_from_fault_data(edges)
elons = [1, 1, 2, 3]
elats = [-20.2, -20., -19.7, -19.5]
numpy.testing.assert_allclose(polygon.lons, elons)
numpy.testing.assert_allclose(polygon.lats, elats)
def test_dip_90_self_intersection(self):
edges = [Line([Point(1, -2, 10), Point(2, -1.9, 10),
Point(3, -2.1, 10), Point(4, -2, 10)]),
Line([Point(1, -2, 20), Point(2, -1.9, 20),
Point(3, -2.1, 20), Point(4, -2, 20)])]
polygon = ComplexFaultSurface.surface_projection_from_fault_data(edges)
elons = [3., 1., 2., 4.]
elats = [-2.1, -2., -1.9, -2.]
numpy.testing.assert_allclose(polygon.lons, elons)
numpy.testing.assert_allclose(polygon.lats, elats)
def test_dip_90_two_points(self):
edges = [Line([Point(2, 2, 10), Point(1, 1, 10)]),
Line([Point(2, 2, 20), Point(1, 1, 20)])]
polygon = ComplexFaultSurface.surface_projection_from_fault_data(edges)
elons = [1.00003181, 0.99996821, 0.99996819, 1.99996819, 2.00003182,
2.00003181]
elats = [0.99996822, 0.99996819, 1.00003178, 2.0000318, 2.0000318,
1.9999682]
numpy.testing.assert_allclose(polygon.lons, elons)
numpy.testing.assert_allclose(polygon.lats, elats)
def test_dip_90_three_points(self):
edges = [Line([Point(1, -20, 30), Point(1, -20.2, 30),
Point(2, -19.7, 30)]),
Line([Point(1, -20, 50), Point(1, -20.2, 50),
Point(2, -19.7, 50)])]
polygon = ComplexFaultSurface.surface_projection_from_fault_data(edges)
elons = [1, 1, 2]
elats = [-20.2, -20., -19.7]
numpy.testing.assert_allclose(polygon.lons, elons)
numpy.testing.assert_allclose(polygon.lats, elats)
def test_1(self):
edge1 = Line([Point(0, 0), Point(0.03, 0)])
edge2 = Line([Point(0, 0, 2.224), Point(0.03, 0, 2.224)])
surface = ComplexFaultSurface.from_fault_data([edge1, edge2],
mesh_spacing=1.112)
self.assertIsInstance(surface, ComplexFaultSurface)
self.assert_mesh_is(surface=surface, expected_mesh=[
[(0, 0, 0), (0.01, 0, 0), (0.02, 0, 0), (0.03, 0, 0)],
[(0, 0, 1.112), (0.01, 0, 1.112),
(0.02, 0, 1.112), (0.03, 0, 1.112)],
[(0, 0, 2.224), (0.01, 0, 2.224),
(0.02, 0, 2.224), (0.03, 0, 2.224)],
])
def __init__(
self,
source_id,
name,
tectonic_region_type,
mfd,
rupture_mesh_spacing,
magnitude_scaling_relationship,
rupture_aspect_ratio,
# complex fault specific parameters
edges,
rake):
super(ComplexFaultSource,
self).__init__(source_id, name, tectonic_region_type, mfd,
rupture_mesh_spacing,
magnitude_scaling_relationship,
rupture_aspect_ratio)
NodalPlane.check_rake(rake)
ComplexFaultSurface.check_fault_data(edges, rupture_mesh_spacing)
self.edges = edges
self.rake = rake
def get_rupture_enclosing_polygon(self, dilation=0):
"""
Uses :meth:`nhlib.geo.surface.complex_fault.ComplexFaultSurface.surface_projection_from_fault_data`
for getting the fault's surface projection and then calls
its :meth:`~nhlib.geo.polygon.Polygon.dilate` method passing
in ``dilation`` parameter.
See :meth:`superclass method
<nhlib.source.base.SeismicSource.get_rupture_enclosing_polygon>`
for parameter and return value definition.
"""
polygon = ComplexFaultSurface.surface_projection_from_fault_data(
self.edges)
if dilation:
return polygon.dilate(dilation)
else:
return polygon
def get_rupture_enclosing_polygon(self, dilation=0):
"""
Uses :meth:`nhlib.geo.surface.complex_fault.ComplexFaultSurface.surface_projection_from_fault_data`
for getting the fault's surface projection and then calls
its :meth:`~nhlib.geo.polygon.Polygon.dilate` method passing
in ``dilation`` parameter.
See :meth:`superclass method
<nhlib.source.base.SeismicSource.get_rupture_enclosing_polygon>`
for parameter and return value definition.
"""
polygon = ComplexFaultSurface.surface_projection_from_fault_data(
self.edges
)
if dilation:
return polygon.dilate(dilation)
else:
return polygon
def __init__(self, identifier, name, tectonic_region, aspect_ratio,
edges, mesh_spacing = 1.0):
'''Instantiate class with just the basic attributes
:param identifier:
Integer ID code for the source
:param name:
Source Name (string)
:param tectonic_region:
Tectonic Region Type (String)
:param aspect_ratio:
Ratio of along-strike length to down-dip width (float)
:param edges:
List of edges([Long, Lat, Depth],[Long,Lat, Depth])
:param mesh_spacing:
Spacing of mesh (km) (Non-negative Float)
'''
self.typology = 'ComplexFault'
self.source_id = identifier
self.name = name
self.tectonic_region_type = tectonic_region
self.aspect_ratio = aspect_ratio
self.mfd = None
self.msr = None
number_edges = len(edges)
self.upper_depth = np.inf
self.lower_depth = -np.inf
self.traces = []
for iloc, edge in enumerate(edges):
# Convert edge to line
if np.min(edge[:, 2]) < self.upper_depth:
self.upper_depth = np.copy(np.min(edge[:, 2]))
self.lower_depth = np.copy(self.upper_depth)
if np.max(edge[:, 2]) > self.lower_depth:
self.lower_depth = np.copy(np.max(edge[:, 2]))
self.traces.append(
Line([Point(node[0], node[1], node[2])
for node in edge]))
self.geometry = ComplexFaultSurface.from_fault_data(
self.traces, mesh_spacing)
self.catalogue = None