^ Наац В. И., Наац И. Э., Рыскаленко Р. А.
Параметризованные модели теории переноса в задачах экологического мониторинга.
ПАРАМЕТРИЗОВАННЫЕ МОДЕЛИ ТЕОРИИ ПЕРЕНОСА В ЗАДАЧАХ ЭКОЛОГИЧЕСКОГО МОНИТОРИНГА АТМОСФЕРЫ И ПРИНЦИП
МИНИМАКСА
В. И. Наац, И. Э. Наац, Р. А. Рыскаленко
PARAMETERIZED MODELS OF TRANSPORT THEORY IN PROBLEMS OF ATMOSPHERE ECOLOGICAL MONITORING AND MINIMAX PRINCIPLE
Naath V. I., Naath I. E.,
Riscalenko R. A.
Parameterized computational models of the transport theory are considered in the paper with regard to the atmosphere ecological monitoring problems and minimax method, which permits to estimate the decision effectiveness in the conditions of information uncertainty.
В статье рассматриваются параметризованные вычислительные модели теории переноса применительно к задачам экологического мониторинга атмосферы и метод ми-нимакса, позволяющий оценивать эффективность решений в условиях информационной неопределенности.
Ключевые слова: экологический мониторинг, пограничный слой атмосферы, концентрация загрязняющих веществ, математические модели теории переноса, принцип ми-нимакса.
УДК 519.5, 519.6
Рассматриваются модели для решения прогнозных задач экологического мониторинга пространственно-временной изменчивости поля концентрации загрязняющих веществ (ЗВ) в пограничном слое атмосферы в ситуациях с так называемыми несанкционированными (аварийными) выбросами в атмосферу. Вычислительные модели предназначены для их реализации в структуре информационно-измерительных систем поддержки принятия решений в критических ситуациях. Эти модели основаны с одной стороны на априорной информации в виде уравнений турбулентного переноса ЗВ и предположениях о постоянстве скорости
ветра V (х, г) и поля коэффициента турбулентной диффузии к (х, г) в пределах локального контролируемого региона. С другой стороны они основаны на оперативных сопутствующих прогнозу изменениях поля концентрации ЗВ в некоторых опорных точках (постах экологического контроля). Авторами предлагается алгоритм, позволяющий прогнозировать пространственно-временную динамику поля концентрации ЗВ в критических ситуациях и оценивать эффективность тех или иных принципиальных решений в условиях информационной неопределенности на основе принципа минимакса.
При решении задачи прогноза пространственно-временной изменчивости поля концентрации ЗВ в атмосфере определяющую роль играют математические модели теории переноса субстанции и методы
численного решения соответствующих уравнений. Простейшим из этих уравнений является так называемое полуэмпирическое уравнение турбулентной диффузии в пограничном слое атмосферы, которое можно записать для одномерного случая в виде [1]:
^ + аи + к аф X). (1)
дx дх дх ^ дх 0
Уравнение (1) связывает концентрацию ЗВ и(х, X) в контролируемом объеме среды О = Ох х Qx = [0, X ]х[0, Т] с полями V (х, X), К(х, X), £(х, X). При этом предполагается, что
скорость ветра V(х, X) и коэффициент турбулентной диффузии К(х, X) в точке (х, X) могут быть связаны друг с другом неким функциональным уравнением, скажем ^ (х, X V, К ) = C(consx) . Ниже кратко рассматриваются вопросы практической реализации подобного способа априорного доопределения модели (1) в задачах экологического прогноза в области О.
Напомним прежде всего, что уравнение (1) при численных оценках поля и(х, X) должно быть доопределено граничными и начальными значениями, а именно, и(х, X = 0) = ио (х), и(0, X) = ), и(X, X) = ¡Ц (X). Приближенные аналитические модели указанных функций ио (х), ) и ¡12 (X) должны быть получены на основе измерений концентрации ЗВ в экологических постах (точках контроля). Подобная информационно-измерительная процедура необходима для так называемой «привязки» модели (1) к контролируемой области О. Ясно, что в этих же точках могут быть осуществлены измерения метеопараметров атмосферы (скорость ветра, температура, влажность и т. д.), позволяющих прогнозировать поле скорости ветра в пределах области О на основе соответствующих априорных моделей. Если считать, что поле V(м, X), где М = М (х, у, z) е О, известно в пределах О, то в рамках предлагаемого метода остается найти функциональную связь между параметрами Vx (М, X), Vy (М, X) , Vz (М, X) и коэффициентами турбулентной диффузии Кх (М, X), Ку (М, X), К2 (М, X) в соответствующих направлениях Ох, Оу и Oz . Для этой цели может быть использована эмпирическая формула вида:
1
К = сЬ2
( дVx Л2 (дVx Л2 (дV Л2 (дК Л2 2
1 + 1 + 1 + —- 1
^ дх 0 { & 0 [ дх 0
(2)
где с и Ь - некие константы (параметры), числовые значения которых согласуются с так называемыми турбулентными состояниями пограничного слоя атмосферы. В частности, в соответствии с принятыми рекомендациями, величина с близка к значению 0.4, а параметр Ь (размерность м2) лежит в диапазоне 20-200 для наблюдаемых состояний пограничного слоя атмосферы. Для практического использования соотношения (2) необходимо предварительно
оценить значение параметра Ь . Это возможно на основе измерений значений поля V (Р, X). Не останавливаясь подробно на обосновании формулы (2), отметим, что главное в ней это констатация наличия зависимости значений коэффициентов турбулентной диффузии от степени пространственной изменчивости компонент скорости ветра, что вполне соответствует основным положениям теории турбулентных уравнений сплошных сред. В качестве примера можно сослаться на качественную модель турбулентной диффузии в пограничном слое атмосферы Лайтхмана [2]. Исходя из этого, при решении вычислительных задач вместо формулы (2) предлагается следующая модель функциональной связи полей V (Р, X) и К (Р, X) (тензорная величина):
Наац В. И., Наац И. Э., Рыскаленко Р. А.
Параметризованные модели теории переноса в задачах экологического мониторинга...
K (X ) = L
3
Е Р.
j=1
f \2
dx, V j 0
(3)
где Ру - элементы некой весовой матрицы Р : ри > 0, X Р у = 1, " У = 1,2,3 . Заметим, что
=1
матрица Р в принципе может быть связана с теорией марковских процессов [3]. Если модель переноса (1) доопределяется уравнением связи (3), то целесообразно общую модель писать в виде системы уравнений:
du * d (1Г \ ^
3d
3d
i
dt
=1 dx
=1 dx
1 V
к du
dx
Л
= S (Q, t),
i 0
1
K (Q, t ) = L
3
Е Р ij
j=1
dV (Q, t)
dx,.
,i = 1,2,3,
(4а)
(4б)
где Q - точка области О с координатами Х}, Х2, Х3 . Для численного решения (4) авторами были разработаны соответствующие алгоритмы на основе схем покоординатного расщепления [1]. Основной вопрос, связанный с обобщенной моделью переноса (4), касается оценки элементов матрицы Р = {р у }3Х3, определяющей структуру поля турбулентности реальной среды. Поскольку вести речь о прямых наблюдениях турбулентного поля в локальных объемах движущейся среды не представляется возможным, то остается прибегнуть к косвенной оценке этой матрицы. Ясно, что если матрица Р принадлежит некоторому семейству матриц Щ, то решая последовательно систему (4), можно построить последовательность распределений I; Р)}
на Щ . Поскольку в работе рассматриваются задачи экологического мониторинга с использованием информационно-измерительной системы, то оценить подходящий вариант структурной матрицы Р можно из решения следующей оптимизационной задачи
штшах|(х,Р)-~(х,^)|, (5)
Рей
где u(x, t) - значения поля концентрации ЗВ в контролируемых точках x области W . Критерий (5) в практических применениях удобно использовать в виде
m r r
F (P)= Е ¿v||u (xv, tv; P)-~ (x, t), (7)
v=1
где 1v - некие весовые коэффициенты, связанные с ранжированием точек наблюдения в об-
m
ласти W, Е 1v = 1, 1v > 0. В этом случае заключительной операцией является решение опти-
v=1
мизационной задачи вида:
min F(Р)=F(Р*). (8)
Рей 4 '
Дальнейшее прогнозирование поля концентрации ЗВ осуществляется по распределению u (q, t; Р ) для точек (Q, t) вне области наблюдений (вне точек (xv, tv )). Теоретический анализ
модели (4) и алгоритмической задачи (5) (существование решений и устойчивость) выходит за рамки настоящей работы. Соответствующие регуляризирующие алгоритмы, обеспечивающие
2
2
устойчивость вычислительных схем, особенно в части усвоения эмпирических данных, частично были рассмотрены в ранее опубликованных работах авторов [4,5]. Использование параметризованных моделей переноса типа (3), определенных на пространстве Ш, естественно связано со значительной долей неопределенности в получаемых прогнозах, которые вряд ли могут быть «точными». Подобную неопределенность в данном классе задач предлагается преодолевать в той или иной системе на основе принципа минимакса, как это и принято в теории исследования операций. Решение задач мониторинга окружающей среды требует использования методов этой теории в полной мере.
Следующим шагом в развитии модели (4) является предложение дополнить её соответствующей аэродинамической моделью пограничного слоя атмосферы, а именно уравнением Навье-Стокса, что позволит выполнять расчеты компонент поля скорости ветра. В результате приходим к следующей обобщенной модели теории переноса примесей в пограничном слое атмосферы:
Я/ *
— + аи + (Иу(уи) - сСХу(КУи) = ^ (0, X), (9а)
дХ
д¥ (' X' 2 * *
д- + (|У| - КУ2У = -^ , (9б)
К (0, X ) = Ь2
Е л
1=1
'дГг (0, X)
дх-
, X = 1,2,3. (9в)
Модель (9), как и всякая модель физического явления, определяется набором неизвестных параметров, которые подлежат оценке на основе экспериментальных данных в каждом конкретном случае. К подобным коэффициентам относятся а , К, Ь2 и матрица Р еШ . Зависимость решения и(0, X) и V(0, X) от неизвестных параметров естественно требует указать способ их оценки. Модель (9) ориентирована на практическое решение в рамках информационно-измерительной системы, которая в штатном режиме обеспечивает контроль полей и(0, X) и
V(0, X) в точках (0У, Xv ).
В этой ситуации в целях преодоления неопределенности естественно прибегнуть к принципу минимакса, т. е. к решению оптимизационной задачи типа
плитах« ~0,XI)-и0,XI,р)|2 + «2 ■ |~0п,X-)- ,X-,р)
(10)
РеШ (д^)
Решением этой задачи является Р еШ. Значение Р позволяет далее решать задачи прогноза поля и (0, X; Р ) вне области контроля QxQX.
Разработке алгоритмического и программного обеспечения подобных задач посвящены некоторые работы авторов [6-10]. Однако полностью данная проблема не решена и требуются дальнейшие разработки и исследования.
ЛИТЕРАТУРА
1. Семенчин Е. А., Наац В. И., Наац И. Э. Математическое моделирование нестационарного переноса примеси в пограничном слое атмосферы. - М.: Физматлит, 2003.
2. Лайтхман Д. Л. Физика пограничного слоя атмосферы. - Л., 1970.
3. Наац И. Э., Семенчин Е. А. Математическое моделирование динамики пограничного слоя в задачах мониторинга окружающей среды. - Ставрополь: СГПУ, 1995.
4. Наац В. И. Обратная коэффициентная задача для уравнения турбулентной диффузии. // Вестник СевКавГТУ. - Серия физико-химическая. - 2003. - С. 89-94.
2
2
2
^ Наац В. И., Наац И. Э., Рыскаленко Р. А.
Параметризованные модели теории переноса в задачах экологического мониторинга.
5. Каргин Н. И., Рыскаленко Р. А. Применение вариационных методов в вычислительной модели уравнения Навье-Стокса //Вестник СевКавГТУ. - 2006. - № 3(7). - С. 22-26.
6. Наац В. И., Рыскаленко Р. А. Построение сеточной модели решения нестационарного уравнения переноса примесей в атмосфере и ее программная реализация // Математическое моделирование и информационные технологии в технических, естественных и гуманитарных науках: Материалы III региональной научно-технической конференции. - Георгиевск, 2003. - С.16-18.
7. Наац В. И. Численное моделирование пространственных задач переноса субстанции в пограничном слое атмосферы // Известия вузов. Северо-Кавказский регион. Естественные науки. - 2005. -№ 2(05). - С. 3-13.
8. Наац В. И. Система информационно - вычислительного обеспечения задач экологического мониторинга // Информационные технологии в науке, проектировании и производстве: Тр. XV Всероссийской научно-технической конференции. - Н. Новгород, 2005. - С. 5-6.
9. Рыскаленко Р. А. Численное моделирование поля концентрации примесей с заданными характеристиками в задачах переноса // Современные проблемы математики и естествознания: Материалы Всероссийской научно-технической конференции. - Н. Новгород, 2003. - С. 15.
10. Каргин Н. И., Рыскаленко Р. А. Многочлены Бернштейна и метод наименьших квадратов в вычислительной модели уравнения Навье-Стокса // Известия вузов. Северо-Кавказский регион Естественные науки. - 2007. - № 6.
Об авторах
Наац Игорь Эдуардович, Ставропольский государственный университет, доктор физико-математических наук, профессор кафедры информационных систем в экономике. Сфера научных интересов - математическое моделирование, численные методы, математическая физика. VINaac@yandex.ru
Наац Виктория Игоревна, Ставропольский государственный университет, доктор физико-математических наук, доцент, заведующая кафедрой математического анализа. Сфера научных интересов - математическое моделирование, численные методы, комплексы программ. VINaac@yandex.ru
Рыскаленко Роман Андреевич, Ставропольский государственный университет, кандидат физико-математических наук, доцент кафедры прикладной математики и информатики. Сфера научных интересов - математическое моделирование, численные методы, программирование, информационные технологии. risc-roman@yandex.ru