СПОСОБ ОПРЕДЕЛЕНИЯ МЕСТА ЗЕМЛЕТРЯСЕНИЯ С БОРТА КОСМИЧЕСКОГО АППАРАТА Российский патент 2005 года по МПК G01V9/00 

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

Проблема предсказания землетрясений имеет многовековую историю.

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

Первые наблюдения предвестников землетрясения на космических аппаратах были проведены около 10 лет назад. На спутниках «Интеркосмос-19», «ОГО-6», «Ореол-3» были обнаружены явления резкого повышения интенсивности низкочастотных электромагнитных излучений (0,1-16 КГц) за несколько часов до землетрясения при пролете спутника в некоторой области, лежащей вблизи эпицентра землетрясения. Аналогичные явления наблюдались и на ИСЗ «Геос-1» и «Геос-2», на которых находилась аппаратура, регистрирующая ОНЧ-излучение (60-7500 Гц). Проводились синхронные наблюдения на спутниках и на наземной станции на острове Кергелен. Всплески ЭМИ отмечались за несколько десятков минут для нескольких исследованных землетрясений, эпицентры которых были расположены вблизи силовой линии геомагнитного поля, проходящего через остров Кергелен.

На ИСЗ «Интеркосмос-Болгария-1300» за десятки минут до землетрясения были зарегистрированы вариации электрического поля (0,1-8 Гц) в двух зонах: над очагом и в магнитосопряженной точке [1].

На ряде ИСЗ проводились измерения ионного и нейтрального состава ионосферы, плотности и температуры нейтронов и положительных ионов, энергетического состава фотоэлектронов, потоков заряженных частиц с энергией до 25 КэВ. В ряде случаев были обнаружены изменения в концентрации плазмы над очагами готовящихся землетрясений. Формировался локальный провал в концентрации плазмы, доходящий до 20%.

При обработке экспериментальных данных от прибора «Мария» со спутника «Интеркосмос- Болгария-1300» было обнаружено 11 резких флуктуации в темпах счета, 5 из которых произошли в отсутствии геомагнитных возмущений и сопровождались через 2,5-3 часа землетрясениями. Спектрометр «Мария» зарегистрировал 11 более чем десятикратных возрастании интенсивности частиц с энергией более нескольких десятков МэВ, около 80% которых составляли протоны. Аппаратура «Мария» была установлена также на ОК «МИР».

Из всего многообразия способов контроля готовящегося землетрясения наиболее отработанным можно назвать способ, заключающийся в регистрации ОНЧ-излучения аппаратурой КА [2]. Данный способ основывается на наличии так называемого «тримпи-эффекта», суть которого заключается в том, что за несколько часов до землетрясения в результате появляющихся в коре напряжений генерируется мощное ОНЧ-излучение с частотами порядка 20 КГц. Возникающее излучение регистрируется за несколько часов до землетрясения при пролетах спутника в некоторой области, вытянутой вдоль параллели и лежащей над эпицентром землетрясения. Недостатками данного способа являются трудности локализации предполагаемого землетрясения и сильное воздействие на приемник ОНЧ-излучения различных помех, в том числе и от самого КА.

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

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

Из известных способов наиболее близким к предлагаемому является «способ определения места и времени землетрясения с борта космического аппарата» [3], который и выбран в качестве прототипа.

Сущность данного способа состоит в том, что регистрируют с борта космического аппарата низкочастотное электромагнитное излучение, в момент превышения его интенсивности фонового уровня дополнительно сканируют участки подстилающей земной поверхности в рентгеновском диапазоне спектра с энергией 2-25 КэВ и по наличию интенсивности излучения, превышающей фоновое значение не менее чем на 20 стандартных отклонений, уточняют местоположение землетрясения.

Недостатки данного способа связаны с трудностями локализации места предполагаемого землетрясения и с сильным воздействием на приемник ОНЧ-излучения различных помех.

Технической задачей изобретения является точная и однозначная пеленгация эпицентра ожидаемого землетрясения и повышение помехоустойчивости приемного ОНЧ-комплекса.

