WWW.SELUK.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

 

Pages:   || 2 | 3 | 4 | 5 |   ...   | 7 |
-- [ Страница 1 ] --

П.Ф. Демченко, А.В. Кислов

СТОХАСТИЧЕСКАЯ ДИНАМИКА

ПРИРОДНЫХ ОБЪЕКТОВ

Броуновское движение

и геофизические

приложения

Москва

ГЕОС

2010

УДК 519.2

ББК 22.171

Д 12

Демченко П.Ф., Кислов А.В. Стохастическая динамика природных объектов.

Броуновское движение и геофизические примеры – М.: ГЕОС, 2010. – 190 с.

ISBN 978-5-89118-533-3

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

и глобальная климатическая система в случае изучения вариаций климата;

это и озера – при анализе динамики уровней воды, и ледники – при исследовании вариаций их размеров;

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

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

Издание осуществляется при финансовой поддержке Российского фонда фундаментальных исследований по проекту № 10-05-07055-д © П.Ф. Демченко, А.В. Кислов © ГЕОС

ВВЕДЕНИЕ

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

Декларированный единый подход к анализу хаотической динамики за ключается в том, что в каждом случае авторы используют один и тот же ап парат стохастических дифференциальных уравнений и неравновесной стати стической механики. Возможность его применения базируется на понимании присущего природным процессам дуализма: с одной стороны, их динамика хаотична, с другой стороны, они подчиняются, как правило, детерминистиче ски сформулированным физическим законам. Как объединить случайное по ведение и детерминированное описание? Один из возможных путей решения этой проблемы заключается в том, чтобы определять параметры статистиче ских моделей случайных процессов исходя из «первых принципов» (детер минированных законов), то есть использовать фундаментальные физические законы сохранения массы, энергии и импульса для расчета вероятностных характеристик поведения природных процессов. Возможно, именно в таком ключе можно интерпретировать известную фразу Альберта Эйнштейна о том, что «Бог не играет в кости» (хотя она была высказана по другому поводу – в дискуссии об основах квантовой механики).

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

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

Временные масштабы, на которых еще не сказываются обратные связи, отно сят к спектральному интервалу, где энергия увеличивается с ростом частоты (соответствующий случайный процесс называют в этом случае «красным шумом», или, более конкретно, «красным шумом» называют ситуацию, когда энергия растет обратно пропорционально квадрату частоты). На тех времен ных масштабах, на которых активизируются стабилизирующие процессы, дальнейший рост энергии флуктуаций прекращается. Мощность флуктуаций выходит на плато (не зависит от частоты), которое называется «белый шум».





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

Совсем не обязательно считать, что широкий диапазон масштабов (вре менных и (или) пространственных) изменчивости какой-либо физической ха рактеристики (температуры, скорости и т.д.) связан с наличием иерархически согласованного набора объектов со своими характерными временами реакции на внешние воздействия, которые дают вклад в результирующий спектр мощности флуктуаций. Широкий диапазон масштабов присущ движению сплошной среды. Так, в случае турбулентных флуктуаций скорости в несжи маемой жидкости существуют два пространственных масштаба. Это, во первых, внешний масштаб L, определяемый размером области течения жид кости. Во-вторых, l – масштаб вихрей, на которые начинают оказывать влияние силы молекулярной вязкости: l ~ v ( – скорость диссипа ции кинетической энергии, – кинематическая вязкость). Для локально одно родной изотропной турбулентности распределение энергии флуктуаций – спектр S – по волновым числам k в инерционном интервале L k l определяется процессом передачи кинетической энергии от больших мас штабов к меньшим (каскад Ричардсона). В этом интервале не содержится бо лее никаких масштабов, связанных с определенными физическими объекта ми, а спектр подчиняется закону Колмогорова–Обухова: S ~ k [Монин, Яглом, 1967]. Еще одним из таких примеров является фликкер – шум с час тотным спектром S ~ ( – круговая частота), широко распространенный в технических устройствах, в частности в электронных лампах [Рытов, 1976;

Климонтович, 1982]. Для описания таких флуктуаций широко применяется ряд математических методов, связанных с теорией броуновского движения (см. далее). Например, спектр КолмогороваОбухова можно получить для движения жидкой частицы под действием белого шума [Голицын, 2004].

Для описания флуктуаций природных объектов будут рассматриваться не все существующие стохастические методы, а только связанные с применени ем теории броуновского движения. Термин броуновское движение обязан своим происхождением шотландскому ботанику Роберту Броуну, который в 1887 г. исследовал беспорядочное движение взвешенной в воде пыльцы рас тений, причем длительное время он посвятил изучению вопроса, не является ли это движение следствием проявления жизнедеятельности. В 1905 г. Аль берт Эйнштейн дал теорию этого процесса и фактически заложил основы анализа флуктуаций методами стохастических дифференциальных уравне ний. Основой для применения концепции броуновского движения к природ ным объектам является возможность разделения совокупности флуктуаций их динамики на «быстрые» и «медленные», согласно принятой в статистиче ской физике терминологии. Важно, что на временах реакции медленных пе ременных на внешнее воздействие быстрые переменные теряют память об их предыдущем состоянии и могут рассматриваться как случайный процесс с заданной статистикой. Предельным случаем такой ситуации с разделением времен является трактовка воздействия быстрых переменных на медленные как воздействие белого шума – случайного процесса с исчезающе малым временем корреляции, так называемого дельта-коррелированного случайного процесса [Гардинер, 1986;

Кляцкин, 1980]. Впервые эту идею для описания флуктуаций медленных объектов окружающей среды высказал Клаус Хас сельманн [Hasselmann, 1976]. Реализация этой идеи, базирующаяся на об ширном эмпирическом материале (раздел 1.1), позволяет применять к расче ту статистических характеристик интересующих нас переменных развитый аппарат стохастических дифференциальных уравнений.

Простейшим (но широко встречающимся в приложениях) примером таких уравнений является линейное уравнение Ланжевена, анализ которого дан в разделе 1.2. Оно содержит все необходимые составляющие для расчета ста тистических характеристик флуктуаций медленной переменной, связывая ее изменение во времени с дельта-коррелированным (или асимптотически при ближающейся к ней) случайным воздействием и систематическим вкладом стабилизирующей обратной связи. Результирующий случайный процесс но сит название процесса Орнштейна–Уленбека. В оригинальной трактовке Альберта Эйнштейна случайная сила соответствует действию частых соуда рений тяжелой частицы с многочисленными окружающими ее более легкими молекулами. Взаимодействие с этим «облаком легких частиц» оказывает на поведение тяжелой частицы двоякое влияние. С одной стороны, эти частые некоррелированные соударения вызывают случайные блуждания тяжелой частицы – броуновское движение. С другой стороны, именно из-за взаимо действия с облаком легких молекул возникает макроскопическое трение – стабилизирующая обратная связь. У природных объектов иной природы, ко торые рассматриваются в данной книге, есть свои стабилизирующие обрат ные связи. Расчет коэффициентов этих обратных связей и параметров слу чайных воздействий для различных природных объектов составляет значи тельную часть главы 3.

Современные методы неравновесной статистической механики позволяют выводить уравнения Ланжевена (в общем случае – нелинейные) для медлен ных переменных, исходя из уравнений для полной системы. Для этого, в ча стности, используется техника проекционных операторов [Mori et al., 1980].

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

