ПРЕПРИНТ N 3 ( 1072 )

БАЗА 1-МИНУТНЫХ ЦИФРОВЫХ ГЕОМАГНИТНЫХ ДАННЫХ ПО ОБСЕРВАТОРИИ МОСКВА

(Май 1990 г. - Декабрь 1994 г.)

А.Н.Зайцев, В.И.Одинцов, В.Г.Петров, С.Б.Скобелев

Введение: ретроспектива цифровых магнитно-вариационных наблюдений в обсерватории Москва.

Магнитная обсерватория Москва ( Красная Пахра) играла роль главной обсерватории Советского Союза со времени переезда НИИЗМ в 1944 году в Подмосковье. По мере становления НИИЗМ и превращения его в ИЗМИРАН роль и значение обсерваторских наблюдений возрастало. Особенно ярко это проявилось в период проведения Международного Геофизического Года 1957 -1958 гг., когда Правительством СССР были выделены значительные средства на проведение планетарных геофизических исследований. К этому времени относится и появление новых магнитных обсерваторий и периферийных отделений ИЗМИРАН по территории СССР в Мурманске, Иркутске, Якутске, Алма-Ате и других городах, впоследствии превратившихся в самостоятельные институты.

Важный шаг после МГГ был сделан в середине 60-х годов, когда была реализована программа создания комплексных магнитно-ионосферных станций. При этом ИЗМИРАН играл роль основного научно-методического центра по всем аспектам исследований, включая наземные магнитные наблюдения. Вместе с тем, затратный характер постановки исследований, отрыв их от Высшей школы и передача Магнитно-ионосферной службы из Академии Наук СССР в Гидрометслужбу в 1975 году привел к дезинтеграции единого научного руководства ИЗМИРАН. В самом ИЗМИРАН также были предприняты шаги объективно сужавшие пространство для развития исследований - передача службы Мирового Центра Данных в Геофизический комитет АН и передача СКБ в ИОФАН. Первое привело к существенному торможению развития информационного обеспечения исследований в ИЗМИРАН, второе подорвало развитие экспериментальных исследований.

Попытки организовать цифровые наблюдения в магнитной обсерватории Москва предпринимались с середины 50-х годов. Одна из них включала установку цифровой вариационной станции конструкции Б.Е.Брюнелли. Эта станция имела электромеханическое сканирующее устройство. Основу цифрового преобразователя составлял двухзаходный винт с движущейся кареткой, на которой были закреплены фотоэлементы, отслеживавшие положение светового луча вариометра. К сожалению, в те времена затраты на организацию регулярных наблюдений на цифровой станции были больше, чем затраты на ручную обработку видимой записи, поэтому станция Брюнелли не была включена в регулярную работу обсерватории Москва.

Следующий период в попытках ввести цифровую регистрацию вариаций в обсерватории Москва начался после изобретения кварцевого вариометра В.Н.Бобровым [1]. Массовое производство вариометров в 60-х годах привело к появлению большого количества различных вариантов построения цифровых магнитно-вариационных станций.

Одной из первых успешных попыток получения геомагнитных данных в цифровой форме была разработка платформы сбора данных для обеспечения советско-американского эксперимента по передаче магнитно-вариационных данных через геостационарный спутник АТС-6 [2]. Основу платформы составляла Цифровая магнитовариационная станция ЦМВС-1 с кварцевыми датчиками [3] конструкции В.И.Одинцова. В течении 1975-1976 годов были организованы регулярные наблюдения вариаций магнитного поля в цифровом виде на обсерваториях Архангельск, Минск и Одесса, которые затем продолжались еще несколько лет. Всего было изготовлено 13 экземпляров ЦМВС-1, при эксплуатации которых на сети магнитных обсерваторий СССР был накоплен важный практический опыт. Важным следствием внедрения цифровых станций стало появление работ по разработке концепции и созданию баз геомагнитных данных, включая систему управления базами данных [4 - 5].

В 70-е годы цифровые магнитовариационные станции появились также на магнитных обсерваториях социалистических стран. В рамках специальной рабочей группы КАПГ были проведены сравнительные испытания различных конструкций магнитовариационных станций [6], которые показали перспективность разработок на базе кварцевых вариометров.

На основе накопленного опыта в 1980 году была разработана станция ЦМВС-2 [7], один из первых экземпляров которой был установлен в обсерватории Москва для опытной эксплуатации в 1986 году. К сожалению, из-за организационных трудностей ввести цифровую регистрацию в регулярную программу работы обсерватории Москва не удалось и архив данных по работе ЦМВС-2 за период 1986-1987 гг. к настоящему времени почти полностью утрачен. Всего было выпущено 57 экземпляров цифровых станций ЦМВС-2, размещенных как на сети магнитных обсерваторий СССР, так и в зарубежных обсерваториях [8]. До сих пор некоторые из них находятся в рабочем состоянии. В 1993 году проведена работа по модернизации электроники станции ЦМВС-2, которая позволит обеспечить ее дальнейшую эксплуатацию.

