Размещено на http://www.allbest.ru/

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

Специальность 25.00.35 - Геоинформатика

Автореферат диссертации на соискание учёной степени

доктора технических наук

Якушев Денис Игоревич

Санкт-Петербург 2008

Работа выполнена в Санкт-Петербургском государственном электротехническом университете «ЛЭТИ» им. В.И. Ульянова (Ленина).

Научный консультант - доктор технических наук, профессор Ковчин И.С.

Официальные оппоненты:

доктор технических наук, профессор Биденко С.И.

доктор технических наук, профессор Гончаров В.К.

доктор физико-математических, старший научный сотрудник Макштас А.П.

Ведущая организация - Государственный научно-исследовательский навигационно-гидрографический институт (ГНИНГИ)

Защита диссертации состоится "_____" ____________ 2008 г. в ______ часов на заседании диссертационного совета Д 212.197.03 Российского государственного гидрометеорологического университета по адресу: 195196, Санкт-Петербург, пр. Металлистов, 3, ауд. 406^б.

С диссертацией можно ознакомиться в библиотеке Российского государственного гидрометеорологического университета.

Автореферат разослан "_____" ____________ 2008 г.

Ученый секретарь

диссертационного совета

д. т. н., проф. Бескид П.П.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы. Компьютеризация всех сторон человеческой деятельности, происшедшая в последней четверти XX века, многократно повысила объёмы фиксируемой информации, в том числе и геоинформации. Накоплены огромные массивы данных, которые требуют анализа с целью последующего использования результатов на практике. Причём процесс накопления данных идёт значительно быстрее, чем их анализ, формулирование выводов и рекомендаций. Таким образом, именно анализ накопленных данных на сегодняшний день является одним из самых узких мест технического прогресса.

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

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

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

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

На основании анализа имеющихся данных о пространственно-временных геофизических процессах с полигармонической структурой был сделан вывод о необходимости разделения объектов исследования на два класса по признаку устойчивости искомых параметров во времени: на стационарные и нестационарные.

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

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

Для некоторых геофизических процессов, обычно описываемых с помощью полигармонических моделей со стационарными параметрами, эмпирически выявляются некоторые отклонения от принятого описания. Например, для ряда солнечной активности (в числах Вольфа) с середины ХIХ века эмпирически установлено отклонение продолжительности периодов от среднего значения (Швабе С.Г., 1851), а с середины ХХ века - некоторая их асимметрия (Эйгенсон М.С., 1963). Учёт этих отклонений может оказаться существенным для понимания, описания и прогнозирования исследуемых геофизических процессов. В то же время модели подобного рода нестационарности параметров практически неизвестны.

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

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

Для достижения поставленной цели требуется решение следующих исследовательских задач.

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

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

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

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

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

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

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

Научную новизну работы определяют следующие положения.

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

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

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

4. Разработана полигармоническая модель ряда солнечной активности, включающие составляющую с периодом 1800 лет, объясняющую современное глобальное потепление климата Земли с позиций изменения солнечной активности. Найдены геоинформационные подтверждения прогноза изменения солнечной активности на протяжении последних 2500 лет на примере климата Европы.

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

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

При этом основную ценность представляют следующие практические результаты диссертационного исследования.

1. Геоинформационный метод и программа для расчёта стационарных параметров полигармонических моделей пространственно-временных геофизических процессов с неравномерной дискретизацией, позволяющие выделять составляющие с периодами, превышающими интервал наблюдения, необходимыми для построения долгосрочных прогнозов. Программа внедрена в ГосНИОРХ.

2. Геоинформационный метод и программа для расчёта параметров моделей взаимозависимости пространственно-временных геофизических процессов, позволяющий объяснять причины взаимозависимости на основе сходства продолжительностей длительных периодов. Программа внедрена в ФГУП "ЦНИИМ" и ФГУП ЦКБ МТ "Рубин".

3. Метод и программа расчёта нестационарных параметров полигармонических моделей пространственно-временных геофизических процессов с полигармонической структурой.

4. Прогноз потепления климата Европы, основанный на полигармонической модели ряда солнечной активности.

5. Прогноз изменения уровня Каспийского моря, основанный на его полигармонической модели.