скорости (зависит от текущих значений медленных переменных), интеграла памяти (зависит от изменений медленных переменных в прошлом) и корот копериодной составляющей – «случайной силы». Вывод обобщенных урав нений Ланжевена с изложением техники проекционных операторов дан в разделе 2.2.

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

В некоторых случаях, наоборот, важно описать качественные изменения в поведении коэффициентов обратных связей. В этом случае линейная теория не может быть использована. В качестве примера можно привести пересы хающие озера, когда система находится у порога применимости процедуры линеаризации, принятой для описания динамики уровня воды озера, напол ненного водой. К аналогичным примерам можно отнести динамику влажно сти почвы в режиме избыточного увлажнения, когда важную роль играют процессы образования стока. В этом случае теория броуновского движения позволяет находить нелинейные уравнения для плотности вероятностей нере гулярных процессов в природных объектах. В разделе 2.3, следуя работам Кляцкина [1980, 2002], приводится один из методов вывода такого уравнения (уравнения Фоккера–Планка) из уравнения Ланжевена методом вариацион ных производных. Приведенный пример расчета времени корреляции флук туаций влагозапаса почвы по теории броуновского движения во всем диапа зоне изменения внешних параметров (включая смену режимов увлажнения) при сравнении с результатами моделирования методом Монте-Карло показы вает эффективность теории.

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

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

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

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

ГЛАВА 1. ХАОТИЧНОСТЬ КАК ТИПИЧНОЕ ЯВЛЕНИЕ

ДИНАМИКИ ОКРУЖАЮЩЕЙ СРЕДЫ

1.1. Эмпирические данные о динамике природных процессов Времення динамика природных процессов крайне редко представляет собой регулярный колебательный процесс. Обычно (но не всегда) такое по ведение является признаком детерминированности влияющих факторов и линейности отклика системы на внешнее воздействие. Даже суточный ход температуры вопреки распространенному мнению далеко не всегда представляет собой кривую, близкую к синусоиде, поскольку на формирование термического режима оказывает воздействие не только закономерный (астрономический) ход высоты Солнца над горизонтом в течение суток, но также неравномерная закрытость горизонта (в горных местностях), нерегулярное воздействие облачности, трансформирующее радиационный баланс, и адвективные изменения температуры за счет динамики воздушных масс разного происхождения. На рис. 1.1.1 показан пример суточных ходов температуры и модуля скорости ветра (задаваемого изменением термического режима, в том числе и суточными изменениями стратификации атмосферы). Из рис. 1.1.1 следует, что амплитуда меняется от суток к суткам, а в некоторых случаях (например, 11.06.2008) вообще нельзя говорить о существовании закономерных суточных изменений.

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

Хорошим примером процессов, близких к периодическим, являются при ливные явления в океанах и морях.

Рис. 1.1.1. Изменение температуры (t) и модуля скорости ветра (U) в приповерх ностном слое (07.06–16.06.2008, Хибины, Кольский полуостров) (по данным измере ний на Хибинской учебно-научной станции географического факультета МГУ) Имеется, наконец, совсем немного явлений, у которых цикличность или, скорее, ритмичность (поскольку «цикличность» все-таки подразумевает гар монический вид динамики) определяется во многом внутренними свойствами системы. Так, регулярный характер имеет так называемая квази-двухлетняя цикличность – удивительно закономерное чередование западных и восточ ных потоков воздуха экваториальной стратосферы с периодом, близким к месяцам [Холтон, 1979]. Однако данная периодичность не является супер гармоникой годового хода – это случайное совпадение. Возникновение ква зидвухлетней цикличности обусловлено нелинейным взаимодействием вер тикально распространяющихся экваториальных волн с стратосферной циркуляцией [Кулямин и др., 2008].

Ритмический характер имеют перестройки общей циркуляции атмосферы, известные как Южное колебание и Северо-Атлантическое колебание – NAO (North Atlantic Oscillation). Последнее отражает согласованную динамику Ис ландского минимума и Азорского максимума атмосферного давления над Атлантическим океаном (рис. 1.1.2), или, более общо, является наиболее от четливо регионально выраженной частью так называемого Арктического ко лебания – последовательного согласованного усиления и ослабления межши ротного градиента давления. Как и другие механизмы дальних связей, NAO Рис. 1.1.2. Межгодовые изменения индекса NAO (North Atlantic Oscillation) наиболее отчетливо выражено при использовании осредненных данных (средних за месяц или сезон). Для диагностирования NAO используют специ альный индекс который рассчитывается как разность между нормированны ми значениями аномалий давления, измеренного на станциях, характеризую щих условия Азорского антициклона (р ) и Исладской депрессии (р ):

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

NAO хорошо проявляется в холодное время года. На рис. 1.1.2 представле ны межгодовые изменения. Видно, в частности, что отрицательные значения, определяющие ослабление вторжений атлантического воздуха на Европу, ха рактерны для 50-х и 60-х годов ХХ столетия, затем наступил период господ ства положительной фазы NAO, продолжавшийся примерно до середины х годов, сменившийся затем этапом перемежающихся положительных и отри цательных значений I, продолжающимся по настоящее время.

Южное колебание представляет собой последовательность переходов состоя ния атмосферы и тропического Тихого океана из стадии Эль-Ниньо в стадию Ла Нинья. На рис. 1.1.3 представлена динамика индекса Южного колебания (South ern Oscillation Index), рассчитываемого как разность нормализованных аномалий атмосферного давления по данным измерений двух станций (обычно это Дарвин (Австралия) и Сант-Яго (Чили)). Однако далеко не все диагностируемые год за годом аномалии настолько отчетливы во всей совокупности признаков, чтобы быть четко отнесены к какому-то определенному классу.

Рис. 1.1.3. Межгодовые изменения индекса Южного колебания SOI (South Oscilla tion Index) Рис. 1.1.4. Динамика климата за последний миллион лет по данным содержания тяжелого изотопа кислорода в отложениях приповерхностного планктона (по мате риалам [Imbrie et al., 1984]) глубоководного бурения донных отложений океана, скважина ODP 677, Рифт Коста-Рика) Более длительные «декадные» колебания характеризуются специальными индексами PDO (Pacific Decadal Oscillation) в Тихоокеанском регионе и AMO (Atlantic Multidecadal Oscillation) в Атлантическом океане.

Переходя к «сверхдлинным» процессом, выделим хорошо выраженный ~100 000-летний цикл изменений климата, который четко виден на рис. 1.1.4.

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

Методы спектрального анализа позволяют определять спектр мощности слу чайного процесса – средний квадрат Фурье-амплитуды разложения процесса по колебаниям различной частоты. Так получается информация о ритмах, проявляющихся «в среднем». Таков, например, спектр колебаний интенсив ности Южного колебания у которого преобладает 3,7 и 4,9-летняя периодич ность (рис. 1.1.5).

