WWW.DISS.SELUK.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА
(Авторефераты, диссертации, методички, учебные программы, монографии)

 

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

Российский научный центр «Курчатовский институт»



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

УДК 519.676


ВОЛКОВА Елена Викторовна


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



05.13.18 – Математическое моделирование,

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


АВТОРЕФЕРАТ

диссертации на соискание ученой степени

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


Москва 2009

Работа выполнена в Российском научном центре «Курчатовский институт»


Научный руководитель: доктор физико-математических наук, профессор

Лебедев Вячеслав Иванович


Официальные оппоненты: доктор технических наук, профессор

Куранов Николай Петрович

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

Мордашев Владимир Михайлович

Ведущая организация: ГУП МосНПО «Радон»



Защита диссертации состоится ____ ______________ 2009 г. в ____ ч ____ мин. на заседании диссертационного совета Д 520.009.06 в Российском научном центре «Курчатовский институт» по адресу 123182 г. Москва, пл. Курчатова, 1.



С диссертацией можно ознакомиться в библиотеке РНЦ «Курчатовский институт»



Автореферат разослан ____ ______________ 2009 г.


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

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

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


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

Актуальность работы

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

Методики анализа чувствительности численных моделей к входным параметрам являются эффективным инструментом для характеристики погрешности модели. Если анализ погрешности позволяет указать численную оценку погрешности модели, то анализ чувствительности дает дополнительную информацию о модели, которая может использоваться для ее эффективной верификации. Анализ чувствительности особенно важен при анализе моделей с большим количеством параметров. Ведущие российские и зарубежные организации, занимающиеся исследованиями по оценке техногенного воздействия на окружающую среду, такие как Европейская Комиссия по Окружающей Среде, Геологическая Служба США (USGS), Агентство по Защите окружающей Среды (EPA) и др., включают этап анализа чувствительности в свои рекомендации при проведении численного моделирования.

Практика «статистического анализа компьютерных экспериментов» применительно к численным моделям физических явлений активно развивается с 80-х годов XX-го века. Среди первых зарубежных обзоров можно выделить работы Ш. Закс, Дж. Клайджнен, А. Сальтелли, С. Тарантола и др. В последнее время проводится большое количество исследований, посвященных анализу чувствительности численных моделей процессов в окружающей среде – в частности, моделей фильтрации и переноса загрязнения в подземных водах. Среди авторов можно выделить Дж. Хелтона, Ф. Камполонго, Р. Маддалена, Х. Ма, Б. Иоосса и других.

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

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

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

Цели и задачи работы

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

В работе были поставлены и решены следующие задачи: (1) изучение теоретической базы статистического анализа моделей и разработка методик анализа чувствительности численных моделей фильтрации и массопереноса в подземной гидросфере, (2) построение численной модели миграции радиоактивного загрязнения для условий площадки ВХРАО и геофильтрационной модели работы системы инженерной защиты от подтопления г. Казани; (3) применение разработанных методик анализа чувствительности и погрешности для указанных численных моделей.

Первая численная модель была создана с целью оценки и прогнозирования распространения радиоактивного загрязнения в подземных водах с площадки ВХРАО на территории РНЦ «Курчатовский институт». Эта модель была реализована в двух вариантах, отличающихся подробностью дискретизации и программной реализацией. Анализ чувствительности этих вариантов численных моделей включал:

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

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

Результаты, выносимые на защиту

  1. Разработка методики анализа чувствительности численной модели переноса радиоактивного загрязнения для условий площадки ВХРАО РНЦ «Курчатовский институт», результаты анализа чувствительности модели и полученные рекомендации по эффективной калибровке ее входных параметров, а также рекомендации по оптимизации сети наблюдательных скважин.
  2. Разработка методики анализа чувствительности численной геофильтрационной модели работы системы инженерной защиты от подтопления г. Казани, результаты анализа чувствительности модели и полученные рекомендации по эффективной калибровке входных параметров модели.

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

Для численной модели применительно площадки ВХРАО задача анализа чувствительности модели была решена путем построения «поверхности отклика», т.е. репрезентативной математической модели изучаемого кода («метамодели»), заменяющей исходную модель в статистических расчетах с малым временем вычислений и сохраняющей статистические свойства начальной модели. Построение поверхности отклика является отдельной задачей в каждом индивидуальном случае применения статистических методов анализа к моделям с большим временем вычислений.