Поставленная задача решается тем, что согласно способу определения места землетрясения с борта космического аппарата, включающего регистрацию низкочастотного электромагнитного излучения и суждение по превышению интенсивностью излучения фонового уровня о местоположении эпицентра землетрясения, при этом в момент превышения интенсивностью низкочастотного излучения фонового уровня дополнительно сканируют с борта космического аппарата участки подстилающей земной поверхности в рентгеновском диапазоне спектра и по наличию и размерам участка земной поверхности, характеризующегося интенсивностью рентгеновского излучения с энергией 2-25 КэВ, превышающей фоновое значение не менее чем на 20 стандартных отклонений, уточняют местоположение эпицентра землетрясения, низкочастотное электромагнитное излучение принимают на три антенны, размещенные на космическом аппарате в виде геометрического прямого угла, в вершине которого помещают антенну первого приемного канала, общую для антенн второго и третьего приемных каналов, размещенных в азимутальной и угломестной плоскостях соответственно, преобразуют принимаемое низкочастотное излучение по частоте с использованием двух гетеродинов, частоты которых разносят на удвоенное значение промежуточной частоты

f_Г2-f_Г1=2f_up

и выбирают симметричными относительно несущей частоты f_c низкочастотного излучения

f_с-f_Г1=f_Г2-f_c=f_up,

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

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

где d₁, d₂ - измерительные базы;

λ - длина волны;

α - азимут эпицентра ожидаемого землетрясения;

β - угол места эпицентра ожидаемого землетрясения.

Однако фазовому методу пеленгации эпицентра ожидаемого землетрясения свойственно противоречие между требованиями повышения точности и устранения неоднозначности фазовых измерений. Действительно, согласно вышеприведенным формулам для повышения точности пеленгации эпицентра ожидаемого землетрясения необходимо увеличить относительные размеры баз d₁/λ и d₂/λ. Однако с увеличением d₁/λ и d₂/λ наступает момент, когда Δϕ₁ и Δϕ₂ превосходят значение 2π, т.е. наступает неоднозначность отсчета угловых координат α и β. Для устранения указанной неоднозначности используют корреляционную обработку принимаемых по трем каналам низкочастотных электромагнитных излучений. При этом для преобразования по частоте принимаемых низкочастотных электромагнитных излучений используют два гетеродина, частоты f_Г1 и f_Г2 которых разносят на удвоенное значение промежуточной частоты

f_Г2-f_Г1=2f_up

и выбирают симметричными относительно несущей частоты f_c низкочастотного излучения (фиг.5)

f_с-f_Г1=f_Г2-f_c=f_up,

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

В предлагаемом способе, так же как и в прототипе, дополнительно сканируют с борта КА в рентгеновском диапазоне спектра участки подстилающей земной поверхности, регистрируют в верхней атмосфере образование фотонов при взаимодействии с атмосферой потоков электронов, высыпающихся из внутреннего радиационного пояса под воздействием ОНЧ-излучения, и по наличию и размерам участка подстилающей земной поверхности, характеризующегося интенсивностью рентгеновского излучения с Е=2-25 КэВ, превышающей фоновое значение не менее чем на 20 стандартных отклонений, уточняют местоположение эпицентра ожидаемого землетрясения.

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

Достигается это следующим образом.

Распространяясь от Земли и достигая L - оболочек, ОНЧ-излучение эффективно рассеивает электроны с энергиями порядка 30 КэВ - 1 МэВ, которые, попадая сверху из радиационного пояса, поглощаются атмосферой на высотах 70-100 км. Энергетический эквивалент потока при этом достигает 0,1 эрг/см²·с. При этом основным процессом взаимодействия с воздухом является ионизация атомов. Но, как следует из рассмотрения сечений взаимодействия, существует вероятность порядка 10^-4 образования фотона тормозным излучением. При данном механизме примерно в половине случаев образуется фотон с энергией, сравнимой с энергией электрона, т.е. в районе 10-30 КэВ, и примерно половина таких фотонов, двигаясь вверх, не поглотится атмосферой.

Энергия электрона в 30 КэВ соответствует 5-10^-8 эрг и поток частиц составляет 2-10⁶ см^-2·с^-1. Если район генерации ОНЧ-излучения ограничен на Земле областью размерами в 100 км, то светимость в рентгеновском диапазоне слоя, где происходит поглощение электронов, составляет

Е=0,1·(5·10^-8)^-1·10⁴ фотонов/с.

На высоте орбиты в 300 км поток фотонов будет составлять

Е=0,25·2·10¹⁶/(4πR)²≅2 фотона/см²,

где R=200 км.

Предположим, что на спутнике установлен пропорциональный счетчик площадью 100 см, регистрирующий рентгеновские фотоны с эффективностью 0,5. Его собственный фон, остающийся после селекции заряженных частиц, равен примерно 2 ими/см²·с. Прибор пролетает область высыпания электронов примерно за 10 с. В таком случае от фона будет накоплено 2000 ими, а от излучающей области 1000 ими, что составляет 20 стандартных отклонений и представляет вполне регистрируемую величину.