Наконец, огромное множество процессов представляет собой чередование экстремумов разных знаков, амплитуд и периодичностей, которые появляются, Рис. 1.1.5. Спектральная плотность вариаций индекса Южного колебания (по строен по данным колебания индекса SOI – см. рис. 1.1.3) исчезают, образуют повторяющиеся группы или следуют изолированно. Не ставя немыслимую задачу перечисления всех квази-ритмических явлений, от метим, например, циклы Дансгора–Оешгера (продолжительностью в несколько сотен лет), четко зафиксированные не только в Гренландии и море Ирмингера (Лабрадорской котловине) (рис. 1.1.6), но проявляющиеся практически повсе местно в мире. Рис. 1.1.6 свидетельствует, что на протяжении холодной позд неплейстоценовой эпохи климатический режим Северо-Атлантического регио на неоднократно испытывал быстрые изменения, во время которых температу ра поднималась практически до современного теплого (последние 10 тыс. лет) уровня, а затем столь же быстро опускалась. Эта последовательность аномалий иногда повторялась несколько раз подряд (стадии 11, 10, 9 и 7, 6, 5), а иногда имели место единичные события. События Хайнриха (выброс необычно боль шого количества айсбергов из ледниковых щитов (Гренландского, Скандинав ского, Лаврентийского, Исландского), диагностируемый по резкому росту в донных отложениях продуктов абразионной деятельности ледников на суше) также происходили нерегулярно во времени.

Рис. 1.1.6. Колебания климата в Северной Атлантике и Гренландии а – динамика изотопа тяжелого кислорода по данным бурения ледникового щита Гренландии (проект GRIP);

б – реконструированные вариации температуры воды на поверхности моря Ир мингера (Лабрадорская котловина);

в – реконструированные вариации числа частиц (размером более 150 мкм) в грамме донных отложений (103 1/г). Цифры – номера теплых стадий Дансгора– Оешгера. Заштрихованы интервалы, отвечающие времени наблюдавшихся событий Хайнриха Спектральный анализ различных временных рядов показывает, что пики на кривых функций спектральной плотности лишь изредка обладают ста тистической значимостью. Если собрать эти случаи, то оказывается, что фак тически можно обнаружить колебания с любыми периодами. Эту точку зре ния подтверждают результаты, представленные в табл. 1.1.1, в которой даны многочисленные примеры обнаруживаемых в динамике природных процес сах периодичностей (для определенности авторы ограничились только меж годовыми вариациями, не превышающими первые сотни лет).

Если пики не достигают уровня значимости (что типично имеет место при обработке данных наблюдений или реконструкций), то отвечающие им флук Таблица 1.1.1. Периоды, соответствующие статистически значимым колебаниям, обнаруженным в спектрах различных индикаторов в разных регионах северного по лушария (по [Stocker, Mysak, 1992;

Добровольский, 2002], с добавлениями авторов) Примечание. Цифры в скобках – значение периодичности колебаний. Буквами обозначен вид данных, использованных для диагноза флуктуаций: B – биологические данные, C – по 14С, G – гляциологические, H – исторические, Hy – гидрологические, I – инструментальные, O – по 18О, T – дендроклиматические, TSS – глобально осредненные данные по температуре морской воды на поверхности.

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

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

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

Плодотворным методом является подбор для описания временных серий стохастических моделей. Среди них наиболее эффективны, по-видимому, так называемые модели авторегрессии. Существуют различного уровня теорети ческие обоснования этого подхода (например, основанные на принятии принципа наибольшей энтропии [Dobrovolski, 1992]), однако наиболее важна, вероятно, эмпирически доказанная эффективность использования моделей такого рода в конкретных приложениях. В рамках данного подхода предпо лагается, что центрированный случайный процесс аппроксимируется рядом где c коэффициенты авторегрессии, индекс «i» обозначает дискретный момент времени, – последовательность нормально распределенных некор релированных величин.

Практическое применение ряда (1.1.1) для аппроксимации случайных процессов различной природы показывает, что крайне редко для описания временной динамики природных процессов требуется M 1 [Добровольский, 2002], поскольку с практической точки зрения одного слагаемого, как прави ло, достаточно для воспроизведения процесса с возможной для конкретного случая точностью. Это существенно упрощает ситуацию, позволяя использо вать широко известные модели теории случайных процессов. В самом деле, в том случае, когда М = 1, модель авторегрессии первого порядка принимает вид марковского процесса первого порядка Его важной разновидностью является случай так называемого винеровско го процесса при c = 1. В рамках модели предусмотрен и случай M = 0, когда =, то есть моделью временной динамики служит белый шум.

Эти принципиальные особенности поведения данного случайного процесса полезно рассмотреть с точки зрения спектрального анализа. Как известно [Мате матическая энциклопедия…, 1977;

Gilman et al., 1963], выражение для спектраль ной плотности процесса авторегрессии первого порядка описывается выражением в котором v – линейная частота. Причем, если c 0, то S(v) = const, соответ ствующий модели белого шума. Обратная зависимость (1.1.3) функции спек тральной плотности от частоты часто называют «спектром красного шума».

Рассмотрим важный случай, когда изучаются низкие по сравнению с час тотой Найквиста частоты. Тогда, раскладывая косинус в степенной ряд и ог раничиваясь первым членом разложения, получим выражение Если дополнительно ограничить частотный диапазон и снизу: v (1– c)/ c ), то спектральная плотность оказывается обратно пропорциональна квадрату частоты. В настоящей работе именно последнюю модель, с законом «2» будем называть красным шумом.

Закономерность S(v) ~ v можно рассматривать и с той точки зрения, что продолжительности аномалий определенного знака (понимаемых как поло вина «периода» колебаний ) и их величины («амплитуды» a ) подобны на разных масштабах. В самом деле, имея в виду, что спектр по определению S ~ a v, а v ~, получаем, что на любом временном масштабе a ~.

Отсюда можно сделать два вывода.

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

Отметим, что получившаяся зависимость a ~ v может быть интерпре тирована в терминах понятий, введенных в разделе 1.2, как закономерность, описывающая величину амплитуды графика функции на отрезке, равном мас штабу частоты (формула (1.2.19)). В этом случае показатель степени представ ляет собой индекс фрактальности и характеризует то, что у любого ряда, спектр которого «красный», он равен 1/2. Поскольку, как показано в разделе 1.2, индекс фрактальности связан с фрактальной размерностью, как = D –1, то получаем, что у авторегрессионного процесса первого порядка D = 3/2. Такая размерность есть типичная особенность броуновского движения.

Рассмотрим некоторые примеры, когда в спектрах, построенных по эмпи рическим рядам, проявляется закон красного шума. На рис. 1.1.7 показаны спектры колебаний некоторых палеоиндикаторов.

Несмотря на то, что эти ряды совершенно различны как по тем величинам, которые они отображают, так и по дискретности и продолжительности, их Рис. 1.1.7. Функции спектральной плотности (в логарифмических координатах) временных рядов, представленных на рис. 1.1.4 и рис. 1.1.6, а [Wunsch, 2003] По оси абсцисс отложены частоты (1/год). Пунктир – прямая уравнения регрессии, прове денная по методу наименьших квадратов. Коэффициент регрессии есть показатель степени для частоты колебаний (формула (1.1.4)), получился равным 2.3 и 1.8, соответственно Рис. 1.1.8. Функция спектральной плотности (в логарифмических координатах) колебаний температуры поверхности океана (средняя для умеренных широт Тихого океана) за 1950–1994 гг., построенная по четырем различным базам данных [Dommenget, Latif, 2002] Прямая линия соответствует закону красного шума спектры обладают сходным поведением – на частотах, существенно превы шающих частоту Найквиста спектр описывается функцией, практически сов падающей c ~. На более низких частотах диагностируется низкочастотное колебание (100 000-летняя периодичность) и спектр выходит на плато. Такие закономерности, конечно, проявляются в спектрах далеко не всегда. Здесь спе циально были подобраны случаи, когда применимость авторегрессионной мо дели первого порядка очевидна и возможно, по крайней мере в принципиаль ном плане, построение физических моделей явлений.

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

