Способ измерения ионосферных предвестников землетрясений

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

 

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

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

По параметрам переходного колебательного процесса (как вариации плотности электронной концентрации в слоях ионосферы) определяют гипоцентр проекции очага на ионосферу и характеристики ожидаемого сейсмического удара [см., например, «Краткосрочный прогноз катастрофических землетрясений с помощью радиофизических наземно-космических методов», Доклады конференции ОИФЗ им. О.Ю. Шмидта, РАН, М, 1998 г., стр. 64-65, 109, 127-129, 138].

Известен «Способ предсказания землетрясений» путем измерения волн плотности электронной концентрации при полете космического аппарата непосредственно на высотах соответствующих слоев ионосферы, Патент Ru №2 205 430, 2003 г - аналог. В способе аналога регистрируют волновой процесс, возникающий в среде накануне удара, вычисляют фазовый центр волнового процесса и его период, рассчитывают характеристики предстоящего удара по их зависимостям от периода, дополнительно получают серию регистрограмм дискретных отсчетов A(Li) электростатического потенциала внешней поверхности космических аппаратов при их полете по орбитам Li непосредственно через область ионосферного образования, находят фазовый центр процесса как точку пересечения траверз восходящих и нисходящих витков космических аппаратов, проведенных к участкам регистрограмм, где допплеровская частота измеряемого процесса равна нулю, отождествляют эту точку с проекцией гипоцентра очага на ионосферу, вычисляют период Т, магнитуду М и время ожидаемого удара tx из соотношений:

где ΔL=L2-L1 - разница пространственных периодов двух симметричных относительно траверзы полуволн регистрограммы;

- отношение пространственных периодов двух симметричных, относительно траверзы, полуволн регистрограммы,

v - скорость акустических волн в ионосфере;

Vr - радиальная скорость движения измерителя относительно фазового центра волнового процесса.

К недостаткам аналога следует отнести:

- все существующие космические аппараты имеют внутренние источники питания, «заземленные» на корпус. Создание дополнительного внешнего корпуса (в качестве обкладки конденсатора) изолированного от «заземленного» представляет технические трудности;

- невысокая чувствительность измерений из-за малой емкости создаваемого конденсатора.

Известны дистанционные методы измерений плотности электронной концентрации ионосферы N [1/м3], путем ее зондирования на частоте ниже критической [см., например, «Космонавтика», Энциклопедия, М, Изд. Сов. энциклопедия, 1986 г., стр. 161, Ионозоид] - ближайший аналог.

В ближайшем аналоге задающий генератор плавно изменяет частоту настройки приемно-передающего устройства в диапазоне от 1 до 20 МГц для получения амплитудно-частотной характеристики (АЧХ) отраженного от ионосферы сигнала. Ионозоид включает импульсный ВЧ передатчик, приемник, электронно-лучевой индикатор, задающий генератор передатчика является гетеродином приемника, чем достигается сопряжение и синхронизация приемника и передатчика. Высоту до отражающей поверхности ионосферы определяют по времени запаздывания отраженного сигнала. Плотность электронной концентрации слоя определяют пересчетом из АЧХ значений критической частоты и времени запаздывания.

Недостатками ближайшего аналога являются:

- невозможность по параметрам регистрируемого сигнала рассчитать ионосферные предвестники землетрясения;

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

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

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

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

рассчитывают время удара tу=4,7T и ожидаемую магнитуду как

где А0 - установившееся значение амплитуды сигнала предвестника, равное

Δt=(t2-t1)=(t3-t2) интервал времени между отсчетами измерений A1, А2, A3.

Изобретение поясняется чертежами, где:

фиг. 1 - плотность электронной концентрации в слоях ионосферы;

фиг. 2 - суточный ход критической частоты (плотности электронной концентрации) в слое F2 а) невозмущенном состоянии, б) возмущенном, накануне сейсмического удара;

фиг. 3 - пеленгация фазового центра волн плотности электронной концентрации двумя пунктами, разнесенными на измерительной базе;