Положения, выносимые на защиту.

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

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

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

4. Полигармоническая модель ряда солнечной активности, включающая составляющую с периодом 1800 лет, объясняющую современное глобальное потепление климата Земли с позиций изменения солнечной активности. Геоинформационные подтверждения прогноза изменения солнечной активности на протяжении последних 2500 лет на примере климата Европы.

5. Полигармоническая модель ряда уровня Каспийского моря, включающая составляющую с периодом 124 года. Геоинформационные подтверждения прогноза изменения уровня Каспийского моря на протяжении последних 100 лет.

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

Якушев Д.И. Алгоритмы математического моделирования/Д.И. Якушев. СПб.: МГП “Поликом”, 2002г.-100с.

Результаты работы опубликованы в журналах из перечня ВАК: "Известия ГЭТУ" и "Акустическом журнале", а также в других журналах.

Основные положения докладывались на Международных и Российских научно-технических конференциях, таких как "Международная конференция по мягким вычислениям и измерениям": SCM-99, SCM-00, SCM-01, SCM-02 (Санкт-Петербург), Международном научном конгрессе "Фундаментальные проблемы естествознания" (Санкт-Петербург, 1998), Международных научных конференциях "Циклы" (Ставрополь, 1997, 1999, 2002), Международной научно-технической конференции "Physcon" (Санкт-Петербург, 2003), Международном симпозиуме "Speech and Computer" SPECOM'2002 (St.Petersburg, 2002), Международной научно-технической конференции "Методы, средства и технологии получения и обработки измерительной информации", (Пенза, 2002), 55-ой конференции профессорско-преподавательского состава (СПбГЭТУ, 2002), научно-технической конференции "Мониторинг и прогнозирование чрезвычайных ситуаций" (Санкт-Петербург, 1998); школе-семинаре "Актуальные вопросы организации и производства судебных экспертиз" (Санкт-Петербург, 1998), докладывались на семинарах Санкт-Петербургского отделения Метрологической академии России.

Публикации по работе. По теме диссертации опубликовано 44 научных работы, из них 1 монография, 5 статей в журналах из перечня ВАК, 7 рукописей, депонированных в ВИНИТИ, 3 статьи и тезисы к 23 докладам на международных и всероссийских научно-технических конференциях и семинарах.

Структура и объём работы. Диссертация состоит из введения, пяти глав, заключения, 1 приложения и списка литературы, включающего 360 наименований. Основная часть работы изложена на 256 страницах машинописного текста. Работа содержит 46 рисунков и 47 таблиц.

периодограмма геоинформационный полигармонический модель

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

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

В первой главе рассмотрены эмпирические основы установления полигармонической структуры геоинформационных моделей пространственно-временных геофизических процессов. На основе анализа априорной информации об этих процессах установлено, что существует необходимость их описания моделями как со стационарными, так и с нестационарными параметрами. Например, гармоническая структура со стационарными параметрами для ряда изменения солнечной активности впервые была установлена Швабе С.Г. (1851) и в настоящее время не оспаривается (Витинский Ю.И., 1986), хотя имеются эмпирические предпосылки для построения моделей этого процесса и с нестационарными параметрами (например, Кливленд, 1993).

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

,

где x(t_i) - моделируемый временной ряд, n - количество гармоник, A_j=const, T_j=const, _j, j=1,n - амплитуда, период и фаза j-ой гармоники, m_x=const - среднее значение модели, u(t_i) - аддитивный шум.

Геомодель априорно не предусматривает:

- требование отсутствия пропусков в исходном временном ряде x(t);

- требование равномерной дискретизации исходного временного ряда x(t);

- количество входящих в сигнал гармонических составляющих;

- критерий близости между моделью и оригиналом;

- требование кратности искомых периодов;

- ограничения на продолжительность искомых периодов T_j.

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

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

В последующих главах разработаны геомодели и методы расчёта их параметров, свободные от ограничений, выявленных в первой главе.

Во второй главе разработан геоинформационный метод расчёта стационарных параметров полигармонических моделей пространственно-временных геофизических процессов, свободный от недостатков, сформулированных в главе 1. Основными преимуществами разработанного метода являются:

