СПОСОБ ПРОГНОЗИРОВАНИЯ ЗЕМЛЕТРЯСЕНИЙ Российский патент 2008 года по МПК G01V1/00 

Изобретение относится к области геофизики, а именно к сейсмологии, и может быть использовано для прогнозирования места, силы и времени тектонических землетрясений.

Известен способ прогнозирования землетрясений по пространственно-временному распределению слабых землетрясений, при котором выявляют область сейсмического затишья, оконтуренную кольцевой сейсмичностью, а по размерам этой области и продолжительности затишья судят о возможности землетрясения [1].

Недостатками этого способа являются отсутствие возможности локализации и определения энергетического потенциала очаговой зоны (ОЗ) - объема, в котором заключена основная доля упругой энергии ожидаемого землетрясения - и связанные с этим трудности организации наблюдений в эпицентральной области ожидаемого землетрясения для осуществления его прогноза в реальном времени.

Наиболее близким к предлагаемому является способ прогнозирования землетрясений, при котором в пределах прогностического полигона создают сеть сейсмических станций с возможностью регистрации слабых землетрясений и локализации их очагов, проводят на полигоне режимные наблюдения и по результатам обработки на ЭВМ наблюдательных данных осуществляют прогноз землетрясений [2].

Недостатками этого способа является недостаточная информативность прогноза, заключающаяся в отсутствии количественных данных о координатах потенциальной очаговой зоны и величине накопленной в ней энергии.

Предлагаемый способ прогнозирования землетрясений основывается на приведенных ниже эмпирических и теоретических положениях.

Распределение упругой составляющей сдвиговой деформации ε_τ(r) в окрестности тангенциально нагруженной шарообразной ОЗ в первом приближении может быть задано выражением

где ε_τ(0) - упругая составляющая сдвиговых деформаций в центре ОЗ, ε_τ(r) - то же на расстоянии r от центра ОЗ, - объем ОЗ, r₀ - ее радиус. Учитывая приближенный характер проводимых оценок, можно принять для всего объема ОЗ , а в области r>r₀ использовать формулу (1).

Для точек, лежащих на поверхности Земли, , где Δ - эпицентральное расстояние, h - глубина очага, формула (1) принимает вид

Для эпицентра имеем

откуда

Подставляя это выражение в (2), получим

Из этой формулы следует, что касательные напряжения в приповерхностном однородном изотропном слое земной коры, вызванных подготовкой землетрясения на глубине h, имеют круговое распределение вокруг его эпицентра. Следовательно, по данным площадных наблюдений ε_τ можно определить эпицентр готовящегося землетрясения, а по измеренным значениям ε_τ(Δ,h), ε_τ(0,h) и известному Δ - глубину очага h, т.е. по этим наблюдениям можно определить координаты потенциальной ОЗ.

Для других предвестников, пропорциональных ε_τ (например, уровня микросейсмической эмиссии), очевидно, имеет место аналогичная формула

где А(0,h) и А(Δ,h) - амплитуды аномалий предвестника соответственно в эпицентре и на расстоянии Δ от эпицентра. В случае микросейсмической эмиссии в роли A(0,h) и A(Δ,h) формуле (5) может быть использован уровень микросейсмической эмиссии (УМЭ) - число микроземлетрясений в единицу времени в единице объема горных пород в заданном энергетическом диапазоне. Как и в предыдущем случае, этих данных достаточно для определения по формуле (5) координат 03 готовящегося землетрясения. Для уменьшения помех, обусловленных процессами, не связанными с подготовкой землетрясения, для наблюдений микросейсмической эмиссии следует выбрать горизонтальный консолидированный слой прочных и жестких (компетентных) горных пород мощностью в несколько километров, расположенный ниже приповерхностного рыхлого и неоднородного слоя, но заведомо выше ОЗ (на глубине 3-10 км, в зависимости от геологических условий). В плане он должен охватывать зону подготовки землетрясения, представляющую собой круг с радиусом Δ_max, равным расстоянию от эпицентра, на котором обнаруженный уровень деформаций составляет 10^-8, т.е. соизмерим с уровнем земноприливных деформаций (радиус зоны подготовки землетрясения [2]). Еще одним условием, предъявляемым к этому слою, которого можно назвать компетентным слоем, является существенное превышение скорости накопления напряжений в нем за счет процессов в ОЗ, над скоростью их релаксации за счет микросейсмической эмиссии. В противном случае регистрируемые параметры не будут адекватно отражать процесс накопления напряжений в ОЗ.

Для оценки энергетических характеристик потенциальной ОЗ используем соотношение (3), переписав его в виде

ε_τ(0)·V≈ε_τ(0,h)·h³.

Умножив обе части этого соотношения на τ_cr/2, где τ_cr - критические касательные напряжения для пород в ОЗ, получим

С учетом ε_τ(0)≈ε_cr (см. выше) и (1/2)·τ_cr·ε_cr=w_cr, где w_cr - плотность упругой энергии в ОЗ, находим, что в левой части равенства (6) имеем полную упругую энергию E_cr, заключенную в ОЗ, так что

Излучаемая в виде сейсмических колебаний доля этой энергии Е_s=η·Е_cr, где η - сейсмический КПД очага, составляющая лишь несколько процентов, так что

(в системе СИ). Прологарифмируя это выражение, получаем оценку энергетического класса ожидаемого землетрясения

с помощью которой, используя также известную формулу Гуттенберга-Рихтера

можно оценить его магнитуду

Учитывая формулу (4) и используя соотношение τ_cr=G·ε_cr, где G и ε_cr - соответственно модуль сдвига и критические значения относительной деформации для горных пород, получим