Достоверность

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

Практическая значимость

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

- численной модели переноса радиоактивного загрязнения, разработанной при реабилитации площадки ВХРАО на территории РНЦ «Курчатовский институт» и ликвидации на ней старых хранилищ радиоактивных отходов;

- численной геофильтрационной модели системы инженерной защиты от подтопления г. Казани, разработанной для обоснования режима работы и оптимизации дренажа;

- численной модели миграции загрязнения тяжелыми металлами и биоочистки, разработанной для площадки Hanford 100H (штат Вашингтон, США).

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

Апробация работы

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

  • площадке ВХРАО на территории РНЦ «Курчатовский институт»;
  • территории действия системы инженерной защиты г. Казани от подтопления;
  • площадке Hanford 100H (штат Вашингтон, США).

Основные результаты работы докладывались на международной конференции «FEM-MODFLOW» (Чехия, г. Карловы Вары, 2004), семинаре CETAMA «Отбор проб и характеризация «От отбора проб до анализа» (Франция, г. Монпелье, 2005), XI всероссийской школе-семинаре «Современные проблемы математического моделирования» (г. Новороссийск, 2005), Третьей Курчатовской молодежной школе (Москва, РНЦ «Курчатовский институт», 2005), международной конференции «Моделирование процессов миграции радионуклидов в окружающей среде и вопросы развития метабаз данных для объектов Советского ядерного комплекса (РАДЛЕГ-РАДИНФО-2005)» (Москва, 2005), Второй всероссийской конференции «Современные проблемы изучения и использования питьевых подземных вод» (г. Звенигород, 2006).

Личный вклад автора

Автор участвовал в постановке и решении основных задач диссертации. С участием автора были разработаны численные модели переноса радиоактивного загрязнения в районе площадки ВХРАО РНЦ «Курчатовский институт» на основе программных модулей MODLOW и MT3DMS и численная геофильтрационная модель работы системы инженерной защиты от подтопления г. Казани. Личный вклад автора состоит в изучении и практическом использовании современных методов анализа чувствительности численных моделей, решении ряда задач по статистическому моделированию и анализу данных, разработке численной модели переноса радиоактивного загрязнения в районе площадки ВХРАО на основе программного пакета MARTHE и предложению нового подхода к анализу и верификации численных моделей фильтрации и массопереноса.

Публикации

По теме диссертации опубликовано 13 печатных работ, в том числе 8 – в соавторстве. Из них: 3 – в материалах и сборниках трудов международных конференций, 4 – в виде тезисов докладов на всероссийских конференциях, 2 – в реферируемых журналах, 2 – в иностранных журналах, 2 – в сборниках трудов научно-исследовательских институтов.

Структура и объем диссертации

Диссертация состоит из введения, 3 глав, заключения и списка литературы. Общий объем диссертации 145 страниц, в том числе 38 рисунков, 34 таблицы. Список литературы включает 109 наименований.

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


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

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

Процессы фильтрации и массопереноса в подземной гидросфере описываются уравнениями математической физики. Для их решения существуют широко используемые гидрогеологами программные модули – MODFLOW, VS2D, GMS, MARTHE, FEFLOW и др., в основу которых положены алгоритмы различных численных методов решения уравнений в частных производных. В разработку теоретических основ и численных методов решения задач фильтрации и массопереноса внесли серьезный вклад Я. Бэр, П.Я. Полубаринова-Кочина, В.М. Шестаков, В. Кинзельбах, Н.Н. Веригин, Б.С. Шержуков, Н.П. Куранов, В.И. Лебедев, Г.И. Марчук, Н.С. Бахвалов, В.И. Агошков, Р.П. Федоренко, А.В. Лапин, Г.М. Кобельков, Ю.В. Василевский, А.А. Самарский и другие.

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