- возможность включения в модель гармонических составляющих с периодом, превышающим длину реализации;

- возможность обработки временных рядов с пропусками и неравномерной дискретизацией;

- высокая разрешающая способность (<2%).

Для расчёта параметров модели был применён метод минимизации остаточной дисперсии:

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

На первом этапе построения модели проводится построение периодограммы - сканирование множества частных моделей, включающих только одну гармоническую составляющую:

При расчётах принимается следующие нулевые гипотезы: A=0, =₀(1+)ⁱ<T, =.0, где T - время наблюдения; ₀ - минимальный рассматриваемый период; - принятая разрешающая способность. Экспериментально показано, что

- рассматривание периодов >T не приводит к изменению получаемых результатов;

- сканирование пространства частных моделей допустимо проводить только по параметру , так как полученный результат не зависит от начальных значений A и . Графическое иллюстрация результатов сканирования приведена на рис. 1.

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

Для расчётов использовался разработанный автором алгоритм поиска локального минимума функции многих переменных. Отличием разработанного алгоритма от схожих, реализованных в известных математических пакетах (например, Matlab, Mathcad), является его программная устойчивость. Так при расчёте параметров полиномиальных моделей в известных пакетах предусмотрены четыре варианта завершения расчётов:

- получение искомого результата;

- возникновение программной ошибки;

- превышение максимально допустимого числа итераций;

- превышение максимально допустимого времени расчётов.

В разработанном алгоритме предусмотрена остановка расчётов только при получении искомого результата. Три последние возможности исключены. Достижение программной устойчивости оказалось возможным благодаря проведённым исследованиям причин возникновения трёх последних вариантов. Проведённый анализ показал, что иногда при производстве расчётов возникает ситуация, когда происходит неограниченное приращение одного из искомых параметров. Причины этого установлены не были. Но оказалось возможным ввести ограничение на последовательное количество приращений параметров, что позволило обеспечить программную устойчивость расчётов. С другой стороны, введение ограничения на последовательное количество приращений могло повлечь за собой преждевременную остановку поиска локального минимума остаточной дисперсии по текущему параметру. Выход из создавшейся ситуации был найден, во-первых, в установке достаточно большого значения начального приращения, а, во-вторых, в наложении ограничений на возможный диапазон исходных данных 0.011000. В противном случае производится их масштабирование, что эквивалентно изменению единицы измерения. Схема программы приведена на рис. 1.

В схеме программы приняты следующие обозначения:

X₀ - значение параметра X до приращения;

N_- - количество отрицательных приращений с шагом _max;

N₊ - количество положительных приращений с шагом _max;

N_max - максимально допустимое количество приращений;

_max - максимально допустимое приращение;

_min - минимально допустимое приращение;

X^* - вычисленное значение параметра X;

да - действие привело к уменьшению дисперсии;

нет - действие не привело к уменьшению дисперсии.

X=X-max

N-=N-+1

2

X*

X0

да

N-=0

нет

N+=0

нет

N->Nmax

да

X=X+max

да

1

3

<min

N+=N++1

N-=0

нет

нет

X=X-

да

да

2

да

нет

X=X+max

нет

N+>Nmax

X=X+2

нет

3

да

нет

X=X

1

=/2

X=X-

Рис. 1 Схема программы приращения текущего параметра X

На втором этапе проводится анализ результатов сканирования путём сравнения полученных значений минимумов остаточной дисперсии с пороговым значением, которое определяется исходя из следующего критерия:

,

где D[y] - дисперсия исходного временного ряда; D[u] - дисперсия шума; N - количество отсчётов исходного временного ряда, k - коэффициент, определяющий уровень значимости принятия решения.

Первоначально критерий был выведен из критерия В.И. Романовского (Видуев Н.Г., Кондра Г.С., 1969). Впоследствии было установлено, что для случаев N₁=N₂=N этот критерий приближает критерий Фишера. При этом соотношение между коэффициентом k и уровнем значимости (по Фишеру) определяется уравнением, полученным путём аппроксимации значений критерия Фишера:

=0.61exp(-0.64k).

