Способ краткосрочного прогноза землетрясений
METHOD OF SHORT-TERM EARTHQUAKE FORECASTING

Изобретение относится к геофизике, а именно к способу краткосрочного прогноза землетрясений с использованием данных наземно-космического мониторинга.
8
Известен способ обнаружения радиационного загрязнения (патент RU 2561305 C2 МПК G01T 1/169, G01W 1/00, Опубликовано: 27.08.2015), заключающийся в дистанционной регистрации зон латентного тепла в приземном слое атмосферы, вызываемым ее нагревом радиоактивным излучением. При прогнозировании в одном масштабе для одной и той же местности создают карты распределения латентного тепла в атмосфере, полученные на основе анализа излучений в инфракрасном диапазоне частот 8-14 мкм, и карты распределения оценочных поправок к химическому потенциалу паров воды в атмосфере, полученные в результате измерений температуры и влажности в ее приземном слое на основе анализа радиоизлучений в диапазоне сантиметровых и миллиметровых волн, затем сравнивают данные по аномалиям к фону латентного тепла и аномалиям к фону оценочных поправок к химическому потенциалу паров воды, и места совпадения аномальных зон по обеим картам постулируют как места радиационных загрязнений.
Известное техническое решение может быть использовано для определения места прогнозируемого землетрясения, в случае, когда радиационное загрязнение обусловливается радоном, выходящим из поверхности Земли в атмосферу при подготовке землетрясения.
8
Известен способ краткосрочного прогнозирования землетрясений (патент RU 2395105 C1 МПК G01V 9/00, Опубликовано: 20.07.2010). Прогноз землетрясений включает измерение радиационного фона атмосферы над зоной подготавливаемого землетрясения, отслеживание динамики изменения сигнала предвестника на интервале его существования, прогноз характеристик сейсмического удара по параметрам регистрируемого сигнала. Для увеличения скорости деления радона используют активное зондирование атмосферы коллимированным пучком элементарных частиц с подвижного носителя, работающего в режиме сканирующего обзора и патрулирования контролируемой зоны. Синхронизируют тракты зондирования и приема путем введения порогового напряжения в тракт приема, осуществляют временную и энергетическую селекцию частиц мгновенной реакции радона на пучок зондирования.
При этом место ожидаемого землетрясения связывается с определяемой зоной аномального радиационного фона атмосферы, размеры которой определяются непосредственно из семейства регистрограмм, по пеленгу размеров зоны.
Недостатком известного технического решения является низкая достоверность прогнозирования места ожидаемого землетрясения, обусловленная свойством выделившегося из поверхности Земли радона, имеющего плотность, превышающую плотность воздуха, растекаться по поверхности в виде «стелющегося тумана».
8
y
y
2
1
3
2
0
1
2
3
1
2
3
Известен способ определения параметров землетрясения (патент РФ №2302020, МПК G01V 9/00, Опубликовано: 10.08.2006. Определение параметров землетрясения включает регистрацию показателей, характеризующих изображения подстилающей поверхности, в виде зависимости амплитуды сигнала A (x, y) от пространственных координат по нескольким каналам приема, совместную обработку изображений каналов приема, прогнозирование землетрясений по числовым характеристикам результирующего сигнала, осуществление визуальными наблюдениями с орбитальной станции обнаружения протяженных гряд облачных линейных аномалий и их одновременную съемку видеокамерой и гиперспектрометром на спектральных линиях атома водорода, отслеживание изменения суммы среднеквадратических значений сигналов в каналах приема на последовательных витках, определение постоянной времени T переходного процесса как , прогнозирование времени удара t≈4,7T и магнитуды удара (M) из соотношения lgt[сут]=0,54M-3,37, отождествление координаты очага землетрясения с координатами снимков облачных линейных аномалий, где Δt=t-t=t-t - интервал времени между последовательными съемками на витках; σ - предельная величина результирующего сигнала, при котором происходит удар; , где σ, σ, σ - сумма среднеквадратических значений сигналов в каналах приема в моменты измерений t, t, t.
Недостатком известного технического решения также является низкая достоверность прогнозирования времени ожидаемого землетрясения, обусловленная использованием статистической обработки данных для нестационарного и, в связи с этим, нестабильного и имеющего низкую воспроизводимость при проведении измерений переходного процесса.
Наиболее близким техническим решением является эмпирическая схема краткосрочного прогноза землетрясений (Дода Л.Н. Эмпирическая схема краткосрочного прогноза землетрясений / Л.Н. Дода, И.В. Степанов, В.Л. Натяганов // Доклады Академии наук: сер. Геофизика. - 2013. - T. 453 .- №5. - стр. 551-552).
Эмпирическая схема краткосрочного прогноза землетрясений заключается в том, что:
1) определяют прогнозную дату землетрясения исходя из того, что землетрясения происходят в среднем на 14-е или 21-е сутки после геоэффективных явлений на Солнце типа корональных выбросов массы или солнечных вспышек, по формуле:
s
где d - дата геоэффективного события, вызывающего на Земле геомагнитные возмущения, а значение n (число кэррингтоновых циклов Солнца) обычно равно 0, иногда 1 и еще реже 2;
2) определяют пространственное положение сейсмомагнитных меридианов, представляющих собой проекции на геоид наиболее возмущенной совокупности силовых трубок геомагнитного поля, спровоцированного геомагнитными возмущениями от геоэффективного явления на Солнце, по формуле:
s
s
*
s
где j=i+1 при k(i+1)-k(i)≥2; i=0…8; λ, ϕ - долгота и широта точки на сейсмомагнитном меридиане, λ - долгота подсолнечной точки из Астрономического ежегодника на дату d геоэффективного явления на Солнце; i=0 для долготы по Гринвичу (λ=0); j - порядковый номер особого трехчасового интервала, на котором выполняется условие возникновения цепочки возможных землетрясений (разность соседних значений k геомагнитных Kp индексов не менее 2); коэффициент γ=tg11°≈0,19, который берется со знаком "плюс" при |γϕ|<λ, где  - долгота сейсмомагнитного меридиана на экваторе; λ=λ-45°⋅j. При этом, так как неравенство k(i+1)-k(i)≥2 может быть выполнено для нескольких значений i из интервала от 0 до 8, то формула (2) может определять несколько сейсмомагнитных меридианов, соответствующих различным j;
3) определяют вероятные места возникновения землетрясений как окрестности пересечений сейсмомагнитных меридианов с литосферными разломами, которые совпадают с границами литосферных плит по модели Берда, накопившими на границах литосферных плит достаточный для сброса упругой энергии потенциал механических напряжений и деформаций.
При этом прогнозная зона определяется радиусом 700 км от геометрического места пересечения сейсмомагнитного меридиана с границами литосферных плит (Дода Л.Н. Результаты сейсмопрогнозного эксперимента по Камчатке / Л.Н. Дода, В.Л. Натяганов, С.А. Шопин // Доклады Академии наук: сер. Геофизика. - 2016. - Т. 469. - №6. - стр. 730);
4) определяют по космическим снимкам окрестность над пересечением сейсмомагнитных меридианов с границами литосферных плит, в которой присутствует облачная сейсмоиндуцированная структура, указывающая на место, где механические напряжения и деформации приближаются к пределу, за которым следует разрыв сплошности земных недр, т.е. происходит землетрясение;
4) определяют магнитуду прогнозируемого землетрясения исходя из определенного по космическому снимку максимального линейного размера облачной сейсмоиндуцированной структуры по формуле:
0
где L - максимальный линейный размер облачной сейсмоиндуцированной структуры, L=1 км.
Недостатками известного технического решения является низкая точность при определении места прогнозируемого землетрясения, ограниченного радиусом 700 км.
Задачей предложенного технического решения является повышение точности при определении места прогнозируемого землетрясения. Это достигается тем, что в отличии от известного технического решения определяют прогнозную дату землетрясения, определяют вероятные места возникновения землетрясения как окрестности радиусом 700 км от пересечений сейсмомагнитных меридианов с литосферными разломами, выявляют на космических снимках окрестность пересечения сейсмомагнитных меридианов с литосферными разломами над которой присутствует облачная сейсмоиндуцированная структура, определяют максимальный линейный размер облачной сейсмоиндуцированной структуры, определяют магнитуду землетрясения, в окрестности пересечения сейсмомагнитных меридианов с литосферными разломами над которой присутствует облачная сейсмоиндуцированная структура в одном масштабе с космическими снимками строят карты распределения поправки атмосферного химического потенциала,
при этом строят карты среднесуточного распределения поправки атмосферного химического потенциала в течение не более 18 суток, отсчитываемых от 5 суток, предшествующих прогнозной дате землетрясения, карты дополнительно покрывают регулярной сеткой, имеющей ячейки со сторонами 2 град. на 2 град., выявляют ячейки, содержащие локальные максимумы поправки атмосферного химического потенциала, суммируют число локальных максимумов в каждой ячейке сетки для каждой карты, строят карты распределения сумм локальных максимумов, строят область, охватывающую локальные максимумы на картах распределения сумм локальных максимумов, полученную область постулируют как наиболее вероятное место возникновения землетрясения.
Примеры выполнения предложенного технического решения поясняются чертежами. На Фиг. 1 представлены линии сейсмомагнитных меридианов на карте в проекции Меркатора и окрестности (круговые зоны) радиусом 700 км с центрами в точках пересечения сейсмомагнитных меридианов с границами литосферных плит, на Фиг. 2 - окрестности над пересечением сейсмомагнитных меридианов с границами литосферных плит, в которых присутствует L-образная облачная сейсмоиндуцированная структура, на Фиг. 3 - карта среднесуточного распределения поправки атмосферного химического потенциала на 27.08.2015 г., на Фиг. 4 - карта среднесуточного распределения поправки атмосферного химического потенциала на 28.08.2015 г., на Фиг.5 - карта среднесуточного распределения поправки атмосферного химического потенциала на 29.08.2015 г., на Фиг. 6 - карта среднесуточного распределения поправки атмосферного химического потенциала на 30.08.2015 г., на Фиг. 7 - карта среднесуточного распределения поправки атмосферного химического потенциала на 31.08.2015 г., на Фиг. 8 - карта среднесуточного распределения поправки атмосферного химического потенциала на 01.09.2015 г., на Фиг. 9 - карта среднесуточного распределения поправки атмосферного химического потенциала на 02.09.2015 г., на Фиг. 10 - карта среднесуточного распределения поправки атмосферного химического потенциала на 03.09.2015 г., на Фиг. 11 - карта среднесуточного распределения поправки атмосферного химического потенциала на 04.09.2015 г., на Фиг. 12 - карта среднесуточного распределения поправки атмосферного химического потенциала на 05.09.2015 г., на Фиг. 13 - карта среднесуточного распределения поправки атмосферного химического потенциала на 06.09.2015 г., на Фиг. 14 - карта среднесуточного распределения поправки атмосферного химического потенциала на 07.09.2015 г., на Фиг. 15 - карта среднесуточного распределения поправки атмосферного химического потенциала на 08.09.2015 г., на Фиг. 16 - карта среднесуточного распределения поправки атмосферного химического потенциала на 09.09.2015 г., на Фиг. 17 - карта среднесуточного распределения поправки атмосферного химического потенциала на 10.09.2015 г., на Фиг. 18 - карта среднесуточного распределения поправки атмосферного химического потенциала на 11.09.2015 г., на Фиг. 19 - карта среднесуточного распределения поправки атмосферного химического потенциала на 12.09.2015 г., на Фиг. 20 - карта среднесуточного распределения поправки атмосферного химического потенциала на 13.09.2015 г., на Фиг. 21 - карта суммированного в ячейках сетки числа локальных максимумов среднесуточного распределения поправки атмосферного химического потенциала за 27.08.2015 г., на Фиг. 22 - карта суммированного в ячейках сетки числа локальных максимумов среднесуточного распределения поправки атмосферного химического потенциала за 27.08.2015-28.08.2015 г.г., на Фиг. 23 - карта суммированного в ячейках сетки числа локальных максимумов среднесуточного распределения поправки атмосферного химического потенциала за 27.08.2015-29.08.2015 г.г., на Фиг. 24 - карта суммированного в ячейках сетки числа локальных максимумов среднесуточного распределения поправки атмосферного химического потенциала за 27.08.2015-30.08.201 г.г., на Фиг. 25 - карта суммированного в ячейках сетки числа локальных максимумов среднесуточного распределения поправки атмосферного химического потенциала за 27.08.2015-31.08.2015 г.г., на Фиг. 26 - карта суммированного в ячейках сетки числа локальных максимумов среднесуточного распределения поправки атмосферного химического потенциала за 27.08.2015-01.09.2015 г.г., на Фиг. 27 - карта суммированного в ячейках сетки числа локальных максимумов среднесуточного распределения поправки атмосферного химического потенциала за 27.08.2015-02.09.2015 г.г., на Фиг. 28 - карта суммированного в ячейках сетки числа локальных максимумов среднесуточного распределения поправки атмосферного химического потенциала за 27.08.2015-03.09.2015 г.г., на Фиг. 29 - карта суммированного в ячейках сетки числа локальных максимумов среднесуточного распределения поправки атмосферного химического потенциала за 27.08.2015-04.09.2015 г.г., на Фиг. 30 - карта суммированного в ячейках сетки числа локальных максимумов среднесуточного распределения поправки атмосферного химического потенциала за 27.08.2015-05.09.2015 г.г., на Фиг. 31 - карта суммированного в ячейках сетки числа локальных максимумов среднесуточного распределения поправки атмосферного химического потенциала за 27.08.2015-06.09.2015 г.г., на Фиг. 32 - карта суммированного в ячейках сетки числа локальных максимумов среднесуточного распределения поправки атмосферного химического потенциала за 27.08.2015-07.09.2015 г.г., на Фиг. 33 - карта суммированного в ячейках сетки числа локальных максимумов среднесуточного распределения поправки атмосферного химического потенциала за 27.08.2015-08.09.2015 г.г., на Фиг. 34 - карта суммированного в ячейках сетки числа локальных максимумов среднесуточного распределения поправки атмосферного химического потенциала за 27.08.2015-09.09.2015 г.г., на Фиг. 35 - карта суммированного в ячейках сетки числа локальных максимумов среднесуточного распределения поправки атмосферного химического потенциала за 27.08.2015-10.09.2015 г.г., на Фиг. 36 - карта суммированного в ячейках сетки числа локальных максимумов среднесуточного распределения поправки атмосферного химического потенциала за 27.08.2015-11.09.2015 г.г., на Фиг. 37 - карта суммированного в ячейках сетки числа локальных максимумов среднесуточного распределения поправки атмосферного химического потенциала за 27.08.2015-12.09.2015 г.г., на Фиг. 38 - карта суммированного в ячейках сетки числа локальных максимумов среднесуточного распределения поправки атмосферного химического потенциала за 27.08.2015-13.09.2015 г.г., на Фиг. 39 - область, охватывающую локальные максимумы на картах распределения сумм локальных максимумов поправки атмосферного химического потенциала, постулируемая как наиболее вероятное место возникновения землетрясения, на Фиг. 40 - положение эпицентра и прогнозные окрестности: 1 - положение эпицентра, на Фиг. 41 - фактическое событие и область, охватывающую локальные максимумы на картах распределения сумм локальных максимумов поправки атмосферного химического потенциала, постулируемая как наиболее вероятное место возникновения землетрясения: 1 - положение эпицентра.
Способ краткосрочного прогноза землетрясений осуществляется следующим образом.
Определяют пространственное положение сейсмомагнитных меридианов, представляющих собой проекции на геоид наиболее возмущенной совокупности силовых трубок геомагнитного поля, спровоцированного геомагнитными возмущениями от геоэффективного явления на Солнце, по формуле (2):
s
λ=±γϕ+(λ-45°⋅j),
s
s
*
*
s
где j=i+1 при k(i+1)-k(i)≥2; i=0…8; λ, ϕ - долгота и широта точки на сейсмомагнитном меридиане, λ - долгота подсолнечной точки из Астрономического ежегодника на дату d геоэффективного явления на Солнце; i=0 для долготы по Гринвичу (λ=0); j - порядковый номер особого трехчасового интервала, на котором выполняется условие возникновения цепочки возможных землетрясений (разность соседних значений k геомагнитных Kp индексов не менее 2); коэффициент γ=tg11°≈0,19, который берется со знаком "плюс" при |γϕ|<λ, где  - долгота сейсмомагнитного меридиана на экваторе; λ=λ-45°⋅j. При этом, так как неравенство k(i+1)-k(i-k)≥2 может быть выполнено для нескольких значений i из интервала от 0 до 8, то формула (2) может определять несколько сейсмомагнитных меридианов, соответствующих различным j.
Определяют вероятные места возникновения землетрясений как окрестностей радиусом 700 км от пересечений сейсмомагнитных меридианов с границами литосферных плит.
Выявляют на космических снимках окрестность над пересечением сейсмомагнитных меридианов с границами литосферных плит, в которой присутствует облачная сейсмоиндуцированная структура, указывающая на место, где механические напряжения и деформации приближаются к пределу, за которым следует разрыв сплошности земных недр, т.е. происходит землетрясение.
8
Облачные предвестники землетрясений, называемые как облачные линейные аномалии (патент РФ №2302020, МПК G0IV 9/00), Опубликовано: 10.08.2006), или облачные линейные структуры (Doda L.N. Earthquakes forecasts following space- and ground-based monitoring / L.N. Doda, V.R. Dushin, V.L. Natyaganov, N.N. Smirnov, I.V. Stepanov // Acta Astronautica. - 2011. - V. 69. - стр. 19), или сейсмоиндикаторы (Дода Л.Н. Космический мониторинг предвестников землетрясений / Л.Н. Дода, Л.А. Пахомов, И.В. Степанов // Наука в России. - 2009. - №6. - стр. 34), или облачные сейсмотектонические индикаторы (Степанов И.В. Подсистема ГИС для решения задач сейсмического мониторинга и прогнозирования // Геоматика. - 2010. - №3. - стр. 99), или облачные сейсмоиндуцированные структуры (Дода Л.Н. Эмпирическая схема краткосрочного прогноза землетрясений / Л.Н. Дода, И.В. Степанов, В.Л., Натяганов // Доклады Академии наук: сер. Геофизика - 2013. - Т. 453. - №5. - стр. 552), кроме линейной формы могут иметь углообразную форму, а также 3-й 4-угольную форму (Дода Л.Н. Космический мониторинг предвестников землетрясений / Л.Н. Дода, Л.А. Пахомов, И.В. Степанов // Наука в России. - 2009. - №6 - стр. 34).
Определяют максимальный линейный размер облачной сейсмоиндуцированной структуры.
Определяют магнитуду прогнозируемого землетрясения по формуле (3):
0
где L - максимальный линейный размер облачной сейсмоиндуцированной структуры, L=1 км.
Определяют прогнозную дату землетрясения исходя из того, что землетрясения происходят в среднем на 14-е или 21-е сутки после геоэффективных явлений на Солнце типа корональных выбросов массы или солнечных вспышек, по формуле (1):
s
где d - дата геоэффективного события, вызывающего на Земле геомагнитные возмущения, а значение n (число кэррингтоновых циклов Солнца) обычно равно 0, иногда 1 и еще реже 2;
Строят карты среднесуточного распределения поправки атмосферного химического потенциала в течение не более 18 суток, отсчитываемых от 5 суток, предшествующих прогнозной дате землетрясения.
Из формулы (1) следует, что наибольший прогнозируемый срок ожидаемого землетрясения составляет 21+2=23 суток. Как эмпирически установлено и показано в примере осуществления предлагаемого технического решения, в интервале времени за 5 суток от прогнозируемого землетрясения до даты землетрясения на картах среднесуточного распределения поправок атмосферного химического потенциала локальные максимумы (количество и расположение) практически не изменяются. Это делает малоинформативным построение карт среднесуточного распределения поправки атмосферного химического потенциала в интервале времени от 5 суток до даты возникновения прогнозируемого землетрясения.
Таким образом, промежуток времени, на котором целесообразно построение карт среднесуточного распределения поправки атмосферного химического потенциала составляет не более 23-5=18 суток, отсчитываемых от 5 суток, предшествующих прогнозной дате землетрясения.
Карты среднесуточного распределения поправки атмосферного химического потенциала дополнительно покрывают регулярной сеткой, имеющей ячейки со сторонами 2 град, на 2 град.
Выявляют ячейки, содержащие локальные максимумы поправки атмосферного химического потенциала.
Суммируют число локальных максимумов в каждой ячейке сетки для каждой карты среднесуточного распределения поправки атмосферного химического потенциала.
Строят карты распределения сумм локальных максимумов поправки атмосферного химического потенциала.
Строят область, охватывающую локальные максимумы на картах распределения сумм локальных максимумов поправки атмосферного химического потенциала.
Полученную область постулируют как наиболее вероятное место возникновения землетрясения.
Пример 1. Прогноз сейсмического события, произошедшего в Чили 16 сентября 2015 г. со следующими параметрами:
- дата, время: 16.09.2015 г., 22:54;
- координаты: 31.55° ю.ш.; 71.58° з.д.;
- магнитуда: 8.3;
- глубина гипоцентра, км: 22.4.
Параметры землетрясения взяты из базы данных Европейского Средиземноморского сейсмологического центра, доступные в Интернет по адресу: http://http://www.emsc-csem.org/.
1.1 Определение пространственного положения сейсмомагнитных меридианов (СММ) производится по формуле (2):
Символическая запись СММ (m1), (m2), (m3) состоит из ряда блоков, разделенных , при этом:
1-й блок - дата геоэффективного солнечного явления, вызывавшего геомагнитные возмущения;
2-й - долгота подсолнечной точки на дату геоэффективного явления, выбрана из Астрономического ежегодника;
3-й - 3-х часовые индексы геомагнитной активности в виде трех подблоков, разделенных точкой с запятой: код вида станции (H - по данным высокоширотных станций, M - среднеширотных, B - станции Боулдере, P - планетарные); номер 3-х часового интервала, в котором наблюдалось геоэффективное геомагнитное возмущение; значения 3-х часовых соседних индексов в шкале от 0 до 8, разделенные запятой;
4-й - долгота СММ на экваторе, рассчитываемая по формуле (2).
Данные Kp-индексов геомагнитной активности (3-х часовых индексов) размещены в открытом доступе в Интернет на сайтах Немецкого исследовательского центра в области наук о Земле (GeoForschungsZentrum Potsdam, GFZ Potsdam) и Центра космической погоды Национального управления океанографических и атмосферных исследований (NOAA Space Weather Prediction Center, SWPC).
1.2 Определение вероятных мест возникновения землетрясений как окрестностей радиусом 700 км от пересечений сейсмомагнитных меридианов с границами литосферных плит производится в следующей последовательности.
Строятся линии СММ на карте в проекции Меркатора и определяются вероятные места возникновения землетрясений как окрестности (круговые зоны) радиусом 700 км с центрами в точках пересечения СММ с границами плит, показанными на Фиг. 1 толстыми линиями. Вероятные места возникновения землетрясений - окрестности с радиусом R=700 км - определены от пересечения СММ (m1) и (m2) с южной границей плиты Альтиплано - окрестности A и B и СММ (m3) с западной границей Южно-Американской плиты - окрестность C.
Центр окрестности A: (21.5° ю.ш., 67.6° з.д.).
Центр окрестности B: (37.6° ю.ш., 74.5° з.д.).
Центр окрестности C: (30.7° ю.ш., 72.6° з.д.).
1.3 Выявление на космических снимках окрестности над пересечением сейсмомагнитных меридианов с границами литосферных плит, в которой присутствует облачная сейсмоиндуцированная структура.
Глобальный мониторинг облачного покрова с использованием метеоспутников позволил обнаружить 28.08.2015 г. на снимках облачности со спутника GOES-13 в районе центрального Чили L-образный ОСТИ (Фиг. 2). Данный участок на Фиг. 2 для удобства анализа обведен овальной линией, а окрестности над пересечением сейсмомагнитных меридианов с границами литосферных плит обозначены A, B и C.
Из Фиг. 2 видно, что облачная сейсмоиндуцированная структура присутствуют в окрестностях B и C, что определяет их как вероятное место возникновения землетрясения.
1.4 Определение максимального линейного размера облачной сейсмоиндуцированной структуры. Максимальный линейный размер облачной сейсмоиндуцированной структуры, определенный по космическому снимку, представленному на Фиг. 2, составляет L=2700 км.
1.5 Определение магнитуды прогнозируемого землетрясения производится по формуле (3), .
1.6 Определение прогнозной даты землетрясения осуществляется в соответствии с формулой (1):
27.08.2015+[(14∨21)±2]=10.09.2015∨17.09.2015±2 сут;
1.7 Производится определение временного интервала, составляющего 18 суток, отсчитываемых от 5 суток, предшествующих прогнозной дате землетрясения:
19.09.2015-5-18=27.08.2015
Определенный временной интервал составляет: 27.08.2015-13.09.2015
8
1.8 Расчет поправки атмосферного химического потенциала выполняется для прямоугольной области, накрывающей окрестности B и C, имеющей координаты: -45° ю.ш… -23° ю.ш' -82° з.д… -65° з.д. осуществляется по формуле (патент RU 2561305 C2 МПК G01T 1/169, G01W 1/00, Опубликовано: 27.08.2015):
g
где ΔU - поправка атмосферного химического потенциала, эВ; T - температура у поверхности Земли, °C; H - влажность у поверхности Земли, %.
Данные о приземной температуре и влажности взяты из базы данных GFS-ANL по состоянию на 0,6,12 и 18 часов каждых суток на регулярной сетке 0,5×0,5 градуса - продукт GFS Analysis, доступный в сети Интернет на сайте Национального Центра Информации об окружающей среде Национального управления океанических и атмосферных исследований США: https://www.ncdc.noaa.gov/data-access/model-data/model-datasets/global-forcast-system-gfs.
Построенные карты среднесуточного распределения поправки атмосферного химического потенциала приведены на Фиг. 3-20. Границы окрестностей B и C на Фиг. 3-21 показаны белыми окружностями.
1.9 Карты распределения поправки атмосферного химического потенциала покрываются регулярной сеткой, имеющей ячейки со сторонами 2 град, на 2 град. Сетка показана пунктирными линиями на Фиг. 3-20.
Сетка, покрывающая окрестности B и C имеет размерность 17×20 ячеек.
1.10 Производится определение положения локальных максимумов на каждой построенной карте поправки атмосферного химического потенциала. Положения максимумов показаны кружочками на Фиг. 3-20.
1.11 Производится подсчет числа локальных максимумов в каждой ячейке сетки на каждой карте распределения поправки атмосферного химического потенциала и строятся карты распределения сумм локальных максимумов, показанные на Фиг. 21-38.
Анализ Фиг. 21-38 показывает, что уже с 01.09.2015 положение основных максимумов на картах Фиг. 21-38 стабилизируется.
1.12 Строится область, охватывающую локальные максимумы на картах распределения сумм локальных максимумов, представленных на Фиг. 21-38. Построенная область показана на Фиг. 39.
1.13 Полученную область постулируется как наиболее вероятное место возникновения ожидаемого землетрясения.
1.14 Фактическое событие произошло 16.09.2015 г в окрестности C, показанной на Фиг. 1, 2 с магнитудой 8,3. На Фиг. 40 показано положение эпицентра землетрясения и прогнозные окрестности A, B и C.
1.15 Положение эпицентра фактического события и выявленного наиболее вероятного места возникновения землетрясения показаны на Фиг. 41. Событие произошло в выявленной зоне.
Вышеизложенное подтверждает повышение точности при определении пространственного положения прогнозируемого землетрясения при использовании предложенного технического решения.
Таким образом, предложенное техническое решение при выбранных расчетных соотношениях, приемах и операциях обеспечивает повышение точности при определении пространственного положения прогнозируемого землетрясения. Это уменьшается неопределенность при планировании и проведении комплекса мероприятий и работ, направленных на минимизацию ущерба для населения и объектов хозяйственной инфраструктуры от наступления прогнозируемого катастрофического события - землетрясения.