фиг. 4 - функция изменения амплитуды регистрируемого сигнала во времени;

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

Техническая сущность изобретения заключается в следующем. Накануне сейсмического удара в атмосфере происходит раскачка очага землетрясения [см., например, Патент Ru №2170446 кл. G.01.V, 9/00, 2001 г.] В приповерхностном слое атмосферы возникают акустолитосферные волны, которые, при их распространении вверх, служат «спусковым крючком» для возникновения плазменных волн электронной концентрации в слоях ионосферы. Исходная электронная концентрация N[1/m3] в слоях ионосферы иллюстрируется фиг. 1 [см., например, «Космонавтика», Энциклопедия, под ред. В.П. Глушко, М, Сов. энциклопедия, 1985 г., стр. 143]. Зарегистрированные в ряде экспериментов волны плотности электронной концентрации в ионосфере иллюстрируются графиками фиг. 2 [см., например, «Краткосрочный прогноз катастрофических землетрясений с помощью радиофизических наземно-космических методов», Доклады конференции ОИФЗ им. О.Ю. Шмидта, РАН, М, 1998 г., стр. 109]. За несколько часов до сейсмического удара, изменение критической частоты слоя F2 из-за изменения плотности электронной концентрации может достигать 40…50%. Для измерения пространственных волн в ионосфере предлагается осуществлять ее зондирование с двух, разнесенных на измерительной базе, пунктов. Признаками предвестниками землетрясения являются: гипоцентр (координаты) очага, ожидаемое время сейсмического удара (ty) и ожидаемая магнитуда (М).

Пеленгация фазового центра ионосферных волн иллюстрируется фиг. 3. Фазовый центр ионосферных волн находят как точку пересечения радиус-векторов. Положение радиус-векторов в пространстве полностью определяется косинус-направляющими. В прямоугольной системе координат, задаваемой ортогональными парами измерителей (13, 14), (15, 16), (17, 18), (19, 20) косинус-направляющая вектора равна отношению его проекции (Пр) на данную ось к длине вектора. Длины векторов R1, R2 фиг. 3 находятся как корень квадратный из суммы проекций:

Проекции радиус-векторов пропорциональны сигналам на входе приемников в ортогональных плоскостях х, у.

Кроме направления в пространстве, ионосферная волна характеризуется периодом (Т) и амплитудой A(t), которая изменяется во времени, фиг. 2. Из математики известно [см., например, Пискунов Н.С., «Дифференциальное и интегральное исчисления для ВТУЗов», учебник т. 1, 5-е издание, М, Наука, 1964 г., стр. 457-458], что сама функция и скорость ее изменения связаны дифференциальным уравнением первого порядка, общим решением которого является экспонента. Экспоненциальная зависимость обладает тем свойством, что по трем ее дискретным отсчетам может быть восстановлена вся функция и определен предел А0, к которому стремится экспонента:

Где A1, А2, А3 - амплитуды сигналов в моменты измерений (отсчетов) соответственно t1, t2, t3 В свою очередь, постоянную времени процесса Т определяют из соотношения:

Δt - интервал времени наблюдений между отсчетами A1, А2.

По постоянной времени переходного процесса прогнозируют характеристики ожидаемого сейсмического удара. Время удара - это интервал времени, за который амплитуда сигнала, с вероятностью близкой к единице, достигает установившегося значения А0, для экспоненты tуст=4,7T (с вероятностью 0,99). Магнитуду удара определяют из соотношения Гутенберга-Рихтера: [см., например, «Краткосрочный прогноз катастрофических землетрясений с помощью радиофизических наземно-космических методов», Доклады конференции ОИФЗ им. О.Ю. Шмидта, РАН, М, 1998 г., стр. 10, стр. 13].

Пример реализации способа