При этом, если N>30, _F<10%, а при N>40 _F<5%, где N - количество отсчётов моделируемого временного ряда.

Во всех расчётах использовался именно этот критерий, который, по сравнению с критерием Фишера, характеризуется наглядностью, удобством восприятия и простотой программной реализации. Так при N=100 и k=1 можно говорить о том, что значимым является уменьшение дисперсии на 10%, если k=2, то, соответственно, на 20%.

В результате моделирования климатических временных рядов и последующей проверке адекватности результатов совместно с доктором географических наук Антоновым А.Е. было установлено, что для большинства исследованных временных рядов коэффициент k должен устанавливаться равным 0.8 (=0.36). В отдельных случаях необходимо его снижение до k=0.7 (=0.39). Здесь и далее для удобства восприятия после значения коэффициента в скобках указан уровень значимости (по Фишеру).

В результате анализа периодограммы определяется количество гармоник в искомой полигармонической модели и нулевые гипотезы для дальнейших расчётов. Графическое представление результатов для ряда уровня Каспийского моря (1891-1985гг.) приведено на рис. 2.

Рис. 2 Периодограмма ряда уровня Каспийского моря

Заметим, что по оси периодов отмечены значения, соответствующие нулевой гипотезе, а максимумы функции (-lg()), соответствуют найденным значениям периодов.

В результате анализа периодограммы установлено, что n=2.

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

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

На четвёртом этапе проводится документирование результатов геоинформационного моделирования, которое включает в себя:

- документирование произведённых расчётов;

- расчёт прогноза;

- расчёт числовых характеристик погрешности параметров модели;

- расчёт гистограммы погрешности;

- расчёт коэффициентов авторегрессионной модели погрешности;

- вывод наложенных графиков исходного временного ряда (сплошная линия) и его полигармонической модели (пунктир).

Результаты моделирования ряда уровня Каспийского моря приведены в табл. 1.

Таблица 1

Параметры полигармонической модели ряда уровня Каспийского моря

Номер гармоники	Амплитуда (см)	Период (годы)	Начало периода (год)
1 2 Среднее модели	24.79.3 148.011. 90.10.	32.32.4 124.11.	1921.54.7 1998.65.1

Рис. 3 Наложенные графики ряда уровня Каспийского моря (сплошная линия) и его полигармонической модели (пунктир), а также график прогноза изменения уровня Каспийского моря до 2027 года

Для уровня Каспийского моря составленный прогноз показывает, что его повышение продолжится до 2027 года до уровня начала XX века. Заметим, что исследовался ряд уровня Каспийского моря с 1891 по 1985 год. Таким образом, тот факт, что за с 1977 по 1995 год уровень Каспийского моря поднялся на 2.5 метра, является подтверждением сделанного прогноза.

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

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

1. Реализована автоматическая работа программы, что обеспечивает: воспроизводимость результатов исследований, позволяет исключить субъективное воздействие на результат вычислений, снижает требования к оператору программы, снижает до минимума время, затрачиваемое оператором на производство вычислений.

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

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

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

5. Длительность выделяемых периодов находится в пределах от 2 интервалов дискретизации до нескольких длин исследуемой реализации. В случае наличия априорных данных о присутствии в исследуемом временном ряде составляющих с периодом длительностью менее двух интервалов дискретизации, программа позволяет производить расчёт их параметров.

6. Программа обеспечивает высокое (настраиваемое) разрешение (<2%) выделяемых периодов, что необходимо в некоторых видах исследований.

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

8. Учитывается интерференция гармонических составляющих.

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

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

Рис. 4 Наложенные графики ряда яркости квазара NCG 4151 с неравномерной частотой дискретизации и его полигармонической модели (пунктир)

В третьей главе получены результаты моделирования ряда солнечной активности. Основным отличием рассчитанной модели от известных является наличие в ней 1800-летней гармонической составляющей, объясняющей современное глобальное потепление климата Земли с позиций изменения солнечной активности. Приведены геоинформационные подтверждения прогноза изменчивости солнечной активности на протяжении последних 2500 лет на примере климата Европы.

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

Рис. 5 Периодограмма ряда солнечной активности

Таблица 2

