Учет процессов массопереноса загрязняющих компонентов в системах автоматизации проектирования водозаборов подземных вод Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

УДК 628.1:658.5

UDC 628.1:658.5

УЧЕТ ПРОЦЕССОВ МАССОПЕРЕНОСА ЗАГРЯЗНЯЮЩИХ КОМПОНЕНТОВ В СИСТЕМАХ АВТОМАТИЗАЦИИ ПРОЕКТИРОВАНИЯ ВОДОЗАБОРОВ ПОДЗЕМНЫХ ВОД

Стародубцев Виктор Сергеевич

д.т.н., к.г.-м.н., профессор

Российский государственный социальный

университет филиал в г.Воронеж, Воронеж,

Россия

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

Ключевые слова: ПРИРОДНО-ТЕХНИЧЕСКАЯ СИСТЕМА, ВОДОЗАБОР ПОДЗЕМНЫХ ВОД, МАССОПЕРЕНОС, СТРУКТУРНАЯ ИДЕНТИФИКАЦИЯ, СИСТЕМЫ АВТОМАТИЗАЦИИ ПРОЕКТИРОВАНИЯ

ACCOUNTING OF PROCESSES OF WEIGHT CARRYING OF POLLUTING COMPONENTS IN SYSTEMS OF AUTOMATION OF DESIGNING OF WATER INTAKES OF UNDERGROUND WATERS

Starodubtsev Victor Sergeevich Dr.Sci.Tech, Cand.Geol-Min.Sci., professor Russian state social university branch in Voronezh, Voronezh, Russia

In the article, application of methods of self-organizing for identification of processes of weight carrying polluting components of underground waters for the purpose of optimization of technogenic loading on environment is considered

Keywords: NATURAL-TECHNICAL SYSTEM, WATER INTAKE OF UNDERGROUND WATERS, WEIGHTCARRYING, STRUCTURAL IDENTIFICATION, SYSTEMS OF AUTOMATION OF DESIGNING

Развитие территориально-производственных комплексов (ТПК) структурно порождает развитие промышленного, транспортного, пищевого и других производств. В свою очередь, промышленное производство в условиях превалирования экономических законов над экологической безопасностью вызывает расширение и интенсификацию загрязнения геосфер Земли. Источниками загрязнения стали практически все промышленные предприятия, транспорт, зоны отдыха, крупные сельскохозяйственные и животноводческие комплексы. Наиболее подвержена антропогенному воздействию гидросфера Земли, где техногенное воздействие может сказываться как за счет прямого загрязнения промышленными, коммунальными, поверхностными стоками, загрязнения нефтью и нефтепродуктами, которое стало обычным явлением для мест их хранения и продажи (нефтебазы, автозаправочные станции и т.д.), так и за счет вторичного комплекса загрязнений гидросферы, формирующегося на основе загрязняющих веществ, которые содержатся в атмосферных осадках. В связи с этим необходим системный подход к http://ej.kubagro.ru/2012/03/pdf/30.pdf

развитию природно-техногенных процессов в природно-технических системах (ПТС) ТПК.

Рассмотри реализацию системного подхода на примере идентификации процессов массопереноса загрязняющих компонентов подземных вод в ПТС инфильтрационного водозабора подземных вод (ВПВ) №11 г.Воронежа.

ВПВ № 11 находится в 16 км севернее промышленного центра на правобережной пойменной террасе. Водозабор состоит из 48 эксплуатационных скважин (ЭС), расположенных вдоль берега водохранилища (рисунок 1). Эксплуатируемый неоген-четвертичный водоносный комплекс представлен мелкозернистыми песками, переходящими в основании в крупнозернистые и гравелистые разности с галькой. По химическому составу воды ВПВ №11 пресные с минерализацией 0,25-0,4 г/л, гидрокарбонатного кальциево-магниевого типа.

Общий дебит ВПВ №11 в настоящее время составляет 201 тыс.м3/сут.

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

о

3

Рисунок 1. ПТС ВПВ №11 с участками погребенных пойменных

отложений

Такая ситуации вызвала необходимость детального изучения процессов массопереноса железа и марганца в системах водозаборов средствами математического моделирования процессов массопереноса загрязняющих компонентов подземных вод. Для этого были разработаны Мр:/^ .kubagro.ru/2012/03/pdf/30.pdf

методика и алгоритм структурной идентификации процесса массопереноса загрязняющих компонентов подземных вод [1], на базе которых было создано программное обеспечение - комплекс программ МАБРЕКЕКОБ [1].

С учетом трудностей решения уравнения геомиграции [3] предлагается использовать методы структурной идентификации и в частности метод группового учета аргументов (МГУА) [2]. Учитывая один из основных принципов теории МГУА - принцип "свободы выбора решений", дифференциальное уравнение конвективно-диффузионного переноса мигранта в двумерном потоке [3], используем полное описание класса структур для идентификации процесса массопереноса загрязняющих компонентов подземных вод вида