Следовательно, по изображению подстилающей поверхности в рентгеновском диапазоне по параметрам яркого пятна (размеры, спектр, поведение во времени) распознается зона готовящегося землетрясения.

Кроме того, в этом случае повышается надежность контроля готовящегося землетрясения с борта КА.

Предлагаемый способ иллюстрируется фиг.1-7. На фиг.1 обозначены Земля 1, область 2 (≈100 км) образования напряжений в земной коре;

образование 3 ОНЧ-излучения, траектория 4 полета КА, поле зрения 5 рентгеновского датчика; КА 6, внутренний радиационный пояс 7; траектория 8 движения электронов по магнитным силовым линиям; высыпание электронов 9 под воздействием ОНЧ-излучения; образование фотонов 10 при взаимодействии потоков электронов с атмосферой на высоте 70-100 км; атмосфера 11.

Возможность идентификации зоны готовящегося землетрясения иллюстрируется фиг.2 и 3. На фиг.2 обозначены: зона 12 готовящегося землетрясения, трасса 13 КА, проходящая над очагом готовящегося землетрясения, трассы 14 и 15 КА, проходящие вблизи очага готовящегося землетрясения.

На фиг.3 показаны в зависимости от времени интенсивности потока, регистрируемого рентгеновским детектором с полем зрения ≈10° в диапазоне 2-25 КэВ при полете КА над очагом готовящегося землетрясения (зависимость 16) и вблизи очага готовящегося землетрясения (зависимости 17 и 18). Время соответствует полету КА на минимальном удалении от эпицентра готовящегося землетрясения.

Предлагаемый способ может быть реализован с помощью известных средств.

Для регистрации ОНЧ-излучения может использоваться приемный ОНЧ-комплекс, включающий в свой состав:

- трехкомпонентные ОНЧ-детекторы магнитного поля (типа МДП-3);

- шестикомпонентный приемник-анализатор (типа АНЧ-6);

- приемник ОНЧ-сигналов, измеряющий амплитуду и фазу СДВ радиопередатчиков (типа ПОНЧ-6);

- микропроцессор-анализатор (типа МА-2);

- приемник-пеленгатор ОНЧ-сигналов.

В качестве КА могут использоваться существующие орбитальные комплексы, перспективные комплексы типа УКП, выводимые на геостационарные орбиты, и другие аппараты.

Для сканирования с борта КА участков земной поверхности могут использоваться системы управления типа «Дельта», «Чайка», устанавливаемые на ОС «Салют» и других КА.

Для регистрации фотонов в рентгеновском диапазоне может использоваться аппаратура типа СКР, установленная на КА «Астрон», либо «Барс», а также радиотелеметрическая система типа БР 9 ЦУ.

На фиг.4 изображена структурная схема приемника-пеленгатора. Частотная диаграмма, поясняющая образование дополнительных каналов приема, представлена на фиг.5. Взаимное расположение приемных антенн показано на фиг.6 геометрическая схема расположения КА и эпицентра 12 ожидаемого землетрясения показана на фиг.7.

Приемник-пеленгатор содержит приемные антенны 19, 20 и 21, гетеродины 22 и 23, смесители 24, 25 и 26, усилители 29, 30 и 31 промежуточной частоты, перемножители 27, 32 и 33, узкополосные фильтры 28, 34 и 35, блоки корреляторов 36 и 37, пороговые блоки 38, 39, 42 и 43, ключи 40, 41, 44 и 45, фазовые детекторы 46 и 47, блок 48 регистрации.

Принимаемые низкочастотные электромагнитные излучения:

i₁(t)=У₁·cos[2π(f_c±Δf)t+ϕ₁]

i₂(t)=У₂·cos[2π(f_c±Δf)t+ϕ₂]

i₃(t)=У₃·cos[2π(f_c±Δf)t+ϕ₂], 0≤t≤T_c

где У₁, У₂, У₃ - интенсивности излучения;

f_c, ϕ₁, ϕ₂, ϕ₃, T_c - несущая частота, начальные фазы и длительность излучения;

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

с выходов приемных антенн 19, 20 и 21 поступают на первые входы смесителей 24, 25 и 26 соответственно, на вторые входы которых подаются напряжения гетеродинов 22 и 23:

u_Г1(t)=υ_Г1·cos(2πf_Г1t+ϕ_Г1),

u_Г2(t)=υ_Г2·cos(2πf_Г2t+ϕ_Г2).

Причем частоты f_Г1, и f_Г2, гетеродинов 22 и 23 разнесены на удвоенное

f_Г2-f_Г1=2f_up