Параметры полигармонической модели ряда солнечной активности

Номер гармоники	Амплитуда (числа Вольфа)	Период (годы)	Начало периода (год)
1 2 3 4 Среднее модели	11322 165 215 275 15030	1800400 985 9.990.01 11.050.01	221090 17406 1706.30.5 1702.20.5

Определение доверительных интервалов параметров полигармонической модели ряда солнечной активности производилось методом Монте-Карло.



Рис. 6 Наложенные графики ряда солнечной активности (сплошная линия) и его полигармонической модели (пунктир)

С целью проверки работоспособности разработанного метода и адекватности параметров полигармонической модели ряда солнечной активности на основании данных табл. 2 был сгенерирован идеальный полигармонический ряд. Периодограмма этого ряда приведена на рис. 7; наложенные графики ряда и его модели - на рис. 8; результаты моделирования - в табл. 3.

Рис. 7 Периодограмма идеального полигармонического ряда

Рис. 8 Наложенные графики идеального полигармонического ряда и его модели

Таблица 3

Параметры полигармонической модели идеального полигармонического ряда

Номер гармоники	Параметры гармоник сгенерированного ряда	Рассчитанные параметры полигармонической модели
1	Амплитуда = 113. Период = 1800. Фаза = 4.4	Амплитуда = 112.77 Период = 1798.14 Фаза = 4.400
2	Амплитуда = 16. Период = 98. Фаза = -2.5	Амплитуда = 16.00 Период = 98.00 Фаза = -2.500
3	Амплитуда = 21. Период = 10. Фаза = 2.3	Амплитуда = 21.00 Период = 10.00 Фаза = 2.300
4	Амплитуда = 27. Период = 11. Фаза = -1.2	Амплитуда = 27.00 Период = 11.00 Фаза = -1.200
Среднее	153.	152.77

Анализ результатов моделирования солнечной активности показывает, что первые две выделенные гармоники относятся к общеизвестному 11-летнему циклу солнечной активности и указывают на то, что он представляет собой биение с периодами в 9.99 и 11.05 лет. Полученные результаты согласуются с данными Тернера и Шустера. В последнее время на двойственную природу 11-летнего цикла указывал, в частности, Сачок Г.И.

Более высокочастотные циклы (с периодом менее 9 лет), а также обсуждаемый географами 22-летний цикл, в моделируемом ряде солнечной активности выявлены не были. Исследованный Э. Брюкнером 33-летний климатический цикл в исследуемом ряде также выявлен не был. Выявленный вековой (98-летний) цикл солнечной активности так же хорошо известен географам (Витинский Ю.И., 1986).

Наиболее важным с прогностической точки зрения является 1800-летний цикл (Эйгенсон М.С., 1963), характеризующийся амплитудой, превышающей суммарную амплитуду всех остальных гармоник.

На 1800-летний цикл солнечной активности впервые было указано Шнитниковым А.В. (1949). Он провёл сопоставление различных природных ритмов Земли и обнаружил их привязанность к астрономическому явлению, называемому Большой сарес, когда Солнце, Земля и Луна оказываются на одной прямой. При этом Земля расположена на наименьшем удалении и от светила, и от спутника. В этом случае достигают набольшего значения приливные силы. Большой сарес повторяется через 1800 лет (с отклонениями) и сопровождается расширением земного шара в экваториальной полосе за счёт приливной волны, в которой участвуют Мировой океан и земная кора. Вследствие этого происходит изменение момента инерции планеты, и она замедляет своё вращение. Изменяется также положение границы полярного ледового покрова, происходит подъём уровня океана. Большой сарес отражается на климате Земли - по-иному начинают чередоваться засушливые и влажные периоды.

Выводы Шнитникова А.В. были подтверждены как в Институте физики Земли АН СССР, так и самыми разнообразными исследованиями в области географических наук, некоторые результаты которых сведены в табл. 4.

Таблица 4

Источники, указывающие на существование и проявление 1800-летнего цикла солнечной активности