Анализ чувствительности тесно связан с решением так называемых обратных задач – задач изучения свойств объектов или процессов, интерпретации результатов экспериментов по каким-либо измеренным значениям «косвенных» величин. Обратные задачи часто являются некорретными. В развитие методов решения обратных и некорректных задач внесли значительный вклад Ж.Адамар, А.Н. Тихонов, В.Я.Арсенин, А.Б. Бакушинский, А.В.Гончарский, В.В. Васин, В.К. Иванов, В.П. Танана, М.М. Лаврентьев, В.А. Морозов, А.М. Денисов, А.И. Прилепко и другие. Анализ чувствительности играет важную роль при решении обратных задач, обеспечивая исследователя вспомогательной информацией о типе зависимости между входными и выходными параметрами модели, наиболее «значимых» параметрах и проч.

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

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

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

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

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

Проблема представления функций в виде комбинации функций меньшего числа переменных берет свое начало с одной из задач, сформулированных Гильбертом еще в 1900 г., исследованием которой в разные годы занимались А.Н. Колмогоров, В.И. Арнольд, А.Г. Витушкин, Г.М. Хенкин и другие.

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

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

,

где и при. При этом первый интеграл берется от 0 до 1 по всем переменным. На основе такого разложения модельной функции f строятся глобальные индексы чувствительности Si1…ik, называемые также индексами Соболя, которые характеризуют вклад дисперсии переменных xi1,…,xik в дисперсию функции f. Термин «дисперсия» происходит из вероятностного подхода к вычислению интегралов, входящих в индексы.

На практике для вычисления индексов чувствительности Соболя используется метод Монте-Карло, требующий значительного числа запусков модели для оценки каждого индекса. Таким образом, глобальный анализ чувствительности попадает в область статистических методов обработки данных. Это быстро развивающаяся в настоящее время область прикладной математики, в которой многими авторами уже создано большое количество подходов и методов. Среди таких авторов, например, А.В. Крянев, Г.В. Лукин, В.Я. Арсенин, В.Н. Вапник, Ф.П. Васильев, Н. Джонсон, Ф. Лион, Г. Крамер, М.М. Лаврентьев, Ш. Закс, Г.И. Марчук, Ю.Н. Тюрин, В.П. Шутяев и другие.

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

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

Вторая глава диссертации описывает практическое исследование по численному моделированию переноса радиоактивного загрязнения на площадке ВХРАО РНЦ «Курчатовский институт» (Москва) и применению методов анализа чувствительности к построенным моделям для оценки погрешности и сравнения моделей.

Аппарат анализа чувствительности, описанный в первой главе, был применен к двум численным моделям распространения радиоактивного загрязнения с площадки ВХРАО. Первая из моделей основана на программных кодах MODFLOW и MT3DMS (модель I), вторая – на программном пакете MARTHE (модель II). Наличие двух моделей объясняется необходимостью верификации результатов прогнозных расчетов. Одной из целей исследования было сравнение результатов анализа чувствительности двух различных моделей одной и той же площадки.

Площадка ВХРАО имеет площадь около 2 га и расположена около внешнего периметра Центра (Рис. 1) в непосредственной близости от жилых кварталов г. Москвы. На этой площадке находились 11 старых хранилищ радиоактивных отходов, сооруженных между 40-ми и 70-ми годами прошлого столетия. Обследование грунтов показало, что радиоактивное загрязнение присутствует не только на поверхности, но и имеет тенденцию к распространению в подземные воды и связано, главным образом, с радионуклидами Cs-137 и Sr-90. Численная модель переноса радиоактивного загрязнения с подземными водами для условий площадки ВХРАО была создана для выполнения прогнозных расчетов распространения загрязнения и сопровождения реабилитационных работ.

Рис. 1. Площадка ВХРАО РНЦ «Курчатовский институт»

Численная модель на основе программного пакета MODFLOW (модель I) была создана группой специалистов ДАР\ВОДГЕО и РНЦ «Курчатовский институт» под руководством А.В. Расторгуева и описывала распространение радионуклидов Sr-90 и Cs-137 в двух верхних водоносных горизонтах. Численная модель на основе MARTHE (модель II) была разработана автором. В сравнении с моделью I, модель II рассматривает меньшую область моделирования и меньший временной интервал и описывает перенос радионуклида Sr-90 только в верхнем водоносном горизонте.

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