Эс

дг

■ = а

+ а2

Э у2

+ а3

( Эс V--

+ а4

( Эс V --

д У

+ (1)

+ а5 ¥ \~к + а6 ¥ 12к + а7 ¥ 1-к + а8 ¥ 4-к + а 9 ¥ 5-к + а 10 ¥ 16к + а11,

где с - концентрация ионов загрязняющего компонента (например, железа или марганца) в подземных водах (прогнозируемая переменная в мг/л); а1 -а4 - соответствующие коэффициенты при производных; ? - время; х, у -пространственные координаты; к - запаздывание по времени, к = 1, 2, 3;

3 0

¥ - водоотбор в тыс.м /сут; ¥2 - температура воздуха в С ; ¥3 - осадки в мм; ¥4 - рН поверхностных вод; ¥5 - содержание О2 в поверхностных водах водохранилища мг/л; ¥6 - содержание ионов хлора в поверхностных водах водохранилища мг/л; а5 - а10 - соответствующие коэффициенты при ¥ 1 -¥6; а11 - свободный член.

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

Блок А

1. Ввод данных.

2. Построение сетки на области моделирования.

Блок В

3. Формирование модели, исходя из полного описания (1).

4. Нормирование и центрирование данных.

хп - Хс,

X ,, = -1---^,

1 в

с-

где X1,, - текущее значение переменной;

- = 1, ... , К, где N - длина выборки;

1 = 1, ... , К, где К - количество переменных;

Хс - среднее значение исходной переменной;

Бс - дисперсия исходной переменной.

5. Определение зависимых величин с помощью коэффициента парной корреляции (Я,).

1 =

м 2

1 Х у,- Х У,}

^-----------, 1 *

м2 -м1 +1

1 1=-,

где М1 - начало выборки, а М2 - конец выборки, 1, ] = 1, 2, ..., п, где п -количество переменных взятых для проведения эксперимента по идентификации целевых функций изучаемых процессов.

Блок С

6. Задание частного описания модели, исходя из полного описания (1).

7. Определение коэффициентов модели по методу наименьших квадратов.

Блок Б

8. Определение значения критерия несмещенности (псм) модели.

I (г а - Гб ) 2

-=1

п см =

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

Блок Е

9. Определение значения критерия сходимости (1) для N1 лучших по критерию несмещенности моделей.

£ (У ш -У,; У

I = -1

П 7 ,

£ у,

, = 1

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

Блок Б

10. Определение значения критерия эпигнозного прогноза (Р) для N2 лучших по критерию сходимости моделей.

£ У -у,, )2

р _ 1=Б+1________

£ 7 .

£ у,,

1=Б+1

Критерий эпигнозного прогноза позволяет отбирать лучшие модели прогностической направленности.

Блок О

11. Определение комбинированного критерия К http://ej.kubagro.ru/2012/03/pdf/30.pdf

К* =

(

+

П

V тах J

где птах и 1тах - максимальное значение критериев минимума смещения и сходимости.

12. Выбор оптимальной модели по совокупности критериев.

В связи с большой протяженностью ВПВ №11 (около 4-х км) была проведена квантификация ПТС ВПВ №11 на подсистемы с целью выявления особенностей процесса массопереноса ионов железа и марганца. При квантификации ПТС ВПВ №11 учитывалось то, что в середине 80-х годов производилось расширение ВПВ №11 за счет намыва берега. Всего было выделено 3 подсистемы: 1 подсистема - ЭС № 14, 17 и 18; 2 подсистема - ЭС № 28, 29 и 45; 3 подсистема - ЭС № 9, 10, 11 и12.

В результате проведения эксперимента по структурной идентификации процесса массопереноса ионов марганца в подземных водах ПТС ВПВ №11 были получены следующие прогностические модели (таблица 1).

Таблица 1 - МОДЕЛИ ПРОЦЕССОВ МАССОПЕРЕНОСА ИОНОВ МАРГАНЦА В ПОДЗЕМНЫХ ВОДАХ ПТС ВПВ №11

П

Система Модель

Общая с^1 = с\_ + 0,071( с\+1]- с\_1} ) - 0,000049Ч]'1 + + 0,008Ч 3 _ 0,129Ч 4-1 _ 0,011 Ч 6 + 6,857

Подсистема №1 с+ = + 0,392( с] _ 2с] + с\-]+1 ) _ 1,527( с11+и_ с\_1} ) + + 0,000003 Ч1-1 _ 0,242

Подсистема №2 с+ = + 0,051( с\+1] _ 2с{} + си,} ) + 0,131( с] _ 2с] + с‘,]+1 ) _ + _ 0,291( с]- с\]} _1 ) _ 0,000003 Ч\-1 + 0,002 Ч 3 _ 0,002 Ч 6 + 0,31