Источник	Длительность цикла, год	Климатические показатели
Петтерсон, 1914	1800	Колебания климатических явлений (V в. до н.э. - XV в. н.э.)
Клаф, Брукс, 1933	1400	Атмосферные осадки
Антевс, 1938	1700	Движение ледников Новой Англии (позднеледниковая эпоха).
Джиллет, 1938 1949	1700 1800	Гидрологические и геологические явления. Отложения осадочных пород.
Предтеченский, 1945	1500	---
Шнитников А.В., 1949 1951, 1957, 1970, 1976	1800-1900	Колебания общей увлажнённости материков (XL в. до н.э. - современность) Слоистость илистых отложений (XXIII в. до н.э. - современность) Уровень Аральского моря.
Эйгенсон М.С., 1963	1800-1900
Карлстрем Т., 1966	1700
Максимов Е.В.	1850	Формирование забронированных глетчеров.
Дроздов О.А., 1976	1800	Колебания климата.
Беркович К.М., Борсук О.А.	1850	Изменение состава аллювия.
Некрасов И.А., Шейнкман В.С.	1850	Создание специфических форм гляциального рельефа, группирующихся в стадиальные комплексы.
Шилик К.К.	2000	- Продолжительность времени между максимумами новочерноморской и нимфейской трансгрессий. - Продолжительность времени между максимумами фанагорийской и корсуньской регрессий.
Алешинская З.В., Гурова В.С.	1850	Изменчивость растительного покрова вокруг озера Неро.
Никитин А.Л.	1850-1900	Комплексное изучение геоморфологических, палеогеографических и археологических данных района Плещеева озера.
Макеев В.М., Бердовская Г.Н.	1600-1800	Слоистость в озёрных осадках полуострова Ямал.
Кюнтцель В.В.	1850	Оползневая активность на европейской территории СССР.
Левковская Г.М.	2000	Радиоуглеродное датирование древесин, обнаруженных в Арктике.
Кащук А.С.	1800	Период обращения магнитного поля вокруг земной оси.

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

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

Рис. 9 Прогноз 1800-летней составляющей солнечной активности с 700 по 2200 год

К 793 году, когда гармоника находилась на своём максимуме, а температура воздуха на Земле также достигла своего максимума, относится первое упоминание о нападении викингов на северо-восточное побережье Британии. В конце IX века началось заселение викингами острова Исландия 8. Здесь и далее в рамках указано обозначение объекта на рис. 10. Примерно через столетие Эйрик Рыжий, бежал из Исландии и добрался до земли, которую назвал Гренландией 1 (англ. "зелёная земля"), покрытой зелёными лугами. Эйрик основывает колонию викингов в Гренландии. Около 1000 года сын Эйрика Рыжего Лейф Счастливый доплыл до новых земель на западе от Гренландии (видимо, побережье современной Канады), но не смог закрепиться. При этом напомним, что в новое время, несмотря на значительный технический прогресс, исследование северного канадского архипелага представилось возможным только в XIX веке.

Приведённый пример является показательным с точки зрения влияния климата на развитие человеческой цивилизации. За достаточно короткий промежуток времени климат изменился настолько, что цветущий остров стал непригоден для проживания.

Далее идёт спад значений солнечной активности и вызванное им похолодание. Подтверждение этому принесли исследования 29 человеческих зубов, раскопанных на юге Гренландии в могилах с датами жизни умерших. По соотношению изотопов кислорода в эмали зубов удалось вычислить среднюю температуру на острове и проследить, как она менялась с 1100 по 1450 годы. Оказалось, что за этот период среднегодовая температура упала примерно на 1.5^oC. Это похолодание не только погубило поля, доставлявшие питание колонистам, но и сделало остров труднодоступным, окружив его айсбергами и сплошными паковыми льдами. Когда в конце XV века корабли викингов снова смогли пробиться к Гренландии, колония оказалась целиком вымершей (см. рис. 9).

В XVII-XVIII веках значения солнечной активности достигают своего минимума, наступает так называемая эпоха маундера, или малого ледникового периода. В XVII веке полностью замерзали каналы в Голландии, устанавливалась устойчивая ледовая переправа между Данией и Швецией 15. До середины XIX века Темза 16 зимой покрывалась льдом. Далее значения солнечной активности начинают возрастать, и происходят следующие изменения климата в западной Атлантике: каналы в Дании и Голландии больше не замерзают, Темза - тоже. С 1976 по 1989 год толщина прибрежного льда на севере Гренландии сократилась с 6.5 до 4.5 м.

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