На основе результатов численного моделирования было получено, что подземные воды с содержанием Sr-90, превышающем «уровень вмешательства», не покинут пределов площадки ВХРАО. Этот вывод важен с точки зрения планирования дальнейших работ по реабилитации площадки и оценки риска распространения радиоактивного загрязнения на окружающие жилые кварталы.

Локальный анализ чувствительности численной модели I был проведен в несколько этапов. Сначала были определены поля чувствительности модельного прогноза к выбранным параметрам в выбранный момент времени. Затем входные параметры были отсортированы по степени влияния на прогнозные модельные концентрации и вычислены ковариационная матрица и корреляции параметров.

Для выбранного момента времени вычислялись поля чувствительности к 7-ми входным параметрам моделей, описывающих интенсивность источников загрязнения, длины стресс-периодов и коэффициент сорбционного распределения по Sr-90. При этом для каждого из анализируемых параметров модель запускалась по одному разу с измененным на 1% значением параметра.

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

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

На основе вычисленных карт чувствительности параметры были отсортированы по степени влияния на модельный прогноз. Наибольшие значения чувствительности прогнозируемых концентраций Sr-90 были выявлены к удельной активности утечек из одного из источников загрязнения (одного из старых хранилищ) и длительности периода его действия. На основе полученных данных для семи исследуемых параметров были рассчитаны корреляции между анализируемыми параметрами.

Рис. 2. Карта суммарной чувствительности модели I для площадки ВХРАО

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

Для глобального анализа чувствительности модели I было выбрано 17 входных параметров модели, описывающих коэффициент фильтрации различных зон и слоев модели, коэффициент сорбционного распределения по Sr-90 в разных слоях модели, интенсивность инфильтрационного питания различных зон площадки и интенсивность источников загрязнения. Указанный набор входных параметров для глобального анализа чувствительности был определен на основе общих соображений потенциального влияния на прогнозируемые концентрации. Кроме того, представлялось целесообразным включить в анализ все входные параметры, рассмотренные на этапе локального анализа чувствительности, чтобы затем сравнить полученные результаты. В качестве выходных параметров анализа были взяты интегральные характеристики прогнозируемого моделью ореола загрязнения – площадь части ореола с концентрацией более 1 Бк/л для верхнего горизонта и более 0,01 Бк/л для нижнего горизонта. Такой выбор выходных параметров анализа объясняется относительной простотой вычисления этих характеристик, их наглядностью, а также тем, что эти характеристики описывают прогнозируемый ореол «в целом», в меньшей степени подвержены ошибкам моделирования, чем локальные характеристики (такие как, например, прогнозируемые значения концентраций в отдельных точках).

Каждый входной параметр был проанализирован с целью выбора интервала его возможного изменения и функции распределения вероятности на этом интервале. На основе выбранных интервалов и распределений вероятности методом Монте-Карло было сгенерировано 230 вариантов наборов значений входных параметров. Для автоматизации циклических вычислений была разработана специальная программа на языке С++, которая считывала на каждой итерации набор значений входных параметров, создавала соответствующие входные файлы для программных вычислительных модулей MODFLOW и MT3DMS, запускала расчет и на выходе сохраняла значения прогнозных концентраций и напоров на конец прогнозного периода. В общей сложности было выполнено 230 итераций, что заняло около 44 часов машинного времени на персональном компьютере с процессором Pentium 4.

Гистограммы полученных значений выходных параметров, определяющих площадь ореола в верхнем водоносном горизонте (SUM1) и площадь ореола загрязнения в нижнем водоносном горизонте (SUM2) приведены на рис. 3. Как видно из рисунка, выходные параметры имеют распределение, близкое к экспоненциальному. Значения SUM1 и SUM2, соответствующие начальным наборам входных параметров, лежат в «хвосте» соответствующих распределений, то есть им соответствуют относительно небольшие вероятности. Такой результат позволяет сделать вывод о том, что ореол распространения загрязнения в среднем может быть меньше, чем прогнозируемый моделью изначально.

Значение коэффициента корреляции, вычисленное для полученных выходных параметров SUM1 и SUM2, составило 0,247, что позволило сделать вывод о низкой корреляции между этими параметрами.

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