Заявленный способ может быть реализован по схеме фиг. 5. Функциональная схема фиг. 5 содержит два пункта измерителей 1, 2, разнесенных в пространстве на расстояние измерительной базы 3. На каждом из пунктов установлено по две антенны 4, 5 и 6, 7 с диаграммами направленности во взаимно ортогональных плоскостях, оси диаграмм направленности антенн 5, 7 ориентируют по направленности базы 3, с подключенными к антеннам приемниками 8, 9 и 10, 11 соответственно. Гетеродины приемников синхронизированы от единого передатчика 12, работающего в импульсном режиме, подключенного к антенне зондирования ионосферы 4. Выход каждого из приемников подключен к тракту обработки из последовательно включенных порогового устройства и аналогово-цифрового преобразователя соответственно (13, 14), (15, 16), (17, 18) и (19, 20). Все тракты обработки подключены к канальному коммутатору 21, имеющему выход на вход компьютера 22 обработки результатов измерений в составе элементов: процессора 23, оперативно-запоминающего устройства 24, винчестера 25, дисплея 26, принтера 27, клавиатуры 28. В компьютер закладывают программу обработки и программу синхронизации работы элементов измерителей, которую пересылают в программируемую схему выборки измерений 29. Результаты обработки выводят на сайт 30 сети Интернет для передачи потребителям.

Взаимодействие элементов устройства при прогнозировании землетрясений состоит в следующем. Известно «Явление раскачки очага землетрясения перед сейсмическим ударом», Научное открытие №365, 2008 г. [см., Потоцкий В.В., Бюллетень, Научные открытия, РАЕН, М, СПб, 2009 г., стр. 66-68].

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

Применительно к одной из возможных реализаций (фиг. 2, 4) установившееся значение сигнала А0 (в шкале квантования 0…255 уровней) составило величину ~250. Значения A1 и А2 в интервале наблюдений Δt=t2-t1=8 час, соответствовали A1=125, А2=175. Откуда постоянная времени

Ожидаемое время удара ty=4,7T=70 час = 2,9 суток.

Ожидаемая магнитуда удара М≈7,1 балла

Направляющие косинусы гипоцентра очага землетрясения (фиг. 3) для первого пункта α=44°, для второго пункта β=62°.

Все элементы устройства представляют существующие технические разработки и средства аналогов. В устройстве использованы новые, по отношению к аналогам, элементы измерительной аппаратуры фирмы Bruel & Kjair, ENDEVCO (Дания) следующих моделей: канальный коммутатор, пороговое устройство, аналогово-цифровой преобразователь - многофункциональный блок, модель 3560-L. «Антенна для зондирования ионосферы», Патент Ru №2504054, 2014 г., широкополосная, работающая во всей полосе критических частот ионосферы, из двух, скрещенных в ортогональных плоскостях ромбов, подвешенных на опорной мачте из композитного материала, высотой 32 м, создающей геометрию главной диагонали ромбов, работающих в режиме бегущей волны.

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

Способ измерения ионосферных предвестников землетрясений, включающий регистрацию волн плотности электронной концентрации зондируемого слоя ионосферы на частоте ниже критической в виде дискретных цифровых отсчетов сигналов во взаимно ортогональных плоскостях Ax(t), Ay(t) в двух разнесенных на измерительной базе пунктах посредством антенн с диаграммами направленности каждого из пунктов в ортогональных плоскостях, при этом оси симметрии диаграмм направленности антенн по координате (x) ориентируют по направлению базы, обработку зарегистрированных выборок измерений сигналов, расчет направляющих косинусов вектора волн плотности электронной концентрации каждого пункта:

,

отождествление координат гипоцентра очага как проекции точки пересечения направляющих на земную поверхность, вычисление постоянной времени Т сейсмического процесса по изменениям амплитуд выборок измерений на интервале Δt:

,

расчет времени удара ty=4,7T и ожидаемой магнитуды как

где А0 - установившееся значение амплитуды сигнала предвестника, равное

Δt=(t2-t1)=(t3-t2) - интервал времени между отсчетами измерений A1, А2, A3.



 

Похожие патенты:

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

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

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

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

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

Изобретение относится к способам литохимических поисков золоторудных месторождений. Сущность: проводят отбор и физико-химический анализ проб.

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

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

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

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

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

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

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

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

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

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

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

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

Наверх