Следующий пример посвящен неравномерности вращения Земли. На рис.

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

Таким образом, налицо широкая применимость модели авторегрессии первого порядка. Как будет показано в разделе 1.2, существует естественное сходство данной модели с уравнением Ланжевена, что позволяет во всех слу чаях, когда реальный процесс описывается авторегрессионной моделью и Рис. 1.1.9. Спектр флуктуаций угловой скорости врашения Земли По оси абсцисс отложены периоды колебаний в сутках, по оси ординат – логарифм спек тральной плотности в относительных единицах (рисунок любезно предоставлен Н.С. Сидоренковым). Белая линия соответствует закону красного шума спектр является «красным», искать его теоретический аналог. В последующих разделах рассмотренные примеры (и другие задачи) обсуждаются именно с данной точки зрения – создания математической модели, опирающейся на «первые принципы» – фундаментальные физические законы сохранения мас сы, энергии и импульса, описывающей наблюдаемое поведение, соответст вующее проявлению броуновского движения в макромасштабных природных процессах.

1.2. Свойства решения уравнения Ланжевена, классическое В предыдущем разделе было показано, что, несмотря на сложность и мно гообразие природных явлений, в их временном поведении и пространствен ной структуре обнаруживаются идентичные закономерности. Это прежде всего авторегрессионное поведение временных рядов, наглядно отражающее ся в спектрах в виде закономерностей белого и красного шума. Данные веро ятностные особенности можно получить путем решения стохастического дифференциального уравнения Ланжевена. Оно может быть применено к различным процессам, однако первоначально было использовано для описа ния одномерной динамики броуновской частицы в «облаке» легких молекул (см. раздел 2.1). Данное уравнение имеет вид Оно описывает «медленную» динамику состояния W = W(t) инерционного объекта под влиянием быстро флуктуирующего внешнего воздействия [Ах манов и др., 1981;

Рытов, 1976]. Как будет продемонстрировано далее, такая математическая модель успешно применима для описания динамики состоя ния различных природных объектов и поведения различных процессов.

В этом уравнении описывает эффективность линейной обратной связи системы эта величина определяется при выводе уравнения (1.2.1) и зависит от конкретных особенностей рассматриваемого «медленного» процесса, оп ределяя характерное время его эволюции. Интенсивность внешнего воз действия = (t) считается очень быстро флуктуирующей случайной величи ной. Выражение «очень быстро» следует понимать с точки зрения заданного масштаба медленных изменений. Внешнее воздействие может создаваться единственным фактором, а может быть представлено суммарным действием нескольких некоррелированных воздействий причем дисперсия выражается как Уравнению (1.2.1) в зависимости от исследуемого объекта может прида ваться различный смысл. Так, если W обозначает координату лагранжевой частицы, то случайным является распределение скоростей. В других случаях, трактуя уравнение (1.2.1) как выражение второго закона Ньютона, имеем картину случайного распределения ускорений. В этом случае можно говорить о случайных блужданиях в пространстве импульсов [Голицын, 2004 а, б].

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

и получим где обозначено c = (1 + ). Выражение (1.2.5) эквивалентно (при соответ ствующих вероятностных определениях) модели марковского процесса, описываемого соотношением (1.1.2).

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

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

Перейдем к нахождению решения уравнения (1.2.1). Это стохастическое уравнение, определяющее реализацию случайного процесса W = W(t), несу щую, именно в силу случайности, мало полезной информации. Поэтому про цедура решения стохастического уравнения должна быть направлена на полу чение вероятностных характеристик описываемого уравнением случайного процесса. В случае уравнения Ланжевена решение этой задачи может быть по лучено следующим образом. Запишем вид точного решения уравнения (1.2.1) (как детерминированного обыкновенного дифференциального уравнения):

В разделе 3.5 рассматривается случай учета случайного характера началь ного состояния. Сейчас будем считать, что оно представляет собой детерми нированную функцию. Возводя обе части выражения (1.2.8) в квадрат, пред положив, что (t) = 0 и выполняя операцию усреднения, получим выраже ние для дисперсии в виде двойного интеграла Для его вычисления используем функцию корреляции вида (1.2.6), под ставим ее в формулу (1.2.9) и, используя свойство четности функции корре ляции по аргументу, получим Вычисление интегралов, с учетом условия (1.2.7), позволяет получить фор мулу, описывающую изменение дисперсии медленного процесса во времени флуктуаций растет пропорционально времени, Последний результат методологически очень важен, поскольку демонстриру ется возможность «саморазвития» природных объектов без какого-либо внешне го влияния, под действием внутреннего источника климатической системы.

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

Рассмотрим спектральное представление случайной функции W(t), описы ваемой уравнением (1.2.1). Ее функция спектральной плотности должна зави сеть от спектра случайной функции (t). Действительно, представляя W ( t ) = exp ( i t ). Подставляя в уравнение (1.2.1), для любого k получим Если считать, что Ф, Г есть коэффициенты Фурье в спектральном пред ставлении случайных функций W(t) и (t), то из-за условия ортогональности отличаться от нуля будут только такие усредненные произведения, у которых индексы одинаковы, то есть,. Поэтому, умножая левую и правую часть (1.2.14) на комплексно сопряженную величину, получим, характеризуют математическое ожидание энергии слу чайного колебания частоты, то есть описывают вклад данного колебания в общую дисперсию. Данная закономерность справедлива для любого номе ра, то есть для любой частоты. Обозначим S (), S () выражение, демонстрирующее взаимозависимость спектров выходного и входного сигналов.

Для разных моделей случайных процессов S () будет различно, однако в рамках данной книги особенный интерес вызывает ситуация, когда случай ный процесс, вызывающий эволюцию инерционной системы, является дель та-коррелированным. В этом случае прямой расчет показывает, что спектр представляет собой постоянную величину, которую обозначим как S () = S (так называемый «белый шум»), и выражение (1.2.16) принимает вид На сравнительно высоких частотах ( ) спектральная функция воз растает пропорционально. Случайный процесс, характеризуемый таким спектром, называется «красный шум» (см. формулу (1.1.7)). С уменьшением частоты скорость изменения постепенно уменьшается, и при ( ) спектр функции отклика имеет характер белого шума.

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

Условие линейного роста дисперсии со временем (выражение (1.2.12)) ти пично для броуновского движения. Дополнительным условием является тре бование гауссового характера приращений и их статистическая независи мость на непересекающихся интервалах. Предполагая эти условия выполнен ными, можно считать, что при малых временах вариации ведут себя как «классическое» броуновское движение, описываемое винеровской моделью dW / dt = (t). На больших временах, когда эффективно проявляется обратная связь, поведение решения меняется.