Рис. 3. Гистограммы значений выходных параметров SUM1 и SUM2 (крестом отмечено значение, соответствующее начальному набору значений параметров)

Для полученных выборок значений входных и выходных параметров были вычислены индексы чувствительности. Было получено, что параметрами, дающими наибольший вклад в погрешность при моделировании SUM1, являются коэффициенты фильтрации и сорбционного распределения 3-го слоя модели, а также интенсивность инфильтрационного питания в зоне сильных утечек из трубопровода и коэффициенты фильтрации 2-го и 4-го слоев модели. Для погрешности в значении SUM2 наиболее влиятельными параметрами являются коэффициент сорбционного распределения по Sr-90, заданный в 4-ом, 5-ом и 6-ом слоях модели, а также коэффициент фильтрации 4-го и 5-го слоев модели.

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

Результаты глобального анализа чувствительности, проведенного для модели II, подтвердили выводы, полученные для модели I.

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

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

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

Анализ чувствительности численной модели геофильтрации в районе Заречной части города был проведен в несколько этапов. Этапы отличались набором рассматриваемых входных параметров и интервалов их возможного изменения. На каждом этапе для выбранных входных параметров и данных об интервалов их возможного изменения по методике Латинских Гиперкубов генерировалась выборка из 400 наборов значений. Для каждого из 400 наборов значений запускалась геофильтрационная модель (на основе пакета MODFLOW), сохранялись значения напоров в наблюдательных скважинах для каждого стресс-периода модели и значения дренажного годового суммарного стока. На основе выборок значений входных и выходных параметров вычислялись индексы чувствительности, на основе выборки значений выходных параметров строились графики минимаксных интервалов возможных значений напоров в скважинах (например, Рис. 4).

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

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

Выводы

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

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

На основе разработанного алгоритма были проанализированы две альтернативные численные модели переноса радиоактивного загрязнения с площадки ВХРАО РНЦ «Курчатовский институт» (Москва), разработанные на основе программных пакетов MARTHE и MODFLOW, соответственно (одна из моделей построена автором).

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

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

Численная модель переноса радиоактивного загрязнения на площадке ВХРАО обладала значительным временем вычислений (порядка 20 минут), статистические вычисления с использованием модели заняли бы нерационально большое время (порядка 7 лет). Задача глобального анализа чувствительности этой модели была решена путем применения в вычислениях поверхности отклика. В таком случае длительность вычислений составила несколько суток.

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

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

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