В 1986 году была закончена разработка станции ЦМВС-6 [9], которая отвечает требованиям длительной полуавтоматической эксплуатации в условиях обсерваторий, включая полярные обсерватории. Успешный пример работы ЦМВС-6 был получен в обсерватории Мирный в Антарктиде за 1987-1990 гг. [10]. В обсерватории Москва станция ЦМВС-6 была установлена для опытной эксплуатации в 1990 году. Все работы по установке, эксплуатации и обработке данных выполнены силами сотрудников сектора полярных геофизических исследований. С самого начала работа была направлена на получение конечного результата - создание базы 1-минутных геомагнитных вариационных данных.

К сожалению, проводимые работы пришлись на исторически трудный период распада СССР и, как следствие, на период отсутствия необходимой поддержки исследований со стороны государства. В ИЗМИРАН продолжают работать специалисты по магнитным наблюдениям, обсерватория Москва не прекратила существования и регулярно ведутся наблюдения с фоторегистрацией. Поэтому наличие базы 1-минутных цифровых данных может облегчить сохранение определенного уровня наблюдений в обсерватории, а доступность данных - стимулировать проведение исследований.

Издание препринта с данными наблюдений на ЦМВС-6 в обсерватории Москва за 1990-1994 год несомненно будет полезно для всех специалистов, использующих геомагнитные данные в своих исследованиях, а также явится стимулирующим шагом на пути восстановления нормальной работы сети магнитных обсерваторий России. Большой проблемой остается вопрос о сохранении сети магнитных обсерваторий в необходимом количестве. По нашему мнению, навряд ли удастся сохранить все обсерватории, так как затраты на их содержание несоразмерны с получаемыми результатами.

По инициативе проф. В.П.Головкова с 1992 года Министерство Науки и Технической Политики РФ выделяет специальное финансирование на поддержание сети магнитных обсерваторий России. Однако, этого ограниченного целевого финансирования недостаточно для того, чтобы в должной мере обеспечить сохранение и развитие сети существующих магнитных обсерваторий. По нашему мнению, опыт цифровой регистрации геомагнитных вариаций, накопленный в процессе эксплуатации ЦМВС-6 в обсерватории Москва (ИЗМИРАН), может служить отправной точкой в поисках оптимальных решений аппаратурного оснащения сети геомагнитных обсерваторий России для обеспечения будущих программ исследований.

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

В качестве основного инструмента для цифровых магнитовариационных наблюдений в обсерватории Москва была использована одна из серийных цифровых магнитовариационных станций ЦМВС-6 1986 года выпуска. Основное предназначение станции ЦМВС-6 - измерение и регистрация вариаций магнитной индукции трех ортогональных компонент полного вектора геомагнитного поля в условиях полярных обсерваторий. В состав станции входят блок датчиков магнитной индукции, измерительный блок, блок контроля и управления и блок сетевого питания.

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

Измерительный блок станции ЦМВС-6 содержит четыре независимых канала: три основных - для измерения магнитной индукции и один дополнительный, с трехканальным коммутатором на входе, - для измерения сигналов от дополнительных датчиков, например, термометра и наклономеров. Для накопления, хранения и регистрации измерительной информации в состав измерительного блока входят оперативное запоминающее устройство (ОЗУ) с модулятором. Модулятор преобразует накопленную в ОЗУ цифровую информацию в двухчастотный тональный сигнал для ее записи на магнитную ленту и передачи по линии связи на вход ЭВМ. Объем запоминающего устройства рассчитан на непрерывное накопление более чем суточного (33 часа) массива данных, поступающих с периодичностью 60 с. Для удобства работы с данными вывод их из ОЗУ и регистрация на магнитную ленту производятся автоматически, один раз в сутки, в 00.00 UT. Время регистрации суточного массива данных на магнитную ленту и/или ввода ее в ЭВМ составляет примерно 36 с, полоса тонального сигнала от 3125 до 6250 Гц. Временная привязка результатов измерений осуществляется с помощью встроенных кварцевых часов с календарным устройством, рассчитанным на непрерывную круглогодичную работу измерительного блока.