Рассмотрим особенности броуновской частицы (или динамику иного по казателя, описываемого уравнением (1.2.1)) с точки зрения фрактальной гео метрии. Определим необходимые понятия. Как известно, фракталом называ ется структура (например, кривая на плоскости), части которой подобны ей самой [Фракталы …, 1998]. Мерой количественного выражения подобия служит так называемая фрактальная размерность. Она вводится естествен ным путем, если возникает задача найти минимальное количество регуляр ных объектов, покрывающих какой-либо сложно устроенный объект. Так, для покрытия отрезка единичной длинны его маленькими копиями (размером ) их потребуется N штук, причем (1/) = N. Для покрытия единичного квадрата его копиями (с размером стороны, равным ) их требуется (1/) штук, а для куба – (1/) копий. В этих случаях показатель степени отражает размерность пространства. Обобщая данные результаты на случай сколь угодно сложной кривой или поверхности, получим определение размерности в виде Не пытаясь рассматривать проблему в общем виде, будем считать, что фракталы располагаются в евклидовом пространстве. В качестве множества, покрывающего фрактал, используется N клеток, размер стороны которых ра вен. Пользуясь этим определением, рассмотрим более сложный по сравне нию с предыдущими пример. Пусть имеется график функции y = y(t). Без по тери общности можно считать, что рассматривается единичный интервал оп ределения, разделенный на n частей, равных 1/t. Отградуируем аналогично и вертикальную ось, то есть выберем масштаб, равный 1/t. Площадь, распо ложенная над одним i-м подинтервалом с размером t, включающая фраг мент функции y = y(t), равняется N tt, где N – число квадратов, покры вающих в данном масштабе отрезок рассматриваемой кривой. С другой сто роны, учитывая, что приращение функции в пределах этого подинтервала равняется y, можно приближенно написать: y t N tt, причем ясно, что это выражение будет тем точнее, чем меньший масштаб t используется.

Отсюда получается, что N = y / t. Всего имеется n подинтервалов, так что полное число квадратов N = N = y ( t )2 (здесь y – среднее прираще ние за время t).

Дальнейшее продвижение в рамках общего подхода невозможно, поскольку динамика во времени произвольного случайного процесса может происходить по-разному. Для конкретизации результата необходимо принять какую-либо модель случайного процесса. Например, в случае случайного гауссового блуж дания (модель броуновского движения, см. выше), приращение на интервале t, в среднем пропорционально t. Поэтому получается, что N = 1/ (t), то есть размерность траектории броуновского движения равняется 3/2.

Наряду с понятием размерности геометрия фрактала может быть охарак теризована так называемым индексом фрактальности [Дубовиков, Старченко, 2003]. Не останавливаясь на его строгом определении (разность метрической и топологической размерности фрактала), отметим его очевидную полез ность, так как он характеризует размах вариаций фрактала. Действительно, рассмотрим поведение функции y = y(t) на определенном отрезке, который разделим на n частей, каждый имеет размер. Составим величину где Ai () – «размах», то есть разность между максимальным и минимальным значениями функции y(t) в пределах [t, t ]. Если при 0, то показатель степени () и называется индексом фрактальности.

Между D и существует однозначная связь. Чтобы ее установить, рас смотрим, как и ранее, наглядный случай клеточной размерности. Пусть опять, как было использовано выше, Ni – число квадратов (размером х ), покрывающих, в данном масштабе, рассматриваемую кривую на [t, t ]. То гда имеет место неравенство N – A () 0. Разделим каждое слагаемое на и просуммируем по всем i. Получим N() – V () 0. Здесь N() – число клеток, покрывающих график функции y(t). При переходе к пределу, исполь зуя определение индекса фрактальности, получаем N(). Следова тельно, D = + 1.

Введем в рассмотрение еще один полезный показатель фрактала Н, кото рый определим, как H = 2 – D. Использование этой величины позволяет ав томатически обобщить закономерности классического броуновского движе ния на случай так называемого фрактального броуновского движения, или фрактального гауссовского процесса – ФГП [Кроновер, 2000]. Для ФГП закон роста дисперсии пропорционально времени формулируется так:

Здесь 0 H 1. Размерность реализации ФГП равняется D = 2 – H, причем случай с H = 1/2 совпадает с классическим гауссовским случайным процес сом. Функция распределения вероятностей описывается выражением [Кроно вер, 2000] Рассмотренные зависимости естественно обобщаются на случай фрак тальной поверхности (двумерное поле) – в этом случае D = 3 – H. Отметим, что индекс H совпадает с так называемым показателем Херста, который был использован как эмпирический индекс, описывающий зависимость от време ни нормированного на стандартное отклонение размаха флуктуаций. Оказа лось, что эта величина пропорциональна длине интервала в степени Н.

Рассмотренные показатели полезны, в частности, потому, что делают воз можным определение «интегральных» особенностей поведения фрактала. Со поставление поведения временных рядов с различными значениями показателя Херста или индекса фрактальности показало, что при 0 H 1/2 ( 1/2, D 3/2) процесс ведет себя знакопеременно, длительные тренды отсутствуют, система возвращается к среднему значению. Ясно, что это может происходить в том слу чае, если в динамике системы активны отрицательные обратные связи, стабили зирующие ее поведение. При 1/2 H 1 ( 1/2, D 3/2) поведение процесса иное – для него характерен небольшой уровень высокочастотного шума в соче тании с сохранением тенденций. Анализ различных эмпирических рядов пока зал, что действительно могут наблюдаться различные ситуации.

Так, в работе Кузнецова [2006] проанализированы ряды среднемесячных значений приземной температуры воздуха за период 1899–2002 гг. в средних широтах Северного полушария (15–75° с.ш.). Предварительно был отфильт рован высокочастотный шум – это было достигнуто тем, что предназначен ные для анализа ряды «собирались» как сумма по трем главным компонентам при разложении рядов по естественным ортогональным векторам. Все значе ния показателя Херста попали в интервал 1/2 H 1, то есть межгодовые из менения температуры оказались трендоустойчивыми с относительно низким уровнем шумов.

Естественно, что полученные теоретические результаты не будут меняться, если в качестве аргумента случайного процесса рассматривать не время, а про странственную координату. Так, авторы исследовали поля различных метеоро логических величин у земной поверхности. Для эмпирической оценки H приме нено выражение (1.2.20), в котором в качестве аргумента использовалась про странственная координата. Реальные поля двумерны, однако, в силу их сущест венной изотропности и однородности, особенно на сравнительно небольших масштабах, можно рассматривать изменения как функцию одного переменного.

Для среднемесячных значений температуры получилось, что H 0,45. Принимая во внимание существенную неопределенность данных наблюдений, можно счи тать, что фактически H = 1/2. Аналогичный результат получился и для простран ственного распределения как месячных, так и суточных сумм осадков. Можно констатировать, что в данном случае модель классического гауссовского слу чайного процесса служит хорошим приближением реальности.

Теперь вернемся к уравнению Ланжевена. Его решением как решением стохастического уравнения является выражение для дисперсии колебаний (1.2.11). Поставим вопрос о том, какие существуют значения H для аппрок симации поведения дисперсии. Ясно, что на очень малых (по сравнению с выбранным масштабом) временах H = 1/2. При рассмотрении более протя женных отрезков времени, определим H из выражения Получается, что Эта формула показывает, что уже на масштабах времени порядка (2) H 0.

