RU2510053C1 - Способ динамической оценки сейсмической опасности - Google Patents
Способ динамической оценки сейсмической опасности Download PDFInfo
- Publication number
- RU2510053C1 RU2510053C1 RU2012144388/28A RU2012144388A RU2510053C1 RU 2510053 C1 RU2510053 C1 RU 2510053C1 RU 2012144388/28 A RU2012144388/28 A RU 2012144388/28A RU 2012144388 A RU2012144388 A RU 2012144388A RU 2510053 C1 RU2510053 C1 RU 2510053C1
- Authority
- RU
- Russia
- Prior art keywords
- active region
- seismic
- network
- earth
- stations
- Prior art date
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Изобретение относится к области геофизики и может быть использовано для прогнозирования места и тренда (увеличения или уменьшения) сейсмической опасности. Сущность: осуществляют мониторинг ситуации, по крайней мере, в одной зоне ожидаемого сейсмического события, принадлежащей исследуемому сейсмоактивному региону. Формируют в сейсмоактивном регионе наблюдательную сеть из «n» пунктов, разнесенных друг от друга. Одновременно и непрерывно измеряют контролируемый параметр, характеризующий процессы в Земной коре, во всех пунктах наблюдательной сети. Определяют область с повышенной сейсмической активностью по результатам сравнения измеренного контролируемого параметра с пороговым значением, определяемым на основе статистического анализа значений контролируемого параметра для предыдущих сейсмических событий в сейсмоактивном регионе. При этом измерение контролируемого параметра на всех «n» пунктах наблюдательной сети осуществляют с постоянным и одинаковым для всех станций шагом дискретизации по времени Δt и регистрируют его в виде электрического сигнала. Формируют для исследуемого сейсмоактивного региона регулярную сеть, причем каждому из узлов сети принадлежит прилегающая к нему зона исследуемого сейсмоактивного региона. Выбирают временное окно, осуществляют обработку электрических сигналов, полученных от указанных «n» пунктов. На основе указанных сигналов вычисляют одновременно во всех пунктах измерения для каждого узла регулярной сетки в указанном временном окне медианы нормализованной энтропии шума по некоторому числу изменений контролируемого параметра. На основе полученных результатов
Description
Изобретение относится к геофизике и может быть использовано для прогнозирования места и тренда (увеличения или уменьшения) сейсмической опасности.
Известен способ прогнозирования землетрясений, основанный на совместной оценке геофизических, сейсмологических и гидрологических данных (см. F.A.McKeown, S.F.Diehl. Evidence of Contemporary and Ancient Excess Fluid Pressure in the New Madrid Seismic Zone of the Reelfoot Rift, Central United States. U.S.Geological Survey Professional Paper 1538-N, Washington, 1994), согласно которому о наличии сейсмически активных зон судят по результатам регистрации избыточного давления флюидов, которое должно превышать гидростатическое давление в горных породах, причем условием однозначной оценки наличия сейсмически активной зоны является характер структурных разрушений горных пород в районе, где зарегистрировано избыточное давление флюидов.
Недостаток данного способа заключается в том, что он выделяет сейсмоактивную область по ретроспективным данным и не обеспечивает возможности достоверного прогноза и локализации очаговой области ввиду того, что не осуществляется непрерывное измерение флюидной динамики и не обеспечивается достоверность получаемых оценок.
Наиболее близким аналогом является способ оперативного прогноза землетрясений, включающий проведение синхронных измерений интенсивности естественных импульсных электромагнитных полей Земли (ЕИЭМПЗ) в нескольких пунктах контролируемого региона, при этом измерения в каждой точке ведут не менее чем в двух различных направлениях приема сигналов, а чувствительность регистрирующих станций выбирают в соответствии с местными геофизическими условиями таким образом, чтобы регистрируемая станциями интенсивность ЕИЭМПЗ была близка по своим значениям к интенсивности типичного суточного хода ЕИЭМПЗ, выделяют аномальную территорию по наличию скачкообразного изменения интенсивности (Патент RU 2238575 С2, МПК8 G01V 9/00, опубл. 20.10.2004).
Недостатком данного способа является низкая точность прогнозирования предстоящих землетрясений в сейсмоактивных зона.
Задачей настоящего изобретения является разработка способа динамического прогнозирования сейсмического события, например землетрясения или извержения вулкана, который даст возможность на основе измерений процессов, происходящих в земной коре в сейсмоопасной зоне, определить область, в которой может произойти землетрясение, определить тренд (уменьшение или увеличение) сейсмической опасности.
Техническим результатом является повышение точности предсказания зоны предстоящего землетрясения и возможность оценки тренд увеличения или уменьшения сейсмической опасности.
Указанный технический результат обеспечивается тем, что в способе динамической оценки сейсмической опасности, в котором осуществляют мониторинг ситуации, по крайней мере, в одной зоне ожидаемого сейсмического события, принадлежащей исследуемому сейсмоактивному региону, формируют в сейсмоактивном регионе наблюдательную сеть из «n» пунктов, разнесенных друг от друга, одновременно и непрерывно измеряют, по крайней мере, один контролируемый параметр во всех пунктах наблюдательной сети, характеризующий процессы в земной коре, определяют область с повышенной сейсмической активностью по результатам сравнения, по крайне мере, одного измеренного контролируемого параметра, характеризующего процессы в Земной коре, с пороговым значением, определяемым на основе статистического анализа значений контролируемого параметра для предыдущих сейсмических событий в сейсмоактивном регионе, измерение, по крайней мере, одного контролируемого параметра, характеризующего процессы в земной коре, на всех «n» пунктах наблюдательной сети осуществляют с постоянным и одинаковым для всех станций шагом дискретизации по времени At и регистрируют их в виде электрического сигнала, формируют для исследуемого сейсмоактивного региона регулярную сеть из N×M узлов, где N - количество узлов в направлении изменения долготы и М - количество узлов в направлении изменения широты, покрывающей исследуемый сейсмоактивный регион, причем каждому из этих узлов принадлежит прилегающая к нему зона исследуемого сейсмоактивного региона, выбирают временное окно, состоящее из заданного числа L, где L≥64, последовательных временных отсчетов с шагом дискретизации At, осуществляют обработку электрических сигналов, полученных от указанных «n» пунктов, на основе которых вычисляют одновременно во всех пунктах измерения для каждого узла (i,j), где ≤i≤N, ≤j≤М, регулярной сетки в указанном временном окне медианы нормализованной энтропии
шума по некоторому числу v, где v≥3 изменений контролируемого, по крайней мере, одного параметра, характеризующего процессы в земной коре, после удаления тренда полиномом заданного порядка, одинакового для всех пунктов измерения, по формуле
где
cj - вейвлет-коэффициенты ортогонального вейвлета из некоторого набора (словаря) базисов, найденного из условия минимума величины En, на основе которых строят матрицу значений медиан нормализованной энтропии, соответствующих указанному текущему временному окну, визуализируют ее как карту, при этом область с повышенной сейсмической активностью определяют как совокупность зон, прилегающих к узлам (i,j) регулярной сети, для которых нормализованная энтропия
превышает пороговое значение En, при этом оценивают силу будущего землетрясения по размеру области с повышенной нормализованной энтропией.
Кроме того, в способе динамической оценки сейсмической опасности в качестве контролируемого параметра используют параметр из списка: микросейсмические колебания, изменение напряженности магнитного поля Земли, изменение напряженности электромагнитного поля Земли, колебания уровня подземных вод, колебания давления подземных вод, концентрацию одного или нескольких растворенных в подземных водах газов, например углекислого газа, метана, водорода, радона.
На фиг.1 представлены карты распределения значений нормализованной энтропии для двух промежутков времени, указанных в подрисуночной подписи.
На фиг.2 - усредненные карты распределения нормализованной энтропии волновых форм сейсмического шума на Японских островах для 4-х временных фрагментов примерно одинаковой длины после японского мегаземлетрясения 11 марта 2011 года.
Способ динамической оценки сейсмической опасности осуществляется следующим образом.
Выбирают, по меньшей мере, один контролируемый параметр, из числа параметров, характеризующих процессы в земной коре, для мониторинга ситуации, по меньшей мере, в одной зоне ожидаемого сейсмического события, принадлежащей исследуемому сейсмоактивному региону. Формируют в исследуемом сейсмоактивном регионе наблюдательную сеть из «n» пунктов измерения, по меньшей мере, этого одного контролируемого параметра, при этом в исследуемом сейсмоактивном регионе пункты измерения контролируемого параметра разнесены относительно друг друга, Измерения контролируемого параметра выполняются непрерывно и одновременно на всех «n» пунктах измерения наблюдательной сети с постоянным и одинаковым для всех станций шагом дискретизации по времени Δt и регистрируют их в виде электрического сигнала. Формируют для исследуемого сейсмоактивного региона регулярную сеть из N×M узлов, где N - количество узлов в направлении изменения долготы и М - количество узлов в направлении изменения широты, покрывающей исследуемый сейсмоактивный регион, причем каждому из этих узлов соответствует прилегающая к нему одна зона исследуемого сейсмоопасного региона. Осуществляют обработку данных, полученных от указанных «n» пунктов измерения, для определения риска появления предстоящего сейсмического события, его местоположения, времени и силы, включающую выбор временного окна, состоящего из заданного числа L последовательных временных отсчетов с шагом дискретизации Δt, причем L≥64 и определение для каждого узла (i,j), 1≤i≤N, 1≤j≤М, регулярной сетки из N×М узлов медианы нормализованной энтропии
шума изменений контролируемых параметров по некоторому числу v (v≥3) ближайших к узлу (i,j) пунктов измерения среди общего числа n пунктов измерения, покрывающих исследуемый сейсмоактивный регион.
Нормализованную энтропию шума En в каждом пункте измерения вычисляют по формуле:
и где
, cj - вейвлет-коэффициенты ортогонального вейвлета из некоторого набора (словаря) базисов, найденного из условия минимума величины En; набор ортогональных вейвлетов (словарь базисов) состоит из 17 вейвлетов Добеши: 10 обычных вейвлетов с числом обнуляемых моментов от 1 до 10 и 7 так называемых симлетов Добеши с числом обнуляемых моментов от 4 до 10.
Величины нормализованной энтропии En вычисляют одновременно во всех пунктах измерения в одном и том же временном окне длиной L отсчетов после удаления тренда полиномом заданного порядка, одного и того же для всех пунктов измерения, что обеспечивает выделение шума сигналов, получаемых от сети геофизического мониторинга.
Таким образом, после вычисления медианы величин En от v ближайших к каждому узлу (i,j) станций получается матрица
значений медиан нормализованных энтропии, соответствующих текущему временному окну длиной L отсчетов, которую визуализируют как карту.
Совокупность зон, прилегающих к узлам (i,j) регулярной сети, для которых нормализованная энтропия En превышает пороговое значение En* En(i,j)≥En*, определяет подобласть исследуемого региона, оцениваемую как область с повышенной сейсмоопасностью в пределах текущего скользящего временного окна, причем пороговое значение определяют на основе статистического анализа значения En для предыдущих сейсмических событий в этом сейсмоактивном регионе, при этом силу будущего землетрясения оценивают по размеру области с повышенной нормализованной энтропии.
Предпочтительно в качестве контролируемых параметров используют, по крайней мере, один или несколько из следующих параметров: микросейсмические колебания, изменение напряженности магнитного поля Земли, изменение напряженности электромагнитного поля Земли, колебания уровня подземных вод, колебания давления подземных вод, концентрацию одного или нескольких растворенных в подземных водах газов, например углекислого газа, метана, водорода, радона, при этом при обработке данных, полученных от указанных пунктов измерения, используют измерения одного из указанных контрольных параметров, либо при обработке данных, полученных от указанных n пунктов измерения, используют измерения разных указанных контрольных параметров.
В качестве примера реализации изобретения рассмотрим выполненный автором мониторинг ситуации в сейсмоопасном регионе - японских островах в регионе от 30° до 46° с.ш. и от 128° до 146° в.д. В данном регионе имеется развитая сеть пунктов измерения контрольных параметров, характеризующих процессы в земной коре. В качестве контрольного параметра использовали микросейсмические колебания, но можно использовать и иные параметры: изменение напряженности магнитного поля Земли, изменение напряженности электромагнитного поля Земли, колебания уровня подземных вод, колебания давления подземных вод, концентрацию одного или нескольких растворенных в подземных водах газов, например углекислого газа, метана, водорода, радона. Применительно к исследуемому сейсмоопасному региону сформировали регулярную сеть, состоящую из 30×30 узлов, каждому из которых соответствует прилегающая к нему зона сейсмоопасного региона, имеющая те же размеры.
Для мониторинга состояния среды в районе Японии использовались низкочастотные микросейсмические колебания, содержащие информацию о процессах в земной коре. Фактически земная кора является средой распространения колебаний при воздействии на нее атмосферных и океанических процессов. Поскольку передаточные свойства коры зависят от ее состояния, можно ожидать, что статистические свойства микросейсмических колебаний отражают изменения свойств литосферы.
Данные широкополосной сейсмической сети F-net свободно доступны в Интернете по адресу http://www.fnet.bosai.go.ip/top.php?LANG=en. Общее число станций равно 83. Наблюдения ведутся с 1997 года по настоящее время. Анализируемые данные - вертикальные компоненты с шагом по времени 1 сек, которые преобразовывались к шагу по времени 1 минута путем вычисления средних значений в последовательных временных фрагментах длиной 60 значений. Далее рассматривались лишь станции, расположенные выше 30° с.ш., что исключает из анализа данные 6 уединенных станций, расположенных на удаленных небольших островах.
Карты строились путем усреднения ежесуточных карт внутри указанных промежутков времени. Каждая суточная карта нормализованной энтропии вычислялась как матрица медианных значений
для каждого узла (i,j), 1≤i≤N, 1≤j≤M, регулярной сетки из N×M узлов, N=30, M=30. Медианы брались по 5 станциям (v=5), ближайшим к каждому узлу (i,j).
Нормализованные энтропии шума вычислялись в последовательных окнах длиной 1 сутки (1440 минут, то есть L=1440) после удаления в каждом окне тренда полиномом 8-го порядка.
При этом силу будущего землетрясения оценивают по размеру области с повышенной нормализованной энтропии.
На фиг. 1 видно, что область подготовки Великого Японского землетрясения 11 марта 2011 года до события 25 сентября 2003 года представляла собой единую область повышенных значений нормализованной энтропии шума. Однако после 25 сентября 2003 года эта область распалась на 2 части, причем Северная часть реализовалась как область мегаземлетрясения.
Из фиг.2 видно, что область мегаземлетрясения 11 марта 2011 года после события стала характеризоваться относительно низкими значениями нормализованной энтропии шума, тогда как южная область (желоб Нанкай) по-прежнему характеризуется как область высоких значений энтропии шума, причем сами значения нормализованной энтропии там прогрессивно увеличиваются. Это говорит в пользу гипотезы, что в этой области готовится повторное мегаземлетрясение, которое может представлять большую опасность для мегаполиса Токио.
Claims (2)
1. Способ динамической оценки сейсмической опасности, заключающийся в том, что осуществляют мониторинг ситуации, по крайней мере, в одной зоне ожидаемого сейсмического события, принадлежащей исследуемому сейсмоактивному региону, формируют в сейсмоактивном регионе наблюдательную сеть из «n» пунктов, разнесенных друг от друга, одновременно и непрерывно измеряют, по крайней мере, один контролируемый параметр во всех пунктах наблюдательной сети, характеризующий процессы в Земной коре, определяют область с повышенной сейсмической активностью по результатам сравнения, по крайне мере, одного измеренного контролируемого параметра, характеризующего процессы в Земной коре, с пороговым значением, определяемым на основе статистического анализа значений контролируемого параметра для предыдущих сейсмических событий в сейсмоактивном регионе, отличающийся тем, что измерение, по крайней мере, одного контролируемого параметра, характеризующего процессы в земной коре, на всех «n» пунктах наблюдательной сети осуществляют с постоянным и одинаковым для всех станций шагом дискретизации по времени Δt и регистрируют его в виде электрического сигнала, формируют для исследуемого сейсмоактивного региона регулярную сеть из N×M узлов, где N - количество узлов в направлении изменения долготы и М - количество узлов в направлении изменения широты, покрывающей исследуемый сейсмоактивный регион, причем каждому из этих узлов принадлежит прилегающая к нему зона исследуемого сейсмоактивного региона, выбирают временное окно, состоящее из заданного числа L, где L≥64, последовательных временных отсчетов с шагом дискретизации Δt, осуществляют обработку электрических сигналов, полученных от указанных «n» пунктов, на основе которых вычисляют одновременно во всех пунктах измерения для каждого узла (i, j), где 1≤i≤N, 1≤j≤М, регулярной сетки в указанном временном окне медианы нормализованной энтропии
шума по некоторому числу v, где v≥3 изменений контролируемого, по крайней мере, одного параметра, характеризующего процессы в Земной коре, после удаления тренда полиномом заданного порядка, одинакового для всех пунктов измерения, по формуле
где
cj - вейвлет-коэффициенты ортогонального вейвлета из некоторого набора (словаря) базисов, найденного из условия минимума величины En, на основе которых строят матрицу значений медиан нормализованной энтропии, соответствующих указанному текущему временному окну, визуализируют ее как карту, при этом область с повышенной сейсмической активностью определяют как совокупность зон, прилегающих к узлам (i, j) регулярной сети, для которых нормализованная энтропия
превышает пороговое значение Еn*.
где
2. Способ динамической оценки сейсмической опасности по п.1, отличающийся тем, что в качестве контролируемого параметра используют параметры из списка: микросейсмические колебания, изменение напряженности магнитного поля Земли, изменение напряженности электромагнитного поля Земли, колебания уровня подземных вод, колебания давления подземных вод, концентрацию одного или нескольких растворенных в подземных водах газов, например углекислого газа, метана, водорода, радона.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2012144388/28A RU2510053C1 (ru) | 2012-10-18 | 2012-10-18 | Способ динамической оценки сейсмической опасности |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2012144388/28A RU2510053C1 (ru) | 2012-10-18 | 2012-10-18 | Способ динамической оценки сейсмической опасности |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2510053C1 true RU2510053C1 (ru) | 2014-03-20 |
Family
ID=50279745
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2012144388/28A RU2510053C1 (ru) | 2012-10-18 | 2012-10-18 | Способ динамической оценки сейсмической опасности |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2510053C1 (ru) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106155512A (zh) * | 2015-04-20 | 2016-11-23 | 中兴通讯股份有限公司 | 结点位置确定方法及装置 |
CN112305606A (zh) * | 2020-10-16 | 2021-02-02 | 宁夏回族自治区地震局 | 一种基于自然正交函数展开的地震活动场分析方法 |
CN116381773A (zh) * | 2023-04-13 | 2023-07-04 | 南京云创大数据科技股份有限公司 | 一种在地震预报中归一化混杂数据的方法及装置 |
RU2805275C1 (ru) * | 2022-09-05 | 2023-10-13 | Федеральное государственное бюджетное учреждение науки Институт земной коры Сибирского отделения Российской академии наук (ИЗК СО РАН) | Способ краткосрочного определения подготовки сильного сейсмического события |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2102780C1 (ru) * | 1996-12-05 | 1998-01-20 | Явдат Накипович Хамидуллин | Способ контроля землетрясений |
RU2238575C2 (ru) * | 2002-06-13 | 2004-10-20 | Институт мониторинга климатических и экологических систем СО РАН | Способ прогноза землетрясений |
RU2279108C1 (ru) * | 2004-11-19 | 2006-06-27 | Игорь Павлович Шестопалов | Способ определения сейсмической опасности |
-
2012
- 2012-10-18 RU RU2012144388/28A patent/RU2510053C1/ru not_active IP Right Cessation
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2102780C1 (ru) * | 1996-12-05 | 1998-01-20 | Явдат Накипович Хамидуллин | Способ контроля землетрясений |
RU2238575C2 (ru) * | 2002-06-13 | 2004-10-20 | Институт мониторинга климатических и экологических систем СО РАН | Способ прогноза землетрясений |
RU2279108C1 (ru) * | 2004-11-19 | 2006-06-27 | Игорь Павлович Шестопалов | Способ определения сейсмической опасности |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106155512A (zh) * | 2015-04-20 | 2016-11-23 | 中兴通讯股份有限公司 | 结点位置确定方法及装置 |
CN112305606A (zh) * | 2020-10-16 | 2021-02-02 | 宁夏回族自治区地震局 | 一种基于自然正交函数展开的地震活动场分析方法 |
US20220120926A1 (en) * | 2020-10-16 | 2022-04-21 | Earthquake Agency Of Ningxia Hui Autonomous Region | Method and system for analyzing seismic active field based on expansion of empirical orthogonal function |
US11650341B2 (en) * | 2020-10-16 | 2023-05-16 | Earthquake Agency Of Ningxia Hui Autonomous Region | Method and system for analyzing seismic active field based on expansion of empirical orthogonal function |
RU2805275C1 (ru) * | 2022-09-05 | 2023-10-13 | Федеральное государственное бюджетное учреждение науки Институт земной коры Сибирского отделения Российской академии наук (ИЗК СО РАН) | Способ краткосрочного определения подготовки сильного сейсмического события |
CN116381773A (zh) * | 2023-04-13 | 2023-07-04 | 南京云创大数据科技股份有限公司 | 一种在地震预报中归一化混杂数据的方法及装置 |
CN116381773B (zh) * | 2023-04-13 | 2024-07-12 | 南京云创大数据科技股份有限公司 | 一种在地震预报中归一化混杂数据的方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Akkar et al. | Reference database for seismic ground-motion in Europe (RESORCE) | |
Wurman et al. | Toward earthquake early warning in northern California | |
Inbal et al. | Sources of long‐range anthropogenic noise in Southern California and implications for tectonic tremor detection | |
RU2581119C2 (ru) | Способ прогнозирования сейсмического события | |
US20210072413A1 (en) | Real-time array-based seismic source location | |
JP5507903B2 (ja) | 震度推定方法及び装置 | |
RU2510053C1 (ru) | Способ динамической оценки сейсмической опасности | |
Amezquita-Sanchez et al. | Current efforts for prediction and assessment of natural disasters: Earthquakes, tsunamis, volcanic eruptions, hurricanes, tornados, and floods | |
Wang et al. | CAV site-effect assessment: A case study of Taipei Basin | |
KR101914657B1 (ko) | 지진배경잡음을 이용한 지진신호의 위상 및 진폭정보 산출방법 | |
Romanowicz et al. | SeaFOAM: A year‐long DAS deployment in Monterey Bay, California | |
RU2346300C1 (ru) | Способ обнаружения возможности наступления катастрофических явлений | |
Lockman et al. | Magnitude-period scaling relations for Japan and the Pacific Northwest: Implications for earthquake early warning | |
CN112180444A (zh) | 一种地层速度结构的探测方法、装置以及存储介质 | |
Liseikin et al. | Monitoring of the natural frequencies of Chirkey arch dam | |
Alizadeh Zakaria et al. | Investigation of the application of geospatial artificial intelligence for integration of earthquake precursors extracted from remotely sensed SAR and thermal images for earthquake prediction | |
Kalimeris et al. | Multi-spectral detection of statistically significant components in pre-seismic electromagnetic emissions related with Athens 1999, M= 5.9 earthquake | |
Liu et al. | Multiscale random field-based shear wave velocity mapping and site classification | |
Warnana et al. | Assessment to the soil-structure resonance using microtremor analysis on pare-east Java, Indonesia | |
Nguyen et al. | Characteristics of earthquake source and ground motions in Northern Vietnam investigated through the 2020 Moc Chau M5. 0 earthquake sequence | |
Rulenko et al. | An analysis of the relationships between high-frequency geoacoustic emissions and the electrical field in the atmosphere near the ground surface | |
RU2466432C1 (ru) | Способ обнаружения возможности наступления катастрофических явлений | |
Anbazhagan | Site characterization and seismic hazard analysis with local site effects for microzonation of Bangalore | |
Partheniu et al. | Monitoring the black sea natural hazards using new technology and equipment | |
Telesca et al. | Investigating the time-correlation properties in self-potential signals recorded in a seismic area of Irpinia, southern Italy |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20181019 |