Основу измерительных каналов станции составляют аналого-цифровые преобразователи со сбалансированным зарядом и разрядом [12 - 13], обеспечивающие непрерывное интегрирование входного сигнала за период измерения. Измерение осуществляется с периодом 1, 10 или 60 с (по выбору оператора). На выходах аналого-цифровых преобразователей формируются усредненные за период измерения значения цифровых отсчетов, которые накапливаются в ОЗУ. Через каждые 60 информационных отсчетов формируется метка времени, содержащая текущее время, дату и служебную информацию о числе и номере регистрируемых дополнительных каналов.

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

Питание станции осуществляется от источника постоянного тока (аккумуляторной батареи) напряжением от 5 до 9 В, либо от радиоизотопного термоэлектрогенератора РИТЭГ (в условиях Антарктиды). Потребляемая мощность не превышает 3-5 ВА при напряжении 6.3 в. Для питания станции в лабораторных условиях от сети переменного тока напряжением 220 В, частотой 50 Гц служит блок сетевого питания.

Основные технические характеристики станции ЦМВС-6
Размещение аппаратуры и организация цифровой регистрации геомагнитных данных в обсерватории Москва.

Цифровая магнитовариационная станция ЦМВС-6 была развернута в вариационном павильоне ИЗМИРАН в конце марта 1990 г., и, после проведения комплексных испытаний и отработки методики ее эксплуатации, переведена в режим непрерывной цифровой регистрации данных с начала мая 1990 г. Наблюдение за работой ЦМВС-6, ввод данных в ЭВМ, их предварительная обработка и формирование базы цифровых геомагнитных данных обсерватории Москва осуществляется силами сотрудников сектора полярных геофизических исследований.

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

Измерительный блок ЦМВС-6 установлен на расстоянии порядка 3 - 5 м от блока датчиков в коридоре вариационного павильона (в связи с чем подвержен влиянию существенных перепадов температур). Там же находится кассетный магнитофон для цифровой регистрации данных и буферная аккумуляторная батарея. Питание ЦМВС-6 осуществляется от сети переменного тока напряжением 220 В. Аккумуляторная батарея напряжением 6.3 В обеспечивает резервное питание станции на время кратковременных отключений питающей сети. Подзаряд батареи осуществляется непрерывно, пока есть сетевое напряжение. Режим подзаряда выбран таким, что зарядный ток не превышает 50 мА для полностью заряженной батареи (компенсируя ее саморазряд) и увеличивается до 1 - 1.5 А при разряде батареи во время пропадания напряжения сети.

Измерительный блок связан двухпроводной линией со зданием, где находится сектор полярных геофизических исследований, на расстоянии 80 м от вариационного павильона. В этом здании установлен блок контроля и управления ЦМВС-6 и ЭВМ типа IBM PC/XT для ввода и предварительной обработки геомагнитных данных. Такое расположение блоков станции обеспечивает возможность непрерывного контроля за работой ЦМВС-6 непосредственно на рабочем месте в лаборатории. Контроль включает в себя проверку точности хода кварцевых часов станции, регистрацию данных на самопишущий прибор и непосредственный их ввод в ЭВМ в реальном масштабе времени с выводом графика на дисплей ЭВМ.

Чтобы обеспечить непосредственный ввод в ЭВМ суточного массива данных, поступающего по линии связи на вход блока контроля и управления, был введен сдвиг показаний кварцевых часов станции на 8 часов, который затем убирается в процессе предварительной обработки данных. Таким образом, момент вывода суточного массива данных из ОЗУ приходится на рабочее время (11 или 12 MSK), и для его ввода в ЭВМ достаточно один раз в сутки включить на несколько минут персональный компьютер. Для исключения потерь данных из-за несвоевременного включения ЭВМ, они регистрируются на кассетном магнитофоне, расположенном рядом с измерительным блоком ЦМВС-6.

Результаты цифровой регистрации геомагнитных вариаций

В Приложении 1 приведены обзорные магнитограммы за 1990 - 1994 гг., сгруппированные по месяцам. Регистрация данных за указанный период велась на магнитную ленту. Надежность регистрации оказалась достаточно высокой, однако наблюдались отдельные случаи потери данных. Причины их носят скорее организационный характер и вызваны, в основном, несвоевременной заменой магнитной ленты в магнитофоне, низким качеством используемой отечественной магнитной ленты и выключениями станции на время проведения профилактических работ. Поэтому, впредь для устранения потерь данных, станция с февраля 1995 г. переведена в режим непосредственного ввода данных в ЭВМ по линии связи.

Предварительный анализ данных и их сравнение с данными геомагнитной обсерватории Москва (ИЗМИРАН) показали хорошую долговременную стабильность станции, высокую надежность ее работы и широкий динамический диапазон. Примером этому служат образцы записи магнитограмм в периоды больших магнитных бурь 24 - 25 марта 1991 г., 28 - 29 октября 1991 г., 08 - 09 ноября 1991 г. и 09 - 10 мая 1992 г., приведенные в Приложении 2.