и выбраны симметричными относительно несущей частоты f_c. (фиг.5)

f_с-f_Г1=f_Г2-f_c=f_up,

На выходах смесителей 24, 25 и 26 образуются напряжения комбинационных частот. Усилителями 29, 30 и 31 выделяются напряжения промежуточной (разностной) частоты:

где

K₁ - коэффициент передачи смесителей;

f_up=f_c-f_Г1=f_rГ2-f_c- промежуточная частота;

ϕ_up1-ϕ₁-ϕ_/1; ϕ_up2=ϕ₂-ϕ; ϕ_up3=ϕ_Г2-ϕ₃,

которые поступают на входы перемножителей 32 и 33. На выходах последних образуются напряжения:

где

K₂- коэффициент передачи перемножителей;

ϕ_Г=ϕ_Г2-ϕ_Г1;

которые выделяются узкополосными фильтрами 34 и 35 соответственно.

Напряжения u_Г1,(t) и u_Г2(t) со вторых выходов гетеродинов 22 и 23 подаются на входы перемножителя 27, на выходе которого образуется напряжение

где

ϕ_Г=ϕ_Г2-ϕ_Г1,

которое выделяется узкополосным фильтром 28.

Напряжения u_пр1(t),u_пр2(t) и u_пр3(t) одновременно поступают на два входа блоков корреляторов 36 и 37, на первом выходе которых формируются напряжения, пропорциональные корреляционным функциям R₁(τ) и R₂(τ). Указанные напряжения сравниваются с пороговым напряжение υ_пор1 в пороговых блоках 38 и 39. Пороговый уровень υ_пор1 превышается только при истинных азимуте α₀ и угле места β₀.. При превышении порогового уровня υ_пор1 в пороговых блоках 38 и 39 формируются постоянные напряжения, которые поступают на управляющие входы ключей 40 и 41, открывая их. В исходном состоянии ключи 40, 41, 44 и 45 всегда закрыты. На вторых выходах корреляторов 36 и 37 формируются напряжения, пропорциональные корреляционным функциям R₃(τ) и R₄(τ). Указанные напряжения сравниваются с пороговым напряжением υ_пор2 в пороговых блоках 42 и 43.

Так как канальные напряжения u_пр1(t) и u_пр2(t), u_пр1(t) и u_пр3(t) образуются одним и тем низкочастотным электромагнитным излучением, принимаемым на частоте f_c то между указанными канальными напряжениями существует сильная корреляционная связь. Выходные напряжения, пропорциональные корреляционным функциям R₂(τ) и R₃(τ), превышают пороговый уровень υ_пор2 в пороговых блоках 42 и 43. В последних формируются постоянные напряжения, которые поступают на управляющие вход ключей 44 и 45, открывая их. При этом напряжения u₁(t) и u₂(t) с выходов узкополосных фильтров 34 и 35 через открытые ключи 40, 41, 44 и 45 поступают на первые входы фазовых детекторов 42 и 43 соответственно, на вторые входы которых подается напряжение u_Г(t) с выхода узкополосного фильтра 28. На выходах фазовых детекторов 42 и 43 образуются низкочастотные напряжения:

u_H1(α)=υ_H1·cosΔϕ₁,

u_H2(α)=υ_H2·cosΔϕ₂,

где

K₃ - коэффициент передачи фазовых детекторов;

пропорциональные фазовым сдвигам Δϕ₁ и Δϕ₂. Эти напряжения фиксируются блоком 48 регистрации.

Описанная выше работа приемника-пеленгатора соответствует случаю приема низкочастотных электромагнитных излучений по основному каналу на частоте f_с, (фиг.5).

Если ложный сигнал (помеха) принимается по первому зеркальному каналу на частоте f₃₁ или по второму зеркальному каналу на частоте f₃₂, то на вторых входах блоков корреляторов 36 и 37 напряжения отсутствуют. Ключи 44 и 45 не открываются, и указанные ложные сигналы (помехи) подавляются.

По аналогичной причине подавляются и ложные сигналы (помехи), принимаемые по первому комбинационному каналу на частоте f_K1 или по второму комбинационному каналу на частоте f_K2 или по любому другому дополнительному каналу.