Далее после подстановки вместо Δ радиуса зоны подготовки землетрясения Δ_max и характерных для горных пород значений G≈30 ГПа=3·10¹⁰ дж/м³ и ε_cr≈10^-4, имеем

Переходя на более привычные единицы измерения эпицентральных расстояний Δ и глубин очагов h - километры, получим после соответствующих преобразований формулу

которая в области Δ_max>>h (т.е. для сильных землетрясений) удовлетворительно согласуется с известной эмпирической формулой [2]

при η≅0.03. Подставляя это значение η в (12), получим окончательную оценку

Очевидно, в пределах достижимой точности оценки магнитуды готовящегося землетрясения можно предположить, что радиус зоны подготовки землетрясения Δ_max, определяемый по деформационным наблюдениям, совпадет с таковым, определенным по наблюдениям микросейсмической эмиссии. Другими словами, в формулу (14) вместо Δ_max можно подставить расстояние от эпицентра, на котором обнаруженная аномалия УМЭ соизмерима с его фоновым значением.

Таким образом, для прогнозирования места, силы и времени ожидаемого землетрясения по предлагаемому способу следует осуществлять следующую последовательность действий: 1) выявление областей накопления упругих тектонических деформаций в земной коре (внутренние участки взаимодействующих геоблоков или область сейсмического затишья, оконтуренная кольцевой сейсмичностью); 2) выделение с использованием геологических данных компетентного слоя, перекрывающего в плане зону подготовки землетрясения; 3) создание наблюдательной сети и режимные наблюдения УМЭ из выделенного слоя; 4) выявление в изолиниях УМЭ математическими методами статистически значимых кольцевых составляющих, оконтуривающих эпицентры потенциальных землетрясений, определение эпицентра главного толчка и глубины залегания его ОЗ по формуле (5) (прогноз места землетрясения); 5) определение радиуса зоны подготовки главного толчка как расстояния от эпицентра, на котором обнаруженная аномалия УМЭ соизмерима с его фоновым значением и определение по формуле (14) его магнитуды (прогноз силы землетрясения); 6) определение времени главного толчка по комплексным наблюдениям за ОЗ с использованием известных предвестников (прогноз времени землетрясения).

Ключевым и в то же время самым трудным звеном в предложенном способе прогнозирования землетрясений остается наблюдения предвестников не вблизи разломных зон (как это обычно принято), где амплитуды их велики, но трудно интерпретируемы, а во внутренних участках геоблоков, где амплитуды аномалий существенно меньше, но их пространственное распределение носит более регулярный, детерминированный характер, что облегчает решение обратной задачи - задачи интерпретации наблюдений и прогноза землетрясений.

Использование предлагаемого способа позволит количественно оценить магнитуды будущих толчков, а также определить эпицентры и глубины их очагов уже на докритической стадии подготовки, соответствующей этапу долго- и среднесрочного прогноза, и тем самым существенно повышает надежность краткосрочного прогноза времени землетрясения, поскольку позволяет заблаговременно организовать комплексные наблюдения за очаговой зоной.

ЛИТЕРАТУРА

1. Моги К. Предсказание землетрясений. М.: Мир, 1988, 382 с.

2. Прогноз землетрясений, №3. Душамбе - Москва: ДОНИШ, 1984а, 216 с.

Изобретение относится к области сейсмологии и может быть использовано для прогнозирования места, силы и времени землетрясений. Сущность: выявляют область накопления упругих тектонических деформаций в земной коре. Выбирают в этой области ниже рыхлого неоднородного осадочного слоя и заведомо выше характерных очаговых зон однородный горизонтальный слой прочных и жестких пород мощностью в несколько километров - компетентный слой. Проводят режимные наблюдения уровня микросейсмической эмиссии из компетентного слоя. Строят изолинии уровней микросейсмической эмиссии и выявляют в них известными методами статистически значимые кольцевые составляющие, определяют по ним эпицентр ожидаемого главного толчка, а затем глубину залегания его очаговой зоны. Определяют амплитуду. По комплексным наблюдениям очаговой зоны судят о времени землетрясения. Технический результат: повышение надежности краткосрочного прогноза.

Способ прогнозирования землетрясений, включающий создание в пределах прогностического полигона сети сейсмических станций с возможностью регистрации слабых землетрясений и локализации их очагов, режимные микросейсмические наблюдения на станциях, обработку на ЭВМ наблюдательных данных и определение по ним параметров ожидаемого землетрясения, отличающийся тем, что выявляют область накопления упругих тектонических деформаций в земной коре (внутренние участки взаимодействующих геоблоков или область сейсмического затишья, оконтуренная кольцевой сейсмичностью), выбирают в этой области ниже рыхлого неоднородного осадочного слоя и заведомо выше характерных глубин очаговых зон однородный горизонтальный слой прочных и жестких (компетентных) пород мощностью в несколько километров (компетентный слой), проводят режимные наблюдения уровня микросейсмической эмиссии (УМЭ) из компетентного слоя (числа микроземлетрясений в единицу времени в единице объема в заданном энергетическом диапазоне), строят изолинии УМЭ и выявляют в них известными методами статистически значимые кольцевые составляющие, определяют по ним эпицентр ожидаемого главного толчка, а затем - глубину h (км) залегания его очаговой зоны по формуле

где A(0,h) и A(Δ,h) - УМЭ соответственно в эпицентре и на расстоянии Δ от эпицентра (прогноз места землетрясения),

определяют радиус зоны подготовки главного толчка Δ_max (км) как расстояния от эпицентра, на котором обнаруженная аномалия УМЭ соизмерима с его фоновым значением, и по формуле определяют его магнитуду (прогноз силы землетрясения), а по комплексным наблюдениям очаговой зоны судят о времени землетрясения.