Анализ долговременной стабильности результатов измерений

Одним из важнейших требований, предъявляемых к обсерваторским станциям, является их долговременная стабильность. В условиях больших искусственных помех в г. Троицке и с учетом экспериментального характера размещения ЦМВС-6, оценка долговременной стабильности работы станции носит предварительный характер. Нами были проанализированы разности среднечасовых значений, вычисленные по одноминутным данным ЦМВС-6 и по аналоговым магнитограммам, полученным в ЦМО за 1991 г. Предполагается выполнить отдельную работу, посвященную этому вопросу, по данным ЦМО за 1992-1994 гг..

На рис. 1 представлены разности среднечасовых значений трех компонент за январь-февраль 1991.

Рис. 1. Разности среднечасовых значений ЦМВС-6 и обсерватории Москва H, D и Z -компонент за январь-февраль 1991 г.

В эти месяцы было меньше всего помех, связанных с работами вблизи датчиков, однако даже на этих графиках скачки присутствуют. Так, на графике разностей Н компонент видно изменение на +5 нТ в конце второго дня и изменения на -5 нТ в 53 день. Просматривается скачок на величину 2-3 нТ одновременно в D и Z компонентах в 49 день. В течение года такие скачки величиной до 20 нТ наблюдались многократно. Просмотр одноминутных значений показывает, что это действительно мгновенные (в течение 1-2 минут) изменения значений поля, и почти всегда удается найти момент, когда значения возвращались к предыдущему уровню. Все это доказывает, что скачки вызваны внешними воздействиями и при анализе долговременной стабильности должны быть исключены.

Скачки величиной более 5-10 нТ действительно удается достаточно надежно убрать, но маленькие скачки на фоне вариаций в одноминутных значениях не заметны, они проявляются только в разностях и надежно исключить их не удается.

По данным за январь-февраль величину долговременного дрейфа можно оценить в -1 нТ/месяц для Н компоненты, около 0 для D, +3 для Z. По данным за весь 1991 г. дрейф оказывается около -30 нТ/год для Н, -15 нТ/год для D и +12 нТ/год для Z. Данные оценки по всей видимости являются завышенными, т.к. в них может входить как температурный дрейф, так и не устраненные скачки, а также и возможный дрейф опорной станции. Кроме того, следует учитывать то обстоятельство, что датчики ЦМВС-6 и опорной станции расположены в разных помещениях, на расстоянии порядка 10 м один от другого, и на разность их показаний могут влиять градиенты поля в местах установки датчиков. Из приведенных в приложении графиков мгновенных значений поля видно, что изменение спокойного уровня за период 1990-1993 гг. не превышает 20 нТ для H и Z компонент и 40 нТ для D, т.е. годовой дрейф порядка, соответственно, 7 и 15 нТ/год. Конечно, это приблизительная оценка, т.к. здесь мы не учитываем вековой ход и техногенные изменения поля.

При внимательном анализе разностей в них можно заметить присутствие суточной периодичности, что свидетельствует о некотором расхождении в ориентации и/или ошибке определения чувствительности. Однако достоверно определить эти параметры по данной выборке не удалось, это предстоит сделать в дальнейшем. В любом случае суммарная ошибка за счет этих факторов не превышает 2-3 %.

Некоторые вопросы использования баз цифровых геомагнитных данных

Полученная база 1-минутных цифровых геомагнитных данных позволяет поставить геомагнитные исследования на качественно новую основу. Аналоговая запись стандартных магнитограмм обсерватории, как правило, мало пригодна для детального анализа геомагнитных вариаций, в то время как цифровые данные позволяют уверенно выделять вариации с периодами от 5-ти минут и анализировать их спектральный состав [14].

В настоящее время на основе адаптивных методов узкополосной фильтрации В.И. Одинцовым разработан для персонального компьютера IBM PC/AT-486 пакет программ спектрально-временного анализа данных. Результаты анализа 1-мин данных позволяют проследить динамику развития регулярных геомагнитных вариаций с периодами от 5 минут и временным разрешением порядка 4 минут.

В качестве примера использования метода спектрально-временного анализа на Рис. 2а приведена магнитограмма обсерватории Москва за 24 марта 1991 г. с 03.00 до 15.00 UT, а на Рис. 2б показаны соответствующие ей спектрально-временные характеристики геомагнитных вариаций по компонентам H, D и Z в диапазоне периодов от 5 до 20 минут (0.8 - 3.3 мГц).