Если ложные сигналы (помехи) одновременно принимаются по первому зеркальному каналу на частоте f₃₁ и по второму зеркальному каналу на частоте f₃₂, то на вторых выходах блоков корреляторов 36 и 37 образуются напряжения, пропорциональные корреляционным функциям R₃(τ) и R₄(τ). Однако ключи 44 и 45 в этом случае не открываются. Это объясняется тем, что канальные напряжения образуются разными ложными сигналами (помехами), принимаемыми на разных частотах f₃₁, и f₃₂. Между ними существует слабая корреляционная связь; напряжения, пропорциональные корреляционным функциям R₃(τ) и R₄(τ), не превышают порогового напряжения υ_пор2 в пороговых блоках 42 и 43. Ключи 44 и 45 не открываются, и ложные сигналы (помехи), принимаемые одновременно по первому зеркальному каналу на частоте f₃₁ и по второму зеркальному каналу на частоте f₃₂, подавляются.

По аналогичной причине подавляются и ложные сигналы (помехи), принимаемые одновременно по двум или более другим дополнительным каналам.

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

При этом прием низкочастотных электромагнитных излучений осуществляют на три антенны, размещенные на борту космического аппарата в виде геометрического прямого угла, в вершине которого помещают антенну первого канала приема, общую для антенн второго и третьего каналов приема, размещенных в азимутальной и угломестной плоскостях соответственно. Точность пеленгации эпицентра ожидаемого землетрясения достигается увеличением относительного размера измерительных баз d₁/λ, d₂/λ. A возникающую при этом неоднозначность отсчета угловых координат α и β устраняют корреляционной обработкой принимаемых излучений. Зная высоту полета КА и измерив углы α и β, точно и однозначно определяют местоположение эпицентра ожидаемого землетрясения.

Причем пеленгация эпицентра ожидаемого землетрясения осуществляется в двух плоскостях на стабильной частоте, равной разности частот гетеродинов f_Г2-f_Г1.

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

Повышение помехоустойчивости приемного ОНЧ-комплекса достигается подавлением ложных сигналов (помех), принимаемых по дополнительным каналам, также корреляционной обработкой.

Литература

1. Чмырев В.М. и др. Электрические поля и гидромагнитные волны в ионосфере над очагом землетрясения. - Геомагнетизм и аэрономия, 1986, т. 26, №6, с.1020-1022.

2. Гохберг М.Б. и др. Сейсмоэлектромагнитные явления. - М.: Наука, 1988, с.25-34,40-145.

3. Беляев М.Ю., Москаленко Е.И. Способ определения места и времени землетрясения с борта космического аппарата. Патент №2045086 (РФ), G 01 V 9/00, 1992.

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

Способ определения места землетрясения с борта космического аппарата, включающий регистрацию низкочастотного электромагнитного излучения и суждение по превышению интенсивностью излучения фонового уровня о местоположении эпицентра землетрясения, при этом в момент превышения интенсивностью низкочастотного излучения фонового уровня дополнительно сканируют с борта космического аппарата участки подстилающей земной поверхности в рентгеновском диапазоне спектра и по наличию и размерам участка земной поверхности, характеризующегося интенсивностью рентгеновского излучения с энергией 2-25 КэВ, превышающей фоновое значение не менее чем на 20 стандартных отклонений, уточняют местоположение эпицентра землетрясения, отличающийся тем, что низкочастотное электромагнитное излучение принимают на три антенны, размещенные на космическом аппарате в виде геометрического прямого угла, в вершине которого помещают антенну первого приемного канала, общую для антенн второго и третьего приемных каналов, размещенных в азимутальной и угломестной плоскостях соответственно, преобразуют принимаемое низкочастотное излучение по частоте с использованием двух гетеродинов, частоты которых разносят на удвоенное значение промежуточной частоты f_г2-f_г1=2f_ир и выбирают симметричными относительно несущей частоты f_c низкочастотного излучения f_с-f_г1=f_г2-f_c=f_пр,выделяют в трех приемных каналах напряжения промежуточной частоты, перемножают между собой напряжения гетеродинов, напряжения промежуточной частоты первого и второго, первого и третьего каналов, выделяют узкополосные напряжения на частоте, равной разности частот гетеродинов, осуществляют корреляционную обработку напряжений промежуточной частоты первого и второго, первого и третьего приемных каналов, сравнивают напряжения, пропорциональные полученным корреляционным функциям, с пороговыми напряжениями и в случае их превышения измеряют разности фаз между полученными узкополосными напряжениями на частоте, равной разности частот гетеродинов, по значению которых определяют направления на эпицентр ожидаемого землетрясения в азимутальной и угломестной плоскостях, причем частоту f_Г1 первого гетеродина используют для преобразования по частоте низкочастотного излучения, принимаемого по первому каналу, а частоту f_Г2 второго гетеродина используют для преобразования по частоте низкочастотных излучений, принимаемых по второму и третьему каналам.