Подсистема №3 с+ = _ 0,8( с+1] _2<] + си,] ) + 1,35( с+1] _2с] + с^ ) _ _ 1,53 ( с*;?] _ 2с] + с]) _ 2,18 ( с]] _ с\-_1 )-0,87( ]_] ) + + 0,958( сс] ) _ 0,000034 Ч1 _ 0,04 Ч 6 + 5,54

Анализ выражения модели общей системы показывает, что массоперенос ионов марганца в ПТС ВПВ №11 связан с миграцией

марганца по площади депрессионной воронки, на что указывает первая производная по Х, присутствующая в модели с нулевым запаздыванием, что соответствует предварительным выводам о местном источнике загрязнения. Результаты моделирования показывают, что процесс массопереноса ионов марганца для различных подсистем ВПВ №11 имеет много похожего. Для всех подсистем характерно наличие производных, как по оси Х, так и по оси У, что свидетельствует о миграции ионов марганца непосредственно во внутренних областях депрессионной воронки ВПВ и подтверждает версию о наличии источника загрязнения непосредственно в этой области (погребенные пойменные отложения в результате намыва береговой зоны). В то же время возможен и массоперенос ионов марганца из водохранилища (наличие в моделях параметров ¥4 и ¥6 поверхностных вод водохранилища). Следует отметить у всех моделей незначительное (-1) запаздывание по времени у параметра водоотбора ( ч{-1), что может также свидетельствовать, что загрязнение находится непосредственно в районе депрессионной воронки.

Для получения прогностической модели массопереноса ионов железа принятые модели были оценены по критерию эпигнозного прогноза и сценарному критерию. С учетом полного описания (1) модели процесса массопереноса ионов железа в ПТС ВПВ №11 представлены в таблице 2.

Анализ выражения для общей системы показывает, что массоперенос ионов железа в ПТС ВПВ №11 связан с миграцией железа по площади депрессионной воронки, на что указывает вторая производная по Х, присутствующая в модели с запаздыванием (-2), что соответствует значительным размерам ПТС ВПВ №11.

Таблица 2 - МОДЕЛИ ПРОЦЕССОВ МАССОПЕРЕНОСА ИОНОВ ЖЕЛЕЗА В ПОДЗЕМНЫХ ВОДАХ ПТС ВПВ №11

с1+1 = с1-1 + 0,0077 гс]]2-с]~2 ) + 0,000031 ¥1 + 0,407 Ч4-1- 6,914.

^ ' Ч+1 1-1~1 7 ’ 1 ’ 4 ’

Система Модель

Общая с+ = с]- + 0,151^ с\~+2,}-2с] + с12 )-0,000021¥ 1 + + 0,071'¥ 2-2 + 0,1 ¥ 5 + 1,82

Подсистема №1 с+ = <1 - 0,834( с\+1]- с1и ) + 4,633( с^- с^ ) + + 0,00008 ¥ 1-1 + 0,1 ¥ 6 -12,18

Подсистема №2 с] = сЩ + 0,0077( ] - с]]]2_1 ) + 0,000031' 1 + 0,407¥ 1-1 - 6,914

Подсистема №3 с+ = с1-] -1,06( с+1,]-2с!,] + с]-,,] ) + 0,84( с] - 2с]] + с\]+1 ) + + 0,64( с];1]- сЦ] ) + 0,37( ]-1 ) - 0,00057¥ 1 --1,16 ¥1]1 - 6,41 ¥ + 126,5

Наличие в модели температурного параметра Ч2-2 подтверждает сделанные предварительные выводы о влиянии температуры воздуха на содержание ионов железа в подземных водах ВПВ №11. Влияние Воронежского водохранилища отражено в модели параметром содержания

О2 в поверхностных водах водохранилища.

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

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

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

С учетом выявленных особенностей массопереноса ионов марганца и железа в ПТС ВПВ №11 рекомендуется при расширении системы ВПВ создать открытые инфильтрационные сооружения (два бассейна шириной по дну 20 м, длиной 300 метров и глубиной бассейна 3 м) для пополнения запасов грунтовых вод на расстоянии 150-200 метров от уреза водохранилища с размещением 12 эксплуатационных скважин, что обеспечит как качество питьевой воды (удаленность от водохранилища и погребенных по,йменных отложений), так и значительное (до 36 тыс.м3/сут) увеличение ее добычи.

Литература

1. Жуков С.А. Моделирование процессов массопереноса загрязняющих компонентов подземных вод //Экология и промышленность России. № 7. 2008. № 7. С. 24 - 27.

2. Ивахненко А.Г. Индуктивный метод самоорганизации моделей сложных систем.-Киев: Наук. думка, 1982. 296 с.

3. Шестаков В.М. Динамика подземных вод.-М.: Изд-во МГУ, 1979.-368 с.