Рис. 2а. Магнитограмма обсерватории МОСКВА за 24 марта 1991 г.

Рис. 2б. Спектрально-временные представления геомагнитных вариаций по компонентам H, D и Z в диапазоне периодов от 5 до 20 минут.

На основе спектрально-временного анализа геомагнитных вариаций можно выделить несколько режимов генерации пульсаций. Первый режим в интервале 07.30 - 09.30 UT характеризуется нерегулярными шумовыми пульсациями. В колебаниях прослеживаются отдельные регулярные всплески, среди которых можно отметить колебания с частотой 2.6 - 3.3 мГц в интервале 08.35 - 08.50 UT. Можно предположить, что эти пульсации являются следствием развития процесса, начало которого приходится на 07.25 UT, и который приводит к генерации пульсаций с частотой 2.1 - 2.6 мГц в интервале 08.00 - 08.25 UT.

Резкое изменение режима колебаний произошло в 10.00 UT. Шумовые колебания сменяются более регулярными. Максимум в спектре колебаний в интервале 10.00 - 10.25 UT приходится на частоты 2.1 - 3.3 мГц. Затем частота плавно уменьшается до 1.7 - 2.2 мГц к 11.00 UT.

В интервале 12.00 - 14.00 UT произошло изменение отношения H/D от 3.4 до 1.3. В то же время наблюдалось изменение во времени частотного состава пульсаций. Пульсации присутствуют только в H и Z -компонентах, плавно изменяясь по частоте от 1.80 - 2.5 мГц в 12.10 UT до 1.50 - 1.80 мГц к 13.10 UT. Особо следует выделить область пульсаций с частотой 1.6 - 1.75 мГц в интервале от 12.55 до 13.25 UT. Наличие пульсаций только в H и Z -компонентах может свидетельствовать о том, что мы имеем дело с поляризованными колебаниями, направление вектора которых не совпадает с направлениями векторов H, D и Z - компонент и меняется с течением времени. Можно также предположить, что эти пульсации являются вынужденными (вторичными) колебаниями, а условия их генерации зависят от места расположения обсерватории.

При общем анализе всего интервала 07.00 - 14.00 UT можно выделить три пакета колебаний:

Эти пакеты следуют с периодом порядка 140 минут, а частота колебаний в них плавно уменьшается от 3.3 до 1.4 мГц.

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

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

Наличие баз цифровых геомагнитных данных и развитие техники обработки данных с использованием персональных компьютеров позволяет автоматизировать определение секторной структуры межпланетного магнитного поля (СС ММП) на основе информации, получаемой с наземных полярных станций в Арктике и Антарктиде. Работа по определению СС ММП много лет велась в лаборатории магнитных возмущений. Возможности использования цифровых баз данных очевидны также для целей вычисления геомагнитных индексов.

В настоящее время в секторе полярных геофизических исследований проводится работа по созданию архива цифровых геомагнитных данных на оптическом диске (CD-ROM). В архив будут включены данные, полученные с помощью цифровых магнитовариационных станций ЦМВС-2 и ЦМВС-6 во время экспедиционных работ на Карском геофизическом полигоне (Ямал), на станции Мирный (Антарктида), на полярных станциях и в геомагнитной обсерватории Москва.

Заключение

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

По нашему мнению, опыт цифровой регистрации геомагнитных вариаций, накопленный в процессе опытной эксплуатации ЦМВС-6 в обсерватории Москва (ИЗМИРАН), может служить отправной точкой в дальнейшем развитии современных методов исследования геомагнитных вариаций и поисках оптимальных решений аппаратурного оснащения сети геомагнитных обсерваторий России для обеспечения будущих программ исследований.

Вся экспериментальная работа, связанная с установкой и эксплуатацией ЦМВС-6 в магнитной обсерватории Москва осуществлялась при непосредственной поддержке заведующего магнитной обсерваторией к.ф.-м.н. Канониди Х.Д. и ведущего инженера Кайнары Л.Н., за что авторы выражают им искреннюю признательность.

Контактные адрес и телефон

База 1-минутных цифровых геомагнитных данных распространяется сектором полярных геофизических исследований ИЗМИРАН.
Данные могут быть поставлены на дискетах в формате IBM PC.
Адрес для справок: 142092, г.Троицк, Московской области, ИЗМИРАН,
сектор полярных геофизических исследований.
Зав. сектором к.ф.-м.н. А.Н.Зайцев.
Телефон: (095) 334-02-93
Телекс: 112405 СПОЛОХ
E-mail: zaitzev@charley.izmiran.rssi.ru

ЛИТЕРАТУРА