Мюонный годоскоп УРАГАН тема автореферата и диссертации по физике, 01.04.01 ВАК РФ
Компаниец, Константин Георгиевич
АВТОР
|
||||
кандидата физико-математических наук
УЧЕНАЯ СТЕПЕНЬ
|
||||
Москва
МЕСТО ЗАЩИТЫ
|
||||
2008
ГОД ЗАЩИТЫ
|
|
01.04.01
КОД ВАК РФ
|
||
|
КОМПАНИЕЦ Константин Георгиевич
МЮОННЫЙ ГОДОСКОП УРАГАН
01.04.01 - приборы и методы экспериментальной физики
Автореферат диссертации на соискание ученой степени кандидата физико-математических наук
0 5 Л Р И 2008
003456448
КОМПАНИЕЦ Константин Георгиевич
МЮОННЫЙ ГОДОСКОП УРАГАН
.04.01 - приборы и методы экспериментальной физики
Автореферат диссертации на соискание ученой степени кандидата физико-математических наук
Автор:
Работа выполнена в Московском инженерно-физическом институте (государственном университете).
Научный руководитель: доктор физико-математических наук,
профессор
Петрухин Анатолий Афанасьевич, МИФИ, г. Москва
Официальные оппоненты: доктор физико-математических наук,
зав. лабораторией Кузьмичев Леонид Александрович, НИИЯФ МГУ, г. Москва
доктор физико-математических наук, зам. руководителя отдела электроники
Сенько Владимир Александрович, ГНЦ "ИФВЭ", г. Протвино
Ведущая организация: Институт земного магнетизма,
ионосферы и распространения радиоволн им. Н.В. Пушкова Российской академии наук, г. Троицк
Защита состоится 17 декабря 2008 г. в 17 час. 00 мин. на заседании диссертационного совета Д212.130.07 в МИФИ по адресу: 115409, Москва, Каширское шоссе, дом 31, телефон 324-84-98.
С диссертацией можно ознакомиться в библиотеке МИФИ.
Автореферат разослан 14 ноября 2008 г.
Просим принять участие в работе совета или прислать отзыв в одном экземпляре, заверенный печатью организации.
Ученый секретарь диссертационного совета, /\ ]л / доктор физико-математических наук / С.Е. Улин
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность
Космические лучи (КЛ) несут информацию как о физических явлениях в самих источниках излучения, так и о процессах в среде, через которую они проникают. В частности, поток мюонов, образующихся в результате взаимодействий КЛ с ядрами атомов воздуха на высотах 15-20 км, чувствителен к различным термодинамическим процессам в атмосфере Земли как естественного, так и техногенного происхождения. Таким образом, характеристики потока мюонов, регистрируемого на уровне Земли, являются эффективным индикатором динамики различных процессов в окружающей среде, включающей магнитосферу и атмосферу Земли.
Для разделения влияния атмосферных процессов и процессов в околоземном пространстве используется то обстоятельство, что гелиосферные и магнитосферные факторы изменяют поток космических лучей в больших пространственных масштабах, а атмосферные эффекты имеют, как правило, достаточно локальный характер. Решение этих задач возможно с помощью широкоапертурного координатного детектора, способного одновременно регистрировать мюоны с различных направлений в широком диапазоне зенитных углов для получения пространственно-временных характеристик потока мюонов. Такая информация позволяет получать полную картину состояния верхних слоев атмосферы и прослеживать динамику ее изменения, в частности, выявлять возмущенные области, определять направления и скорости их перемещения и оценивать время появления в заданной точке.
Поэтому создание аппаратуры и развитие методов дистанционного мониторинга состояния гелиосферы и магнитосферы, а также термодинамических процессов в атмосфере Земли, приводящих к образованию таких опасных явлений как ураганы, является весьма актуальной задачей.
Цель работы
Разработка и создание Установки для РАспознования Грозовых АНомалий (УРАГАН) - трекового мюонного годоскопа для исследования динамических процессов в атмосфере Земли и околоземном пространстве методом мюонной диагностики.
Научная новизна
Созданная установка впервые позволяет получать с высокой достоверностью непрерывные распределения характеристик потока мюонов из верхней полусферы в широком диапазоне зенитных углов вплоть до 80°. Аппаратура представляет собой новый тип наземного модульного детектора космических лучей, предназначенного для исследования связей между пространственно-временными вариациями потока мюонов космических лучей и различными динамическими процессами в гелиосфере, магнитосфере и атмосфере Земли. Оригинальность разработанной системы быстрого считывания информации мюонного годоскопа подтверждена полученным патентом.
Практическая значимость
Созданная установка позволяет вести непрерывный мониторинг состояния атмосферы и магнитосферы Земли и развивать методы мюонной диагностики для дистанционного и заблаговременного обнаружения опасных явлений (ураганы, магнитные бури и т.п.). Универсальность отдельных элементов и всей системы триггирования и сбора данных позволяет использовать их в системах регистрации данных других многоканальных космофизических установок.
Личный вклад
Автор участвовал во всех этапах разработки и создания мюонного годоскопа УРАГАН. Непосредственно автором была разработана, собрана и налажена аппаратура детектора, созданы алгоритмы и программы, обеспечивающие работу системы сбора данных и триггирования, а также специализированные стенды и методики для тестирования и наладки элементов детектора. Автор принимал непосредственное участие в проведении экспериментов и в получении методических и физических результатов.
Автор защищает
1. Установку модульного типа УРАГАН, обеспечивающую годоскопический режим регистрации мюонов в полном объеме.
2. Общую структуру и основные элементы системы сбора данных и триггирования.
3. Алгоритмы и программы, обеспечивающие работу системы сбора данных и триггирования.
4. Стенды и методики для тестирования и наладки элементов детектора.
5. Полученные характеристики и оценки возможностей мюонного годоскопа УРАГАН.
Апробация
Результаты работы были доложены на российских и международных конференциях: Европейском симпозиуме по космическим лучам (ECRS 2004), Международной конференции по космическим лучам (ICRC 2005), Всероссийской конференции по космическим лучам (BKKJI 2006), Баксанской молодежной шкоде экспериментальной и теоретической физики (2006), научной сессии МИФИ (2007), опубликованы в их трудах, а также в четырех статьях в журналах "Приборы и техника эксперимента" и "Известия РАН. Серия физическая".
Объем и структура диссертации
Диссертация состоит из введения, четырех глав, заключения и списка литературы. Объем диссертации: 120 страниц, 72 рисунка, 7 таблиц, 46 наименований цитируемой литературы.
ОСНОВНОЕ СОДЕРЖАНИЕ ДИССЕРТАЦИИ
Во введении приводится краткий обзор аппаратуры, предназначенной для изучения вариаций космических лучей как источника информации о процессах в среде, через которую они проникают.
Вариации интенсивности космических лучей в месте наблюдения могут быть вызваны изменениями мощности источников (включая ускорительные процессы на Солнце и в Галактике), а также различными модуляционными эффектами при распространении KJI в гелиосфере, магнитосфере и атмосфере Земли.
Из-за малой амплитуды вариаций к аппаратуре, предназначенной для их изучения, предъявляются жесткие требования. Она должна обеспечивать непрерывность измерений, высокую стабильность, а также достаточную статистическую точность. Средства мониторинга вариаций KJI можно разделить на два класса: космические и наземные. Детекторы на борту космических аппаратов позволяют проводить прямые измерения
потока первичных частиц, но, к сожалению, их количество мало, а риск необратимого выхода аппаратуры из строя и их стоимость велики. Кроме того, приемлемую статистическую точность измерений вариаций космических лучей на космических аппаратах, из-за ограниченности их размеров, можно получить только в области малых энергий.
При наземных измерениях вариаций космических лучей обычно регистрируется либо нейтронная, либо мюонная компоненты. Нейтронные мониторы регистрируют нейтроны, которые образуются в веществе установки адронами, достигшими уровня наблюдения. В настоящее время поверхность Земли покрыта сетью нейтронных мониторов, включающей более 50 станций. Нейтронные мониторы являются приборами интегрального типа, их пространственное разрешение определяется лишь светосилой, а результаты относятся к потоку КЛ с невысокой средней энергией (2-10 ГэВ).
Вариации космических лучей исследуются также путем измерений интенсивности мюонной компоненты. Для ее измерений обычно используют стандартные телескопы, ориентированные определенным образом. Для наблюдения анизотропии космических лучей мюонные телескопы приспособлены значительно лучше, чем нейтронные мониторы. Это связано, в первую очередь, с их более узкой диаграммой направленности. Первый многонаправленный мюонный телескоп для непрерывной регистрации потока мюонов был создан в Нагое и работает с октября 1970 г. В настоящее время работают несколько наземных многонаправленных мюонных телескопов: в Нагое (Япония), Хобарте (Австралия), Сан Мартино (Бразилия), Кувейте, регистрирующих поток мюонов с различных направлений в реальном времени и в едином формате. Мюонные телескопы дают информацию, которая соответствует КЛ со средней энергией около 50 ГэВ, и позволяют детально изучать их вариации. Телескопы с крупной мозаичной структурой регистрируют КЛ из небольшого числа (4 - 20) пространственных направлений с угловым разрешением 20-30 градусов. В последние годы было создано два счетчиковых многоканальных телескопа - на горе Норикура и большой счетчиковый телескоп индийско-японского проекта в Ооти. Но даже лучшие мюонные телескопы имеют угловое разрешение около 10 градусов.
Новым шагом в развитии методов наземного мониторинга процессов, модулирующих потоки космического излучения, является создание следующего поколения мюонных детекторов - мюонных годоскопов. Под годоскопическим режимом работы мюонного детектора понимается возможность измерения практически непрерывного пространственного распределения потока мюонов, т.е. когда
дискретность установки не вносит искажений в отображение исследуемых физических процессов. Первой в мире попыткой реализации широкоапертурного мюонного годоскопа для исследования вариаций КЛ на поверхности Земли является установка ТЕМП, созданная в МИФИ в 1995 г., угловое разрешение которой составляет 1-2 градуса.
К сожалению, реконструкция треков по двум точкам, конструктивно заложенная как во всех существующих мюонных телескопах, так и в установке ТЕМП, не может дать надежной идентификации направления пролета мюона. В результате пространственное распределение искажается из-за ложных реконструкций, связанных с дельта-электронами, группами мюонов, каскадными ливнями и случайными шумовыми срабатываниями.
В разделе 1 описываются конструкция трекового мюонного годоскопа УРАГАН, системы, обеспечивающие его работу, а также структура системы триггирования и сбора данных установки и ее основные элементы.
В качестве базового регистрирующего элемента мюонного годоскопа применяются камеры стримерных трубок с пространственной локализацией трека около 1 см х 1 см, с двухкоординатным съемом информации, ранее использовавшиеся в эксперименте МШЕХ. Размер камеры 357x16.7x1.5 см3. Основной элемент камеры - 8 стримерных трубок внутренним сечением 9x9 мм, объединённых в один открытый профиль, изготовленный из поливинилхлорида толщиной 1 мм. На профиль нанесено графитовое покрытие с сопротивлением 50 кОм/квадрат, которое выполняет роль катода. По центру каждой трубки расположена анодная нить, изготовленная из бериллий-медного сплава и посеребрённая. Диаметр нити 100 мкм. Внешние считывающие Х- и У-стрипы изготовлены из алюминиевой фольги толщиной 40 мкм. Стрипы имеют разные размеры: ширина Х-стрипов 4 мм и наклеены они с 6 мм интервалом, а У-стрипы 10 мм ширины с 2 мм интервалом. Характерное время распространения импульса по стрипу - 6 нс/м, импеданс - 50 Ом для Х-стрипов и 25 Ом для У-стрипов. Режим ограниченного стримера обеспечивается трехкомпонентной газовой смесью (аргон и углекислый газ в пропорции 1/2 и насыщенные пары нормального пентана при стабилизированной температуре), а также выбором необходимого высоковольтного напряжения (4100-4300 В).
В основу конструкции широкоапертурного мюонного годоскопа УРАГАН заложен модульный принцип его построения. Основой годоскопа является автономный супермодуль (СМ), расположенный на мобильной платформе (рис.1). Число и взаимное расположение
супермодулей в мюонном годоскопе определяется характером решаемых задач и может изменяться в процессе проведения исследований. Минимальный состав мюонного годоскопа - один СМ.
Каждый автономный супермодуль состоит из восьми плоскостей камер стримерных трубок вместе с системой внешних считывающих пластин (стрипов) (всего 320Х+288У каналов с шагом 1,0 и 1,2 см соответственно). Регистрирующие плоскости располагаются послойно, 5 см зазор по вертикали между плоскостями заполнен сплошными пенопластовыми пластинами. Полная толщина установки по вертикали составляет — 1.1 г/см2 и определяется в основном веществом камер. Каждая плоскость детектора состоит из 20 камер, уложенных с одинаковым зазором 3 мм по горизонтали и имеет эффективную площадь -11,5 м2. Конструкция супермодуля обеспечивает высокую пространственную и угловую точность регистрации мюонов (соответственно 1 см и 1 градус) в широком диапазоне зенитных углов от 0 до 80 градусов.
Система газоподготовки мюонного годоскопа УРАГАН аналогична системе, применявшейся в эксперименте ИШЕХ, с расширением некоторых возможностей, и работает по следующему принципу. Скорость потоков аргона и углекислого газа регулируется с помощью прецизионных контроллеров потока с точностью ± 1 % от полной шкалы и повторяемостью ± 0.2 %. Газовая смесь насыщается парами п-пентана при комнатной температуре (около 20°С) и затем «пробулькивается» через п-пентан с контролируемой рабочей температурой (+7°С). Таким образом, излишек паров п-пентана конденсируется и на выходе системы образуется газовая смесь с требуемыми пропорциями.
Необходимое для режима ограниченного стримера высокое напряжение обеспечивается 40-канальным высоковольтным источником CAEN HV SYSTEM, Model SY-127. Низковольтное питание (5.25 В) считывающей электроники осуществляется блоками БНН-152.
Система триггирования и сбора данных установки УРАГАН реализована как совокупность нескольких идентичных, независимых друг от друга подсистем супермодулей с четкой общей синхронизацией по календарному времени. Таким образом, скорость сбора данных с установки определяется быстродействием системы одного супермодуля. Расширение установки легко достигается добавлением новых супермодулей, при этом возрастает статистическая обеспеченность суммарных результатов, а стабильность и идентичность работы супермодулей контролируется по величине отношения регистрируемых интенсивностей.
Управление потоком данных супермодуля и их первичная обработка осуществляются с помощью периферийной ЭВМ -промышленного компьютера с шинами расширения ISA и PCI. Каждая периферийная ЭВМ супермодуля может работать как автономно, так и в составе системы. Для сбора данных и генерации тестовых последовательностей используется плата цифрового ввода/вывода PCI 7200 ADLINK. Выделение событий и точный подсчет временных интервалов производится специально разработанным триггерным блоком.
Эффективность регистрации заряженной частицы одной плоскостью в основном обусловлена геометрией трубок (толщина стенок и т.п.) и составляет около 90 %. Наличие восьми двухкоординатных плоскостей позволяет, с одной стороны, организовать достаточно «мягкий» триггер, повышая тем самым эффективность регистрации, а с другой стороны, при реконструкции выбирать события только с прямыми треками, действительно соответствующими прохождению мюонов. Так, совпадение сигналов не менее четырех сработавших плоскостей из восьми с 99% эффективностью идентифицирует прохождение заряженной частицы через супермодуль.
Триггерные сигналы от восьми плоскостей каждого СМ приходят на триггерный блок, в котором осуществляется выработка триггера супермодуля, сигнала загрузки регистров плат считывания, прерывания в ЭВМ, флагов состояния и последовательности тактовых импульсов считывания в режиме "запрос - подтверждение", а также подсчет живого времени и измерение темпа счета сигналов плоскостей. Структура триггерного блока включает следующие функциональные узлы: интерфейс с шиной ISA; приемник и формирователь по длительности
триггерных сигналов от плоскостей, он же являющийся сканирующим ОЗУ, в которое с частотой 40 МГц записывается состояние триггеров; логическую матрицу совпадений; схему мониторирования темпа счета триггеров и часы «живого» времени; регистры состояния триггеров и всего блока; схему выработки общесистемного триггера и прерывания; модуль выработки тактовых импульсов.
Условием выработки триггера супермодуля является совпадение не менее четырех сигналов сработавших плоскостей в течение 300 не. Расчетное количество случайных совпадений при одинаковом уровне шумов плоскостей можно получить по следующей формуле:
п4(1..8) = С>4Г\
где V - частота триггерных сигналов плоскостей, т - длительность совпадения сигналов. С учетом измеренных шумов плоскостей (~ 3 кГц) и длительности импульсов 300 не количество случайных совпадений составляет ~1.5 события в сутки (при темпе счета частиц 1.8 кГц).
Один супермодуль включает в себя 4864 стрипа -информационных канала. Вследствие такого большого числа каналов, одним из самых важных и дорогостоящих элементов системы считывания является устройство преобразования аналогового сигнала со стримерной камеры до уровня цифровой логики.
Базовым элементом системы сбора данных установки является специально разработанная плата быстрого считывания (ПБС), обеспечивающая усиление, дискриминацию, формирование, хранение сигналов, предварительное триггирование и последовательную передачу данных с 16 стрипов. ПБС включает в себя восемь специализированных двухканальных усилителей-компараторов, два восьмиразрядных сдвиговых регистра и схему управления. 16-канальная компоновка платы обусловлена технологическими особенностями поверхностного монтажа компонентов. Кроме того, наличие схемы управления на каждых 16-ти каналах повышает гранулированность и гибкость системы в целом.
Общая функциональная схема платы изображена на рисунке 2. Усилители-компараторы реализованы на специализированных двухканальных микросхемах Р4040, которые были разработаны на основе БМК Б1451ХНЗ-4 ("Фархад-2") совместно со специалистами НПО «ФОРМ». Микросхема обеспечивает усиление, дискриминацию по заданному порогу и формирование выходных импульсов фиксированной длительности по двум независимым каналам. Для формирования триггерного сигнала по схеме «монтажное ИЛИ» реализован выход логического ИЛИ по двум каналам с открытым коллектором. На выходе микросхем формируются сигналы ТТЛ уровня длительностью 1 мке, которые поступают на входы параллельной загрузки двух
восьмиразрядных сдвиговых регистров. Схема управления платой включает в себя два О-триггера, мультиплексор прохождения данных через плату, формирователь триггерного сигнала платы и схему взвода О-триггеров по триггерному сигналу платы. Платы включаются последовательно между собой, образуя каналы последовательного чтения данных.
Рис. 2. Плата быстрого считывания данных с камер стримерных трубок.
Отличительной особенностью разработанной системы является наличие двух цепочек последовательного считывания данных. Первая цепочка формируется из О-триггеров схемы управления платой и представляет информацию о сработавших платах. Вторая - формируется из сдвиговых регистров и позволяет получить информацию с каждого канала (стрипа). Кроме того, схема управления включает сдвиговые регистры во вторую цепь лишь при наличии сработавших каналов на данной плате. Таким образом, информация считывается только со сработавших плат, что позволяет обеспечить необходимую скорость сбора данных.
В разделе 2 приводится описание аппаратуры, методик проверки, отбора и наладки элементов годоскопа, а также программного обеспечения установки.
Обеспечение годоскопического режима требует решения следующих основных задач:
• тщательного отбора камер стримерных трубок и плат считывания на специализированных стендах;
• создания адекватного программного обеспечения для организации устойчивого сбора и обработки данных в режиме реального времени и их представления в виде, удобном для дальнейшей физической обработки.
Для обеспечения надежной работы установки все камеры перед сборкой проходили полный цикл подготовки, измерения характеристик и отбора на специализированном стенде, который обеспечивает автоматизированное управление системой подачи высокого напряжения и одновременное тестирование восьми камер. В процессе тестирования измеряется эффективность срабатывания камеры в целом и уровень счета каждой трубки камеры при различных значениях высокого напряжения. Была проведена проверка более 1500 камер стримерных трубок.
Методика проверки платы быстрого считывания включает в себя поканальное определение порога 50% эффективности, вычисление среднего значения и разброса порога, проверку перекрестных наводок с канала на канал и поканальную проверку формирования триггерного сигнала. После изготовления все платы прошли наладку на специализированном стенде. По результатам проверки плат средний порог срабатывания канала составил 3.8 ± 0.4 мВ. С учетом точности установки среднего порога срабатывания плат, разброс порогов срабатывания отдельных каналов микросхем Р4040 не превысил ±10%. При нахождении среднего порога срабатывания установки на плато эффективности, 10% разброс отдельных каналов укладывается в пределах данного плато. При этом неоднородность эффективности различных каналов установки составляет величину не более 2%.
Для обеспечения работы наладочных стендов установки УРАГАН был разработан многофункциональный контроллер. Объединение нескольких функций в одной плате позволило минимизировать количество слотов, требуемых для организации системы, и оптимизировать трафик прохождения данных. Контроллер включает в себя четыре идентичных, параллельно функционирующих канала, каждый из которых обеспечивает прием данных с плат считывания, передачу их в ЭВМ, выработку тестовых последовательностей, подсчет количества триггерных сигналов и участвует в выработке общего мажоритарного триггера. Контроллер дает менее 1 сбоя на 1 млн. событий при частоте триггерных сигналов 10-15 кГц.
Отладка и настройка совместной работы платы сбора и триггерного блока в различных режимах производилась сначала отдельно от установки с имитацией событий генератором. Общая сборка измерительной системы проводилась путем последовательного наращивания цепочек плат быстрого считывания на каждой плоскости с
одновременным контролем темпа счета триггерных сигналов плоскостей и проверкой прохождения тестовых последовательностей данных. По результатам контроля определялось наличие сбоев и устранялись такие неисправности, как шумящие каналы, замыкания, обрывы, аномальные паразитные емкости, разница потенциалов стрипов. После сборки и наладки удалось достигнуть устойчивой работы системы в режиме «запрос - подтверждение» на тактовой частоте 2.5 МГц с фиксированным считыванием четырех ПБС по каждой координате всех плоскостей. Доля событий со срабатыванием не более четырех ПБС составила 99%, а время считывания каждого события при этом составило 38 мкс.
Для аппаратной части установки УРАГАН было разработано специальное программное обеспечение, во многом определяющее быстродействие и точность измерительной системы в целом. Аналогично аппаратуре, программное обеспечение системы такой сложности является модульным и многоуровневым. Общая блок-схема программного пакета для полного цикла обслуживания модульного годоскопа, представлена на рис. 3.
Рис. 3. Структура программного обеспечения установки УРАГАН.
В задачу программы управления супермодуля входит: инициализация и управление основными режимами работы системы супермодуля;
декодирование данных и вычисление параметров трека каждого события в режиме реального времени;
• формирование экспозиционных кадров и их передача вместе с
информацией мониторинга супермодуля программе управления
установкой.
Работа программы управления одного СМ может осуществляться как в автономном режиме, так и под контролем программы управления всей установкой. Процесс управления потоком данных и взаимодействия основных модулей программы с аппаратной частью установки происходит следующим образом.
При инициализации системы происходит определение базовых адресов и номеров прерываний триггерной платы и адаптера ввода-вывода. Задаются начальные значения переменных. Определяются адреса и размеры буферов данных. Устанавливаются векторы обработчиков прерываний. Далее управление передается подпрограммам измерения шумов и тестирования. По окончании мониторинга триггерная плата переводится в режим экспозиции, на адаптере ввода-вывода инициализируется режим прямого доступа к памяти и устанавливается разрешение прерываний, а управление передается подпрограмме обработки событий и формирования кадра. При фиксировании события триггерная плата выставляет сигнал прерывания. Обработчик прерываний считывает информацию о сработавших плоскостях и «живом» времени, инкрементирует счетчик событий и производит запись в кольцевой буфер. Адаптер ввода-вывода считывает данные с ПБС и в режиме прямого доступа к памяти записывает их в кольцевой буфер. Подпрограмма обработки событий постоянно анализирует разницу между указателем чтения и указателем записи кольцевого буфера и, при поступлении данных, производит декодировку и вычисление номеров сработавших стрипов, реконструкцию и определение параметров трека, суммирование «живого» времени и формирование матриц направлений треков в экспозиционном кадре. Данные формируются программой по кадрам, в которых производится суммирование статистики в заданном в конфигурационном файле временном интервале (обычно 1 минута). Ход работы отображается на экране монитора.
Программно отбираются так называемые bound треки - треки, полностью проходящие через всю регистрирующую область супермодуля. Определяются вероятности срабатываний плоскостей при прохождении через них частицы на основе анализа сработавших на треке стрипов. Подсчитывается количество срабатываний каждого стрипа в зарегистрированных за время экспозиции событиях. В кадре записываются матрицы направлений восстановленных треков, живое время экспозиции, календарное время и данные мониторирования.
Для хранения углового распределения потока частиц используются три формата матриц:
1) Линейное разбиение по зенитному (в) и азимутальному (ф) углам, Мы{0, <р) - ячейки Д<р = 4°, Ав = 1°, {¡>=0-360°, в = 0-90°;
2) Линейное разбиение по проекционным углам. М(0у, 0Х) - ячейки
3) Разбиение по тангенсам проекционных углов. М^ву, tg0^•) — ячейки ЬЩвхг= Х%^°!Л5.5, вху= 0-80;
В разделе 3 представлены основные характеристики, влияющие на эффективность работы установки и точность получаемых данных, изложена методика их экспресс-анализа.
Исследована зависимость эффективности регистрации детектора от режима высоковольтного питания, атмосферного давления и температуры в помещении расположения установки. В качестве параметра оценки влияния внешних факторов на работу годоскопа была взята величина эффективности срабатывания каждой плоскости установки. Эффективность срабатывания для /-ой плоскости можно выразить следующей формулой:
где - эффективность срабатывания г-ой плоскости; ТУ, - число срабатываний г-ой плоскости в составе реконструированных треков; Нрек - число реконструированных треков в эффективном объёме детектора.
Исследование зависимости основных параметров регистрации частиц супермодулем от изменения напряжения проводилось при незначительном изменении давления и в спокойной геомагнитной обстановке. Корреляционный анализ напряжения питания и эффективности выявил нелинейную зависимость эффективности и темпа счета СМ от напряжения питания с выходом на характерное плато, в пределах которого и была выбрана рабочая точка.
По результатам линейной аппроксимации изменений эффективности плоскостей супермодулей от давления и температуры были получены соответствующие коэффициенты. Большинство значений барометрических коэффициентов лежит в диапазоне 0.01 - 0.02 %/мбар. Температурные коэффициенты варьируются от 0.02 до 0.12 %/"С. Разброс коэффициентов вызван тем, что продув каждой плоскости специальной газовой смесью осуществляется независимо, и, следовательно, условия формирования стримерного разряда в каждой плоскости (и, более того, в каждой камере) являются индивидуальными. Поправки позволяют учесть такое поведение независимым образом для каждой плоскости.
Зависимость интегрального темпа счета всего супермодуля от атмосферного давления может быть описана стандартным выражением для учета барометрического эффекта.
/го = /-/?(Р-Р0)/0хЮО%,
где ¡5 — барометрический коэффициент. Для расчета барометрического коэффициента необходимо построить корреляционную зависимость темпа счета от давления за выбранный период наблюдений. Далее методом линейной регрессии определяется коэффициент В уравнения аппроксимирующей прямой. Барометрический коэффициент связан с параметром В следующим образом:
В = —х100% • h
Зная барометрический коэффициент /? и значение константы 10, можно вводить поправку на барометрический эффект в данные темпа счета, то есть избавиться от изменений темпа счета, вызванных изменением давления. Полученные значения барометрических коэффициентов для СМ8, СМ10 и СМИ равны 0.195 %/мбар, 0.184 %/мбар и 0,189 %/мбар, соответственно. Различия в значениях барометрического коэффициента обусловлены разными эффективными пороговыми энергиями для СМ8, СМ10 и СМ11, которые определяются местом их расположения. При введении поправок на барометрический эффект для мюонных годоскопов необходимо учитывать, что величина барометрического коэффициента зависит от зенитного угла, так как при увеличении зенитного угла возрастает средняя энергия мюона и средний путь, который мюон пролетает от точки генерации. Данная зависимость представлена на рис. 4.
•0.05
•«.10
(0
IO -0.15
J
-0.25
0 10 20 30 40 50 60 70 зенитный угол, градусы
Рис. 4. Зависимость барометрического коэффициента от зенитного угла для трех супермодулей координатного детектора
■ ■ а '»и
■ СМ8 • СМ10 •*. СМ11
■ И и1 I»1
■ _ ■ »ж
Влияние температуры атмосферы на интегральный поток мюонов имеет более сложный характер, т.к. температурный эффект обусловлен изменением температуры по всей высоте атмосферы. Для его точного учета необходимо знать дифференциальные температурные коэффициенты (ДТК), позволяющие вносить поправку в темп счета установки с учетом изменения температуры на всех высотах.
К сожалению, высотное распределение температуры в каждый конкретный момент времени, особенно в периоды изменений состояния атмосферы, неизвестно. Поэтому оценку температурного эффекта обычно проводят по темпу счета детектора, усредненному за одни сутки, и по температуре на поверхности Земли, также усредненной за эти сутки. Метод оценки такой же, как и для определения параметров барометрического эффекта. Полученный таким образом среднемассовый температурный коэффициент для годоскопа УРАГАН составляет величину - 0,15 %ЛС
Стабильность и идентичность работы двух супермодулей установки была проверена путем сравнения их интегральных темпов счета за период всего 2007 г. (рис. 5).
^ 1,010-,
о
л
с
1,008 1,006 1,004 1,002 1,000 0,998 0,996
1 .Лап 1 Маг 1 Мау
1 Бар 1 Моу
2007 год
Рис. 5. Отношение темпов счета срабатывания двух супермодулей УРАГАН.
Отношение частоты срабатываний двух супермодулей мюонного годоскопа очень близко к единице, среднее значение составляет 1.003 ±0.005, что подтверждает идентичность их отклика на прохождение мюонов. Исследования вариаций полного потока из верхней полусферы одновременно двумя идентичными широкоапертурными мюонными годоскопами проводятся впервые и
позволяют значительно повысить достоверность экспериментальной информации, так как на основе сравнительного анализа работоспособности двух супермодулей установки УРАГАН минимизируются систематические неопределенности, связанные с возможными отклонениями работы регистрирующей системы каждого СМ от штатного режима.
Корректность работы супермодулей мюонного годоскопа УРАГАН может быть проверена по азимутальной анизотропии потока мюонов на поверхности Земли. Хорошо известно, что из-за влияния магнитного поля поток частиц с западного направления немного превышает поток с восточного направления. Но влияние перекрытий и стен здания, в котором размещается УРАГАН, может изменить это отношение, что и было обнаружено на супермодулях СМ 10 и СМ11, которые находились в углу экспериментального зала. С целью проверки симметрии собственной светосилы супермодулей был проведен эксперимент, в ходе которого СМ 10 и СМ 11 были передвинуты в центр зала и размещены рядом с СМ8. При таком размещении СМ 10 и СМ11 в центре зала, симметрично по отношению к стенам здания, регистрируемые ими асимметрии потока стали практически идентичными СМ8 и выявили характерную восточно-западную анизотропию, обусловленную магнитным полем Земли.
В процессе работы мюонного годоскопа постоянно проводится экспресс-анализ всех мониторируемых характеристик: атмосферного давления; температуры в помещении установки; количества триггеров за время набора кадра; количества реконструированных событий в кадре; количества bound-событий в кадре; живого времени экспозиции; шумов плоскостей; количества срабатываний плоскостей на треке; таблиц срабатывания стрипов в плоскостях Х- и Y-проекций. Минутные матрицы угловых распределений не всегда могут использоваться для анализа, из-за больших статистических флуктуаций, но просуммированные часовые матрицы могут быть использованы для контроля направления интегрального вектора потока частиц. Этот вектор характеризует светосилу и эффективность самого супермодуля, динамика изменений этого вектора характеризует изменение характеристик как самого супермодуля, так и потока регистрируемых частиц.
Получаемые данные ежемесячно заносятся в паспорт работы супермодуля, в котором отражено поведение основных характеристик и приводится список интервалов времени, которые необходимо исключить при анализе данных супермодуля.
В разделе 4 обсуждаются исследовательские возможности мюонного годоскопа УРАГАН.
Возможности установки по исследованию зависимости динамики потока мюонов во время атмосферных и внеатмосферных процессов были изучены в ходе проведения непрерывных измерений пространственно-временных вариаций интенсивности потока мюонов в 2005-2007 гг. Аппаратура работала стабильно с малой долей времени нештатной работы (не более 3% от полного времени регистрации). В ходе эксперимента были получены первые данные по пространственно-угловым вариациям потока мюонов во время различных нестационарных процессов в атмосфере Земли и околоземном пространстве
Например, в периоды развития форбуш-понижений наряду с резким падением интенсивности потока мюонов наблюдается значительное искажение азимутального распределения потока мюонов, максимум которого за 4 часа поворачивается примерно на 60°. Этот эффект обусловлен вращением Земли и движением конуса чувствительности детектора относительно облака солнечной плазмы. После прохождения максимальной фазы форбуш-понижения такие искажения постепенно исчезают.
Одним из наиболее ярких результатов использования мюонных годоскопов для сканирования небесной сферы в режиме мониторинга различных процессов модуляции потока KJI является возможность получения "мюонных снимков" верхней полусферы. В период до конца суток 9 июля 2006 г. состояние магнитосферы оставалось невозмущенным, однако СМ УРАГАН зарегистрировали ночью 7 июля около 2:00 UTC изменение темпа счета мюонов в определенном угловом секторе небесной полусферы. Это явление исчезло через 4 часа вследствие вращения Земли. Ситуация повторилась ровно через сутки 8 июля в 2:05 UTC. Однако положение возмущения сместилось на небесной полусфере, а амплитуда изменения темпа счета резко увеличилась. Через сутки, когда в зону апертуры СМ10 и СМ11 опять попал данный сектор небесной сферы, понижение дифференциального углового темпа счета мюонов было зарегистрировано значительно раньше - около 22:00 8 июля (рис. 6). Причем отчетливо прослеживалось движение возмущенной области в течение экспозиции данного сектора небесной сферы. Зенитно-азимутальная анизотропия потока мюонов отчетливо проявилась вновь во время начала магнитной бури и регистрировалась в течение всего ее активного периода.
Рис. 6. Мюонные снимки неба 8 и 9 июля 2006 года (усредненные за 1 час экспозиции).
Полученные данные указывают на то, что установка УРАГАН зарегистрировала движение в гелиосфере сгустка магнитной плазмы, образованной во время КВМ на Солнце 6 июля 2006 г., за двое суток до его прихода к орбите Земли.
Одной из основных задач установки УРАГАН является мониторинг вариаций потока мюонов в периоды атмосферных возмущений, в частности, во время мощных гроз. Прохождение грозы над установкой сопровождается своеобразной "впадиной" в темпе счета мюонов. Следует отметить, что наблюдаемое "отражение" поведения давления отличается от обычного "статического" барометрического эффекта. Во время грозы проявляется как бы "динамический барометрический эффект", который приводит как к изменению коэффициента корреляции между интенсивностью мюонов и давлением. Величина "динамического" барометрического коэффициента составляет = -0.33 %/мбар для СМ10 и Рв = -0.32 %/мбар для СМ11. Эти значения почти в два раза превышают обычный "статический" барометрический коэффициент для мюонов.
Особый интерес представляет поиск скрытых периодичностей в потоке мюонов. С использованием разработанных алгоритмов вейвлет-
(Ши1
анализа была проведена обработка данных установки УРАГАН за период 22 - 29 июня 2005 г. с целью выявить в темпе счета квази-периодические сигналы, которые могли возникнуть вследствие волновых процессов в верхних слоях атмосферы (внутренние гравитационные, а также ударные волны). Следует подчеркнуть, что в ходе анализа было выявлен не только сам факт наличия периодического сигнала в потоке мюонов, но и моменты его возникновения (ночь с 24 на 25 июня) и исчезновения (ночь с 26 на 27 июня). Период колебания составляет около пяти часов. Этот сигнал совпадает по времени с прохождением мощного атмосферного фронта через Московскую область. Необходимо отметить, что данная гармоническая составляющая не была обнаружена по измерениям приземного атмосферного давления.
Во время спокойной геомагнитной обстановки зенитно-азимутальные вариации потока мюонов отражают развитие различных атмосферных неоднородностей, например грозовых ячеек. Для проверки этой гипотезы был проведен предварительный анализ накопленного материала мюонных снимков, полученных суммированием зарегистрированных данных тремя супермодулями годоскопа УРАГАН в пятиминутных интервалах. Для исследования было выбрано две грозы: 13 мая и 13 июля 2007.года. Первый снимок 13 мая в 12:54 показывает спокойный фон мюонов над московским регионом. Небольшие флуктуации составляют величину порядка одного среднеквадратичного отклонения. На следующем снимке в 12:59 заметно появление небольшого возмущения на северо-востоке в районе г. Мытищи. Далее видно, что возмущение продвигается на юг и достигает максимума в 13:14, находясь в районе г. Люберцы. Далее возмущение перемещается на юг по направлению к г. Подольск. Примерно через час после начала грозы в 13:49 возмущение практически рассеивается.
В заключении сформулированы основные результаты работы. 1) Разработка аппаратуры.
• Разработана структура супермодуля УРАГАН на базе камер стримерных трубок, обеспечивающая регистрацию мюонов в диапазоне зенитных углов 0° - 84° с высоким пространственным (~ 1 см) и угловым (~ 0.8°) разрешением.
• Разработаны системы триггирования и сбора данных, как для одного супермодуля, так и для совокупности нескольких идентичных супермодулей с общей синхронизацией по календарному времени.
• Разработаны основные элементы измерительной системы супермодуля - плата быстрого считывания и триггерный модуль, обеспечивающие сбор данных с мертвым временем около 10%, подсчет
живого времени с дискретностью до 125 не, мониторинг шумов плоскостей и тестирование цепей считывания.
• Совместно со специалистами НПО ФОРМ разработана специализированная микросхема двухканального усилителя-дискриминатора-формирователя, предназначенная для преобразования сигналов со стрипов детектора в сигналы ТТЛ уровней заданной длительности.
2) Обеспечение годоскопического режима работы.
• Разработаны методики и созданы стенды для тестирования и наладки элементов детектора. Разработан многофункциональный контроллер, алгоритмы и программное обеспечение для тестирования плат быстрого считывания в автоматическом режиме.
• Проверено и налажено более 1000 плат быстрого считывания, обеспечено их функционирование в области плато эффективности с разбросом порогов срабатывания в узком интервале (3.8 ± 0.3 мВ).
• Последовательно в течение 2005 - 2007 гг. собраны, налажены и введены в эксплуатацию три супермодуля мюонного годоскопа УРАГАН.
• Разработано программное обеспечение для обслуживания установки УРАГАН: программы регистрации и управления супермодулями, обслуживания датчиков атмосферного давления и температуры, on-line контроля работы супермодулей.
• Разработан комплекс программ, обеспечивающий непрерывное накопление мюонных снимков-матриц и получение непрерывных пространственно-угловых распределений регистрируемого потока мюонов.
3) Исследование характеристик мюонного годоскопа.
• Проведено несколько методических серий по изучению характеристик супермодулей.
• Исследовано изменение отклика детектора при вариациях атмосферного давления, температуры, рабочего напряжения и состава газовой смеси. Полученные результаты позволили определить основные факторы, влияющие на стабильность работы супермодулей, получить зависимость эффективности плоскостей от различных параметров.
• Измерены барометрический (-0,18 %/мбар) и температурный (-0,15 %/°С) коэффициенты для различных супермодулей. Измерена восточно-западная асимметрия потока мюонов.
• Разработаны методика экспресс-анализа экспериментальных данных и паспорт супермодуля, которые позволяют формировать базы данных для физического анализа.
4) Исследовательский потенциал созданной установки.
• В 2005 - 2007 гг. обеспечено проведение длительных измерений пространственно-временных вариаций потока мюонов. Аппаратура работала стабильно с малым числом сбоев - отношение времени нештатной работы к суммарному времени не превышало 3%. Отношение частоты срабатываний двух супермодулей мюонного годоскопа близко к единице, среднее значение и разброс составили 1.002 ± 0.005.
• Продемонстрирована возможность обнаружения гелиосферных возмущений, вызывающих форбуш-понижепия, за двое суток до их появления в околоземном пространстве.
• Продемонстрирована возможность наблюдения грозовых явлений в потоке мюонов и подтверждено существование волновых процессов в периоды возникновения и развития ураганов.
Основные результаты диссертации опубликованы в следующих работах
1) Н.Н.Вонсовский, К.Г.Компаниец, Ю.Н.Родин, В.В.Шутенко Многофункциональный контроллер для сбора данных со стримерных камер координатного детектора, Приборы и техника эксперимента, 2000, № 2, с. 58-62.
2) Н.С.Барбашина, А.А.Езубченко, Р.П.Кокоулин, К.Г.Компаниец, А.А.Коновалов, А.А.Петрухин, Д.В.Чернов, В.В.Шутенко, Э.Е.Янсон. Координатный детектор для исследования горизонтального потока космических лучей, Приборы и техника эксперимента, №6, 2000, с. 20-24.
3) D.V. Chernov, N.S. Barbashina, G. Mannocchi, R.P. Kokoulin, K.G. Kompaniets, A.A. Petrakhin, O. Saavedra, V.V. Shutenko, D.A. Timashkov, G. Trinchero, I.I. Yashin. Experimental setup for muon diagnostics of the Earth's atmosphere and magnetosphere (the URAGAN project). Proc. 29th ICRC, Pune, India, 2005, v. 2, pp. 457-460.
4) К.Г.Компаниец. Широкоапертурный мюонный годоскоп большой площади УРАГАН. Труды 7-й Баксанской молодежной школы экспериментальной и теоретической физики, т. 2, М., МИФИ, 2007, с. 127-138.
5) Н.В.Ампилогов, К.Г.Компаниец, Учёт влияния внешних факторов на эффективность работы установки УРАГАН. Научная сессия МИФИ-2007. Сб. научн. трудов. М., МИФИ, 2007, т. 7, с. 89-90.
6) Н.С.Барбашина, Р.П.Кокоулин, К.Г.Компаниец, Дж.Маннокки, А.А.Петрухин, О.Сааведра, Д.А.Тимашков, Дж.Тринкеро, Д.В.Чернов, В.В.Шутенко, И.И.Яшин. Широкоапертурный мюонный
детектор для диагностики атмосферы и магнитосферы Земли. Известия РАН. Серия физическая, Т. 71, № 7,2007, с. 1072-1074.
7) Н.С.Барбашина, Р.П.Кокоулин, К. Г. Компанией, Дж.Маннокки, А.А.Петрухин, О.Сааведра, Д.А.Тимашков, Дж.Тринкеро, Д.В.Чернов, В.В.Шутенко, И.ИЛшин. Широкоапертурный мюонный годоскоп большой площади УРАГАН. Приборы и техника эксперимента, 2008, № 2, с. 26-32.
8) К.Г.Компаниец, Д.А.Тимашков, В.В.Шутенко, И.И.Яшин. Система считывания информации со стримерных камер. Патент РФ на изобретение №2327209, регистр. 20 июня 2008 г.
Подписано в печать 12 ,/,•/. Заказ 3 & О Тираж •/ ОО
Типография МИФИ. Москва, Каширское шоссе, 31
Введение.
Раздел 1.
Мюонный годоскоп на стримерных трубках.
1.1 Камера стримерных трубок.
1.2 Конструкция годоскопа.
1.3 Системы обеспечения работы годоскопа.
1.4 Система сбора данных.
1.5 Система триггирования.
1.6 Платы быстрого считывания.
1.7 Устройство преобразования аналогового сигнала.
Раздел 2.
Обеспечение годоскопического режима работы установки.
2.1 Проверка, отбор и наладка элементов годоскопа.
2.1.1 Стенд для проверки камер стримерных трубок.
2.1.2 Стенд для проверки плат считывания.
2.1.3 Многофункциональный контроллер сбора данных.
2.1.4 Наладка и тестирование плат быстрого считывания.
2.1.5 Сборка и наладка измерительной системы супермодулей.
2.2 Программное обеспечение.
2.2.1 Общие принципы построения программного обеспечения установки.
2.2.2 Алгоритмы работы систем СМ.
2.2.3 Программа управления супермодулем.
2.3 Структура данных супермодуля.
Раздел 3.
Характеристики монного годоскопа УРАГАН и методика экспресс - анализа данных
3.1. Зависимость эффективности регистрации от внешних факторов.
3.1.1 Исследование зависимости эффективности регистрации от напряжения высоковольтного питания.
3.1.2 Исследование зависимости эффективности регистрации от атмосферного давления и температуры.
3.2. Интегральный темп счета и барометрический эффект.
3.3 Температурный эффект.
3.4. Корреляционный анализ работы супермодулей.
3.5 Измерение восточно-западной анизотропии потока мюонов.
3.6. Методика экспресс-анализа.
3.7. Паспорт работы супермодулей.
Раздел 4.
Исследовательские возможности мюонного годоскопа.
4.1 Исследование корреляций между вариациями потока мюонов на поверхности Земли и магнитосферными явлениями.
4.1.1 Динамика относительной азимутальной анизотропии мюонного потока во время форбуш-понижений.
4.1.2 Мюонные снимки динамики возмущений в магнитосфере Земли и гелиосфере.
4.2 Вариации потока мюонов во время атмосферных аномалий.
4.2.1 Антикорреляция темпа счета и атмосферного давления во время гроз
4.2.2 Поиск скрытых периодичностей в потоке мюонов.
4.2.3 Мюонные снимки грозовых ячеек.
Космические лучи (КЛ) несут информацию как о физических явлениях в самих источниках излучения, так и о процессах в среде, через которую они проникают. Временные вариации интенсивности космических лучей в месте наблюдения могут быть вызваны изменениями мощности источников, включая ускорительные процессы на Солнце и в Галактике, а также различными
1 I модуляционными эффектами при распространении космических лучей в гелиосфере, магнитосфере и атмосфере Земли.
За пределами гелиосферы первичные космические лучи высокой энергии характеризуются высокой степенью изотропии и в пределах ошибок измерений их интенсивность практически постоянна во времени. Но в гелиосфере и на поверхности Земли интенсивность космических лучей может сильно изменяться. По характеру изменения интенсивности наблюдаемые вариации можно разделить на регулярные (периодические) и нерегулярные (спорадические или стохастические). К периодическим вариациям относятся сезонные изменения интенсивности, связанные с вращением Земли (солнечно-суточные, звездно-суточные), вращением Солнца (27-дневные) и циклами солнечной активности (11 и 22-летние). К нерегулярным вариациям относятся изменения интенсивности космических лучей, обусловленные метеорологическими эффектами, а также вариации, вызванные солнечными вспышками и другими нестационарными явлениями на Солнце, например, выбросами коронального вещества в межпланетную среду. Как правило, большинство типов вариаций характеризуются небольшим изменением интенсивности, амплитуды которых могут изменяться от долей процента до нескольких процентов. В области небольших энергий потоки частиц от солнечных вспышек могут превышать средний уровень в несколько десятков и сотни раз [1].
Из-за малой амплитуды вариаций к аппаратуре, предназначенной для их изучения, предъявляются жесткие требования. Она должна обеспечивать непрерывность измерений, высокую стабильность, а также достаточную статистическую точность.
Средства мониторинга вариаций КЛ можно разделить на два класса: космические и наземные.
В настоящее время измерения потока частиц непосредственно в точке своего местонахождения проводят космические аппараты АСЕ, WIND, STEREO и GOES и др. Находящиеся на их борту детекторы позволяют в реальном времени измерять величину и направление магнитного поля, плотность, температуру и скорость плазмы солнечного ветра, потоки электронов с энергиями от нескольких эВ до десятков МэВ, протонов- от нескольких кэВ до сотен МэВ, ионов (вплоть до железа) - до энергий 30 МэВ / нуклон и другие характеристики.
Детекторы на борту космических аппаратов проводят прямые измерения потока первичных частиц, но, к сожалению, их количество мало, а риск необратимого выхода аппаратуры из строя и их стоимость велики, кроме того, приемлемую статистическую точность измерений вариаций космических лучей на космических аппаратах из- за ограниченности их размеров можно получить только в области малых энергий.
При наземных измерениях регистрируется вторичное излучение: мягкая, жесткая и ядерно-активная компоненты.
На протяжении полувека большой круг физических явлений, связанных с вариациями интенсивности КЛ, изучается с помощью разветвленной мировой сети нейтронных мониторов. Нейтронные мониторы регистрируют нейтроны, которые образуются в веществе установки адронами, достигшими уровня наблюдения. Поскольку нейтронные мониторы регистрируют вторичные частицы, возникающие в результате многократных актов взаимодействий в глубине атмосферы, их показания соответствуют интегральной интенсивности ядерно-активной компоненты космических лучей в атмосфере, которая слабо зависит от температурного профиля атмосферы. Поэтому нейтронные мониторы удобны для исследования вариаций внеатмосферного происхождения. Необходимо лишь учесть барометрический эффект. В настоящее время поверхность Земли покрыта сетью нейтронных мониторов, включающей более 50 станций [2].
Однако нейтронные мониторы чувствительны только к ядерно-активной компоненте космических лучей и не позволяют определять направления прихода частиц. Нейтронные мониторы являются приборами интегрального типа, их пространственное разрешение определяется лишь светосилой, а результаты относятся к потоку КЛ с невысокой средней энергией (2-10 ГэВ).
Вариации космических лучей исследуются также путем измерений интенсивности жесткой и мягкой компонент. Жесткая компонента космических лучей в атмосфере - поток мюонов, образующихся при распаде заряженных пионов, которые в свою очередь генерируются в процессах множественного рождения при взаимодействии космических нуклонов высоких энергий с ядрами атомов атмосферы. Поскольку энергичные мюоны характеризуются большой проникающей способностью, частицы, образованные в первых актах взаимодействий, дают существенный вклад в суммарный счет детекторов, находящихся на уровне моря или под землей. Эффективная глубина образования мюонов и, соответственно, доля мюонов, которые могут достичь наземных или подземных установок, зависят от состояния атмосферы. Интенсивность мягкой -электронно-фотонной компоненты также подвержена сильному влиянию различных атмосферных возмущений. Для их измерений обычно используют стандартные телескопы, ориентированные определенным образом. Для наблюдения анизотропии космических лучей мюонные телескопы приспособлены значительно лучше, чем нейтронные мониторы. Это связано, в первую очередь, с их более узкой диаграммой направленности.
Первый многонаправленный мюонный телескоп для непрерывной регистрации потока мюонов был создан в Нагое в рамках программы Международного Года Активного Солнца (IASY) и работает с октября 1970 г. по настоящее время [3]. Телескоп смонтирован внутри термостатированной камеры (7.2 х 7.2 х 4.3 м3), установленной в кондиционируемом помещении (117 м2). Эффективная площадь телескопа составляет 36 м2. Мюонный телескоп состоит из двух слоев по 36 детекторов в каждом. Слои разделены по высоте 1.73 м стальной фермой и свинцовыми блоками толщиной 5 см, расположенными под верхней плоскостью для подавления мягкой компоненты космических лучей. Каждый детектор состоит из пластикового сцинтиллятора толщиной 5 см и площадью 1 м2, который устанавливается в корпус из 1.6 мм стали и просматривается 5" фотоумножителем (Hamamatsu Photonics R877). Детекторы установлены по слоям в квадраты (6 х 6м2), ориентированные по сторонам света. Темп счета одного детектора составляет примерно 1.1 х 104 событий/мин в верхнем слое и 0.9 х 104 событий/мин в нижнем. Температурная зависимость эффективности срабатывания детектора составляет -0.035 %/°С но, благодаря термостабилизации, колебания общей эффективности установки не превышают 0.01 %. Импульсы от 72 детекторов верхнего и нижнего слоев поступают на 992 ячейки двойных совпадений, которые объединяются в группы, выделяющие события по 17 направлениям телескопа. Количество событий в 10 мин интервале с каждого направления вместе с меткой времени, показаниями давления и температуры непрерывно записывается, формируя ряд 10 мин кадров.
В настоящее время работают несколько наземных многонаправленных мюонных телескопов (таблица 1): в Нагое (Япония), Хобарте (Австралия), Сан Мартино (Бразилия), Кувейте, регистрирующих поток мюонов с различных направлений в реальном времени и в едином формате [4].
Таблица 1. Мюонные телескопы.
Станция Площадь, m2 Количество направлений Количество слоев
Nagoya, Japan 36 17 2
Hobart, Australia 9 13 2
Sa~o Martinho, Brazilb 28 17 2
Kuwait University 9 13 2
Действующие мюонные телескопы дают информацию, которая соответствует К/1 со средней энергией около 50 ГэВ, и позволяют детально изучать их вариации. Телескопы с крупной мозаичной структурой регистрируют КЛ из небольшого числа (4 - 20) пространственных направлений с угловым разрешением 20-30 градусов.
Следующим этапом развития методов исследования вариаций потока мюонов на поверхности Земли стало создание мюонных телескопов с высокой степенью грануляции. Это дает возможность одновременной регистрации интенсивности мюонов сразу по многим направлениям. В последние годы было создано два счетчиковых многоканальных телескопа - на горе Норикура [5] и большой счетчиковый телескоп индийско-японского проекта в Ооти [6].
Многоканальный мюонный телескоп на вершине горы Норикура (Япония) У начал свою работу с мая 1998 года (географические координаты 36.1° северной широты и 136.6° восточной долготы на высоте 2770 м над уровнем моря). Телескоп состоит из четырех слоев пропорциональных счетчиков, по 44 счетчика в слое. Каждый счетчик представляет собой пятиметровый цилиндр с диаметром 10 см вдоль оси которого натянута 50 мкм вольфрамовая анодная нить. Для подавления мягкой компоненты космических лучей над верхним слоем размещен слой свинца толщиной 5 см. Оси счетчиков в верхнем и в третьем слоях ориентированы в направлении восток - запад, а во втором и в нижнем слоях - в направлении север - юг. Первый и второй слой составляют верхнюю координатную пару, а третий и четвертый - нижнюю. Пары разнесены друг от друга на высоту 80 см. Запись события прохождения мюона происходит по совпадению срабатывания счетчиков в четырех слоях. Направление мюона определяется по положению X-Y сработавших счетчиков. Структура установки приблизительно эквивалентна двум 44 х 44 квадратным сборкам 10 х 10 см детекторов, вертикально разделенных друг от друга на 80 см. Диапазон измерения углов составляет 360° по азимуту и от 0° до 55° по зенитному углу. Угловое разрешение - около 7°.
Эксперимент GRAPES-3 индийско-японского проекта в Ооти (11.4° северной широты 76.7° восточной долготы, 2200 м над уровнем моря) первоначально планировался для исследования космических лучей высоких энергий (выше 3 х 1013эВ). Установка включает в себя 380 сцинциляционных счетчиков ШАЛ, размещенных на площади более 20 000 м2, а также большой мюонный телескоп общей площадью 560 м2. Базовым элементом телескопа является пропорциональный счетчик (PRC) размером 600 х 10 х 10 см, работающий на смеси 90% аргона и 10% метана. Сборка из 58 PRC, размещенных рядом в одной плоскости, представляют собой один слой детектора. Две пары ортогональных слоев счетчиков проложенных бетонным поглотителем и экранированных сверху бетонными плитами общей толщиной 2,4 м называются модулем (рис 1). Пороговая энергия установки определяется слоем поглотителя и составляет 1,1 ГэВ для вертикальных мюонов. Четыре модуля, установленных рядом в одном здании, называется мюонной станцией. Весь телескоп состоит из четырех таких станций.
45° !
Рис. 1. Модуль мюонного телескопа установки GRAPES 3. Триггерным условием срабатывания модуля является прохождение мюона через все четыре слоя PRC. Темп счета одного модуля составляет примерно 3200 8
Гц. Угловое разрешение детектора составляет 9°. В 1998 году возникла идея использовать возможности гигантского мюонного телескопа для исследования вариаций солнечных космических лучей, обусловленных такими явлениями, как вспышки, выбросы коронального вещества, форбуш-понижения. Для решения этой задачи была разработана система быстрой записи мюонных событий, позволяющая записывать треки частиц из 225 направлений, и в 2001 г. вся установка была оснащена новой электроникой. Система записывает около 3 х 106 событий в минуту, что составляет 9 Гбайт данных в сутки. В 2003 г. этой системой сбора данных был оборудован мюонный телескоп в Акено.
Однако описанные выше многонаправленные мюонные телескопы обладают рядом существенных недостатков. В этих установках используются круглые и квадратные пропорциональные камеры, сложенные в виде детектирующих плоскостей. Внутри апертуры имеется много зон пониженной эффективности на границах вдоль каждого счетчика, обусловленных щелевыми зазорами между ними. Даже лучшие мюонные телескопы имеют угловое разрешение около 10 градусов.
Новым шагом в развитии методов наземного мониторинга процессов, модулирующих потоки космического излучения, является создание следующего поколения мюонных детекторов - мюонных годоскопов. Под годоскопическим режимом работы мюонного детектора понимается возможность измерения практически непрерывного пространственного распределения потока мюонов, т.е. когда дискретность установки не вносит искажений в отображение исследуемых физических процессов.
Первой в мире попыткой реализации широкоапертурного мюонного годоскопа для исследования вариаций КЛ на поверхности Земли является установка ТЕМП [7], созданная в МИФИ в 1995г. Годоскоп представляет собой сборку из двух пар координатных плоскостей на основе узких трехметровых сцинтилляционных счетчиков-стрипов, имеет площадь 3 х 3 м2 и содержит 512 регистрирующих каналов.
Конструктивно, отдельный счетчик состоит из обернутой в светоотражающую майларовую пленку и светоизолированной снаружи сцинтилляционной полосы-стрипа (2,5 см х 1,0 см х 290 см), механически соединенной с дюралевым корпусом, где размещается ФЭУ-85 с делителем высокого напряжения. ФЭУ-85 регистрирует вспышки света, образующиеся в результате взаимодействия падающих заряженных релятивистских частиц с веществом сцинтиллятора.
Два слоя счетчиков по 128 штук образуют координатную плоскость. Чтобы не было щелей между счетчиками, соседние пластины в слое размещены с перекрытием 3 мм (рис. 2), что и позволило обеспечить годоско пи чески й режим работы этой установки. Сцинтилляционные стрипы одного слоя уложены перпендикулярно стрипам другого. Триггерным усповием срабатывания установки является прохождение мюона через все четыре слоя стрипов. Если мюон проходит через плоскость, то срабатывают два перпендикулярных стрипа (рис. 3). Их пересечение дает точку мюонного следа. Вторая точка находится во второй плоскости. Через эти точки проводится прямая в пространстве, которая дает положение трека мюона.
Установка снабжена поворотной рамой, позволяющей наклонять детектор до зенитного угла 45°.
Рис. 2. Размещение пластин сцинтилляторов в одной плоскости детектора
Мюоны регистрируются практически с любого направления верхней полусферы. Темп счета детектора составляет около 600 событий в секунду, угловое разрешение - порядка 1-2 градусов. Трек мюона, проходящего через две пары взаимно пересекающихся стрипов, попадает в некоторую ячейку телесного угла.
Y1
Х1
Рис. 3. Схема определения направления движения мюона
На мюонном годоскопе ТЕМП, был получен ряд новых физических результатов. В частности, был успешно апробирован новый способ наблюдения опережающих внутренних гравитационных волн, связанных с отдаленной грозовой активностью в атмосфере Земли [8]. Были получены данные по вариациям потока мюонов, предшествующим геомагнитным возмущениям, которые свидетельствуют о существовании корреляций между динамическими процессами в гелиосфере и изменениями интенсивности мюонов и могут быть использованы в качестве предикторов возмущений магнитосферы Земли [9].
К сожалению, реконструкция треков по двум точкам, конструктивно заложенная как во всех существующих мюонных телескопах, так и в установке ТЕМП, не может дать надежной идентификации направления пролета мюона. В результате прстранственные и угловые распределения искажаются из-за ложных реконструкций, связанных с дельта-электронами, группами мюонов, каскадными ливнями и случайными шумовыми срабатываниями.
Вместе с тем в настоящее время внимание исследователей приковано к изучению гелио- и геофизических явлений, предшествующих и сопровождающих катастрофические явления [10, 11], в том числе магнитные бури и ураганы. Эти исследования имеют особое значение для своевременного обнаружения предвестников подобных катастроф. Поэтому создание новой высокочувствительной аппаратуры и развитие методов дистанционного мониторинга состояния гелиосферы, а также термодинамических процессов в атмосфере Земли, сопровождающих формирование и перемещение ураганов, является весьма актуальной задачей.
Цель работы
Разработка и создание Установки для РАспознования Грозовых АНомалий (УРАГАН) - трекового мюонного годоскопа для исследования динамических процессов в атмосфере Земли и околоземном пространстве методом мюонной диагностики.
Научная новизна
Созданная установка впервые позволяет получать с высокой достоверностью непрерывные распределения характеристик потока мюонов из верхней полусферы в широком диапазоне зенитных углов вплоть до 80°. Аппаратура представляет собой новый тип наземного модульного детектора космических лучей, предназначенного для исследования связей между пространственно-временными вариациями потока мюонов космических лучей и различными динамическими процессами в гелиосфере, магнитосфере и атмосфере Земли.
Оригинальность разработанной системы быстрого считывания информации мюонного годоскопа подтверждена полученным патентом.
Практическая значимость
Созданная установка позволяет вести непрерывный мониторинг состояния атмосферы и магнитосферы Земли и развивать методы мюонной диагностики для дистанционного и заблаговременного обнаружения опасных явлений (ураганы, магнитные бури и т. д.).
Универсальность отдельных элементов и всей системы триггирования и сбора данных позволяет использовать их в системах регистрации данных других многоканальных космофизических установок.
Личный вклад
Автор участвовал во всех этапах разработки и создания мюонного годоскопа УРАГАН. Непосредственно автором была разработана, собрана и налажена аппаратура детектора, созданы алгоритмы и программы, обеспечивающие работу системы сбора данных и триггирования, а также специализированные стенды и методики для тестирования и наладки элементов детектора. Автор принимал непосредственное участие в проведении экспериментов и в получении методических и физических результатов.
Автор защищает
1. Установку модульного типа УРАГАН, обеспечивающую годоскопический режим регистрации мюонов в полном объеме.
2. Общую структуру и основные элементы системы сбора данных и триггирования.
3. Алгоритмы и программы, обеспечивающие работу системы сбора данных и триггирования.
4. Стенды и методики для тестирования и наладки элементов детектора.
5. Полученные характеристики и оценки возможностей мюонного годоскопа УРАГАН.
Апробация
Результаты работы были доложены на российских и международных конференциях: Европейском симпозиуме по космическим лучам (ECRS 2004), Международной конференции по космическим лучам (ICRC 2005), Всероссийской конференции по космическим лучам (BKKJ1 2006), Баксанской молодежной школе экспериментальной и теоретической физики (2006), научной сессии МИФИ (2007), опубликованы в их трудах, а также в четырех статьях в журналах "Приборы и техника эксперимента" и "Известия РАН. Серия физическая".
Объем и структура диссертации
Диссертация состоит из введения, четырех глав, заключения и списка литературы. Объем диссертации: 120 страниц, 72 рисунка, 7 таблиц, 46 наименований цитируемой литературы.
Основные результаты работы можно сформулировать следующим образом.
1) Разработка аппаратуры.
• Разработана структура супермодуля УРАГАН на базе камер стримерных трубок, обеспечивающая регистрацию мюонов в диапазоне зенитных углов 0° - 84° с высоким пространственным 1 см) и угловым 0.8°) разрешением.
• Разработаны системы триггирования и сбора данных, как для одного супермодуля, так и для совокупности нескольких идентичных супермодулей с общей синхронизацией по календарному времени.
• Разработаны основные элементы измерительной системы супермодуля -плата быстрого считывания и триггерный модуль, обеспечивающие сбор данных с мертвым временем около 10%, подсчет живого времени с дискретностью до 125 не, мониторинг шумов плоскостей и тестирование цепей считывания.
• Совместно со специалистами НПО ФОРМ разработана специализированная микросхема двухканального усилителя-дискриминатора-формирователя, предназначенная для преобразования сигналов со стрипов детектора в сигналы ТТЛ уровней заданной длительности.
2) Обеспечение годоскопического режима работы.
• Разработаны методики и созданы стенды для тестирования и наладки элементов детектора. Разработан многофункциональный контроллер, алгоритмы и программное обеспечение для тестирования плат быстрого считывания в автоматическом режиме.
• Проверено и налажено более 1000 плат быстрого считывания, обеспечено их функционирование в области плато эффективности с разбросом порогов срабатывания в узком интервале (3.8 ± 0.3 мВ).
• Последовательно в течение 2005 - 2007 гг. собраны, налажены и введены в эксплуатацию три супермодуля мюонного годоскопа УРАГАН.
• Разработано программное обеспечение для обслуживания установки УРАГАН: программы регистрации и управления супермодулями,
114 обслуживания датчиков атмосферного давления и температуры, on-line контроля работы супермодулей.
• Разработан комплекс программ, обеспечивающий непрерывное накопление мюонных снимков-матриц и получение непрерывных пространственно-угловых распределений регистрируемого потока мюонов.
3) Исследование характеристик мюонного годоскопа.
• Проведено несколько методических серий по изучению характеристик супермодулей.
• Исследовано изменение отклика детектора при вариациях атмосферного давления, температуры, рабочего напряжения и состава газовой смеси. Полученные результаты позволили определить основные факторы, влияющие на стабильность работы супермодулей, получить зависимость эффективности плоскостей от различных параметров.
• Измерены барометрический (- 0,18 %/мбар) и температурный (- 0,15%/°С) коэффициенты для различных супермодулей. Измерена восточно-западная асимметрия потока мюонов.
• Разработаны методика экспресс-анализа экспериментальных данных и паспорт супермодуля, которые позволяют формировать базы данных для физического анализа. S
4) Исследовательский потенциал созданной установки.
• В 2005 - 2007 гг. обеспечено проведение длительных измерений пространственно-временных вариаций потока мюонов. Аппаратура работала стабильно с малым числом сбоев - отношение времени нештатной работы к суммарному времени не превышало 3%. Отношение частоты срабатываний двух супермодулей мюонного годоскопа близко к единице, среднее значение и разброс составили 1.002 ± 0.005.
• Продемонстрирована возможность обнаружения гелиосферных возмущений, вызывающих форбуш-понижения, за двое суток до их появления в околоземном пространстве.
• Продемонстрирована возможность наблюдения грозовых явлений в потоке мюонов и подтверждено существование волновых процессов в периоды возникновения и развития ураганов.
В заключение считаю своим приятным долгом выразить глубокую благодарность моему научному руководителю, профессору А.А.Петрухину за постоянную поддержку в работе и многочисленные обсуждения возникающих проблем.
Автор благодарен всем сотрудникам ЭК НЕВОД, участвовавшимм в создании установки УРАГАН и в проведении долговременных серий измерений. Особенно я признателен Н.С.Барбашиной, В.В.Борогу, А.Н.Дмитриевой, В.В.Киндину, Р.П.Кокоулину, Д.А.Тимашкову, Д.В.Чернову, В.В.Шутенко, И.И.Яшину за помощь в процессе проведения данной работы, плодотворные обсуждения результатов и критические замечания.
Автор благодарит Ю.Н.Родина (СНИИП) за помощь в разработке многофункционального контроллера и С.А.Коршунова (НПО "ФОРМ") за участие в совместной разработке микросхемы Р4040, а также А.Н.Сытина (ИФВЭ) за оказанное содействие в изготовлении электронных плат.
Заключение.
1. Л.И.Дорман, И.Я.Либин. Короткопериодические вариации интенсивности космических лучей. // УФН. Т. 145.1985. вып.З. 403-440.
2. В.Г.Янке. Наземные наблюдения космических лучей и представление информации в интернет в реальном времени. // Солнечно-земная физика. Иркутск. 2002. Вып.2(115). с. 99-103.
3. Z.Fuji, S.Sakakibara, K.Fujimoto, H.Ueno. Multi directional cosmic ray intensities, Nagoya, 1986-1988. // Report of cosmic-ray research laboratory. Nagoya university. Nagoya. Japan. 1990.
4. T.Kuwabara, et. al. (2006), Real-time cosmic ray monitoring system for space weather. //Space Weather, 4, S08001, doi:10.1029/2005SW000204.
5. Y.Ohashi, A.Okada, T.Aoki et. al. New narrow angle muon telescope at Mt. Norikura. // Proc. 25th ICRC. Durban. 1997. v.1. p. 441- 444.
6. S.Kawakami, K.Fujimoto, S.K.Gupta et al. Search of muons in association with large solar flares with the Grapes-3 MMT at Ooty. // Proc. 28th ICRC. 2003. Tsukuba. p. 3405-3408.
7. В.В.Борог, А.Ю.Буринский, В.В.Дронов. Мюонный годоскоп для исследования солнечно-земных связей в области энергий больше 10 ГэВ. // Изв.РАН. сер.физ.,1995, т.59, вып.4, с.191-194.
8. В.В.Борог, В.В.Дронов. Изучение короткопериодических колебаний интенсивности мюонов, связанных с конвективно-грозовыми явлениями в атмосфере Земли. // Изв РАН^ Сер. физ.,1999, т.63, вып.8, с. 1675-1677.
9. В.В.Борог, В.В.Дронов. Изучение динамических процессов в гелиосфере с помощью мюонного годоскопа. // Астрономический вестник, Наука, Москва, 2000, т.34, № 2, стр. 126-130.
10. L.I.Dorman. Space weather and dangerous phenomena on the Earth: principlesiof great geomagnetic storms forecasting by online cosmic ray data. // Annales Geophysicae, 23, 2997-3002, 2005.
11. Г.С.Голицын, А.А.Васильев (ред.). Гидрометеорологические опасности. // M. Крук. 2001.296 с.
12. G. Battistoni, E. Bellotti, C. Bloise et al. The NUSEX detector. // NIM 176(1980) 297-300.
13. Battistoni, E.Belotti, C.BIoise et al. The NUSEX detector. // Nucl. Instrum. and Meth. A. 1986. v.245. p.277.
14. Е. larocci, Plastic Strimer Tubes and Their Applicattions in High Energy. // Physics,NIM 217 (1983) 30.
15. G. Battistoni et al., Operation of Limited Streamer Tubes. // NIM 164 (1979) 57.
16. D.V. Chernov, N.S. Barbashina; G. Mannocchi, R.P. Kokoulin, K.G. Kompaniets, A.A. Petrukhin, O. Saavedra, V.V. Shutenko, D.A. Timashkov, G. Trinchero, I.I.
17. Yashin. Experimental setup for muon diagnostics of the Earth's atmosphere andimagnetosphere (the URAGAN project). // Proc. 29th ICRC, Pune, India, 2005, v. 2, pp. 457-460.
18. В.М.Айнутдинов, М.Б.Амельчаков, Н.С.Барбашина, В.В.Киндин, К.Г.Компаниец, А.А.Петрухин, Д.А.Роом, В.В.Шутенко, Э.Е.Янсон, И.И.Яшин. НЕВОД многоцелевой нейтринный детектор на поверхности. Земли. // Инженерная физика, №4, 2000, с.71 - 80.
19. CAEN HIGH VOLTAGE SYSTEM SY-127, User manual
20. ADLINK Technology Inc.: User's Guide, PCI-7200 / cPCI-7200 12MB/S High Speed Digital Input/ Output Card, .Manual Rev. 2.20: October 14, 2000.
21. AMCC: S5933 32-Bit PCI "MatchMaker", February 12, 1997, Revised October 1998.
22. А.И.Кобзарь. Прикладная математическая статистика. // М.: ФИЗМАТЛИТ, 2006. с. 84-88.
23. М. Гук. Аппаратные интерфейсы ПК. Энциклопедия. // СПб.: Питер. 2002. с. 146-174.
24. ALTERA, Data Sheet: FLEX 10K Embedded Programmable Logic Device Family, January 2003, ver. 4.2.
25. ALTERA, Application Note 116 (Configuring APEX 20K, FLEX 10K & FLEX 6000 Devices)
26. Д.А.Комолов, З.А.Мяльк, А.А.Зобенко, А.С.Филиппов. Системы автоматизированного проектирования фирмы ALTERA Max+Plusll и Quartus II. //изд. «РадиоСофт», г.Москва, 2002г.
27. К.Г.Компаниец, Д.А.Тимашков, В.В.Шутенко, И.И.Яшин. Система считывания информации со стримерных камер. // Патент РФ на изобретение № 2327209, зарегистр. 20. июня. 2008 г.
28. Операционные усилители и компараторы. // М.: Изд. Дом « Додека-XXI», 2002. с. 269, 340.
29. В.Г.Домрачев, П.П.Мальцев, И.В.Новаченко, С.Н.Пономарев Базовые матричные кристаллы и матричные БИС. // М.: Энергоатомиздат. 1992. С.7,66.
30. A.Dyachkov. Amplifiers for radiation monitoring system // Proceedings of the lllrd NEXUSPAN Workshop on Microsystems in Envronmental Monitoring. Moscow. 1997. p. 107-108.
31. В.А.Шамис. C++Builder3 Техника визуального прграммирования. // изд. "Нолидж" Москва 1998
32. Н.Н.Вонсовский, К.Г.Компаниец, Ю.Н.Родин, В.В.Шутенко. Многофункциональный контроллер для сбора данных со стримерных камер координатного детектора. // Приборы и техника эксперимента, №2, 2000, с. 58-62.
33. В.Б.Стешенко. ПЛИС фирмы ALTERA: проектирование устройств обработки сигналов. // изд "ДОДЕКА", г. Москва, 2000
34. В.В.Шутенко. Диссертация: Система сбора и обработки данных экспериментального комплекса НЕВОД-ДЕКОР. // Москва, 2007. с. 63-69.
35. Л.И.Дорман. Метеорологические эффекты космических лучей. // М.: Наука, 1972.
36. Л.И.Дорман, В.Г.Янке. К теории метеорологических эффектов космических лучей //. Изв. АН СССР. сер. физ. 1971. т. 35. №12. с.2571-2586.
37. ЗЭ.Л.И.Дорман, Ю.Я.Крестьянников. Среднемассовая температура и спектрографический метод исследования вариаций космических лучей. // Гемагнетизм и аэрономия. 1977. т.17. № 4. с. 622-627.
38. Дж.Вильсон (ред.). Физика космических лучей. // М. Изд. Иностр. Лит. 1954. т.1. 437 с.
39. Н.Ш.Кремер. Теория вероятностей и математическая статистика. // М.: ЮНИТИ-ДАНА, 2007. с. 246-250.
40. А.В.Белов, Е.А.Ерошенко, К.Г.Иванов, В.Г.Янке. Причины и предвестники форбуш-понижений. // Изв. РАН. сер. физ. 1997. т.61. №6. с.1100-1103.
41. С.Д.Данилов, А.И.Свертилов. Внутренние гравитационные волны, генерируемые при прохождении гроз. // Изв. АН ФАО. 1991. т.27, №3. с. 234
42. Р.Отнес, Л.Эноксон. Прикладной анализ временных рядов. М. Мир. 1982. 428 с.
43. Н.М.Астафьева. Вейвлет-анализ: обзор теории и методы применения. // Успехи физических наук, том 166, № 11, 1996, с. 1145241.