Литература по теме диссертации

  1. Volkov V.G., Tsurikov D.F., Zverkov Yu.A., Rastorguev I.A., Volkova E.V., Rastorguev A.V., Buharin K.L. “Development of models of radionuclide migration from the RRC “Kurchatov Institute” waste disposal site,” Proceedings of the International conference on Finite Element Models, MODFLOW, and More: Solving Groundwater Problems, pp. 95-98, Karlovy Vary, Czech Republic, September 13-16, 2004
  2. Волкова Е.В. «Статистический анализ поля коэффициента фильтрации», Сборник тезисов докладов к научно-практической конференции, посвященной 70-летию ФГУП «НИИ ВОДГЕО» «Водоснабжение, водоотведение, гидротехника, инженерная гидроэкология», с. 144-145, Москва, 2004
  3. Волкова Е.В. «О глобальном анализе чувствительности для численной модели распространения радионуклидного загрязнения с подземными водами для площадки ВХРАО РНЦ «Курчатовский институт», Сборник трудов «Проблемы Инженерной Геоэкологии», стр. 51-62, Изд-во ЗАО «ДАР\ВОДГЕО», 2005
  4. Волкова Е.В. «О статистическом анализе модели переноса радиоактивного загрязнения для площадки ВХРАО РНЦ «Курчатовский институт», Москва», - В сб.: Труды XI Всероссийской школы-семинара “Современные проблемы математического моделирования”, Абрау-Дюрсо, 2005, с. 87 – 96.
  5. Volkova E., Van Dorpe F., Iooss B. «Application de l’analyse de sensibilit globale un modle hydrogologique (site de stockage de dchets radioactifs)», Sminaire CETAMA «Echantillonnage et Caractrisation «Du Prlvement l’Analyse», Montpellier, France, October 10-12, 2005
  6. Волков В.Г., Волкова Е.В., Расторгуев И.А., Зверков Ю.А., Расторгуев А.В. «Разработка численной модели распространения подземными водами радиоактивного загрязнения с площадки ВХРАО РНЦ «Курчатовский институт», идентификация параметров и анализ модели», Международная конференция «Моделирование процессов переноса радионуклидов в окружающей среде и вопросы разработки баз метаданных по радиационным объектам Советского ядерного комплекса (РАДЛЕГ-РАДИНФО-2005)», 31 октября – 2 ноября, Москва, 2005
  7. Царев А.Г., Волкова Е.В. «Локальный анализ чувствительности для модели распространения радиоактивного загрязнения площадки ВХРАО РНЦ «Курчатовский институт», 3-я Молодежная Курчатовская Школа, 14-16 ноября, Москва, 2005
  8. Волкова Е.В. «Об анализе чувствительности модели переноса загрязнения подземными водами для площадки ВХРАО РНЦ «Курчатовский институт», 2-я Всероссийская Конференция «Современные проблемы изучения и использования питьевых подземных вод (памяти Л.С. Язвина)», 7-10 февраля, Звенигород, 2006
  9. Volkova, E., B. Iooss, F. Van Dorpe, “Global sensitivity analysis for a numerical model of radionuclide migration from the RRC “Kurchatov Institute” radwaste disposal site”, Stochastic Environmental Research and Risk Assessment, vol. 22, N. 1, pp. 17-31, 2008
  10. Волков В.Г., Волкова Е.В., Расторгуев И.А., Зверков Ю.А., Расторгуев А.В. «Моделирование переноса радиоактивного загрязнения подземными водами при реабилитации площадки старых хранилищ РНЦ «Курчатовский институт», Атомная Энергия, т. 104, вып. 2, стр. 115-118, 2008
  11. Marrel, A., B. Iooss, F. Van Dorpe, E. Volkova, “An efficient methodology for modeling complex computer codes with Gaussian processes”, Computational Statistics and Data Analysis, vol. 52, N.10, pp. 4731-4744, 2008
  12. Волкова Е.В., Расторгуев И.А., Расторгуев А.В., Копылов А.С., Куранов П.Н. «Численное моделирование для обоснования регламента работы объектов инженерной защиты г. Казани», Сборник трудов НИИ ВОДГЕО, Изд-во ЗАО «ДАР\ВОДГЕО», 2009, в печати
  13. Волкова Е.В. «Анализ чувствительности численной модели переноса радиоактивного загрязнения подземными водами для площадки старых хранилищ РНЦ «Курчатовский институт», принято к публикации в журнале Атомная Энергия, 2009



Автор диссертации выражает благодарность своему научному руководителю, д.ф.-м.н. профессору Лебедеву В.И. за руководство написанием диссертации. Автор глубоко признателен к.т.н. Расторгуеву А.В. за помощь в постановке задач и постоянное внимание к работе. Автор благодарит коллективы НИИ ВОДГЕО, Института проблем безопасного использования ядерной энергии и Научно-технического комплекса «Реабилитация» РНЦ «Курчатовский институт» за оказанную помощь и полезные советы.



 
Похожие работы:

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

«СИДЕЛЬНИКОВ Михаил Владимирович Разработка системЫ ПОДДЕРЖКИ ПРИНЯТИЯ ОРГАНИЗАЦИОННО-управленЧЕСКИХ РЕШЕНИЙ В ПРОИЗВОДСТВЕННОЙ ДЕЯТЕЛЬНОСТИ многоуро в невЫХ комп а ниЙ 05.13.10 – Управление в социальных и экономических системах (экономические науки) АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата экономических наук Ростов-на-Дону – 2007 Работа выполнена в ФГОУ ВПО ПИ ЮФУ на кафедре экономики и прикладной математики Научный руководитель: кандидат...»

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

«Шульга Татьяна Эриковна МАТЕМАТИЧЕСКИЕ МОДЕЛИ ФУНКЦИОНАЛЬНО ИЗБЫТОЧНЫХ ДИСКРЕТНЫХ СИСТЕМ Специальность 05.13.18 Математическое моделирование, численные методы и комплексы программ Автореферат диссертации на соискание ученой степени доктора физико-математических наук Саратов-2010 Работа выполнена в ГОУ ВПО Саратовский государственный технический университет Научный консультант: доктор технических наук, профессор Сытник Александр Александрович