Рис. 10 Пространственное расположение исследуемых временных геофизических процессов на территории Европы

Таблица 5

Геоинформационные подтверждения прогноза изменчивости солнечной активности на основании 1800-летнего цикла

Датировка	События
500 год до н. э. (197)	Начало интенсивного похолодания по схеме расчленения голоцена М.И. Нейштадта.
340 год до н. э. (132)	Вероятно, последнее путешествие в античные времена к Гренландии (остров Туле) 1.
последние 500 лет до н. э. (197 41)	Отмечается особенная суровость климата Скандинавии 2.
8-10 гг. (41)	Замерзал Дунай 18.
Начало н. э. ( 42)	Высокий уровень Каспийского моря 3.
II - III века (42 86)	"В середине II века путь циклонов сдвинулся в лесную зону. В III веке засуха усилилась, количество осадков упало до 100-200 мм в год. Прокормиться больше не удавалось, и некоторые обитатели Центральной Азии вновь, как и две тысячи лет назад, двинулись на запад".
500 год (163)	Потепление климата. "Арктические льды вступили в стадию полуустойчивого существования" (Брукс) 4.
VII-XI века (202 266)	Уровень Каспийского моря на 2-4 м ниже уровня 1975г. (68) 3
VIII век (235 257)	Максимальное распространение дуба в Усманском бору (Воронежская область) 20.
VIII-XII века (235 209)	- Повышенное количество засух в Восточной Европе при повышенном уровне солнечной активности. Полярные сияния нередко наблюдаются жителями Среднего Поволжья 5. - Потепление в горах Европы 6. - Земледелие на Кавказе было возможно на 200-300 м выше современного уровня 7. - Отсутствие морских льдов в районе Исландии 8 и южной Гренландии 1. - Увлажнение было достаточным, чтобы занять земледельческой культурой почти всё междуречье Аму-дарьи и Сыр-Дарьи 17.
X век (266 260)	Возделывание пшеницы в Волжской Болгарии 9.
940±120 гг. (266)	На острове Виктория (72) 19, находящемся между Землёй Франца-Иосифа и Шпицбергеном, из-под отступившего края ледника вытаял ствол дерева, возраст которого (см. первый столбец) был определён радиоуглеродным методом.
X-XI века (266 240)	Поселения норманнов в Северной Америке.
XII-XIII века (240 170)	Похолодание, вызвавшее снижение границы лесов в Хибинах и на Кавказе 7. Сокращение возобновления деревьев широколиственных пород на Карельском перешейке 10.
XIII-XIV века (209 130)	Похолодание в Гренландии 1.
XIV век (170 130)	Европейский климат соответствовал современному.
XVI- XVIII века (92 40)	- Полуостров Канин на всех картах изображался островом 11. - Пшеницу на территории Волжской Болгарии не сеяли из-за суровости климата 9. - Разрастание горных ледников Скандинавии 2, Альп 6, Исландии 8, Аляски
XVII век (62 44)	Период наибольшего продвижения ледников.
конец XVIII - начало XIX в. (40 42)	Максимум ледовитости в Исландии 8.
Начало XIX в. (42)	Высокий уровень Каспийского моря 3.
1840-1940 гг. (43 59)	- Граница земледелия в Канаде переместилась на 100-200 км на север. - Увеличение количества перелётных птиц в Исландии 8 и Гренландии 1. - В Скандинавии 2 начали культивировать овощи, ранее недоступные для климатических условий этих стран. -Повышение границы распространения бука в горах Италии 12. - Хорёк и заяц в Европе распространились на 600 км в сторону полюса. - Увеличение срока навигации: на Западной Двине 13 на 17 дней, на Неве 14 - на три недели. - Лена близ Якутска стала вскрываться на 4 дня раньше, чем обычно. - Площадь оледенения Альп с 1890 по 1940 гг. (50 59) уменьшилась на 25% 6. - В Гренландии ледяной покров отступил настолько, что обнажил земли, которые в XII веке (226) были заняты поселениями и могильниками норманнских колонистов 1. - В Исландии освободились от надвига льдов земли, которые возделывались 600 лет назад (166). Та же картина наблюдается и на Скандинавском полуострове 8. - Повсеместно отмечено повышение снеговой линии. Так, в горах северной части Перу она поднялась на 900 метров.

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