Таким образом, на малых временах поведение решения уравнения Ланже вена представляет собой классический броуновский процесс. Однако с ростом времени начинает все активнее проявляться отрицательная обратная связь, пе реводящая решение к стабилизации поведения системы, в динамике которой не может быть ни положительных, ни отрицательных «сверхдлительных» трен дов. Вместо этого наблюдается появление локальных экстремумов определен ного масштаба, обусловленных сильным возвращающим воздействием отрица тельной обратной связи. В соответствии с введенными определениями при рас смотрении поведения решения на больших масштабах получающийся график зависимости от времени является фракталом с D = 2 ( = 1).

ГЛАВА 2. ФИЗИЧЕСКИЕ ОСНОВЫ ПОСТРОЕНИЯ

СТОХАСТИЧЕСКИХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

ДЛЯ ПРИРОДНЫХ ОБЪЕКТОВ

2.1. Разделение переменных и движений на быстрые и медленные Формальные математические модели, рассмотренные в предыдущей главе, показали свою способность давать статистический прогноз эволюции раз личных природных объектов при условии, что статистика процессов, опреде ляющих их эволюцию, не меняется в течении длительного времени – време ни, достаточного для того, чтобы накопить объем данных, требуемых для формирования статистической модели (определения коэффициентов регрес сии и интенсивности шумов). В действительности внешние по отношению к исследуемому объекту условия изменяются. К примеру, планетарные изме нения климата сопровождаются изменениями статистики природных процес сов регионального масштаба (засух, наводнений, лесных и степных пожаров, лавинной опасности и т.д.). Поэтому не существует гарантии от того, что ста тистические регрессионные модели с эмпирически определенными для со временного периода коэффициентами и уровнем шума будут надежными в дальнейшем. Скорее можно уверенно ожидать обратной картины. Действи тельно, целый ряд эмпирических методик, успешно применявшихся в практике гидрометеорологических прогнозов в 1960–1970-х годах, стал давать неверные результаты в 1980-х годах, когда глобальное потепление климата нарушило ка завшиеся незыблемыми связи. В этом случае методически простая задача по строения физических моделей природных процессов с включением стохасти ческой составляющей трансформируется в задачу статистической физики – теории, позволяющей вычислять параметры стохастических моделей исходя из исходных уравнений баланса энергии, импульса и вещества. В данной главе изложены методы, которые позволяют определять коэффициенты регрессии прогностических моделей исходя из исходных динамических уравнений.

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

Mazo, 1978]. Рассмотрим этот подход на примере явлений, происходящих в земной климатической системе. Согласно определению, данному Всемирной метеорологической организацией (ВМО), это есть «система, состоящая из взаимодействующих физических элементов атмосферы, океана, криосферы, поверхности суши и биомассы, которые из меняются в масштабах времени, превышающих время жизни индивидуаль ных возмущений синоптического масштаба» [Физические основы теории климата и его моделирования, 1977]. Процессы, протекающие в каждой из названных подсистем и при взаимодействии между ними имеют различные временные масштабы. Некоторые из этих масштабов существенно разнесены.

Так, время релаксации аномалий температуры поверхности океана за счет контактного теплообмена с атмосферой, составляет, при учете только верхнего квазиоднородного слоя, несколько месяцев. На этих временах атмосферные воздействия за счет потоков тепла, влаги и импульса, происходящие с харак терными временами корреляции, составляющим несколько суток, воспри нимаются океаном как «шум». При этом статистические моменты этого «по годного» шума (средние, дисперсии и т.д.) в свою очередь могут зависеть от состояния более инерционного объекта (в данном случае – океана). Разделение временных масштабов позволяет применять для расчета статистики инерцион ных объектов математический аппарат теории броуновского движения.

Начало этой теории было положено в пионерных работах А. Эйнштейна, М. фон Смолуховского и П. Ланжевена [см., например, Кайзер, 1990;

Mazo, 1978]. В первой пионерной работе А. Эйнштейна была рассмотрена задача определения среднего квадрата отклонения координаты тяжелой броунов ской частицы массы m от ее первоначального положения r. Для того, чтобы определить (r (t ) r (0))2 (угловые скобки здесь и далее означают статисти ческое осреднение), можно пойти по пути, который будет использоваться ниже – вычислить эту величину, если известна плотность распределения ве роятностей траекторий системы вблизи точки r (2.1.10). Однако Эйнштейн заменил осреднение по множеству траекторий на вычисление плотности рас пределения числа N броуновских частиц (N1) в пространстве n(r,t)=(r,t)N, считая, что в первоначальный момент все броуновские частицы сосредоточе ны в узкой окрестности r (t ) = r (0) r = r (0), а далее облако частиц распро страняется в окружающее пространство по законам молекулярной физики с коэффициентом диффузии D. Таким образом решается классическая задача диффузии c сохранением общего числа частиц N, так что при N = const ре шение для концентрации n(r,t) можно связать с вероятностной характери стикой (r,t) где – оператор Лапласа. Начальное условие для (2.1.1) – это условие сохране ния общего числа частиц N для сосредоточенного первоначально в бесконечно малой окрестности r распределения, заданного через дельта-функцию Дирака:

(r) = (r ). Оно задает симметричное в пространстве распределение вероятно стей для r = r–r, выраженное решением уравнения диффузии r = rr В силу пространственной симметрии плотности вероятности смещения вдоль отдельной оси совпадают (с заменой обозначения пространственного аргумента). Тогда средний квадрат смещения броуновской частицы будет ра вен утроенному среднему квадрату ее смещения вдоль одной из осей (в силу той же пространственной симметрии). Определяя с помощью (2.1.2) средний квадрат смещения броуновской частицы в направлении x, получаем одно из самых известных соотношений статистической физики – пропорциональ ность квадрата смещения времени процесса Следующий шаг в развитии теории сделал П. Ланжевен, который рассмот рел уравнения движения отдельной тяжелой частицы под действием соударе ний с газом легких молекул как движение отдельной частицы, движение ко торой подчиняется второму закону Ньютона (например, по координате x) с учетом линейного сопротивления трения по закону Стокса Здесь a – радиус броуновской частицы, – коэффициент трения легкого газа, обтекающего массивную броуновскую частицу.

После умножения на x каждого члена уравнения и последующего осред нения получим В (2.1.4) скорость броуновской частицы v = (v, v, v ) можно определить из постулированной А. Эйнштейном и П. Ланжевеном выполнимости теоремы статистической физики о равном распределении энергии по степеням свободы В (2.1.5) T – абсолютная температура, k – постоянная Больцмана – внеш ние параметры.

Далее П. Ланжевен предположил, что при осреднении по большому числу соударений легких молекул о поверхность массивной частицы последний член справа в (2.1.4) должен исчезать. Тогда из (2.1.4) с учетом (2.1.5) полу чаем уравнение для x Тогда находим общее решение, как общее решение однородного плюс частное решение неоднородного:

где C и C – произвольные постоянные, из чего следует где C – произвольная постоянная. Асимптотически:

