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

РОССИЙСКАЯ АКАДЕМИЯ НАУК

СИБИРСКОЕ. ОТДЕЛЕНИЕ Институт вычислительного моделирования

На правах рукописи

ЯКУБАЙЛИК ОЛЕГ ЭДУАРДОВИЧ

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

05.13.16 - применение вычислительной техники, математического моделирования и математических методов в научных исследованиях

(в экологии)

Диссертация на соискание ученой степени кандидата физико-математических наук

Научные руководители:

доктор физико-математических наук

Денисенко В.В.

кандидат физико-математических наук

Замай С.С.

Красноярск 1998

Содержание

ВВЕДЕНИЕ.........................................................................................4

ГЛАВА 1. ЗАДАЧИ МОДЕЛИРОВАНИЯ В ЭКОЛОГИИ........30

1.1. Проблемы управления качеством окружающей среды....................30

1.2. Экологическая экспертиза и системы оценки качества

окружающей среды...............................................................................33

1.3. программное обеспечение экологического мониторинга

и вычислительный эксперимент.........................................................36

1.4. Объект исследований данной работы.................................................41

ГЛАВА 2. МАТЕМАТИЧЕСКИЕ МОДЕЛИ РАСПРОСТРАНЕНИЯ ПРИМЕСЕЙ В АТМОСФЕРЕ............ 44

2.1. Общие сведения о моделях.................................................................44

2.2. Классификация существующих моделей........................................57

2.3. Стандартные модели природоохранных служб................................61

2.4. Модель Паскуилла-Гиффорда:............................................................63

2.5. Модель Института экспериментальной метеорологии..................66

2.6. Трехмерные модели переноса и диффузии примесей.......................69

2.7. Аэродинамическое моделирование...................................................72

2.8. районирование зоны загрязнения по степени опасности..............75

2.9. Заключение............................................................................................79

ГЛАВА 3. ЧИСЛЕННЫЕ И АЭРОДИНАМИЧЕСКИЕ ЭКСПЕРИМЕНТЫ И АНАЛИЗ РЕЗУЛЬТАТОВ......................83

3.1. Общая постановка задачи для уравнения переноса-диффузии примеси..................................................................................................83

3.2. Двумерная стационарная аналитическая модель...........................86

3.3. Двумерная численная модель..............................................................92

3.4. Сравнительный анализ диффузионных и эмпирических моделей 114

3.5. Двухслойная численная модель.....................................................:. 129

3.6. Результаты аэродинамического моделирования...........................134

3.7. Выводы и рекомендации.....................................................................151

г

ГЛАВА 4. ТЕХНОЛОГИИ В СОЗДАНИИ ИНФОРМАЦИОННО-ЭКСПЕРТНЫХ СИСТЕМ ПРИРОДООХРАННЫХ СЛУЖБ МЕГАПОЛИСА..................154

4.1. Проблемы реализации прикладных моделей в программном

обеспечении для решения экологических задач............................154

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

4.3. Программное обеспечение и технологии геоинформационных систем в экологическом моделировании. 161

4.4. ГИС-Web-cepbep как эффективный способ построения прикладных распределенных информационных систем

на основе ГИС- и Internet-технологий...........................................165

ВЫВОДЫ.........................................................................................169

ПРИЛОЖЕНИЕ: ПРИКЛАДНЫЕ ПРОГРАММЫ..................171

1. Архитектура программного обеспечения и

учет особенностей операционной системы...................................171

2. Комплекс программ МОНИТОР Версия 2.4....................................175

3. Программное обеспечение проекта ЭКОСФЕРА............................180

4. Опытно-конструкторская разработка МАКЕТ...............................184

5. ГИС-приложение МОНИТОР Версия 3.0..........................................187

СПИСОК ЛИТЕРАТУРЫ.............................................................191

Введение

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

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

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

Для реализации проблемно-ориентированных компьютерных систем с таким функциональным наполнением необходимы проведение серьёзного предварительного анализа требований к системе, разработка концепции и проектирование организации комплекса программ, архитектуры программного обеспечения. Речь идет о создании крупной информационной системы, в основе которой лежат подсистемы строго формализованной отчетности и моделирования распространения примесей (для проведения экологической экспертизы). Задача численного моделирования процессов распространения примесей является одной из центральных в обсуждаемом подходе. К настоящему моменту в прикладной математике сформировалась концепция решения подобных задач - это метод вычислительного эксперимента, основные принципы которого подробно раскрываются в работах H.H. Моисеева, A.A. Самарского и H.H. Яненко [71, 72].

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

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

Создание отвечающей современным стандартам информационно-аналитической системы природоохранных служб немыслимо без привлечения и практического применения результатов научных исследований и технологических решений, полученных за последние годы в информатике. Основополагающие принципы теории программирования и проектирования программного обеспечения раскрываются в работах Э. Дейкстры, Ч. Хоара, Д. Кнута, Н. Вирта, Г. Буча [73-78]. Ряд принципиально важных результатов был получен в исследовательских лабораториях компаний, занимающихся разработкой программного обеспечения (IBM, Sun Microsystems, Microsoft, Oracle, Intergraph и другие); эти результаты воплощены в созданные ими высокотехнологичные программные продукты.

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

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

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

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

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

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

Актуальность исследований

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

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

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

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

• отсутствует или затруднен прогноз состояния экосферы города в зависимости от действий субъектов и состояния объектов управления;

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

Цель и задачи исследований

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

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

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

• Анализ существующих и создание новых математических моделей распространения примесей в атмосфере с точки зрения их включения в систему экологической экспертизы, т.е. обеспеченных необходимыми начальными данными и выполнимых за разумное время на стандартном персональном компьютере;

• Выработка рекомендаций по использованию моделей в составе экологической экспертной системы и по проведению дополнительных исследований и натурных измерений - с целью повышения точности оценки и прогноза распространения примесей в атмосфере;

• Проектирование архитектуры и состава комплекса программ экологической экспертизы, и на этой основе - создание конечных программных продуктов для экологических служб;

• Апробация разработанного программного обеспечения и методов интегрированной обработки разнородных данных при решении практических задач экологической тематики.

Научная новизна

Проведенные исследования позволили получить ряд новых результатов:

• Разработана концепция поэтапного создания программного обеспечения информационно-аналитической системы природоохранных служб на примере Красноярского края. Детально исследованы вопросы реализации первого этапа (промышленная подсистема). Выработаны требования к модельному и программному обеспечению.

• Показана целесообразность использования простых моделей распространения примесей в составе программного обеспечения экологической экспертизы. При этом двумерные численные модели могут дополнить официальные модели МАГАТЭ, ОНД-86 в части учета пространственных неодно-родностей подстилающей поверхности

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

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

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

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