Способ и система для краткосрочного многопараметрического прогноза землетрясения

(51) МПК НАЦИОНАЛЬНЫЙ ЦЕНТР ИНТЕЛЛЕКТУАЛЬНОЙ СОБСТВЕННОСТИ СПОСОБ И СИСТЕМА ДЛЯ КРАТКОСРОЧНОГО МНОГОПАРАМЕТРИЧЕСКОГО ПРОГНОЗА ЗЕМЛЕТРЯСЕНИЯ(71) Заявители Ледовской Иван СергеевичЛедовской Александр Иванович(72) Авторы Ледовской Иван Сергеевич Ледовской Александр Иванович(73) Патентообладатели Ледовской Иван СергеевичЛедовской Александр Иванович(57) 1. Способ краткосрочного многопараметрического прогноза землетрясения путем определения времени, оставшегося до начала землетрясения, а также его магнитуды и координат его эпицентра, в котором с радиолокационных станций последовательно через заданные интервалы времени 1 излучают электромагнитные колебания заданной частоты в область атмосферы над районом предполагаемого землетрясения, каждый раз принимают колебания, отраженные от ионизированного образования в указанной области,измеряют время их распространения до ионизированного образования и обратно, а также приращение частоты и принятого колебания по сравнению с излученным, а затем определяют координаты указанного образования с учетом последнего из измеренных времен распространения и предварительно принимают за зону эпицентра участок земной поверхности под этим образованием, определяют первое значение ост 1 искомого времени в соответствии с выражением 16281 1 2012.08.30 ост 1( пи ) /( /1 ) ,где п - заранее найденное пороговое критическое приращение частоты принятого колебания по сравнению с излученным- приращение указанной частоты за интервал времени 1,и определяют первое промежуточное значение магнитуды землетрясения по заранее построенному графику зависимости приращения частоты и от магнитуды, далее с помощью прибора, установленного на космическом аппарате, известная орбита которого проходит за пределами ионосферы, фиксируют свечение в ионосфере газовых линий 5577 и 6300 и определяют второе положение эпицентра и второе промежуточное значение магнитуды на основе накопленных статистических данных по положению светящегося участка ионосферы и по его площади соответственно, и одновременно направляют с космического аппарата на поверхность Земли электромагнитные колебания, на поверхности Земли измеряют дальность до аппарата посредством измерения времени распространения этих колебаний, определяют величину ошибки измерения указанной дальности путем сравнения его измеренного значения с известным его значением, определенным по известным параметрам указанной орбиты, и при превышении этой ошибкой заранее заданного эталонного значения находят точку, лежащую на прямой распространения колебаний на высоте 100 км от земной поверхности, проецируют ее на земную поверхность и принимают точку проекции за третье положение эпицентра, затем уточняют положение эпицентра землетрясения путем сопоставления друг с другом всех его полученных значений и перемещаются в район эпицентра с уточненными координатами, где измеряют мощность инфранизкочастотных флуктуаций горизонтальной составляющей электрического поля Земли и определяют направление, вдоль которого указанная мощность максимальна, перемещаются вдоль этого направления с шагом 500-700 м в направлении возрастания инфранизкочастотных флуктуаций вертикальной составляющей электрического поля Земли и окончательно определяют как эпицентр землетрясения точку, в которой отношение второй из указанных мощностей к первой максимально, затем выбирают последовательность опорных точек измерения, удаленных от точки эпицентра на расстояния, кратные 500-700 м,измеряют мощность вертикальной составляющей инфранизкочастотных флуктуаций одновременно в эпицентре и первой опорной точке, аналогично повторяют измерения для каждой из опорных точек и находят радиусвзаимной корреляции, равный расстоянию между эпицентром и опорной точкой, в которой значение автокорреляционной функции мощности вертикальной составляющей указанных флуктуаций уменьшается враз по сравнению с эпицентром, и вычисляют третье значение М 3 магнитуды землетрясения в соответствии с выражением 3 э(/э ) ,где Мэ и э - соответственно эталонная магнитуда и эталонный радиус взаимной корреляции для данного района в спокойном состоянии,затем измеряют в точке эпицентра мощности инфранизкочастотных шумов в земной коре 1(1) и 1(2) на фиксированных частотах 1 и 2 соответственно, через заданный интервал времени 2 повторно измеряют указанные мощности и находят их новые значения 2(1) и 2(2), и затем определяют второе значение ост 2 искомого времени в соответствии с выражениемост 2 фп ( )ф/ф ,фп (1 ) где фп ( ) фп(1) и фп(2) - заранее найденные пороговые значения мощностей инфранизкочастотных шумов для данного района,а затем находят итоговое значение магнитуды землетрясения и времени, оставшегося до его начала, путем сопоставления друг с другом всех указанных выше промежуточных значений магнитуды и величин ост 1 и ост 2. 2. Система для краткосрочного многопараметрического прогноза землетрясения, содержащая наземные стационарные пункты прогноза, каждый из которых связан с радиолокационными станциями для определения параметров готовящегося землетрясения и положения его эпицентра, выполненными с возможностью определения координат ионизированных образований в атмосфере над районом предполагаемого землетрясения и частоты отраженных от них электромагнитных колебаний, спутниковые пункты прогноза,выполненные с возможностью наблюдения за свечением ионосферы над районом предполагаемого землетрясения, а также излучения электромагнитных колебаний на наземные пункты их приема с целью исследования параметров ионосферы и передачи на Землю полученной информации, и мобильные пункты прогноза, выполненные с возможностью определения параметров инфранизкочастотных шумов в земной коре, а также связанный со всеми наземными пунктами прогноза центральный пункт сбора и обработки информации,выполненный с возможностью определения времени, оставшегося до начала землетрясения, а также его магнитуды и координат его эпицентра путем сопоставления друг с другом данных, полученных указанными наземными пунктами. Изобретение относится к геофизике, конкретно к области исследования процессов,происходящих в земной коре (ЗК), и может быть использовано для оценки состояния сейсмической активности Земли и краткосрочного многопараметрического прогноза землетрясения (ЗМТ). Известен способ и система прогноза ЗМТ, заключающийся в установлении за пределами опасной зоны не менее двух синхронизированных по фазе СДВ-ДВ передатчиков и приемников. Трассы передатчик-приемник перекрывают опасные зоны и смежные районы сеткой с линейными размерами каждой из ячеек не менее 50-150 км. Осуществляют каждым приемным устройством регистрацию фазы сигналов всех передатчиков, измеряют по всем трассам абсолютное отклонение текущей фазы от средних значений на данное время суток, сравнивают отклонения текущей фазы с заданными пороговыми значениями и при превышении пороговых значений выделяют район готовящегося землетрясения по пересечениям воздушных трасс. Определяют время до ЗМТ на период до пяти суток 1, 2. Недостатками известного способа и системы являются низкая вероятность прогноза,обусловленная использованием одного признака (изменение фазы), сложностью выделения этого прогнозного признака из мешающих возмущений на огромных трассах распространения (например, на трассах Либерия-Омск, Реюньон-Омск, Норвегия-Омск и др.) сложность охвата трассами приемник-передатчик всех потенциально опасных регионов низкая глубина прогноза - до пяти суток переменная вероятность прогноза в течение суток, выше в ночное время и ниже в дневное время. Кроме того, конструкции СДВ сигналов фазовой радионавигационной системы (ФРНС) Омега в настоящее время практически выведена из эксплуатации. Наиболее близким по совокупности признаков к заявляемому является способ 3, где предлагается автоматизированная система мониторинга краткосрочных прогнозных предвестников землетрясений (ЗМТ) на базе интегрированных информационных полей станций ДВ ИФРНС, СНРС и аппаратуры потребления этих систем, суждении о готовящемся 16281 1 2012.08.30 ЗМТ по аномалиям и превышениям пороговых значений на трассах распространения на центральном пункте мониторинга. Недостатками известного способа являются низкая вероятность достоверного прогноза предвестников ЗМТ и отсутствие прогноза силы, места и времени наступления ЗМТ. Кроме того, радионавигационные системы ДВ и СДВ, на базе которых планировалась реализация предлагаемого способа, в настоящее время выводятся из эксплуатации. Задачей настоящего изобретения является расширение функциональных возможностей путем одновременного прогноза всех трех параметров ЗМТ место, время и сила, а также повышение достоверности правильного прогноза. Поставленная цель достигается тем, что в способе и системе прогноза землетрясения,основанном на измерении вариаций параметров прогнозных признаков и суждении о подготавливаемом землетрясении (ЗМТ) по изменениям параметров этих признаков о месте,силе и времени наступления события. Это достигается 1. В способе краткосрочного многопараметрического прогноза землетрясения (ЗМТ) путем определения времени, оставшегося до момента наступления ЗМТ, а также его магнитуды и координат его эпицентра, в котором с радиолокационных станций (РЛС) последовательно через заданные интервалы времени 1 излучают электромагнитные колебания заданной частоты в область атмосферы над районом предполагаемого ЗМТ, каждый раз принимают колебания, отраженные от ионизированного образования в указанной области,измеряют время их распространения до ионизированного образования и обратно, а также приращение частоты и принятого колебания по сравнению с излученным, а затем определяют координаты указанного образования с учетом последнего из указанных времен распространения и, предварительно принимая за зону эпицентра участок земной поверхности под этим образованием, определяют первое значение ост 1 искомого времени в соответствии с выражениемост 1( пи ) /( / 1 ) ,где п - заранее найденное пороговое критическое приращение частоты принятого колебания по сравнению с излученным- приращение указанной частоты за интервал времени 1,и определяют первое промежуточное значение магнитуды землетрясения по заранее построенному графику зависимости приращения частоты и от магнитуды, далее с помощью прибора, установленного на космическом аппарате, известная орбита которого проходит за пределами ионосферы, фиксируют свечение в ионосфере газовых линий 5577 и 6300 и определяют второе положение эпицентра и второе промежуточное значение магнитуды на основе накопленных статистических данных по положению светящегося участка ионосферы и по его площади соответственно, и одновременно направляют с космического аппарата на поверхность Земли электромагнитные колебания, на поверхности Земли измеряют дальность посредством измерения времени распространения этих колебаний, определяют величину ошибки измерения указанной дальности путем сравнения его измеренного значения с известным его значением, определенным по известным параметрам указанной орбиты, и при превышении этой ошибкой заранее заданного эталонного значения находят точку, лежащую на прямой распространения колебаний на высоте 100 км от земной поверхности, проецируют ее на земную поверхность и принимают точку проекции за третье положение эпицентра, затем уточняют положения эпицентра землетрясения путем сопоставления друг с другом всех его полученных значений и перемещаются в район эпицентра с уточненными координатами, где измеряют мощность инфранизкочастотных флуктуаций горизонтальной составляющей электрического поля Земли и определяют направление, вдоль которого указанная мощность максимальна, перемещаются вдоль этого направления с шагом 500-700 м в направлении возрастания инфранизкочастотных флуктуаций вертикальной составляющей электрического поля Земли и окончательно определяют как эпицентр землетрясения точку, в которой отношение второй из указанных 4 16281 1 2012.08.30 мощностей к первой максимально, затем выбирают последовательность опорных точек измерения, удаленных от точки эпицентра на расстояния, кратные 500-700 м, измеряют мощность вертикальной составляющей инфранизкочастотных флуктуаций одновременно в эпицентре и первой опорной точке, аналогично повторяют измерения для каждой из опорных точек и находят радиусвзаимной корреляции, равный расстоянию между эпицентром и опорной точкой, в которой значение автокорреляционной функции мощности вертикальной составляющей указанных флуктуаций уменьшается в е раз по сравнению с эпицентром, и вычисляют третье значение М 3 магнитуды землетрясения в соответствии с выражением 3 э 1 ( к /э ) , где Мэ иэ - соответственно эталонная магнитуда и эталонный радиус взаимной корреляции для данного района в спокойном состоянии,затем измеряют в точке эпицентра мощности инфранизкочастотных шумов (ИНЧШ) в земной коре 1(1) и 1(2) на фиксированных частотах 1 и 2 соответственно, через заданный интервал времени 2 повторно измеряют указанные мощности и находят их новые значения 2(1) и 2(2) и затем определяют второе значение ост 2 искомого времени в соответствии с выражениемост 2 фп ( )ф/ф ,где фп ( )ф фп(1) и фп(2) - заранее найденные пороговые значения мощностей инфранизкочастотных шумов для данного района,а затем находят итоговое значение магнитуды землетрясения и времени, оставшегося до его начала, путем сопоставления друг с другом всех указанных выше промежуточных значений магнитуды и величин ост 1 и т 2. 2. В системе краткосрочного многопараметрического прогноза землетрясения, содержащей наземные стационарные пункты прогноза, каждый из которых связан с радиолокационными станциями для определения параметров готовящегося землетрясения и положения его эпицентра, выполненными с возможностью определения координат ионизированных образований в атмосфере над районом предполагаемого землетрясения и частоты отраженных от них электромагнитных колебаний, спутниковые пункты прогноза,выполненные с возможностью наблюдению за свечением ионосферы над районом предполагаемого землетрясения, а также излучения электромагнитных колебаний на наземные пункты их приема с целью исследования параметров ионосферы и передачи на Землю полученной информации, и мобильные пункты прогноза, выполненные с возможностью определения параметров инфранизкочастотных шумов в земной коре, а также связанный со всеми пунктами прогноза центральный пункт сбора и обработки информации, выполненный с возможностью определения времени, оставшегося до начала землетрясения, а также его магнитуды и координат его эпицентра путем сопоставления друг с другом данных, полученных указанными наземными пунктами. Это позволяет, во-первых, одновременно прогнозировать все три параметра подготавливаемого землетрясения (ЗМТ) место, время, сила - с глубиной прогноза от 20 до 90 суток во-вторых, повысить вероятность правильного прогноза по всем параметрам до значений 0,850,99. 5 16281 1 2012.08.30 На фиг. 1 представлена структурная схема (блок-схема) практической реализации способа и системы, где приняты обозначения 1 - наземные пункты прогноза, реализующие прогноз всех трех параметров ЗМТ по параметрам и характеристикам инфранизкочастотных шумов (ИНЧШ) 2 - наземные радиолокационные пункты прогноза, реализующие прогноз параметров ЗМТ при работе их как в штатном, так и специальных режимах 3 - спутниковые системы прогноза, реализующие прогноз ЗМТ в режиме обмена информацией с наземными устройствами 4 - спутниковые системы прогноза, реализующие прогноз ЗМТ по свечению газовых линий в возбужденной ионосфере 5 - центральный пункт сбора и обработки информации (ЦПСОИ) 6 - линии связи. Сущность практической реализации способа и системы базируется на разработанной авторами теории и практике о триединой взаимосвязанной многопараметрической физической системе Земля-Тропосфера-Ионосфера, согласно которой при подготовке ЗМТ в окрестностях (зоне) эпицентрального (ЭПЦ) района всегда одновременно происходят аномальные изменения параметров и характеристик во всех трех физических средах диэлектрической и магнитной проницаемостей, удельной электрической проводимости,концентрации носителей зарядов, интенсивности газовыделения и газового взаимодействия и аэрозольных образований, напряженности электрических и магнитных полей происходят образования ионизированных структур (облаков-сгустков) в тропосфере эпицентрального района (ЭПЦР) и имеющих тенденции перемещения к ионосфере (ИОН) свечение в ИОН определенных газовых линий излучение ИОН как открытого колебательного контура на определенных резонансных частотах искажение полезной информации между спутниковыми радионавигационными системами (СРНС) отражения и модуляция электромагнитного колебания (ЭМК), излучаемого радиолокационными станциями (РЛС) различного диапазона частот и т.д. 4-10. Вариации названных и других параметров и характеристик в триединой системе при подготовке ЗМТ проявляются достаточно ярко, могут дистанционно и задолго до события быть обнаружены и измерены приборно после идентификации и приведения к единому времени и единым координатам могут быть использованы в целях краткосрочного и высоковероятностного прогноза всех трех параметров ЗМТ место, время, сила. Справедливость указанного следует из краткого рассмотрения вклада в прогнозные признаки каждой из составляющих впервые предложенной авторами триединой многопараметрической взаимосвязанной системы. Земля - третья от Солнца планета 11 со средним диаметром в 6371,032 км. Главными составляющими Земли, вносящими вклад в формирование прогнозных признаков подготавливаемого ЗМТ, являются земная кора, мантия (верхняя и нижняя), внешнее и внутреннее ядро в верхней части мантия содержит литосферу и астеносферу, процессы в которых связаны с тектоническими движениями, вулканизмом и др. Согласно 12-16, Земля не является сплошной (монолитной) структурой и в общем случае представима семью литосферными плитами Евразийская, Индийская, Тихоокеаническая, Австралийская, Африканская, Антарктическая, Американская. Размеры плит от сотен до нескольких тысяч километров и могут перемещаться 12-14. Причины глубинных сил, перемещающих плиты, не совсем ясны, но предполагается, что они возникают в результате медленных течений горячего пластического вещества в недрах Земли,которые вызываются тепловой конвекцией в сочетании с динамическими эффектами вращения Земли 14, 15. Во всех случаях это вещество поднимается наверх из недр Земли,оттесняя плиты в стороны, а в других местах плиты проскальзывают краями вдоль других есть также области, где одна плита заталкивается под другую (подобное происходило, например, при ЗМТ в Индонезии и в Японии 11 марта 2011 года). 16281 1 2012.08.30 Согласованности в движении плит не существует, поэтому они перемещаются в произвольном направлении, и это заставляет каменную толщу растрескиваться, создавая мелкие разрывы 17. При этом по обе стороны трещин отдельные слои не соответствуют друг другу как по физическим параметрам, так и характеристикам. Геофизики, геологи и сейсмологи такие нарушения сплошности пород, возникающие в результате движения плит, называют разломами. На местности визуально разломы практически не определяемы, их определение производят аппаратурными измерениями по тепловому контрасту, по электромагнитному излучению, по выходу радона и гелия из недр Земли и т.д. Согласно сложившейся классификации 12-17, крупные трещины между плитами называют глубокофокусными (глубинными) или региональными (субрегиональными), простирающимися по глубине от первых километров до десятков-сотен, а по ширине (от борта к борту) - от сотен до нескольких тысяч метров. Мелкофокусные разломы, иногда называемые локальными, имеют глубину до единиц километров, а ширину - до сотен метров. На первый взгляд все разломы на земной поверхности как бы срослись с окружающей поверхностью, но это не так. Имеют право на существование так называемые живущие (активные) разломы, которые находясь под воздействием флуктуационных энергетических,механических и др. явлений в земной коре (ЗК), взаимодействий с космосом и планетами,всегда находятся в режиме изменяющихся параметров и характеристик 12-17. Этому способствует то, что на один квадратный метр земной поверхности ежесекундно притекает энергетическая подпитка за счет солнечной энергии - 10-2103 Дж/м 2 атмосферного давления - 102 Дж/м 2 грозовых разрядов - 102 Дж/м 2 атмосферного электричества 10-5 Дж/м 2 гравитации Солнца и Луны - 10-1 Дж/м 2 и др. 14, 19. Естественно, что земная поверхность (ЗП) часть этой энергии принимает и растрескивается, заставляя разломы жить часть этой энергии рассеивается, а часть, накопившись и преобразовавшись в другие виды энергии (например, в электрическую и тепловую), через разломы возвращается на поверхность 14, 19. Укажем также, что изменение параметров разломов и накапливаемой в них, например,электрической энергии происходит и за счет воздействия на них упругих колебаний как за счет сейсмических волн проходящих, так и локальных ЗМТ, происходящих на Земле практически ежеминутно 13. Продолжая поиск наиболее информативных прогнозных признаков подготовки ЗМТ,авторы пришли к необходимости рассмотрения вариаций параметров и характеристик генетически связанных с Землей тропосферы и ионосферы через физическую среду - атмосферу. Атмосфера - газовая среда (оболочка), окружающая Землю и вращающаяся вместе с ней 11, 20-22. Эта среда разделена на отдельные оболочки, которые вследствие диффузии не имеют границ как с наружной стороны (Космос), так и с внутренней стороны - поверхность ЗК. Оболочки разделяют так Тропосфера. Определена до высот в 17 км в зависимости от континентов. Представляет собой азотную среду - 78 содержит 21 кислорода 0,9 аргона и в незначительных долях водород, гелий, углекислый газ и др. Стратосфера - азотно-кислородная среда - до высот в 80 км. Ионосфера - среда, в которой за счет космических (солнечных) лучей имеет место газовый фотоэффект, порождающий огромное количество заряженных частиц (носителей зарядов) - ионов. Здесь концентрация заряженных частиц в зависимости от слоев , ,с границами высот 70-80 110-140 км и более 200 км достигает значений 105 эл/см 3 106 эл/см 3 и более 19, 20. Все газовые оболочки ионизированы как от воздействия солнечных излучений, так и от погранично-электрических процессов конвекционных токов 19, 20 и грозовых разрядов 22. Средняя и нижняя часть атмосферы, кроме того, подвергается интенсивной ионизации за счет контактных напряжений, термоэлектрических, химико-электрических 16281 1 2012.08.30 процессов и других слоев в верхних слоях ЗК, а также радиации радиоактивности вод 23,что вытекает из того, что все естественные воды рек, озер, морей содержат растворимые соли радия, тория, распада радона и др. 23. Таким образом, в силу наличия носителей зарядов и взаимодействия атомов и молекул, содержащихся в газовой среде, она может быть охарактеризована электрической проводимостью, диэлектрической и магнитной проницаемостями, напряженностью электрического и магнитного полей и др. Следовательно, атмосфера за счет постоянной ионизации и наличия носителей зарядов является проводящей средой. При этом между положительно заряженной ионосферой (ИОН) и отрицательно заряженной поверхностью Земли всегда существует напряженность электрического поля (ЭП),которое всегда направлено нормально к Земле 19-21, т.е. под действием ЭП атмосферы интенсивностью в 1 В/см через слой газовой среды в 80 км от положительно заряженных слоев ионосферы (ИОН) происходит движение зарядов к отрицательно заряженной поверхности Земли. Возникающая при этом напряженность ЭП у поверхности Земли в зависимости от времени суток и регионов земного шара определяется значением в 100-130 В/м, а значение тока, стекающего на поверхность, равно 10-12 А/м 2 21. Фундаментальной основой Земли является ее незамкнутость 19, что означает ее постоянную открытость для воздействия, например, грозовых разрядов, несущих электрический заряд до 50 Кулон, создавая так называемое стороннее напряжение величиной до 103106 кВ 22. Поскольку в общем объеме земного шара ежесекундно происходит более 300 грозовых разрядов молний 22, понятен их вклад в так называемые земные токи электротеллурические токи (ЭТТ). Огромный вклад в эти токи, естественно, вносят экзогенные процессы (процессы на поверхности и в верхних слоях ЗК, обусловленные в основном энергией Солнца и деятельностью живых организмов) и эндогенные процессы(процессы, протекающие в недрах Земли, обусловленные энергией, выделяемой при развитии физико-химических процессов в них, и проявляющиеся в виде суммы сил температуры суммы сил давления суммы сил электрического происхождения и, конечно, в виде тектонических движений и т.д.) 19, 20, 22, 24, 28. Практические исследования и измерения плотности ЭТТ 19, 20 дают значения в пределах 10-6104 А/м 2. Становится понятным, что совокупное воздействие перечисленных сил приводит к возбуждению физико-химических процессов в Земле и (что очень важно) к образованию в ней огромных энергетических источников и, в первую очередь, в зонах литосферы, в ЗК и в верхней мантии 19. Этому способствует также и то, что Земля как геологическая структура, состоящая из массивов пород, не имеет сплошности и имеет огромное количество разнофокусных разломов и трещин, а свойства ЗК определяются минералогическим составом с разнообразием по полидисперсности, химическому составу, многокомпонентности, раздробленности 12-20, что усиливает образование и развитие электрических процессов за счет естественного насыщения трещин влагой и естественные контакты с потенциало-образующими системами 19, 20 контакты двух пород, обладающих электронной проводимостью контакты водонасыщенных пород контактные напряжения типа пятно и дырка 20 ЭП контактов горных пород (ГП) фильтрационно-электрическое напряжение водопадов ЭП жидкого потока речное ЭП диффузионно-электрическое поле кинетическое ЭП речного потока и др. Важно, что подземные воды и флюиды (жидкие и газообразные легкоподвижные компоненты магмы с содержанием хлора, водорода, фтора, углеводорода) 23, 24, являясь легкоподвижными компонентами ЗК, фильтруясь в ее тонких порах и трещинах, вносят огромный вклад в формирование заряженных конструкций в ЗК. Обмен зарядами в подобных структурах ведет к разделению зарядов и образованию двойных электрических слоев (ДЭС) или областей разделенных зарядов (ОРЗ) 19. Поскольку толщина ДЭС 16281 1 2012.08.30 скачке потенциала на границе раздела фаз в один вольт напряженность ЭП составляет десятки миллионов В/см 19, 25. Заполнение влагой и воздухом полостей на стыках твердых минералов, а также трещин и пор повсеместно соприкасающихся компонентов с различным химическим составом приводит к тому, что породы оказываются сложными и многокомпонентными физико-химическими системами с огромным количеством (ДЭС). Другими словами, наличие любого контакта двух сред вызывает диффузию носителей зарядов, а также адсорбцию ионов 19, 24 (адсорбция наблюдается при контакте твердой среды как с жидкой, так и газообразной средой), что приводит к образованию на контактах систем устойчивых ДЭС, внутри которых сосредотачиваются мощные ЭП 19, 24. Очевидно, что группирование ДЭС в локальных пространствах (объемах участков ЗК) превращает ДЭС в термодинамически неустойчивую многослойную область, а всю систему в энергоактивное образование, на которое, согласно указанному ранее, воздействуют как экзогенные, так и эндогенные силы и силы планетарного масштаба. На континентах в общем случае выделяют 1, 20, 24 два основных типа энергоактивных зон (ЭАЗ) электрокинетической и полигенной природы. ЭАЗ электрокинетической природы представляют собой первый от поверхности капиллярный слой ЗК мощностью 0,1-5,0 м, ограниченный по площади областью высокой водонасыщенности. Здесь наблюдается аномальная суточная активность удельного электрического сопротивления, диэлектрической и магнитной проницаемостей, электрохимической активности и др., что деформирует ДЭС и интенсифицирует процесс электротермокапиллярного эффекта и вариации напряженности ЭП при изменении температуры 19. Электробарокапиллярные потенциалы в приповерхностном слое ЗК также варьируются в результате смещения диффузной части ДЭС, и при изменении давления напряженность ЭП также существенно изменяется. Действие экзогенных, эндогенных сил и сил планетарного масштаба подобно действию, что и в зонах электрокинетической природы, но приуроченных к природным разломам, вызывает более интенсивные явления 19, 20, 24. Так, под действием только одной из составляющих планетарных сил (например, сил давления) внутритрещинное и внутрипоровое пространство изменяется минимум на величину толщины ОРЗ (ДЭС) 10-8 м и этого становится вполне достаточно, чтобы в тонких трещинах и порах во время сжатия(например) произошло массовое перекрытие ДЗС. Это приводит к возникновению нелинейного эффекта взаимной поляризации и других свойств в ЗК и, что главное - ведет к образованию конденсаторов с огромной емкостью и каналов высокой проводимости,концентрирующих электротеллурические и другие токи в Земле 19. Размеры ЭАЗ здесь определяются величинами в 0,1-20,0 м, а диэлектрическая проницаемость может достигать величин более 104 (для сравнения - для воды она равна 80), превращая энергоактивные образования в конденсаторы с огромной емкостью для накопления геоэлектрической энергии огромных величин. Локализованная здесь напряженность ЭП достигает пробойных для диэлектриков значений - 105108 В/см 19, 25. В таких условиях развития физико-химических процессов по образованию и развитию ОРЗ ЭАЗ происходит 1. При воздействии планетарных и экзоэндогенных сил, а также процессов сжатия и разряжения (например, при прохождении сейсмических (упругих) волн от проходящих ЗМТ) ОРЗ (ДЭС) стремятся к объединению, образуя мощные энергетические источники. 2. Объединяющиеся ОРЗ формируют конструкции (зоны) с повышенными значениями напряженности ЭП, вплоть до 108 В/см. 3. Под действием высоких параметров ЭП образуются каналы высокой проводимости,в которых диэлектрическая проницаемость растет до значений 104 и более 25. 4. Образующиеся каналы высокой проводимости (малого сопротивления и затухания),во-первых, создают возможность передачи (трансляции) на большие расстояния как по го 9 16281 1 2012.08.30 ризонтальному, так и вертикальному направлениям как тепловой (Кормянский котел. Советская Беларуссия. 20.03.2010), так и электрической энергии (Курильские огни. Труд. 13.07.1974) во-вторых, при определенных естественных формах и размерах разломов излучают электромагнитные колебания (ЭМК) 19 в-третьих, если проследить за развитием физико-химических процессов в разрезе от поверхности Земли до ее ядра, то становится очевидным, что в результате взаимодействия между внешним и внутренним ядрами происходит подток жидкого вещества 19, 20, 24, 25 с образованием высокоэнергетической субстанции, выходящей через так называемый мантийный канал на поверхность Земли. Такой процесс Э.В. Бороздич 27 назвал короткопериодным локальным возмущением(КПЛВ). Выходящая при этом из недр Земли высокоэнергетическая субстанция представляет собой сгусток электрических зарядов (плазма), которые, распространяясь в атмосфере у земной поверхности, вызывает свечения, образования в форме шаров, гантелей и др., достигая ИОН-свечения кислородных линий. Подобные процессы происходят нечасто и, перемещаясь в атмосфере, сгустки быстро компенсируются зарядами противоположного знака. Согласно изложенному, в высокоэнергетических областях диэлектрическая проницаемость, достигая величин более 104, способствует глубоко залегающим астеносферным зонам преобразовываться в конденсаторы с огромной емкостью, которые, накапливая электрическую энергию, создают плотность тока не менее 100 А/см 2, что при такой величине диэлектрической проницаемости среды и воздуха в 1/3610 Ф/м плотность геоэлектрической энергии достигает значений примерно в 100 Дж/м 2. Выход на поверхность Земли энергии такого порядка несет беду - Чернобыль, Армения (Спитак) и др. 5. Землетрясение - результат разрыва или нарушения сплошности как коровых, так и мантийных конструкций (слоев) ЗК с образованием так называемых дислокаций - участков ЗК с нарушением первичного залегания ГП 15, вызываемых как тектоническими, так и магматическими процессами, в которых огромную роль играют электрические процессы. Простейший пример свершения разрыва сплошности участка ЗК можно проиллюстрировать на резиновом жгуте. Представим себе резиновый бинт (жгут) шириной 5-6 см,толщиной 1-2 мм с незначительным дефектом (пусть сделанным специально) и станем руками растягивать этот жгут перед грудью. Процесс растягивания можно довести до такого состояния, когда эластичность жгута достигнет почти предела, а растягивающая сила продолжает прикладываться. Возникнет ситуация, когда жгут исчерпает пределы сопротивляемости (прочности) и начнет растрескиваться (разрываться) по границам дефекта,лавинообразно создавая новые дефекты и при дальнейшем приложении растягивающей силы - удар по рукам, наступает полный разрыв жгута. Применительно к ЗМТ аналогом резинового жгута может служить участок земных недр, например мантии, с дефектами - разломами, расчленяющими его на отдельные спаянные блоки. Очевидно, что для разрыва спаек между блоками энергия должна быть огромной и, конечно, процесс накопления энергии для разрыва и поднятия (сдвига) блоков является очень длительным и зависит от многих факторов. При прогнозе ЗМТ важно определить механизм, который определяет спусковой механизм - время освобождения накопленной энергии, ее величины и места освобождения. Если таких механизмов в применении к прогнозу ЗМТ нет, возможность предсказания ЗМТ становится весьма проблематичной 13-15. В настоящее время определено огромное количество прогнозных признаков и предвестников ЗМТ 1, 5, 6-10, 12-16, 23, но до настоящего времени нет ни одного надежного способа (метода) прогноза. Примером являются последние ЗМТ в Японии, в Туве (Россия), Индонезии, которые не были спрогнозированы. Это объясняется тем, что предлагаемые способы базируются, как правило, на основе одного предвестника, что априори предполагает исключительно малую вероятность правильного прогноза. В 13 указывается, что для прогноза параметров подготавливаемого события с высокой вероятностью необходимо знать точные размеры сейсмогенерирующего, например, разлома, прочность 16281 1 2012.08.30 вмещающих пород, шероховатость стенок разлома, силу, с которой они прижаты друг к другу, и множество других количественных характеристик, которые практически определить невозможно. В соответствии с этим авторы пришли к тому, что необходимы обобщенные показатели,характеризующие процессы подготовки ЗМТ, а прогноз должен быть построен на основе многопараметричности. Именно такой подход предлагается реализовать в настоящем изобретении, где на основе предложенной многопараметрической взаимосвязанной триединой физической системы реализуются измерения вариаций параметров, изменяющихся одновременно при подготовке ЗМТ, что впервые задолго до момента наступления ЗМТ позволяет прогнозировать все три параметра место, время и сила - с высокой вероятностью. Кратко покажем механизм развития процесса подготовки ЗМТ с учетом того, что точкой начала подготовки и развития события является гипоцентр, а ее проекция на земную поверхность Земли - эпицентр. Допустим, что в некоторой точке разломной зоны земных недр начинает формироваться гипоцентр, возникновение которого обуславливается проявлением в некоторой ослабленной области действием эффекта Ребиндера 26 (ослабление прочности ГП более чем в два раза по отношению к окружающей среде в результате действия адсорбции при контакте пород с влагой). Это в районах глубинных разломов способствует формированию энергоактивных зон (ЭАЗ), накапливающих высокоэнергетические субстанции. Здесь в ОРЗ (ДЭС), как в межфазовых поверхностях, происходит локализация значительной доли всей массы и свободной энергии полидисперсных многофазных систем. Это превращает ОРЗ в термодинамически неустойчивую многослойную область, а всю систему в мощное энергетическое образование. Подобные энергообразования способны к быстрому изменению электромагнитных свойств под действием экзо- и эндогенных сил, а также сил планетарного масштаба, которые деформируют ОРЗ ЭАЗ, вызывая при этом различные виды поляризации 15, 19, 25, 26. Важно, что действие только одной из сил планетарного масштаба (например, давления) оказывается достаточным, чтобы в тонких трещинах и порах Земли во время сжатия, например, произошло перекрытие ОРЗ. В результате в течение короткого времени может возникнуть (и возникает) эффект активной взаимной поляризации 19, который в свою очередь приводит к аномальной динамике электрических и других свойств в Земле, а также к образованию каналов высокой проводимости,концентрирующих все токи, действующие в окружающей среде 19. Здесь при определенной объемной доле концентрации ОРЗ диэлектрическая проницаемость (это оговаривалось ранее) нелинейно достигает величин 104, превращая энергообразования в огромные емкостные накопители геоэлектрической энергии 25 локализованная в ОРЗ напряженность ЭП достигает пробойных для диэлектриков величин в 106108 В/см как вдоль, так и вкрест границ раздела фаз. Огромное значение в этом случае приобретают флюиды, которые в высокотемпературном поле вздымаются вверх по раскрывающимся разломам, устремляясь по пути наименьших потерь, производя структурно-вещественные преобразования слоев ЗК и мантии с образованием заряженных дислокаций. Здесь в ЭАЗ с вихревым развитием огромных значений токов возникают огромные температуры, приводящие к формированию сверхвысоких давлений с образованием газовых фракций, стремящихся к выходу на поверхность Земли. Возникает своеобразный взаимосвязанный и практически лавинообразный процесс с активизацией геоэлектрических процессов в очаге, приводящих также к лавинообразному возрастанию напряженности ЭП, которое порождает здесь возникновение искрений, свечение предметов и формирование плотных электрических сгустков, которые в последующем, объединяясь, перемещаются в атмосфере, повышение температуры, давления и др. Развиваясь, очаг в процессе формирования выделяет среднюю мощность, определяемую 40 выражением,4 р- разность потенциалов ДЭС оч - емкость очага р - время релаксации. Примерный размер очага можно определить из полуэмпирического выражения оч 103,51,8,где- магнитуда. Очаг долго развивается до некоторого критического состояния, при котором трещины начинают объединяться и расширяться, достигая критических значений температуры,давления, накопленной электрической энергии и др. Максимальное время для достижения, например, порогового значения давления в очаге также можно определить 40 через магнитуду ожидаемого ЗМТп 3,710 9,51,8 . Ожидаемые значения напряженности ЭП, плотности токов и напряжений можно также определить через параметры очаговой зоны. Например, ток утечки находится 34 из выражения 4 о 10 7, 00,3612,о 2 12 где о - глубина залегания очага о - разность потенциалов между очагом и поверхностью к - проводимость- коэффициент. В общем случае на больших глубинах, соответствующих области залегания магматических очагов в поле высоких значений диэлектрической проницаемости, т.е. при 104 19,25, плотность накопленного заряда составляет нз 1 ттнз 10 7 К / см 2 , где 1 тт 1011 А / см 2 . Такая плотность заряда создает напряженность поля накопленного заряданз 1 Ф / м 10 4. 36 109 Соответствующая этому плотность тока достигает величины 10 2 А / см 2 , где мо 10 Ом/м - удельное сопротивление магматического очага. Естественно, напряженность поля в 100 В/см создает геоэлектрическую энергию,среднее значение которой равно Значение этой энергии велико. Одна часть этой энергии в очаге переходит в тепловую и механическую энергию, а другая часть по глубинным разломам выносится на поверхность Земли и далее в Атмосферу. При этом величина плотности накопленного заряда ограничивается напряженностью пробоя конденсатора ОРЗ - Еп 106 В/см, где п 1, 2,3, , и постоянно возрастает, увеличивая свое значение в зависимости от параметров и характеристик конкретного геофизического энергообразования. Таким образом, совокупное взаимодействие всех приведенных факторов ведет к активизации как всех физико-химических процессов в очаге, так и процессов от воздействия экзо- и эндогенных сил и сил планетарного масштаба, а также упругих (сейсмических) волн проходящих ЗМТ. Результатом этих взаимодействий становятся процессы повыше 12 16281 1 2012.08.30 ния напряженности ЭП, температуры, давления, которое в очаге может достигать нескольких десятков тысяч бар, и др. При развитии названных выше процессов в районах ЭАЗ разломов интенсивно возникают и проявляются в недрах Земли каналы проводимости, во-первых способные передавать электрическую энергию на достаточно большие расстояния, аккумулируя ее до огромных значений во-вторых способствующие превращению ее в тепловую энергию нагрева, например, воды, а также плотные электрические сгустки и выводу их на поверхность. При этом за счет объединения ОРЗ в локальные ЭАЗ формируются конструкции с высокой напряженностью ЭП вплоть до 108 В/см в Земле образуются участки с отрицательным сопротивлением - эффект Зенера проявляется пинч-эффект, т.е. эффект сжатия (фокусирования) 19 носителей зарядов(например, электронов) в узкий шнур (пучок) сечением в 10-20 м (подобное наблюдалось при ЗМТ в Армении 5-7 - наблюдения авторов и при аварии на ЧАЭС (Невский А.М. Газета Республика. 2.09.1999 г., что намного меньше сечения вмещающих природных разломов и т.д. Выход на поверхность высоких значений вертикальных токов и повышение напряженности ЭП в приповерхностной атмосфере способствует ускорению движения свободных носителей зарядов, повышая их кинетическую энергию до величины, достаточной для ударной ионизации молекул окружающего воздуха. Вследствие этого появляются вторичные носители зарядов, что приводит к ионизационному или газовому усилению 29. Очевидно, что совокупные механизмы выхода на поверхность газов, механизм газового усиления и интенсивный выход на поверхность Земли мощных электрических образований в приповерхностной атмосфере формируют плотные ионизированные сгустки, которые, стремясь к объединению, образуют зарядовые образования различной формы шары различного размера, эллипсы, дугообразные образования и др. Названные сгустки плазмоиды - имеют огромную плотность заряженных частиц (в том числе и заряженных частиц пыли, захватываемых с поверхности Земли) с концентрацией от 1010 до 1016 эл/см 3. Сформировавшиеся в районе эпицентра (ЭПЦ) плотные ионизированные образования(плазмоиды) перемещаются в сторону ионосферы, вступая с ней во взаимодействие 5, 6,7, 8, 10, 55, 56, возбуждая в ней свечение газовых линий кислорода, а также излучение электромагнитных колебаний открытом колебательным контуром Земля - ионосфера. Кроме того, ионизированные образования, имея высокие концентрации носителей зарядов, отражают ЭМК, излучаемые РЛС и распространяющиеся в эпицентральном районе с ионизированной средой 5, 6, 7-10, 55, 56. Все изложенное и дальнейшие представления авторами как визуально, так и приборно наблюдались многократно и в различных регионах бывшего СССР и показываются ими реальными результатами на конкретных ЗМТ. В 11 часов 41 минуту (по местному времени) 7 декабря 1988 года на территории Армении (район Спитака) произошло катастрофическое ЗМТ. Никто из специалистовсейсмологов и различного рода прогнозистов наступления его не предсказал, хотя оно сигнализировало о процессе подготовки различными признаками в течение нескольких месяцев. И это несмотря на то, что помимо сейсмопунктов на территории Армении и Грузии располагались сейсмообсерватории Тбилиси, Ереван сейсмостанции Ленинакан,Кировакан, Эчмиадзин и др. Авторы впервые в 1984 году теоретически, в лабораторных условиях и на местном полигоне с помощью малогабаритной РЛС 1 РЛ 336 А практически отработали на локальных тротиловых взрывах и электрофорных машинах и показали возможности прогноза ЗМТ с помощью РЛС. Но природа опередила нас. Переждав ажиотаж после ЗМТ, авторы отправились в Армению, где проработали с 15 по 28 января 1989 года. Исследования проводились по трассе Тбилиси, Марнеули, Болниси, Казарети, Калинино, Степанаван, Спитак, Ленинакан, Алагяз и обратно. Никаких 16281 1 2012.08.30 задач априори воинским формированиям, дислоцированным в этом регионе, на наблюдения за процессами подготовки ЗМТ нами не ставилось. Исследования первых же документов объективного контроля (фотографий) воздушной обстановки на выносных индикаторах кругового отображения (ВИКО) показали правильность наших наработок. На ВИКО РЛС признаки подготовки ЗМТ начали фиксироваться более чем за три месяца до момента ЗМТ, т.е. с 11 августа 1988 года. ЗМТ произошло 7.12.1988 г. (более ранних наблюдений авторам предоставлено не было) и продолжало наблюдаться (отображаться на ВИКО) как до момента свершения события(7.12.1988 г.), так и при последующих афтершоках. Процессы были исследованы на 16-ти РЛС, работающих на различных частотах, а также на радиорелейных и связных радиостанциях. Подтверждено актом, подписанным командованием войскового соединения 27 января 1989 года. Анализ результатов наблюдений и представленных материалов наблюдений на РЛС,работающих как в активном, так и пассивном режимах, показал 1. Процесс развития подготовки ЗМТ в Спитаке начал регистрироваться на ИКО РЛС почти за 90 суток (11.08.1988 г.) до свершения события в виде искровых разрядов (светлые точки-блестки) и плотных энергетических засветок (ионизированные светлые образования - плазмоиды) - засветки по всему экрану - фиг. 2. 2. По мере развития процессов подготовки ЗМТ обстановка существенно изменилась по состоянию на 11 часов 15.08.1988 г. на экране ВИКО отображались как искрения(центр экрана), так и плазменные образования (вблизи центра экрана), а также светлые полосы по всему диаметру экрана - фиг. 3 по состоянию примерно на 10 часов 26.08.1988 г. все указанные процессы значительно усилились - фиг. 4 а на фиг. 5 по состоянию на 10 час. 35 мин. 27.08.1988 г. отображаемые процессы подобны, указанным на фиг. 2, но отличаются, во-первых, плотностью ионизированных образований, во-вторых, движением заряженных плазмоидов к ИОН - граница экрана ВИКО. На фиг. 6 показана картина отображения процесса развития ЗМТ по состоянию на 13 час. 07 мин. 30.08.1988 г., где на ВИКО отображены светлые полосы по всему диаметру экрана и при работе РЛС в пассивном режиме работы РЛС, т.е. при работе только на прием без включения передатчика на излучение и приема сигналов от отражающих объектов - плазменных образований. Отображаемые на экране светлые полосы - активные помехи, т.е. электромагнитные излучения открытого колебательного контура Земля - ионосфера, образующегося в процессе подготовки ЗМТ 5-8, 10 - работы авторов . 3. Результаты развития процессов подготовки ЗМТ отображались не только на ВИКО РЛС светились горные вершины перемещались светящиеся объекты в виде шаров, цилиндров, эллипсов и др. в атмосфере эпицентрального района. Плазмоиды, достигая ИОН,возбуждали ее, а отклик отображался на экранах РЛС в виде светлых полос по всей дальности, принимаемых РЛС как помехи как по прямому, так и по боковым лепесткам диаграммы направленности антенн РЛС - фиг. 3, 4, 6 5, 6, 7, 9, 10. 4. Образующиеся в процессе широкомасштабной подготовки ЭМТ плотные пучки электронов с определенной энергией, достигая ИОН, вызывают ее возбуждение и интенсивное свечение на частотных спектральных линиях кислородных 5577 и 6300 Ангстрем. Авторами с помощью разработанной ими оптико-электронной аппаратуры, установленной на орбитальной станции Салют-4, исследовано и экспериментально показано, что изменяющаяся концентрация зарядов в ИОН, в процессе подготовки ЗМТ, приводит к нарушению распространения электромагнитных волн (ЭМВ) и потере информации, передаваемой с космических аппаратов. Авторы и это явление предложили использовать как прогностический признак в качестве оценки вероятности проявления наступающего ЗМТ. 5. В процессе развития события в период с 11 августа по 7 декабря 1988 года наблюдалась своеобразная периодичность в формировании и проявлении на индикаторах круго 14 16281 1 2012.08.30 вого обзора (ИКО) с интервалом в 2,5-3,5 часа. При этом существенно изменялись приращения как по фазе, так и по частоте и амплитуде 5, 8. Изложенное подтверждено наблюдениями и исследованиями на конкретных событиях в Грузии, Киргизии, Камчатке и др. На фиг. 7 и 8 в качестве примера показано отображение на ИКО процесса развития ЗМТ на территории Грузии, состоявшегося 19.04.1989 г. Из изложенного следует предлагаемый авторами реальный способ прогноза всех трех параметров ЗМТ место, время, сила - с помощью РЛС. Радиолокационный способ прогноза ЗМТ Прогноз эпицентра ЗМТ. Операция прогноза эпицентра ЗМТ с помощью РЛС заключается в определении координат дальности и угловых координат при работе РЛС в штатном режиме классическими методами 30, 31, 32. Это означает, что как при реализации этой операции, так и при прогнозе других параметров ЗМТ РЛС, выполняя свои штатные задачи по наблюдению за воздушным пространством во вверенном секторе ответственности, обнаруживают отображаемые на ИКО аномальные проявления в виде искрений,плотных засветок и светлых радиальных полос - активные помехи (фиг. 2-8). Принимаются решения на распознавание и идентификацию (такие алгоритмы имеются в каждой РЛС) и после обработки соответствующих образов принимаются решения на отнесение засветок к признакам, характеризующим начало или течение процесса подготовки ЗМТ. В зависимости от магнитуды подготавливаемого ЗМТ и глубины залегания очага проявление признаков подготовки могут обнаруживаться и отображаться на экранах ИКО РЛС по срокам от нескольких десятков суток до нескольких месяцев. Наблюдения авторов показывают ЗМТ в Армении (Спитак) с магнитудой 9,0-9,5 баллов по шкале Рихтера начали отображаться на ИКО РЛС с августа 1988 года и продолжались до 7 декабря 1988 года в Грузии процесс подготовки ЗМТ фиксировалея с 1 апреля по 19 апреля 1989 года (ЗМТ 4,6 балла) в Киргизии при магнитуде в 4,5 баллов (август 1992 года) процесс подготовки наблюдался за 12 суток. Подобные наблюдения имеются по Камчатке и др. Преимущество способа, впервые разработанного и предложенного 6, 7, 8, 55, 56 авторами, состоит в том, что РЛС, например, ПВО не отвлекаются от своих штатных задач,прогноз проводится одновременно. Действительно, в режиме работы РЛС имеет место так называемый режим контроля функционирования станции (КФС), который предполагает проверку работоспособности всего комплекса аппаратуры. Чаще всего этот процесс реализуется в 9.00 и 21.00 час и длится менее 30,0 минут. Именно в эти периоды могут быть обнаружены аномальные процессы в атмосфере контролируемого сектора. В случае положительного результата начинают определять координаты (дальность и азимут) эпицентра ЗМТ и др. Дальность до эпицентра (ЭПЦ) - до скопления в атмосфере ионизированных образований, концентрация которых здесь может достигать значений от 1010 до 1016 эл/см 3, определяется временем распространения излученного РЛС ЭМК до ионизированного образования и обратно, т.е. Д/2. Азимут (угловые координаты) ЭПЦ (ионизированного образования) определяется по максимуму сигнала, принимаемого антенной системой РЛС, ширина диаграммы направленности антенны РЛС обычно не превышает 2-10 градусов 31-32. Поскольку ионизированные образования не сосредоточены в точке, а распределены по пространству, то применительно к конкретным отображениям на ИКО РЛС (фиг. 3, 4, 5) азимут определяется при просмотре образования с левой и с правой стороны, где амплитуда принимаемого сигнала спадает до значения 0,7 от максимума, суммарный результат делят на 2 - это и есть азимут ЗПЦ ЗМТ. С целью облегчения определения как дальности, так и азимута на экранах ИКО РЛС наносят координатные сетки по дальности - это концентрические окружности с определенным шагом по дальности (например, в 10, 20 км и т.д.) для азимута - это азимутальные секущие, наносимые на экране, например, через 10 градусов, фиг. 2-8. Долговременность 16281 1 2012.08.30 процесса подготовки ЗМТ позволяет многократно производить обзор пространства и однозначно отселектировать отображение на экранах, например, от молниевых разрядов (если они имеют место в периоды измерений), которые кратковременны 22, от энергетических проявлений ЗМТ и путем уточнения определиться в правильности обнаружения координат ЭПЦ. Все эти действия операторы РЛС производят практически автоматически. Определение времени наступления ЗМТ. В общем случае как радиолокационную цель 30-32 ионизированное образование следует рассматривать как протяженную цель, т.е. как объект, состоящий из набора элементарных отражателей. Для таких целей отраженное ЭМК при приеме РЛС можно записать в виде выражения где- характеристика диаграммы направленности РЛС 2 - волновое числодлина волны- расстояние- амплитуда отраженного ЭМК. Очевидно, что в этом выражении имеет место частота, которую можно выделить как математически, так и инструментально согласно следующему 2 где- амплитуда сигнала на -й частоте- ширина полосы пропускания -го фильтра на -й частоте- количество фильтров. Общая полоса пропускания устройства фильтрации обычно выбирается из условий необходимой точности измерений, связанной с вероятностью правильного прогноза, т.е.в-н 31, 32. Поскольку колебание в каждом узкополосном фильтре является квазигармоническим,то закон его изменения через напряжение можно описать выражением,где- начальная фаза. Беря производную от этого выражения, получаем,где- амплитуда производной напряжения на выходе каждого -го фильтра. Учитывая обозначение, получаем Согласно этому выражению авторами синтезирована и практически реализована схема, приведенная на фиг. 9, которая при необходимости подключается к выходу усилителя промежуточной частоты приемника РЛС и работает согласно алгоритму. Отраженные от ионизированных образований ЭМК через антенну 1 поступают на устройство 2, преобразуются с помощью смесителя 4 и опорного генератора 5 и через усилитель 6 поступают на систему узкополосных фильтров 7.17 п и с последних на измерители амплитуд напряжений 8.18 п и далее на первый сумматор 9, а также на устройства дифференцирования 10.110 п, после которых - на устройство измерения 11.111 п и далее на сумматор 13. Далее с сумматоров 9 и 13 сигналы поступают на схему отношений 12, где формируется напряжение, пропорциональное значению частоты. Это напряжение поступает на устройства 14 и 15. В устройстве 14 производится сравнение полученного результата с соответствующим пороговым значением, проградуированным в размерности частота, и делается вывод о времени, оставшемся до момента наступления ЗМТ согласно выражению 16281 1 2012.08.30 Очевидно, что при организации прогноза с применением РЛС предварительно решается задача по накоплению статистического материала и установлению величины п. Эта задача решается для каждого сейсмоактивного региона на основе наблюдений как локальных, так и проходящих ЗМТ, а также по афтершокам их. Такая задача достаточно просто решается, например, для районов Камчатки, Сахалина, Кавказа, где ЗМТ с магнитудой 44,5 балла происходят очень часто. Для накопления статистики пригодны также афтершоки произошедших ЗМТ, поскольку они после свершения ЗМТ проявляются с большой частотой и большой магнитудой. Так, после свершения основного события в Спитаке в период с декабря 1988 г. по конец января 1989 г. сейсмостанциями Кавказа было зафиксировано более 170 афтершоков при ЗМТ в Японии в период с 11 марта по 11 апреля 2011 г. отмечено более 700 афтершоков, среди которых наблюдались и опасные ЗМТ с магнитудой 6,0 и 7,5 баллов. Естественно, что подобный рой позволяет достаточно точно определиться с законом изменения приращения частоты отраженного ЭМК в период подготовки ЗМТ и произвести запись в память вычислителя необходимого порогового значения. При необходимости определения приращения частоты для больших магнитуд производится экстраполяция. Аналогично можно реализовать прогноз времени, оставшегося до наступления ЗМТ, с помощью РЛС, работающих в импульсном режиме, т.е.(2)ост 1 ( ииип ) / ( и /1 ) ,где ии, ип - измеренное и пороговое значения длительности импульса и приращение длительности импульса ЗМТ. Последнее объясняется тем, что при подготовке ЗМТ среда распространения ЭМК оказывается ионизированной диспергирующей, а следовательно, приводящей к изменению коэффициента преломления и зависимости его от частоты 31. Это означает, что при прохождении импульсного колебания через возбужденную среду эпицентрального района подготавливаемого ЗМТ будет происходить искажение формы импульсов за счет различного влияния среды на частоты, формирующие форму импульсов 31. Прогнозирование магнитуды ЗМТ. В настоящее время как и все параметры ЗМТ определяются после свершения события, методов такого определения магнитуды наработано достаточно много 12-16, 23, но методов прогноза магнитуды априори (до свершения события) практически нет. Это доказало и ЗМТ с афтершоками, произошедшее на территории Тувы в декабре 2011 г. - январе 2012 г. Авторы впервые предлагают практически работающие способы прогноза магнитуды ЗМТ до его свершения. Это вытекает из конкретных результатов наработок по реальным событиям, на основе которых показаны возможности прогноза параметров подготавливаемого ЗМТ, в том числе и магнитуды, с использованием РЛС 5-8. При этом опираясь на собственные работы как по исследованию образцов горных пород 43-А.с.19408286 и исследования в безэховой камере с созданием ионизированной среды совокупностью электрофорных машин и контроля изменения параметров ЭМК (в том числе и частоты),распространяющихся в ней, с помощью мини-РЛС типа 1 РЛ 136 А. Наблюдения за реальными событиями в Армении, Грузии и др. позволили однозначно подтвердить, что частота ЭМК, отраженных от ионизированных образований в ЭПЦ подготавливаемого ЗМТ, связана с магнитудой. Это следует также и из патента авторов 6565 РБ, где приведено конкретное выражение, связывающее зависимость диэлектрической проницаемости среды через концентрацию носителей зарядов в ЭПЦ с магнитудой 4 2(3)22 где ,и- заряд, концентрация, масса носителей зарядов- частота столкновения носителей зарядов, равная 2,4107 / (в 1/сек)- давление в мм рт. ст.- абсолютная температура. 16281 1 2012.08.30 После подстановки конкретных значений, входящих в приведенное выражение величин, получаем простое выражение 8, 10, связывающее частоту с концентрацией носителей зарядов, зависящую от магнитуды подготавливаемого ЗМТ(4)0,506103 . На основании изложенного очевиден алгоритм определения магнитуды подготавливаемого ЗМТ расчеты РЛС, выполняя свои штатные задачи, производят обзор сектора наблюдения на несущей (рабочей) частоте, которая термостатирована и не может изменяться в условиях нормальной (невозбужденной) обстановки в атмосфере процесс подготовки ЗМТ приводит к изменению в пространстве ЭПЦ параметров атмосферы, формируя здесь ионизированную среду, концентрация зарядов в которой может изменяться от 1010 до 1016 5-9 в зависимости от магнитуды подготавливаемого ЗМТ, что однозначно удовлетворяет условиям отражения излученного РЛС ЭМК, принимаемого с приращением частоты запуская ранее рассмотренный механизм набора статистических данных по приращению частоты в зависимости от магнитуды локальных или проходящих ЗМТ, а также афтершоков их (при необходимости экстраполируя), определяют пороговое значение частоты для возможных опасных магнитуд и записывают их в ячейку 15 устройства (фиг. 9) производя наблюдения за сейсмообстановкой, отмечают приращения текущего значения частоты, сравнивают его с опорным значением и принимают решения об опасности. При необходимости текущие значения приращения частоты можно контролировать, например, с помощью частотомера 3-34. Результаты наблюдений передают на вышестоящий пункт приема решений. Кроме того, магнитуду подготавливаемого ЗМТ с помощью РЛС можно определять путем определения площади стабилизировавшегося ионизированного образования в районе ЭПЦ при определении его азимута. Ионосферный способ прогноза ЗМТ Авторами показано 8, 9, что на поверхности Земли в процессе подготовки ЗМТ в эпицентральном районе (зоне) формируются две области пульсирующих зарядов противоположного знака, что с учетом процессов в ИОН, генетически связанных с этими процессами, образуется эквивалентный колебательный контур, излучающий электромагнитные волны на собственных частотах, которые регистрируются РЛС наземного базирования в виде активных помех (фиг. 3, 4, 6-9). Показано, что в процессе широкомасштабной подготовки ЗМТ образующиеся пучки носителей зарядов, достигая ИОН, на длине свободного пробега в ней развивают энергию более 750 эВ, что вызывает ионизацию молекул ионами и интенсивное свечение спектральных линий кислородных. Эксперименты авторами были выполнены с помощью разработанной оптико-электронной аппаратуры, установленной на орбитальной станции Салют-4. В ходе экспериментов было установлено также, что изменяющаяся концентрация зарядов в ИОН в процессе подготовки ЗМТ приводит к нарушению распространения электромагнитных волн и потере информации, поступающей на земные объекты с космических аппаратов, орбита которых располагается выше -слоя ИОН. Результаты исследований с уточнением фактов ЗМТ были уточнены по каталогу таковых в Институте физики Земли АН СССР, г. Москва. Уточнение координат (широта 59 градусов, долгота 57,6 градуса) позволило постфактум определиться, что сейсмическое событие произошло в районе г. Кизил Пермской области 28.08.1874 г. Исследования в последующем были продолжены с преломлением на прогноз и по признакам искажения (иногда и потере) информации. Это следует из того, что при излучении ЭМК с космического аппарата (КА), находящегося выше -слоя и движущегося по известной орбите, принимают это колебание на земной поверхности аппаратурой с известными координатами, определяют измеренное 16281 1 2012.08.30 значение дальности между источником излучения ЭМК и точкой приема его путем измерения времени прохождения колебаний между ними, сравнивают его с известным значением дальности по разности полученных значений определяют величину ошибки измерения дальности и при превышении этой величиной эталонного значения судят о подготовке ЗМТ. На основе изложенного в работах 8, 10 - под руководством И.С. Ледовского было показано, что над участком подготавливаемого ЗМТ в ИОН возникает пространство с существенно повышенной концентрацией носителей зарядов по отношению с соседними участками, более чем в 2,5 раза, что изменяет время прохождения сигнала через этот участок, а следовательно, измеренная разность расстояний представляет собой ошибку, так называемую дифференциальную погрешность, величина которой в значительной степени зависит от состояния ИОН и может быть использована в интересах прогноза ЗМТ. Алгоритм работы достаточно прост. Допустим, что с источника, находящегося на КА, входящего, например, в состав орбитальной группировки типа ГЛОНАСС, излучается радионавигационный сигнал, трасса распространения которого проходит через обычную среду ИОН без процесса подготовки ЗМТ. Этот сигнал принимается станцией выработки дифференциальных поправок, координаты точки приема и местоположение КА в любой момент времени известны. На земной поверхности определяется время распространения сигнала между объектами и разность расстояния между ними. Это и является одной из основных функций станции выработки дифференциальных поправок, которая используется при реализации настоящего способа. Если трасса распространения сигнала проходит через возмущенную область ИОН, то эта разность превышает пороговую величину, задаваемую в блоке сравнения на Земле, индицируется на индикаторе, сигнализируя о начале подготовки ЗМТ. Результат передается на ЦПСОИ, где производится привязка к единым координатам и единому времени и выработка соответствующих решений. Для определения координат предполагаемого ЭПЦ достаточно опустить перпендикуляр от возбужденной зоны на земную поверхность, например, с 100 км. Поскольку процесс подготовки ЗМТ длится достаточно долго, то операция прогноза многократно уточняется путем облета КА возбужденной зоны района ЭПЦ с уточнением контура ее. Разработанный авторами оптико-электронный прибор ДИСК-1 (дистанционный измеритель свечения кислорода) может быть адаптирован к ГЛОНАСС. Способ прогноза позволяет определить начало процесса подготовки ЗМТ определять ЭПЦ подготавливаемого ЗМТ прогнозировать примерную магнитуду подготавливаемого ЗМТ по размеру возбужденной зоны как в ИОН, так и на Земле. Глубина прогноза подготавливаемого ЗМТ составляет от 20 до 30 суток. Инфранизкочастотный способ прогноза ЗМТ Авторы в своих работах неоднократно указывали, что процесс подготовки ЗМТ следует рассматривать как термодинамическую флуктуацию в некоторой ограниченной области объема земных недр 5, 7, 8, 10. При этом установлено, что эта термодинамическая флуктуация всегда сопровождается изменением энергетических процессов в возбужденной части объема с последующим восстановлением равновесного значения их либо активизацией их. Поскольку установлено, что в земных недрах протекают различного вида токи 5,7, 8, 19, 20 и др., то согласно постулату Ван дер Зила А 48 Там где протекает электрический ток, там существуют его флуктуации. В соответствии с этим, рассматривая участок ГП как среду, характеризующуюся некоторым удельным сопротивлением, и развивая наработанные результаты для микросистем с целью определения их надежности и продолжительности жизни 45, 46, 47, 50-54, авторы распространили свои подходы к процессам подготовки ЗМТ как адекватным, но проявляющимися в макросистемах. 16281 1 2012.08.30 При этом подходы авторов базируются на предложенной ими резистивно-емкостной и дислокационной 43, 45-47, 50, 52-54 шумовой модели кристаллических пород ЗК, связанных с процессами образования и эволюции заряженных дислокаций при трещинообразовании 43 и миграции флюидов (подобно показанному в 53) в ходе подготовки и развития сейсмологического события. В рамках этой концепции авторы теоретически и практически 43, 52-54 показали возможность использования нового устойчивого прогнозного признака ЗМТ - инфранизкочастотного шума (ИНЧШ). Такой предвестник оказывается существенно более эффективным, чем известный признак вариаций удельной проводимости среды 49, и позволяет достоверно прогнозировать начало сейсмической активности, местоположение очага и время ожидаемого события, что подтверждено 17 авторскими свидетельствами и патентами. Согласно классической теории шумовых процессов 44, 48 установлено, что в резистивно-емкостных и полупроводниковых структурах имеют место тепловые, дробовые, импульсные, мерцательные шумы. Особое место в ряду шумов занимают так называемые избыточные (токовые), или фликкер-шумы. Установлено, что мощность тепловых, дробовых шумов в широком диапазоне частот распределяется равномерно 44 и поэтому в области инфранизких частот имеет малое значение и информации не несет. Причинами фликкершумов (назовем их инфранизкочастотными шумами - ИНЧШ) являются дефекты в структурах 43-54 и различные нестационарные процессы, протекающие, например, в ГП. Этот вид шумов всегда связан с протеканием тока, поэтому его иногда называют токовым 44, 48. Дисперсия э.д.с. этого вида шума определяется согласно выражению шф/гдеи- постоянные спектра,- текущее значение частоты. На практике производят измерение мощности ИНЧШ, поэтому для расчета их количественных значений удобнее пользоваться выражением 44, 48 вида шф/ ,-13-16 где 510 210 и 0,82,2. Очевидно, что мощность ИНЧШ обратно пропорциональна частоте и возрастает с ее уменьшением. На практике диапазон наблюдаемых ИНЧШ определен диапазоном частот от 10-6 Гц до 106 Гц 44. Впервые авторами применение параметров ИНЧШ в интересах реального прогноза ЗМТ было предложено в 1984 году 9, что явилось развитием и продолжением работ 43,45-47 и более 75 авторских свидетельств, и адекватности физических процессов, происходящих как в элементах и системах промышленных устройств, так и в ЗК при их деградации и разрушении. Понимание процессов, происходящих в названных элементах и ЗК, следует из подобия представления их, например, в виде композиционного резистора, т.е. резистора, состоящего из материалов различной проводимости с наличием разделительных слоев(пленок) как твердого, так и жидкого состояния, например, для ЗК 35. Структура такого резистора представима в виде матрицы параллельно-последовательно соединенных контактирующих зерен 47, 50, 51 и характеризующихся длиной 1, ширинойи толщиной . При этом количество зерновых слоев равно 1, ,соответственно. С целью упрощения понимания и анализа, положим ток протекает вдоль зерен в направлении длины образца частицы зерен являются одинаковыми с контактной площадьюфлуктуации на контактах являются некоррелированными. Сопротивление одного слоя в этом случае в образце определяется согласно простому выражению/ . Среднеквадратическое значение флуктуации (шумов) этого слоя вычисляется/( /) / (1 /) . 16281 1 2012.08.30 Переходя в приведенных выражениях для зерен/ - число контактов ( - длина участка) 1, для спектральной плотности мощности ИНЧШ можно записать( )2 / 104. Допуская, что в ГП возникли нарушения сплошности, приводящие к изменению величини(при 1), спектральная плотность мощности (СПМ) существенно изменяется 9, 43 и может быть определена в виде выражения( )10002 / 10 4. Очевидно, что при возникновении дефектов в структуре участков ГП, например, при трещинообразовании в моменты подготовки ЗМТ уровень спектральной плотности мощности ИНЧШ возрастает в 1000 раз 43, что говорит о значительно большей чувствительности способа по отношению к способу предложенному О. Барсуковым 49. Это первое. Второе. Огромный вклад в развитие процесса подготовки ЗМТ вносят заряженные дислокации (об этом авторы говорили), под которыми применительно к процессам подготовки ЗМТ следует понимать дефекты проводящей среды любого происхождения и размера, окруженные энергетическим барьером или содержащие потенциальные ямы 15, 47, 51. При этом в процессе деформации и трещинообразования заряженные дислокации перемещаются, размножаются и объединяются, обуславливая дальнейшие деформации, которые по мере роста могут привести (и приводят) к лавинообразным процессам разрушения. Влияние дислокаций на структуру, качество и другие характеристики твердотельных устройств и систем резистивного и полупроводникового типа обсуждены и показаны в работах 15, 43-47, 50-54 и др., что позволило установить закономерности, связанные с перемещением зарядов, скоростью их перемещения, концентрацией носителей зарядов, уровнем ИНЧШ, их спектром. Показано, что, поскольку ИНЧШ возникает и развивается в ограниченных областяхдислокациях, движения носителей зарядов связаны с ними и порождаются ими и это справедливо как для микросистем, так и для макросистем. Поскольку флуктуации заряда дислокации являются источником флуктуирующего токового диполя, то они порождают и флуктуации напряжения на любом наблюдаемом участке земной коры при подготовке ЗМТ. С учетом кратко изложенного можно записать выражение для СПМ, например, по напряжению( ) / (/) ,где- выбираемая средняя частота наблюдения. Практические наблюдения 43 показали, что наличие любого дефекта в образце ГП при деформациях происходят размножения деформаций и резкое возрастание ИНЧШ в широком интервале частот, что хорошо оценивается математически через эмпирическое выражение для СПМ по напряжению( )2 /21 /. Полагая здесь 6 (практически 5-20), это выражение можно записать в виде 4,910-13(1/), В/Гц. Допуская, что два образца ГП с СПМ по напряжению 14,910-11(1/), В/Гц (бездефектный образец) и 21,5610-10(1/), В/Гц (образец с дефектом) дают сравнение,оцениваемое величиной в 320 раз. Таким образом, при развитии активности дислакационного процесса в период подготовки ЗМТ и развитии трещин в ЗК существенно повышается общий уровень ИНЧШ и СПМ их, что расширяет возможности прогноза параметров ЗМТ подтверждено 19 авторскими свидетельствами и патентами на изобретения. Предсказание времени наступления ЗМТ Способ предполагает простую аппаратурную реализацию, финансово малозатратную и надежную реализацию. 16281 1 2012.08.30 В целях повышения достоверности принятия решения при прогнозировании соответствующих параметров ЗМТ способ предполагает возможности накопления статистики по маломощным ЗМТ с целью уточнения закона изменения СПМ. Согласно указанному ранее,процесс накопления статистики является достаточно простым, особенно в сейсмоактивных регионах. Схема практической реализации достаточно проста и содержит керамические неполяризующиеся электроды с медными вставками с заполнением медным купоросом разработки ВИРГ (известны в применении геологами и геофизиками), устанавливаемые в выбранном участке земной поверхности прецезионный регулируемый проволочный резистор,необходимый для согласования сопротивления электродов отведения со входным сопротивлением измерителя малошумящий усилитель на полевых транзисторах с включенными в него полосовыми декадными фильтрами усилитель с регулируемым коэффициентом усиления квадратичный детектор на основе микросхемы 140 МА 1 Б индикатор. Устройство выполнено авторами 57, 58. Измерительное устройство выполнено компактно и может быть использовано как совместно со стационарными сейсмостанциями 18, так и в автономном режиме выездными бригадами наблюдения за сейсмоактивностью. С целью понимания предлагаемого, приведем пример. 1. С помощью источника э.д.с. (аккумулятор емкостью 190 Ампер-часов и напряжением 12 Вольт) в межэлектродном пространстве ЗК возбуждается электрический ток с наложением на него флуктуации ИНЧШ исследуемого участка. 2. Изменением величины сопротивления проволочного резистора добиваются максимума показаний мощности на индикаторе с учетом устанавливаемых частот. Допустим,что на выбранных частотах отношение мощностей достигло величины фо 5,87. 3. Далее проводят наблюдение за возрастанием величины фо в течение некоторого промежутка времени, например 250 часов (выбирается произвольно). Естественно,что выдержка такого времени приводит к разряду аккумулятора, искажению результатов эксперимента, поэтому краткосрочный разряд и изменение уровня и спектра ИНЧШ имитировалось использованием шунтов. Допускаем, при вторичном измерении отношение мощностей возросло до значения фо 241,62. Априори также допускаем, что для наблюдаемого сейсмически активного района пороговое отношение мощностей ИНЧШ равно фоп 163,364, которое определелено согласно выражению фопфо 2 фо ост 2/,где фо/(фо 2-фо)/(41,62-5,87)/2500,143 часа. Время, оставшееся до момента наступления ЗМТ, определяется из выражения ост 2(фоп-фо 2)/(фо/)(163,364-41,62)/0,143851,0 час. Настоящий алгоритм прогноза как на этапе подготовки реализации, так и при практическом прогнозе включает как расчетные операции, так и операции по набору статистики. При этом этапы его выполнения таковы 1. Согласно выражению для спектральной плотности ИНЧШ в параметрах мощности производят расчет таковой для различных значений постоянных, т.е.-16 где 510 2100,82,2- частота. Для различных значенийистроят графики (создают таблицы) зависимостей мощности в координатах мощность - время и вносят в память вычислителя с указанием пороговых значений для каждого сейсмоактивного региона. Постояннаяхарактеризует свойства ГП каждого из регионов, она определяется и задается геологами или геофизиками. В частности, для Беларуси она задается НПО Белгеология или Институтом геологии и геофизики АН РБ 18 для Казахстана, например, Институтом сейсмологии АН КА и др. Показательв общем случае задается с интервалом в 0,2 частота задается в декадах 0,1-1,0 Гц 1,0-10,0 Гц 10,0-100,0 Гц через интервалы, не допускающие больших погрешностей с учетом свойства инвариантности мощности 44. 16281 1 2012.08.30 2. В процессе повседневных наблюдений за сейсмообстановкой в регионах с использованием ранее приведенного устройства, устанавливаемого на существующих сейсмостанциях, производят набор статистики. Строят графики в координатах мощность - время для различных магнитуд с указанием пороговых значений - фиг. 10. 3. Производят сравнение расчетных и наблюдаемых результатов, отдавая предпочтение наблюдаемым результатам. 4. В зависимости от наступающей сейсмической активности с выбранным временным интервалом производят измерения мощности шумов ЗК, например, 1(1) и (2) на фиксированных частотах 1 и 2 (в одной из декад) соответственно, через выбранный интервал времени 2 повторно измеряют указанные мощности и находят их новые значения 2(1) и 2(2) и затем определяют второе значение ост 2 искомого времени в соответствии с выражением(6) ост 2 фп-ф/ф,где фпфп(1)/фп(2) ф 2(1)/2(2) /фф-1(1)/Р(2)/,где фп(1) и фп(2) - заранее определенные пороговые значения мощностей ИНЧШ для данного региона. 5. При отработанной программе определения времени, оставшегося до момента наступления ЗМТ, это возможно производить и с помощью графического представления мощность - время, построенного ранее на основе расчетов по п. 1 и сравнения с набранной статистикой. В этом случае входят в график мощность - время (фиг. 10) и, например, по кривой 2 для заданного порогового значения мощности определяют соответствующее время. На графике обозначено фп и фи (и) - пороговое значение мощности ИНЧШ и мощность их на момент измерения и. 6. Если возникает необходимость уточнения значения времени, оставшегося до момента наступления ЗМТ, определенных от нескольких источников прогноза, то это возможно, например, сравнения такового через ост 1 и 2. Прогнозирование положения эпицентра ЗМТ В основу способа прогноза ЭПЦ ЗМТ впервые положено измерение флуктуаций горизонтальной и вертикальной составляющих напряженности ЭП Земли, которые однозначно обусловлены и связаны с процессами подготовки ЗМТ. Это следует из того, что, согласно изложенному, землетрясению предшествует накопление деформаций, объединение и перемещение заряженных дислокаций, возникновение и размножение зарядо-разрядных процессов и др. Названное в совокупности приводит к возникновению и распространению продольных и поперечных упругих (сейсмических) волн, которые деформируют не только глубинные, но и приповерхностные слои ЗК, приводя к вариациям их проводимости, группированию и флуктуациям различных параметров в районе очага подготавливаемого ЗМТ и в том числе к проявлению флуктуаций напряженности электрического поля Земли. При этом очевидно, что за флуктуации горизонтальной составляющей напряженности ЭП Земли (ЭПЗ) ответственность несут продольные сейсмические волны, а за флуктуации вертикальной составляющей - поперечные волны. Понятно, что определяющий вес в ЭПЦ подготавливаемого ЗМТ принадлежит флуктуациям вертикальной составляющей напряженности ЭПЗ, амплитуда которых по мере удаления от ЭПЦ уменьшается. Таким образом, очевидно, что процесс реализации способа предполагает измерение флуктуаций вертикальной и горизонтальной составляющих напряженности ЭПЗ, сравнение их и суждение по этим измерениям и сравнениям о положении ЭПЦ подготавливаемого ЗМТ. Впервые предложено авторами 5-10. Схема практической реализации способа такова, что при реализации операции определения времени, оставшегося до момента наступления ЗМТ, предложенной авторами (22 описания), где вместо неполяризующихся электродов отведения к усилителю подключается преобразователь (датчик) флуктуации напряженности ЭПЗ - флюксометр, конструкция которого приведена в 21. Датчики флуктуаций как вертикальной, так и 23 16281 1 2012.08.30 горизонтальной составляющих ЭПЗ конструктивно одинаковы, отличаясь лишь расположением приемных пластин на 90 градусов. Диаграмма направленности этих датчиков сформирована по закону квадрата косинуса и имеет относительно выраженный максимум. Алгоритм реализации способа 1. Эпицентр ЗМТ при реализации настоящего способа является по сути уточняющим и проводится выездной бригадой после предварительного определения его, например, радиолокационным способом. В этом случае в начальной точке операции производят измерение мощности флуктуаций вертикальной и горизонтальной составляющих ЭПЗ, определяют величину их отношения, а также азимутальный угол 1, соответствующий максимуму мощности флуктуаций горизонтальной составляющей напряженности ЭПЗ в этой точке. Азимут определяют как угол между плоскостью местного меридиана и осью измерителя, совпадающей с максимумом диаграммы направленности входного преобразователя. При измерении мощности флуктуаций горизонтальной составляющей напряженности ЭПЗ второй координатный угол - полярный угол между местной вертикалью и осью измерительного устройства - равен 90 градусам, полярный угол при измерении вертикальной составляющей равен нулю. При реализации способа могут использоваться два или один канал флюксометров в последнем случае выполняется переключение их, что увеличивает время переходного процесса при измерениях, но снижает стоимость реализации и требует разворота датчика на 90 градусов при измерении горизонтальной составляющей напряженности ЭПЗ. 2. Далее изменяют азимут (направление) оси измерителя, проводя подобные измерения и отсчет через каждые 510 градусов до тех пор, пока не будет достигаться максимум мощности соответствующих флуктуаций ЭПЗ. 3. Регистрируют азимут 1 и рассчитывают отношение мощностей флуктуации ИНЧШ вертикальной и горизонтальной составляющих ЭПЗ, допуская, что в точке измеренное значение угла 130 градусов (при отсчете угла по часовой стрелке), а отношение в/г 0,12 перемещают аппаратуру из точки 1 по направлению измеренного азимута 1. Шаг измерений зависит от условий местности и на больших расстояниях от ЭПЦ предполагаемого ЗМТ может составлять километры. 4. Далее, действуя аналогично, производят измерения в точке 2, где азимутальный угол 212,5 градусов, а отношение в/г 0,35, что является показателем приближения к предполагаемому ЭПЦ подготавливаемого ЗМТ. Если направление перемещения выбрано ошибочно (например, противоположно), то отношение мощностей флуктуаций будет иметь меньшее значение. Ошибка может возникнуть по причине формы диаграммы направленности флюксометра - гантельная, косинусквадратная. В этом случае возвращаются в точку 1 и повторяют измерения при перемещении по азимуту 1 в направлении точки 2. 5. Перемещают аппаратуру в точки 3, 4 и т.д., последовательно измеряя отношение мощностей флуктуаций вертикальной и горизонтальной составляющих напряженности ЭПЗ. Показателем близости к ЭПЦ подготавливаемого ЗМТ является размытость или отсутствие максимума указанного отношения мощностей флуктуаций. 6. С целью уточнения положения ЭПЦ после операций, например, в точке 6 производят еще один шаг к точке 7 (последней) по азимуту точки 5 и соответствующие измерения. Дальнейшее уточнение положения ЭПЦ и его предполагаемого района можно провести путем обследования с шагом в 0,51,0 км, производя измерения по изложенной методике. За ЭПЦ подготавливаемого ЗМТ принимают некоторый центр области, внутри которой отношение мощностей флуктуаций составляющих напряженности ЭПЗ остается примерно постоянным. 16281 1 2012.08.30 Поскольку процесс подготовки ЗМТ длится достаточно долго, то координаты ЭПЦ могут уточняться путем повторных измерений расчетом выездной бригады, сформированной, например, на какой-то сейсмообсерватории и укомплектованной аппаратурой,размещенной на автомобиле типа Нива. Прогнозирование магнитуды ЗМТ Сущность способа вытекает из того, что флуктуации рассмотренных составляющих напряженности ЭПЗ в ЭПЦ подготавливаемого ЗМТ и в некоторой опорной точке (некоторых точках) порождены одной и той же причиной - распространением сейсмических волн (СМВ), источником которых является район очага подготавливаемого ЗМТ. Флуктуации составляющих напряженности ЭП в указанных точках, естественно, оказываются коррелированными, и их связь количественно может быть охарактеризована через автокорреляционную функцию (АКФ), имеющую максимум в некоторой точке. Применительно к предлагаемому способу прогноза значение АКФ в максимуме определяется расположением датчиков флуктуаций составляющих напряженности ЭПЗ, и конкретно для ЭПЦ оно будет равно единице, поскольку регистрируемые флуктуации напряженности ЭП в ЭПЦ порождены одним и тем же типом СМВ, например продольными. Следовательно, при реализации способа для достижения максимума в некоторой опорной точке датчик флуктуаций необходимо разворачивать до совпадения максимума его диаграммы направленности с направлением распространения, например, продольной СМВ, приходящей в опорную точку, т.е. плоскость раскрыва датчика (его апертура) должна совпадать с плоскостью распространения СМВ, т.е. они должны быть параллельны. При отсутствии реальных методов определения сил, действующих в районе подготавливаемого ЗМТ 12-17 до момента его свершения, используются косвенные методики определения магнитуды триангуляционный метод засечек 12, 13, методики наложения карт ожидаемых землетрясений (КОЗ) 15-17, а также методики КОЗ и алгоритмов ,реализуемых после свершения события. Впервые предлагается практически работающий способ прогноза магнитуды подготавливаемого ЗМТ до момента его свершения. При этом в качестве параметра, характеризующего магнитуду подготавливаемого ЗМТ, предложен (и используется) радиус взаимной корреляции, т.е. расстояние, на котором значение ВКФ в максимуме убывает в е раз по сравнению со значением ВКФ в ЭПЦ ЗМТ. Это связано с тем, что амплитуда СМВ,распространяющихся из очага подготавливаемого ЗМТ, уменьшается, уменьшается и значение ВКФ в опорной точке (на это указывалось выше). Очевидно, что радиус ВКФ практически не зависит от глубины залегания очага подготавливаемого ЗМТ, поскольку расстояние относа точек измерения и глубина залегания очага - величины не соизмеримые. Структурная схема реализации представлена на фиг. 11, где обозначено 1 и 5 - датчики флуктуаций ЭП 2 и 6 - широкополосные усилители 3 и 7 - анализаторы спектра (декадные фильтры) 4 - двухканальный самописец 8 - коррелометр с квадратичным детектором. Кроме указанных устройств, используются специальные высокостабильные серебряно-цинковые аккумуляторы и малогабаритные средства связи между исследователями в ЭПЦ и в опорных точках. Комплект аппаратуры с исследователями (выездная бригада) размещаются на автомобилях типа Нива, ГАЗ-69 и др. Работы проводятся после установления местоположения ЭПЦ подготавливаемого ЗМТ и ведутся с повтором при необходимости. Алгоритм реализации 1. Устройства измерения флуктуаций (с флюксометрами в качестве датчиков ЭП 51) составляющих напряженности ЭПЗ устанавливают в ЭПЦ подготавливаемого ЗМТ и в одной из опорных точек на расстоянии в 500-700 м от ЭПЦ. 2. Производят измерения вертикальной составляющей напряженности ЭПЗ в обеих точках с последующим преобразованием в электрическое колебание, усилением широко 25 16281 1 2012.08.30 полосными усилителями в каналах до одинакового уровня (амплитуды) путем регулирования коэффициента усиления одного из усилителей. 3. Усиленные и выравненные по амплитуде колебания через анализаторы спектра подают на коррелометр, по индикатору которого путем изменения времени задержки колебания (сигнала), пришедшего от ЭПЦ, регистрируют зависимость ВКФ (выраженную через значения напряжения) от времени задержки. Определяют время задержки, при котором достигается максимум ВКФ, фиксируя его. 4. Производят разворот флюксометра в опорной точке относительно горизонтальной оси и по индикатору определяют максимум ВКФ. Подобные измерения производят последовательно в выбранных опорных точках, каждый раз определяя зависимость максимума ВКФ флуктуаций вертикальной составляющей напряженности ЭПЗ от местоположения опорной точки до ЭПЦ подготавливаемого ЗМТ. 5. Перемещение измерительной аппаратуры производят до тех пор, пока не будет определена точка, где значение ВКФ не достигнет значения 0,37 от ее максимального значения в ЭПЦ подготавливаемого ЗМТ. 6. Последовательно повторяют измерения для каждой опорной точки и определяют зависимость максимального значения ВКФ флуктуаций ЭПЗ от расстояния опорной точки до предполагаемого ЭПЦ, определяя величину радиуса взаимной корреляции к. 7. Магнитуду подготавливаемого ЗМТ в прогнозируемом ЭПЦ вычисляют в соответствии с выражением МзМэ(к/кэ),где э и кэ - соответственно эталонная магнитуда и эталонный радиус взаимной корреляции для наблюдаемого района в спокойном состоянии. Таким образом, инфранизкочастотный способ по измеряемым и вычисляемым параметрам и характеристикам ИНЧШ, развивающихся в ЗК при подготовке ЗМТ, позволяет с минимальными затратами прогнозировать одновременно все три параметра ЗМТ место,время, сила. Минимум затрат на реализацию способа объясняется как простотой конструктивного решения измерительной аппаратуры, так и тем, что она располагается на существующих сейсмостанциях и для передачи информации используются те же каналы передачи. Кроме того, аппаратура измерения может выполняться в выездном варианте. Следует иметь в виду, что каждый из способов может быть реализован самостоятельно (автономно), но вероятность прогноза параметров подготавливаемого ЗМТ каждым из способов различна, различно и время обнаружения начала процесса подготовки 1. Наиболее чувствительным является способ с применением РЛС. Это объясняется, вопервых, тем, что чувствительность их по мощности принимаемого ЭМК превышает 10-14-10-15 Вт/см 2, что позволяет вести процесс обнаружения начала подготовки для ЗМТ в 3,54,0 баллов за 10-12 суток, для ЗМТ от 8,0 и более - за 85-90 суток до свершения события вовторых, вероятность правильного прогноза параметров ЗМТ составляет не менее 0,75-0,8. 2. Ионосферный и наземно-космический способы начинают уверенно работать при магнитудах подготавливаемого ЗМТ с 4,0-4,5 баллов с вероятностью в 0,3-0,4 и по срокам в 10-12 суток. 3. Прогноз по ИНЧШ по параметрам обнаружения приближается к способу прогноза с помощью РЛС, но уступает ему по вероятности правильного прогноза. Очевидно, что совокупное применение, т.е. комплексное и многопараметрическое прогнозирование с одновременным использованием предлагаемых способов, позволяет вести уверенный прогноз всех трех параметров подготавливаемого ЗМТ с вероятностью правильного прогноза 0,90-0,99 и сроком обнаружения начала подготовки до 90 суток. Законченность решения краткосрочного многопараметрического прогноза ЗМТ обуславливается сведением разработанных способов прогноза в единую систему прогноза,включающую наземные стационарные пункты прогноза, каждый из которых связан с РЛС для определения всех трех параметров подготавливаемого ЗМТ спутниковые пункты про 26 16281 1 2012.08.30 гноза с оптико-электронной аппаратурой прогноза по параметрам свечения кислородных линий в ИОН над районом ЭПЦ, а также аппаратуру излучения ЭМК на наземные пункты через возбужденную ионосферу как стационарные, так и мобильные пункты прогноза с аппаратурой прогноза эпицентра, магнитуды и времени наступления подготавливаемого ЗМТ по параметрам и характеристикам ИНЧШ в ЗК. Информация с названных пунктов прогноза по своим каналам связи передается на центральный пункт сбора и обработки информации (ЦПСОИ), где после идентификации всех поступивших сведений, приведения к единому времени и единым координатам принимаются решения на дальнейшие действия либо по месту, либо по передаче результатов в вышестоящие инстанции, например на центральный пункт МЧС. Алгоритм реализации изложенного поясняется 1. Допустим, что воздушное пространство Беларуси контролируют, например, 5 (пять) РЛС кругового или секторного наблюдения. При обнаружении всеми или некоторыми из них начала процесса подготовки ЗМТ в вверенном районе ответственности их операторы по изложенным методикам определяют параметры ЗМТ (предварительно уточнив, что наблюдается именно процесс подготовки ЗМТ) и по существующим каналам связи передают определенные параметры на вышестоящий командный пункт (КП), на который также поступают данные с других РЛС. На КП производится сравнение поступивших данных, приведение их к единым координатам, единому времени (при необходимости). Каждому параметру присваивается код (шифр), а также шифруется и код соответствующего КП. Зашифрованная информация передается на ЦПСОИ, где она дешифруется и по ней принимается решение на дальнейшие действия. 2. На территории Беларуси в настоящее время имеет место шесть сейсмических станций 18 Нарочь, Солигорск, Гомель, Брест, Могилев и Минск, а также несколько пунктов сейсмического наблюдения. На указанных сейсмообъектах согласно реализации способа прогноза располагается аппаратура прогноза по параметрам ИНЧШ, дубляж которой выполнен также и в мобильном (подвижном) варианте. С помощью этой аппаратуры (устройств) параллельно со стационарной аппаратурой сейсмостанций ведутся наблюдения за сейсмической обстановкой. Получаемые результаты начала и процесса протекания подготовки ЗМТ по существующим каналам связи системы сейсмонаблюдения передаются на центральный пункт, находящийся в Минске в Институте геологии и геофизики АН РБ. Здесь производится обработка поступающей информации как по времени, так и по координатам, обобщается, шифруется и передается на ЦПСОИ. 3. Ионосферный и наземно-космический способы начинают уверенное обнаружение процесса подготовки ЗМТ при магнитудах в 4,0-4,5 баллов. Информация, получаемая оптико-электронной аппаратурой о начале подготовки ЗМТ по существующим каналам связи либо непосредственно, либо через пункт космической связи, находящийся в НИИ информатики АН РБ, передается на ЦПСОИ. Из изложенного следует, что достаточно простые, относительно малозатратные и надежные способы прогноза, совмещаемые с существующими как комплексами наблюдения, так и со средствами связи, объединяемые в единую систему наблюдения за развитием сейсмических событий, дают высоковероятностный прогноз всех трех параметров ЗМТ место, время и сила - за 10-90 суток до свершения ЗМТ с вероятностью правильного прогноза 0,85-0,99. В относительно наглядном представлении структурная схема способа и система в развернутом виде представлена на фиг. 12, где (по аналогии с фиг. 1) обозначено 1 - наземные пункты прогноза по параметрам и характеристикам ИНЧШ как стационарного типа, совмещенные с существующими сейсмостанциями, так и мобильные в виде выездных бригад 2 - радиолокационные пункты наблюдения, реализующие прогноз при совмещении функции прогноза с действующими РЛС, например, ПВО 16281 1 2012.08.30 3 и 4 - космические аппараты как специального ионосферного и наземно-космического прогноза, так сопутствующие этому прогнозу, т.е. выполняющие какие-то штатные функции, например, в системе ГЛОНАСС 5 - центральный пункт сбора и обработки информации (ЦПСОИ), выполняющий операции сбора и обработки информации с приведением ее к единым координатам, единому времени и принятия решений на действия по подготавливаемому ЗМТ и др. 6 и 7 - линии связи и пункт ретрансляции при необходимости увеличения расстояния разноса между пунктами прогноза. Примерный алгоритм прогноза процесса параметров и характеристик подготавливаемого ЗМТ реализуется так 1. Радиолокационные станции ПВО (например), выполняя свои штатные задачи, согласно возложенным на них дополнительным функциям, обнаруживают, например, изменение координат местника (телевышка, вышкаи др.) или какого-либо отражающего объекта, отличающегося по параметрам от стандартной цели 2. На РЛС в штатном режиме производится так называемая первичная обработка обнаруженного объекта классификация как стационарного (неподвижного в пространстве) или движущегося, отнесение по известным признакам к цели самолет, например, или ионизированное образование, характеризующее признак начала подготовки ЗМТ, например,по объему засветки экрана ВИКО. Принимается решение о типе объекта начинается работа по нему, получаемая информация через центральный пункт группировки ПВО в выбранном коде (шифре) передается на ЦПСОИ. При поступлении информации от нескольких РЛС, например минуя КП, производится ее обработка с целью отнесения ее к ранее поступившей информации, т.е. к определению идентичности или отнесению к другому объекту и к неотнесению ее к прогнозу. В случае обнаружения начала процесса подготовки ЗМТ посредством информационного обмена между объектами прогноза начинается их активизация по прогнозу. Поскольку чувствительность каждого из способов прогноза различна во-первых, а, во-вторых,способы прогноза, работающие по прогнозу из космоса, могут включаться в активную работу только тогда, когда орбита объекта прогноза проходит над эпицентральным районом подготавливаемого ЗМТ, то существует и некоторое запаздывание в передаче информации на ЦПСОИ. Однако это не является недостатком разработанного способа, поскольку процесс подготовки ЗМТ длится десятки суток и обновление информации будет происходить многократно, что позволяет получать устойчивые совокупные результаты прогноза от всех объектов, участвующих в процессе прогноза, и тем самым повысить вероятность правильного прогноза всей системы до величины, сколь угодно приближающейся к единице. Способ, включающий несколько локальных способов прогноза, является исключительно наглядным, поскольку процесс начала подготовки ЗМТ начинает отображаться на индикаторах кругового обзора на ранних стадиях его зарождения свечение газовых линий кислорода также начинает фиксироваться за десятки суток до свершения события прогноз по ИНЧШ по времени адекватен прогнозу с помощью РЛС и наглядно отображается как по возрастанию уровня ИНЧШ, так и по контрастному изменению спектральной плотности шумов и др. Ни один из существующих и применяемых в настоящее время способов прогноза не позволяет этого делать. Способ может быть легко совмещен с аппаратурой радионавигационной системы ГЛОНАСС. Применение настоящего изобретения позволит предотвратить многочисленные человеческие жертвы и сохранить огромные материальные ценности. 1. Маренко В.Ф. Исследование связи сейсмотектонических процессов с возмущениями нижней ионосферы методом радиопросвечивания на сверхдлинных волнах / Кандидатская диссертация. - Омск, 1989. - 155 с. 2. А.с. СССР 1382222, МПК 01 3/12. 3. Балов А.В., Жолнеров В.С., Зарубин С.П. и др. Радионавигационные системы дальнего действия с наземным и космическим базированием как система электромагнитного мониторинга краткосрочных предвестников землетрясений // Новости навигации НТЦ. Интернавигация. - 2004. -3. - С. 8-20. 4. Воинов В.В., Гуфельд И.Л., Ледовской И.С. и др. Исследование аномальных электромагнитных излучений в диапазоне СДВ и КВ в условиях сейсмической активности // Известия АН СССР. Физика Земли. - 1990. -4. - С. 93-98. 5. Воинов В.В., Гуфельд И.Л., Ледовской И.С. и др. Эффекты в ионосфере и атмосфере перед Спитакским землетрясением 7 декабря 1998 г. // Известия АН СССР. Физика Земли. - 1992. -3. - С. 96-101. 6. Ледовской И.С. Землетрясение произойдет. Жертв не будет // Газета Советская Белоруссия. - Минск. - 2.09.1989. 7. Ледовской И.С., Омельчук А.П., Шеломенцев Т.И. и др. Комплексная система прогноза землетрясений с использованием радиотехнических методов // Спасение, защита,безопасность в науке, технике, технологии. - М., 1995.- С. 148. 8. Воинов В.В., Демьяненко А.В., Ледовской И.С. и др. Применение сети контрольнокорректирующих станций для сейсмического мониторинга Земли Сб. трудов 111-й Международной конференции. Планирование глобальной радионавигции. - М., 9-11 октября 2000. 9. Ледовской И.С. и др. Способ предсказания времени наступления землетрясения А.с. СССР 1182462, МПК 01 3/00. 10. Воинов В.В., Демьяненко А.В., Ледовской И.С. и др. Способ прогноза землетрясения Патент РБ 6565, 2004. 11. Советский энциклопедический словарь. - М. Советская энциклопедия, 1979. - 1600 с. 12. Рикитаке Т. Предсказание землетрясений. - М., 1979. - 388 с. 13. Эйби Дж.А. Землетрясения / Пер. с англ. - М. Недра, 1982. - 264 с. 14. Болт Б. В глубинах Земли. О чем рассказывают землетрясения / Пер.с англ. - М. Мир, 1984. - 184 с. 15. Курскеев А.К. Проблемы прогноза землетрясений. - Алма-Ата Наука, 2003. - 264 с. 16. Соболев Г.А., Пономарев А.В. Физика землетрясений и предвестники. - М. Наука,2003. - 270 с. 17. Садовский М.П., Болховитинов Л.Г., Писаренко В.Ф. Деформирование геофизической среды и сейсмический процесс. - М. Наука, 1987. - 100 с. 18. Белорусский сейсмологический бюллетень. Вып. 2. - Минск, 1992. - 235 с. 19. Баласанян С.Ю. Динамическая геоэлектрика. - Новосибирск Наука, 1990. - 235 с. 20. Краев А.П. Основы геоэлектрики. - Л. Недра, 1965. - 587 с. 21. Имянитов И.М. Приборы и методы для изучения электричества атмосферы. - М. Гостехиздат, 1957. - 483 с. 22. Юман М. Молнии. - М. Мир, 1972. - 327 с. 23. Киссин И.Г. Землетрясения и подземные воды. - М. Наука, 1982. - 176 с. 24. Рысс Ю.С. Геоэлектрические методы разведки. - Л. Недра, 1983. - 283 с. 25. Духин С.С., Шилов В.Н. Динамические явления и двойной слой в дисперсных системах и полиэлектролитах. - Киев Наукова думка, 1972. - 128 с. 26. Ребиндер П.А., Щукин А.Д. Поверхностные явления в твердых телах в процессе деформации и разрушения // Успехи физических наук. - М., 1978. Т. 105. Вып. 1. - С. 3-42. 27. Бороздич Э.В. // Знание Сила. -6. - 1991. 29 16281 1 2012.08.30 28. Прикладная геофизика / Пер. с англ. Под ред. М.Б.Лейбова. - М. Недра, 1980. 29. Грановский В.А. Электрический ток в газе / Под ред. А.А.Сена и В.Е.Голана. - М. Наука,1871. - 544 с. 30. Степаненко В.Д. Радиолокация в метрологии. - Л. Гидрометеоиздат, 1978. - 343 с. 31. Современная радиолокация (анализ, расчет и проектирование систем) / Пер. с англ. Под ред. Ю.В.Кобзарева. - М. Советское радио, 1968. - 496 с. 32. Радиотехнические системы / Под ред. Ю.М.Казаринова. - М. Советское радио,1068. - 496 с. 33. Электромагнитные предвестники землетрясений. - М. Наука, 1982. - 88 с. 34. Поиск электромагнитных предвестников землетрясений // ИФЗ АН СССР. - 1988. 35. Царв В.П. Особенности нефтегазообразования в зонах тектоно-сейсмической активизации. - Новосибирск. Наука, 1988. - 322. 36. Гефан Г.Д., Иванов В.Б., Хазанов Об одном механизме усиления свечения верхней ионосферы при воздействии мощных радиоволн // Геомагнетизм и аэрономия. - Том ХХ.4. - 1981. - С. 727-729. 37. Григорьев А.И., Гершензон Н.И., Гохберг М.Б. О природе свечения атмосферы при землетрясениях // Докл. АН СССР. - 1988. - Т. 300. -5. - С. 1087-1090. 38. Явление генерации электромагнитных волн ионосферными токами под воздействием на ионосферу модулированного радиоизлучения. Эффект Гетманцева Г.Г. Открытие 231 // Открытия СССР. - 1980. - М. ВНИИПИ, 1981. 39. Зеленова Т.И., Леченька А.Д. Ионосферные эффекты, связанные с монероновским землетрясением // Изв. АН СССР. Физика Земли. - 1989. -10. - С. 6-12. 40. Жуков С.Н., Куксенков В.А., Петров В.А. Физические основы прогнозирования механического разрушения // Докл. АН СССР. - 1981.- Т. 259. -6. - С. 1350-1559. 41. Гуфельд И.Л., Маренко В.Ф., Пономарев Е.А., Ямпольский В.С. Исследование Д-области ионосферы методом наклонного зондирования на сверхдлинных волнах // Поиск электромагнитных предвестников землетрясений. - М. ИФЗ АН СССР, 1988. - С. 150-169. 42. Флуктуации электромагнитного поля Земли в диапазоне СНЧ / Под. ред. М.С.Александрова. - М. Наука, 1972. - 194 с. 43. Ледовской И.С., Воинов В.В. и др. Способ определения предела длительной прочности горных пород А.с. СССР 19408286. 44. Букингем М. Шумы в электронных приборах и системах / Пер. с англ. - М. Мир,1986. - 399 с. 45. Ледовской И.С. Низкочастотные шумы полевых транзисторов Сб. Радиоэлектроника в приборостроении. Вып. 1. - Минск, 1973. - С. 17-21. 46. Ледовской И.С. Оценка качества тонкопленочных резисторов по параметрам инфранизкочастотных шумов Сб. Электронная техника. Серия 8. - Минск. - Вып. 4 (115). 1985. - С. 23-25. 47. Ледовской И.С., Воинов В.В. Влияние дислокаций на спектральную плотность инфранизкочастотных шумов полупроводниковых приборов // НТШ. Шумовые и деградационные процессы в полупроводниковых приборах. - М., 1986. 48. Ван дер Зил А. Шум. (Источники, описание, измерение). Пер. с англ. / Под ред. А.К.Нарышкина. - М. Сов. радио, 1973. - 288 с. 49. Барсуков О.Б. Явление изменения эффективного удельного сопротивления горных массивов перед местными землетрясениями // Открытия в СССР 1979 г. - М.ЦНИИПИ,1980. - С. 29. 50. Воинов В.В., Ледовской И.С. Связь дислокаций с параметрами инфранизкочастотного шума в полупроводниковых и резистивных структурах // НТШ. Шумовые и деградационные процессы в полупроводниковых приборах. - М., 1987.