Это означает (с учетом симметрии r = 3x ) Сравнение (2.1.6) и (2.1.3) приводит к известному соотношению (соотно шение Эйнштейна) молекулярной теории броуновского движения между ко эффициентом диффузии D и коэффициентом линейного трения (6 a) Соотношение Эйнштейна (2.1.7) является первым примером так называе мых флуктуационно-диссипативных соотношений, оно связывает флуктуаци онную характеристику (коэффициент диффузии) с диссипативной (коэффи циент трения).

Эти первые работы положили начало теории броуновского движения, они содержали первое стохастическое дифференциальное уравнение – уравнение Ланжевена – и уравнение для плотности распределения вероятностей – про стейший аналог уравнения ФоккераПланка. Однако в работе П. Ланжевена содержалось интуитивное предположение об обращении в ноль последнего члена в (2.1.4). Для создания методов расчета таких членов потребовалось не сколько десятилетий, в течение которых была создана теория стохастических интегралов [Гардинер, 1986].



Pages:   || 2 | 3 | 4 | 5 |   ...   | 7 |
 


Похожие материалы:

«Федеральное агентство по образованию РФ Владивостокский государственный университет экономики и сервиса _ Н.Г. МИЗЬ А.А. БРЕСЛАВЕЦ КОРЕЯ – РОССИЙСКОЕ ПРИМОРЬЕ: ПУТЬ К ВЗАИМОПОНИМАНИЮ Монография Владивосток Издательство ВГУЭС 2009 ББК 63 М 57 Ответственный редактор: Т.И. Бреславец, канд. фил. наук, профессор Дальневосточного государ ственного университета Рецензенты: С.К. Песцов, д-р полит. наук, профессор Дальневосточного государ ственного университета; И.А. Толстокулаков, канн. ист. наук, ...»

«Министерство сельского хозяйства РФ Российская академия сельскохозяйственных наук Федеральное агентство по образованию Администрация Воронежской области ГОУВПО Воронежская государственная технологическая академия ГОУВПО Московский государственный университет прикладной биотехнологии ГОУВПО Московский государственный университет пищевых производств ГОУВПО Санкт-Петербургский государственный университет низкотемпературных и пищевых технологий Ассоциация Объединенный университет имени В.И. ...»

«ВЫСШАЯ ШКОЛА ЭКОНОМИКИ НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ УНИВЕРСИТЕТ ПРИ УЧАСТИИ ВСЕМИРНОГО БАНКА И МЕЖДУНАРОДНОГО ВАЛЮТНОГО ФОНДА XI МЕЖДУНАРОДНАЯ НАУЧНАЯ КОНФЕРЕНЦИЯ ПО ПРОБЛЕМАМ РАЗВИТИЯ ЭКОНОМИКИ И ОБЩЕСТВА В трех книгах Ответственный редактор Е.Г. Ясин Издательский дом Высшей школы экономики Москва, 2011 УДК 330.101.5(063) ББК 65.012 О-42 Идеи и выводы авторов не обязательно отражают позиции представляемых ими организаций © Оформление. Издательский дом ISBN 978-5-7598-0861-9 (кн. 3) ISBN ...»

«Министерство сельского хозяйства Российской Федерации ФГБОУ ВПО Уральская государственная академия ветеринарной медицины Разработка и внедрение новых технологий получения и переработки продукции животноводства 20 марта 2013 г. Материалы международной научно – практической конференции Троицк-2013 УДК: 631.145 ББК: 65 Р - 17 Разработка и внедрение новых технологий получения и переработки продукции Р - 17 животноводства20 марта 2013 г.,. / Мат-лы междунар. науч.-практ. конф.: сб. науч. тр.– ...»

«Министерство сельского хозяйства Российской Федерации ФГОУ ВПО Уральская государственная академия ветеринарной медицины Инновационные подходы к повышению качества продукции АПК 21 марта 2012 г. Материалы международной научно-практической конференции Троицк-2012 УДК: 631.145 И-66 ББК: 65 Инновационные подходы к повышению качества продукции АПК, И-66 21 марта 2012 г. г: материалы междунар. науч.- практ. конф. / Урал. гос. академия вет. медицины. – Троицк: УГАВМ, 2012. – 148 с. Редакционная ...»

«Чернышев В.Б. Экология насекомых Москва 1996 ББК 28.68 Ч47 УДК 574.001; 595.7.15 Рецензенты: кафедра энтомологии Санкт–Петербургского университета, чл.– кор. РАН, профессор Ю.И.Чернов, профессор Г.А.Мазохин–Поршняков Издание финансируется Российским фондом фундаментальных исследований Чернышев В.Б. Экология насекомых. Учебник. – М.: Изд–во МГУ, 1996 – 304 с.: ил. ISBN 5–211–03545–3 В учебнике рассмотрены основные принципы экологии насекомых, показаны особенности образа жизни насекомых, ...»

«Т.А.Работнов ИСТОРИЯ ФИТОЦЕНОЛОГИИ Москва Аргус 1995 ББК 28.58. Р13 УДК 581.55 Научный редактор д.б.н., профессор В.Н.Павлов Р13 Работнов Т.А. История фитоценологии: Учебное пособие. - М.: Аргус, 1995. - 158 с. ISBN 5-85549-074-2 В учебном пособии рассмотрены основные этапы развития фитоценологии, включая современный период, детально охарактеризовано совершенствование методических подходов к исследованию растительности, сделан обзор важнейших направлений этой науки в настоящее время. Автор, в ...»

«МИНИСТЕРСТВО СЕЛЬСКОГО ХОЗЯЙСТВА РОССИЙСКОЙ ФЕДЕРАЦИИ Федеральное государственное образовательное учреждение высшего профиссионального образования Алтайский государственный аграрный университет Н.Е. Борисенко, О.В. Кроневальд ВЕТЕРИНАРНО-САНИТАРНАЯ ЭКСПЕРТИЗА ПРОДУКТОВ ВЫНУЖДЕННОГО УБОЯ ЖИВОТНЫХ, ПРИ ВЫЯВЛЕНИИ БОЛЕЗНЕЙ И ПРИ ИЗМЕНЕНИЯХ, ВОЗНИКАЮЩИХ В ПРОЦЕССЕ ХРАНЕНИЯ МЯСА Учебно-методическое пособие для лабораторно-практических занятий и самостоятельной работы для студентов и слушателей отдела ...»

«Федеральное государственное образовательное учреждение высшего профессионального образования Российский государственный аграрный университет – МСХА имени К.А. Тимирязева _ Студенческое научное общество имени Н.И. Вавилова 61-я СТУДЕНЧЕСКАЯ НАУЧНАЯ КОНФЕРЕНЦИЯ Секция ГЕНЕТИКА, СЕЛЕКЦИЯ И БИОТЕХНОЛОГИЯ 19 марта 2008 г. Сборник тезисов Москва, 2008 УДК 575:573.6:631.524 Сборник тезисов участников 61 студенческой научной конференции секции Генетика, селекция и биотехнология, состоявшейся 19 марта ...»

«ТЕХНОГЕННЫЕ ПОВЕРХНОСТНЫЕ ОБРАЗОВАНИЯ ЗОНЫ СОЛЕОТВАЛОВ И АДАПТАЦИЯ К НИМ РАСТЕНИЙ Пермь, 2013 МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования ПЕРМСКИЙ ГОСУДАРСТВЕННЫЙ НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ УНИВЕРСИТЕТ О.З. Ерёмченко, О.А. Четина, М.Г. Кусакина, И.Е. Шестаков ТЕХНОГЕННЫЕ ПОВЕРХНОСТНЫЕ ОБРАЗОВАНИЯ ЗОНЫ СОЛЕОТВАЛОВ И АДАПТАЦИЯ К НИМ РАСТЕНИЙ Монография УДК 631.4+502.211:582 ББК ...»