Согласно выдвинутой гипотезе, потепление продолжится ещё приблизительно 600 лет, что может существенно сказаться на сложившейся экономической и социальной структуре многих регионов. Направленность капиталовложений, призванных смягчить и предупредить негативные последствия изменений климата будет определяться исходя из того, какая гипотеза о причинах всемирного потепления климата будет принята. В качестве альтернативных можно рассматривать выдвинутую гипотезу, указывающую на природную обусловленность происходящих изменений, и "парниковый эффект", указывающий на их антропогенный характер. Если принимается последняя гипотеза, то средства должны быть направлены на уменьшение выброса в атмосферу углекислого газа. Если же рассматривается предложенная гипотеза, которая предполагает длительное природно обусловленное потепление климата, то необходима соответственная перестройка инфраструктуры регионов в масштабе всего государства.

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

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

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

Математическое определение периода :

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

Исследования проводились на примере ряда солнечной активности (в числах Вольфа) (см. рис. 6). Для каждого из выделенных 26-ти периодов эмпирическим путём была выбрана зависимость амплитуды от времени, описывающая асимметричность исследуемых периодов (нарастание значений идёт быстрее, чем их убывание):

где a = const, b - параметр асимметричности амплитуды. В отношении гласных звуков было установлено, что величина параметра b коррелирует с такими перцептивными (слуховыми) характеристиками гласных как сила и ударение.

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

,

где y(ti) - отсчёты исследуемого ряда; myj - среднее значение модели j-го периода; aj=const - постоянная составляющая амплитуды; bj - параметр асимметричности амплитуды; j - продолжительность моделируемого периода; i0 - первый отсчёт периода.

Расчёты производились с помощью метода минимизации остаточной дисперсии, разработанного в главе 2.

Близость параметров b, рассчитанных для каждого из выделенных периодов позволила допустить постоянство этой характеристики на протяжении исследуемой реализации: b = 0.4.

Установить зависимость изменения периодов от времени не представилось возможным. Таким образом, было установлено, что в принятом описании исследуемого временного ряда от периода к периоду изменяются два параметра: a_j и m_yj.

Зависимости этих характеристик от времени были описаны полигармоническими моделями, рассчитанными методом, разработанным в главе 2. Результаты моделирования представлены ряда a(t) приведены в табл. 6.

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

- Биение со значениями периодов близкими к определённым для ряда солнечной активности (9.99 и 11.05 лет): 10.2 и 11.1 года. Заметим, что подобное биение было выявлено и в результате моделирования ряда экстремумов солнечной активности: 10 и 10.9 лет.

- Вековую составляющую с периодом 98 лет, совпадающим с аналогичным периодом, полученным для ряда солнечной активности (98 лет).

- Составляющую с периодом 157 лет. Заметим, что значение этой составляющей, приведённое в табл. 6, получено на третьем этапе - расчёте параметров итоговой полигармонической модели (см. главу 2). По окончании же второго этапа - построения периодограммы - период этой составляющей был равен 2171 году.

Таблица 6

Результаты моделирования ряда a(t)

Номер гармоники	Амплитуда (числа Вольфа)	Период (год)	Начало периода (год)
1 2 3 4 Среднее модели	8. 8. 9. 17. 28.	157. 98. 10.2 11.1	1825. 1720. 1713. 1711.

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

Аналогичные результаты были получены для ряда m_y(t).

На основании полученных результатов представилось возможным осуществить прогноз амплитуды и среднего значения за период на ближайшие три периода. Результаты приведены в табл. 7.

Таблица 7

Прогноз значений амплитуды a(t) и среднего m_y(t)

j (№ периода)	27	28	29
start - end	1988 - 1998	1998 - 2009	2009-2019
a(t)	17.	17.	15.
my(t)	32.	32.	27.

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

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

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