«АФАНАСЬЕВ Александр Николаевич ГРАФО-АНАЛИТИЧЕСКИЙ ПОДХОД К АНАЛИЗУ И КОНТРОЛЮ ПОТОКОВ РАБОТ В АВТОМАТИЗИРОВАННОМ ПРОЕКТИРОВАНИИ СЛОЖНЫХ КОМПЬЮТЕРИЗОВАННЫХ СИСТЕМ Специальность: 05.13.12 – Системы автоматизации проектирования (промышленность) АВТОРЕФЕРАТ диссертации на соискание ученой степени доктора технических наук Ульяновск – 2011 Работа выполнена на кафедре Вычислительной техники Ульяновского государственного технического университета Научный консультант: доктор...»

«Кружилов Иван Сергеевич МЕТОДЫ И ПРОГРАММНЫЕ СРЕДСТВА ПОВЫШЕНИЯ ЭФФЕКТИВНОСТИ РАСПОЗНАВАНИЯ ГРУПП ЗВЕЗД В АВТОНОМНОЙ АСТРОНАВИГАЦИИ Специальность 05.13.11 – Математическое и программное обеспечение вычислительных машин, комплексов и компьютерных сетей АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Москва – 2010 Работа выполнена в ГОУ ВПО Московском энергетическом институте...»

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

«Майоров Николай Николаевич Разработка и исследование математического, лингвистического и программного обеспечения подсистемы САПР построения математических моделей гетерогенных объектов Специальность: 05.13.12 Системы автоматизации проектирования ( приборостроение) АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Санкт-Петербург 2008 Работа выполнена в государственном образовательном учреждении высшего профессионального образования...»

«И сламов Ринат Робертович МОДЕЛИРОВАНИЕ ДИНАМИЧЕСКИХ СИСТЕМ С ГИРОСКОПИЧЕСКОЙ СТРУКТУРОЙ ПРИ ПАРАМЕТРИЧЕСКИХ ВОЗМУЩЕНИЯХ Специальность 05.13.18 - Математическое моделирование, численные методы и комплексы программ АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук Уфа 2007 Работа выполнена на кафедре вычислительной техники и защиты информации в ГОУ ВПО Уфимский государственный авиационный технический университет Научный руководитель...»

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

«Коротиков Сергей Викторович ПРИМЕНЕНИЕ СЕТЕЙ ПЕТРИ В РАЗРАБОТКЕ МНОГ О ПОТОЧНОГО ПРОГРАММНОГО ОБЕСПЕЧЕНИЯ С ОГРАНИЧЕННЫМИ РАЗД Е ЛЯЕМЫМИ РЕСУРСАМИ НА ПРИМЕРЕ ЦЕНТРОВ ДИСТАНЦИО Н НОГО УПРАВЛЕНИЯ И КОНТРОЛЯ Специальность 05.13.11 – Математическое и программное обеспечение вычислительных машин, комплексов и компьютерных сетей АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Новосибирск – 2008 Работа выполнена в Государственном образовательном...»

«Жегуло Ольга Анатольевна ИССЛЕДОВАНИЕ И РЕАЛИЗАЦИЯ Непроцедурных преобразований программ для построения расширяемой системы распараллеливания 05.13.11  — Математическое и программное обеспечение вычислительных машин, комплексов, систем и сетей Автореферат диссертации на соискание учёной степени кандидата технических наук Ростов-на-Дону – 2007 Работа выполнена на кафедре информатики и вычислительного эксперимента факультета математики, механики и компьютерных наук Южного...»

«Волков Александр Валерьянович ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ПРЕОБРАЗОВАНИЯ ЛАЗЕРНЫХ ИМПУЛЬСОВ ПРИ ОДНОФОТОННОМ И ДВОЙНОМ РЕЗОНАНСАХ Специальности: 05.13.18 – Математическое моделирование, численные методы и комплексы программ 01.04.21 – Лазерная физика АВТОРЕФЕРАТ диссертации на соискание учёной степени кандидата физико-математических наук Саратов 2009 Работа выполнена в ГОУ ВПО Саратовский государственный технический университет Научный руководитель : доктор физико-математических...»