«МИНИСТЕРСТВО СЕЛЬСКОГО ХОЗЯЙСТВА РОССИЙСКОЙ ФЕДЕРАЦИИ ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ САРАТОВСКИЙ ГОСУДАРСТВЕННЫЙ АГРАРНЫЙ УНИВЕРСИТЕТ ИМЕНИ Н.И. ВАВИЛОВА БИОТЕХНОЛОГИЯ: РЕАЛЬНОСТЬ И ПЕРСПЕКТИВЫ В СЕЛЬСКОМ ХОЗЯЙСТВЕ Материалы Международной научно-практической конференции К 100-летию СГАУ имени Н.И. Вавилова САРАТОВ 2013 УДК 579.64:60 ББК 30:40.5 Биотехнология: реальность и перспективы в сельском хозяйстве: Материалы ...»

«ФГБОУ ВПО Ульяновская государственная сельскохозяйственная академия Научно-исследовательский инновационный центр микробиологии и биотехнологии Ульяновская МОО Ассоциация практикующих ветеринарных врачей АКТУАЛЬНЫЕ ПРОБЛЕМЫ ИНФЕКЦИОННОЙ ПАТОЛОГИИ И БИОТЕХНОЛОГИИ Материалы V-й Всероссийской (с международным участием) студенческой научной конференции 25 – 26 апреля 2012 года Ульяновск – 2012 Актуальные проблемы инфекционной патологии и биотехнологии УДК 631 Актуальные проблемы инфекционной ...»

«РЕСПУБЛИКА АРМЕНИЯ МИНИСТЕРСТВО ОХРАНЫ ПРИРОДЫ НАЦИОНАЛЬНАЯ ПРОГРАММА ДЕЙСТВИЙ ПО БОРЬБЕ С ОПУСТЫНИВАНИЕМ В АРМЕНИИ ЕРЕВАН 2002 НАЦИОНАЛЬНАЯ ПРОГРАММА ДЕЙСТВИЙ ПО БОРЬБЕ С ОПУСТЫНИВАНИЕМ В АРМЕНИИ Руководитель Программы: Вардеванян Ашот Ответственный редактор: Балоян Самвел Консультант: Дарбинян Нуне Министерство охраны природы Республики Армения выражает глубокую благодарность Программе окружающей среды Организации Объединенных Наций (UNEP), Секретариату Конвенции ООН “По борьбе с ...»

«МИНИСТЕРСТВО СЕЛЬСКОГО ХОЗЯЙСТВА И ПРОДОВОЛЬСТВИЯ РЕСПУБЛИКИ БЕЛАРУСЬ ГЛАВНОЕ УПРАВЛЕНИЕ ОБРАЗОВАНИЯ, НАУКИ И КАДРОВ Учреждение образования БЕЛОРУССКАЯ ГОСУДАРСТВЕННАЯ СЕЛЬСКОХОЗЯЙСТВЕННАЯ АКАДЕМИЯ АКТУАЛЬНЫЕ ПРОБЛЕМЫ ИНТЕНСИВНОГО РАЗВИТИЯ ЖИВОТНОВОДСТВА Сборник научных трудов Выпуск 15 В двух частях Часть 1 Горки БГСХА 2012 УДК 631.151.2:636 ББК 65.325.2 А43 Редакционная коллегия: А. П. Курдеко (гл. редактор), Н. И. Гавриченко (зам. гл. редактора), Е. Л. Микулич (зам. гл. редактора), Р. П. ...»

«МИНИСТЕРСТВО СЕЛЬСКОГО ХОЗЯЙСТВА РОССИЙСКОЙ ФЕДЕРАЦИИ ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ САРАТОВСКИЙ ГОСУДАРСТВЕННЫЙ АГРАРНЫЙ УНИВЕРСИТЕТ ИМЕНИ Н.И. ВАВИЛОВА Факультет электрификации и энергообеспечения АКТУАЛЬНЫЕ ПРОБЛЕМЫ ЭНЕРГЕТИКИ АПК Материалы III Международной научно-практической конференции САРАТОВ 2012 УДК 338.436.33:620.9 ББК 31:65.32 Актуальные проблемы энергетики АПК: Материалы III Международной научно практической ...»

«А.Я. Ала РОЛЬ ГОРИЗОНТАЛЬНОГО ПЕРЕНОСА ГЕНОВ В СЕЛЕКЦИИ A. Ya. Ala ROLE OF HORISONTAL TRANSFER OF GENES IN SELECTION Российская академия сельскохозяйственных наук Russian academy of agricultural sciences Всероссийский научно-исследовательский институт сои All-Russian Soybean Research Institute А.Я. Ала A. Ya. Ala РОЛЬ ГОРИЗОНТАЛЬНОГО ПЕРЕНОСА ГЕНОВ В СЕЛЕКЦИИ ROLE OF HORISONTAL TRANSFER OF GENES IN SELECTION Благовещенск, ПКИ Зея, Blagoveshchensk Zeya, УДК 633.853.52:631. ББК 41. А Ала А.Я. ...»

«Алтайский государственный технический университет им. И.И. Ползунова Сельский туристский бизнес в Алтайском крае Учебное пособие Барнаул • 2009 УДК 379.85 ББК 65.9(2Рос– 4Алт) 497.58 С 279 Авторы: А.Н. Дунец, В.В. Исаев, Н.В. Биттер, Л.И. Донскова, В.С. Ревякин, В.С. Бовтун, Т.Г. Петракова, О.Ю. Герасимова, Е.Л. Панин, А.В. Косицына Рецензент кандидат педагогических наук, доцент С.А. Гокк С 279 Сельский туристский бизнес в Алтайском крае : учебное пособие / под ред. А.Н. Дунца. – Барнаул : ...»

«Василий Скакун ВСЁ, ЧТО БЫЛО НЕ СО МНОЙ, ПОМНЮ. Ставрополь АГРУС 2013 УДК 82-3 ББК 84(2Рос=Рус)6 С42 Скакун, В. Всё, что было не со мной, помню. / Василий Ска- С42 кун. – Ставрополь : АГРУС Ставропольского гос. аграрного ун-та, 2013. – 224 с. ISBN 978-5-9596-0870-5 Каждый здравомыслящий человек, обозревая вокруг себя людей с абсолютно разными чертами характера, рано или поздно просто обязан прийти к разгадке этой тайны. Мы жи вём множество жизней, накапливая в каждой тот или иной опыт – ...»

«Федеральное агентство по образованию Российской Федерации Государственное образовательное учреждение высшего профессионального образования Ивановский государственный химико-технологический университет Т.К. Акаева, С.Н. Петрова Основы химии и технологии получения и переработки жиров Часть 1. Технология получения растительных масел Учебное пособие Иваново 2007 1 УДК 664.34.002(075) Акаева Т.К., Петрова С.Н. Основы химии и технологии получения и переработки жиров. Ч.1. Технология получения ...»






 
© 2013 www.seluk.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.