Интенсивное развитие различных отраслей промышленности, энергетики, автомобильного и железнодорожного транспорта, рост численности населения и другие факторы оказывают неблагоприятное воздействие на окружающую среду, в частности на состояние водных экосистем. В настоящее время для изучения масштаба загрязнений водных объектов все чаще опираются на широкий круг математических методов. Эффективность использования того или иного метода возрастает при правильном понимании сущности явлений и процессов, протекающих в исследуемом объекте, а также в ходе применения системного анализа. Однако, проблема точного математического описания динамики распространения загрязнений из локальных источников по объему ограниченных резервуаров по сей день слабо разработана [1-3]. Это объясняется тем, что данная задача сопряжена с определенными математическими трудностями как в постановке граничных условий, так и в методах аналитического решения.
На основе метода локальных источников, предложенного в работе [4], разработана методика моделирования динамики распространения локализованных примесей по ограниченной области. Обоснование идеи и адаптация этого метода к данному процессу заключается, на наш взгляд, в математической аналогии диффузионного процесса и процесса растекания жидкости по насадочному слою [5].
При использовании метода локальных источников граничные условия формируются на основании тех или иных физических соображений. В реальных условиях размеры источника загрязнения намного меньше, чем размеры исследуемых объектов, поэтому считаем, что на бесконечности концентрация загрязняющего вещества равна нулю. Начальная концентрация задана так: вне источника равна нулю, а внутри источника равна единице.
В данной статье приведены результаты решения двумерной диффузионной задачи, характеризующей распространения загрязнений поверхностного слоя водного объекта с применением метода локальных источников.
Были введены следующие обозначения: 2В -длина области, м; 2А -ширина области, м; 2k, 2n - размеры области начальной локализации примеси по длине и ширине соответственно, м; s - расстояние от центра источника загрязнения до середины области по длине, м; h - расстояние от центра источника загрязнения до середины области по ширине, м; t - период наблюдения; D - коэффициент диффузии.
Графически данная задача может быть представлена в виде рис. I.
Отметим, что оси координат совпадают с серединой источника загрязнения.
В результате расчетных операций в среде Mathcad 13 с учетом вышеприведенных исходных данных получили, что при t = 1605 с неравномерность распределения не превышала 5%.
Далее задачу усложняли, добавляя еще один источник загрязнения. Была изучена методика расчета концентрации примеси в зависимости от расположения этих двух источников загрязнения. Источники загрязнения локализованы в разных областях водной поверхности, как показано на рисунке 3 (а,б).
С добавлением второго источника загрязнения возникает необходимость введения следующих обозначений: k, k1, k2, k3, n, n1, n2, n3 , определяющие размеры области начальной локализации примеси обоих источников по длине и ширине соответственно, м; S1 - расстояние от центра второго источника загрязнения до середины области по длине, м; Һ - расстояние от центра второго источника загрязнения до середины области по ширине, м.
Каждый из случаев, представленных на рисунке 3 имеют свои различия и схожести. Первоначально была рассмотрена задача представленная на рисунке 3 а. При расчете общей концентрации примеси для второго источника загрязнения следует учесть его размеры, расстояние, отделяющее эти два источника, их расположение относительно середины исследуемого водного пространства и т.д.
И вновь, используя метод локальных источников [5, с.125-128], найдем концентрацию примеси для второго источника. Вычислительный процесс для двух источников проводится в одной системе координат.
Определить границы интегралов, характеризующих размеры второго источника по длине (k2, k3) и ширине (n2, n3) достаточно несложно. Трудности возникают при нахождении этих границ по оси абсцисс, то есть значений k2, k3.
Таким образом, получили картину рассеяния примеси на поверхностном слое водного объекта при временах: 10 с, 15 с, 100 с, 200 с, 300 с, 400 с и следующих входных данных: В=15, А=8, k=3, n=2, k2=-17, n2=-2, k3=-9, n3=2, s=7, s1=6, h= h1=5,D=0,3 (см. рисунок 4).
Концентрация примеси для исследуемой водной области складывается из концентраций примеси первого и второго источников. В результате проведенных вычислений с вышеприведенными исходными данными при времени распространения t = 4365с (что составляет 1ч 12 мин 45 с), неравномерность распределения не превышала 5 %.
Затем рассмотрен случай, когда источники загрязнения были локализованы в различных областях исследуемой поверхности, что и представлено на рисунке 3б. Отметим, что как и в предыдущем случае, оси координат совпадают с серединой источника загрязнения.
После нахождения суммарной концентрации загрязняющего вещества необходимо вновь обратиться к вышеприведенным заменам, что позволяет вернуться к начальному условию, связанному с расположением источника примеси.
На рисунке 5 представлены графики, характеризующие распространение примеси в водной среде для двумерного пространства при наличии двух источников загрязнения, расположенных на разных уровнях, при временах:
10 с, 15 с, 100 с, 200 с, 300 с, 400 с и следующих исходных данных: В=15, А=8, k=3, n=2, k2=-17, n2=-11, k3=-9, n3=-7, s=7, s1=6, h=5, h1=4,D=0,3.
Границы интеграла при расчете S2* подобраны таким образом, что A + Һ определяет границу с верхней стороной длины области, а -(A - характеризуют нижнюю границу длины области по оси ординат. -(B-s1) и B+s1 определяют соответственно левую и правую боковые границы объекта относительно оси абсцисс.
Поправка была рассчитана при временах: 10 с, 15 с, 100 с, 200 с, 300 с, 400 с и исходных данных: В=15, А=8, k=3, n=2, k2=-17, n2=-11, k3=-9, n3=-7, s=7, si=6, h=5, hi=4,D=0,3. Результаты приведены в таблице 2.
расчетных операций, с учетом вышеприведенных исходных данных, получили: при t = 9090с, что соответствует 2 часам 31 минуте 30 секундам, неравномерность распределения не превышала 5 %.
Результаты численных экспериментов, а также построенные на их основе графики, физически обосновывают исследуемый процесс и подтверждают целесообразность применения вышеизложенной методики при наличии одного или двух источников загрязнения для анализа состояния поверхностного слоя водного объекта. Предложенный в работе подход позволяет рассчитать концентрацию загрязняющего вещества за определенный период времени, а также наглядно представить картину распространения загрязняющего вещества при определенных исходных данных.
ЛИТЕРАТУРА
- Полубаринова -Кочина П.Я. Теория движения грунтовых вод. - М.: Наука, 1977. - 664с.
- Шакиров Б.С., Куракбаева С.Д. Основные задачи и особенности математического моделирования катастрофических явлений в экологии // Наука и образование Южного Казахстана. - 2006. - №3 (52). - С. 74- ТТ.
- Никаноров А.М., Тепляков Ю.В., Уразаев К.К. Подобие моделей экосистем водных объектов: критерии подобия и способы их оценки // Водные ресурсы. - 1996. - Т.23, №5. - С.578-582.
- Brener A.M. Adaptation of random walk methods to the modeling of liquid distribution in packed columns // Advances in Fluid Mechanics IV, 2002, WIT Press, P. 291-300.
- Болгов Н.П., Майлибаева Л.И., Соколов Н.М., Бренер А.М., Тарат Э.Я., Еримбетов А.К. Распределение жидкости в регулярной полочной насадке, ограниченной стенками // Теоретические основы химической технологии. - 1979, С. 125-128.
- Тихонов А.Н., Самарский АА. Уравнения математической физики., Изд. 5-е. - М.: изд. «Наука», 1977. -736с.
- Куракбаева С.Д. Методика расчета распространения загрязнений в водном объекте // Сборник статей XIV международной научно-практической конференции «Экология и жизнь». - Пенза, 2008. - С.168- 171.