«Кетова Каролина Вячеславовна УДК 517.977.5 + 519.866 : 330.14 (06) разработка методов исследования и оптимизация стратегии развития экономической системы региона Специальности: 08.00.13 – “Математические и инструментальные методы экономики”; 05.13.18 – “Математическое моделирование, численные методы и комплексы программ” Автореферат диссертации на соискание ученой степени доктора физико-математических наук Ижевск 2008 Работа выполнена в Ижевском государственном техническом...»

«Никольский Святослав Олегович МОДЕЛИРОВАНИЕ ОЦЕНКИ ХАРАКТЕРИСТИК НАДЕЖНОСТИ БАНКОВСКИХ ТИРАЖНЫХ ПРОГРАММНЫХ СИСТЕМ НА ОСНОВЕ НЕЙРОСЕТЕВЫХ ТЕХНОЛОГИЙ Специальность 05.13.18 – Математическое моделирование, численные методы и комплексы программ АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Брянск 2006 Работа выполнена в ГОУВПО Брянский государственный технический университет. Научный руководитель Официальные оппоненты Ведущее предприятие Защита...»

«Баринов Сергей Владимирович Разработка и исследование комплекса генетических алгоритмов разбиения схем с учетом временных задержек Специальность: 05.13.12 – системы автоматизации проектирования АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Таганрог – 2008 Работа выполнена в Технологическом институте Южного федерального университета в г. Таганроге. Научный руководитель: заслуженный деятель науки РФ, доктор технических наук, профессор Курейчик...»

«УДК 621.372.54.061 ФИЛИППОВ АНДРЕЙ ВЛАДИМИРОВИЧ Структурно-параметрический СИНТЕЗ РЕЗИСТИВНО-ЕМКОСТНЫХ ЭЛЕМЕНТОВ С РАСПРЕДЕЛЕННЫМИ ПАРАМЕТРАМИ СО СТРУКТУРОЙ СЛОЕВ ВИДА R-CG-NR Специальность: 05.13.01 – Системный анализ, управление и обработка информации (в науке и технике) АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Ижевск 2010 Работа выполнена в государственном образовательном учреждении высшего профессионального образования Ижевский...»

«Котельникова МАрия николаевна Моделирование процессов пенообразования в растворах Специальность 05.13.18 – Математическое моделирование, численные методы и комплексы программ Автореферат диссертации на соискание ученой степени кандидата технических наук АСТРАХАНЬ – 2010 Работа выполнена в ГОУ ВПО Астраханский государственный университет Научный руководитель: доктор химических наук, профессор Алыков Нариман Мирзаевич Официальные оппоненты: доктор химических наук, профессор...»

«ЯКИМЕНКО МАРИАННА ВЛАДИМИРОВНА ИНФОРМАЦИОННО-ИНСТРУМЕНТАРНОЕ ОБЕСПЕЧЕНИЕ ПРОЦЕССА УПРАВЛЕНИЯ РЕГИОНАЛЬНЫМ ТУРИСТСКО-РЕКРЕАЦИОННЫМ КОМПЛЕКСОМ Специальность: 05.13.10 – Управление в социальных и экономических системах (экономические науки) АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата экономических наук Ростов-на-Дону - 2008 Диссертация выполнена на кафедре государственного и муниципального права и управления Технологического института Южного федерального...»

«МАНДРА Андрей Геннадьевич СТРУКТУРНОЕ МОДЕЛИРОВАНИЕ И АВТОМАТИЧЕСКОЕ УПРАВЛЕНИЕ ДИФФУЗИОННЫМИ ПРОЦЕССАМИ ХИМИЧЕСКИХ ТЕХНОЛОГИЙ Специальность 05.13.06 – Автоматизация и управление технологическими процессами и производствами (промышленность) АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Самара 2011 Работа выполнена на кафедре Автоматика и управление в технических системах Государственного образовательного учреждения высшего профессионального...»








 
2014 www.avtoreferat.seluk.ru - «Бесплатная электронная библиотека - Авторефераты диссертаций»

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