Моделирование и оптимизация использования топливного потенциала бытовых отходов на примере Коломенского района

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

Введение

Более десяти лет в нашей стране существует «Концепция перехода РФ на модель устойчивого развития» [1] в которой зафиксированы основные направления и принципы охраны окружающей среды и развития общества. Ключевым понятием концепции является «устойчивое развитие». Такое развитие общества удовлетворяет его насущные нужды без нанесения ущерба удовлетворению потребностей будущих поколений. Во главу угла ставятся вопросы охраны окружающей среды и сбережения ресурсов. При этом энергетическая стратегия России особое внимание отводит местным топливноэнергетическим ресурсам (ТЭР) регионов, к которым может быть отнесены определенные виды отходов. Таким образом, управление отходами можно позиционировать как сложную комплексную задачу, так как они серьезно ухудшают экологию и, в тоже время, требуют при переработке значительных затрат полезных дефицитных ресурсов (свободной площади, энергии и т.д.), хотя сами могут использоваться как ТЭР.

Современные особенности управления отходами

Задача организации эффективного управления отходами для больших территориальных образований в России стоит весьма остро и широко обсуждается в научно-техническом сообществе. При этом концептуальные основы управления отходами в отечественной практике проработаны достаточно полно [2]. Всё же необходимо отметить пробел в части обеспечения специалистов и экспертов программно-вычислительными средствами, разработанными с учетом особенностей самих систем, среди которых:

1. Пространственная сложность, обусловленная географической уникальностью каждого территориального образования, индивидуальными особенностями транспортной и энергетической инфраструктуры, а также сложившимся размещением источников отходов, мест их переработки и захоронения.
2.  Системная и динамическая сложность обусловленная тем, что объекты системы работают в непосредственной связи с многочисленными сферами жизнедеятельности территориального образования (промышленность, ЖКХ, автотранспорт, сельское хозяйство, энергетика), а факторы, влияющие на сферу управления отходами, подвержены временным колебаниям (суточным, месячным, сезонным) и, накладываясь, могут создавать экстремальные условия, которые не очевидны при неинструментальном анализе системы.
3.  Многовариантность путей развития и их сочетаний при ограниченном количестве ресурсов. Данная сложность обусловлена широким набором альтернатив (размещения объектов, выбора способов переработки и производства сопутствующих вторичных ресурсов и товарных продуктов, организации транспорта) из которых для практического воплощения необходимо выбрать наилучшее сочетание, при котором при минимуме затрат будут обеспечиваться экологическая безопасность и эффективность работы системы управления отходами.

Постановка задачи научного управлением отходами

В связи с указанными сложностями для эффективного решения задач развития систем и обеспечения требуемого научного планирования встала необходимость разработки специальной методики и реализации ее в инструментально-вычислительных средствах.
Сфера обращения отходов производства и потребления объединяет в себе несколько подсистем, обусловленных особенностями групп отходов, на которых они специализируются. К таким группам относятся: строительные, сельскохозяйственные, медицинские, твердые бытовые отходы (ТБО) и др. В настоящее время для большинства территориальных образований именно ТБО (или «муниципальные отходы») являются общей актуальной проблемой, требующей эффективного решения. Представленная в данной статье работа затрагивает именно ТБО [3].
Отдельные компоненты ТБО являются ценным вторичным сырьем, которое может быть при современном уровне технологии вовлечено (возвращено) в хозяйственный оборот. Важное направление использования - это реализация их топливного потенциала для генерации тепловой и электрической энергии, выработки квалифицированного топлива, а также комбинация этих методов с неэнергетическими способами переработки отходов. Так как важной стороной вопроса является согласование реализации полученной энергетической продукции (т. е. вовлечения потенциала отходов в энергетические балансы потребителей рассматриваемого территориального образования), было решено использовать применяемые в энергетике способы моделирования и оптимизации систем [4, 5], адаптировав их к особенностям управления отходами. При этом, обязательным рассматривалось дополнение их актуальными разработками зарубежных специалистов, работающих на поприще инструментальновычислительного обеспечения сферы управления отходами [6,7,8].

Методика управления отходами

Итогом представленной работы стали методический подход и реализованная на его основе информационно-аналитическая система [9,10,11], позволяющие на базе балансовой оптимизации, сценарного анализа и имитационных экспериментов создать эффективную информационную основу для научного обоснования среднесрочных (до 5 лет) и долгосрочных (до 20 лет) программ управления отходами для территориальных образований различного уровня. Основа предложенного подхода лежит на принципах системного анализа и заключается в создании адекватной модели исследуемой системы обращения отходов (математической, и далее - компьютерной), изучения ее свойств на сериях сценарных вычислительных экспериментов, а затем переносе результатов на реальные объекты территории.

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

Для реализации эффективного управления таким ТСОТБО был предложен алгоритм представленный на схеме (Рис. 1. ) [11].

Рис. 1. Схема алгоритма разработки научнообоснованных программ обращения отходов

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

Пространственная модель ТСОТБО - реализует подходы пространственного анализа и геомаркетинга [12] и учитывает географическое расположение и взаимное влияние объектов самой системы и ее окружения. Инструментально реализуется в геоинформационной системе ArcView 3.2.

Экономико-математическая модель ТСОТБО - служит для оптимизации структуры системы и выбора из множества альтернатив развития решений, инвариантных к изменяющимся внешним условиям, а также обеспечивает эффективное вовлечение энергетической продукции системы в топливно-энергетические балансы объектов. Модель реализована в систему управления базами данных (СУБД) Microsoft Office Access 2003 в купе с оптимизационным приложением (QSopt или LINDO).

Имитационная модель ТСОТБО - реализует подходы системной динамики [13,14], служит для учета системной и динамической сложности функционирования системы управления отходами, выявления её «узких» мест и контринтуитивного [20] поведения. Реализуется в пакете имитационного моделирования Powersim Studio 2005.

Пространственная модель формируется на этапах 1, 2 и 3 (Рис. 1. ), как географическая основа системы обращения отходов. Отмечаются объекты образования отходов. Далее, согласно существующих методик прогнозирования, определяется их мощность по образованию отходов (этап 2) и анализируются транспортные подсистемы: сеть автомобильных и железных дорог, речной транспорт. При построении модели учитывается территориальное административное деление (в городе - кварталы, в областном районе - сельские округа). Особенно подробно фиксируются в географической модели существующие объекты ТСОТБО, такие как: перегрузочные станции, предприятия рециклинга, компостирования, комплексные и мусоросжигательные заводы, полигоны. На основе перечисленной информации формируется расчетный граф геоинфор- мационной (ГИС) модели. Узлы графа представляют собой точки на территории к которым привязываются близлежащие объекты где образуются отходы, а также места размещения существующих или рассматриваемых к строительству объектов

ТСОТБО. Рёбрами графа являются транспортные связи, пригодные для перемещения рассматриваемых материалов и ТЭР.

Параллельно с формированием ГИС модели на этапах 2 и 3 происходит информационное наполнение БД, где отражается каждый источник отхода (этап 2) и каждый объект ТСОТБО (этап 3). Причем в БД описываются как текущие показатели, так и прогнозные на расчетный срок планирования. Прогнозы строятся на основе существующих методик [15, 16].

Таким образом, перед 4-м этапом имеется ГИС модель и реализованная в базах данных (БД) балансовая модель которые отражают: а) текущее и перспективное образование отходов, б) сложившуюся на момент планирования структуру ТСОТБО, в) перспективные возможности (альтернативы) развития системы.

На 4-м этапе с помощью разработанных в БД процедур проводится оптимизация структуры ТСОТБО на принятом промежутке планирования. Из всех возможных (допустимых) путей развития выбирается такой, при котором совокупные дисконтированные затраты из бюджета территориального образования на сферу обращения ТБО будут минимальны.

Структура БД и используемые в ней процедуры представляют задачу развития ТСОТБО в виде задачи линейного программирования [17]. Основным преимуществом такого подхода к поиску оптимального решения является: а) возможность решения задач большой размерности (до 100 тыс. переменных), б) возможность анализа чувствительности полученного оптимального решения к изменению внешних условий функционирования ТСОТБО, в) возможность анализа двойственной задачи линейного программирования [18]. Последнее наиболее важно, так как именно анализ двойственных оценок позволяет отвечать на вопрос: почему именно такая схема является оптимальной, а также позволяет отслеживать и анализировать «жизненный цикл отходов» [19] в рамках оптимальной структуры ТСОТБО.

В результате, в начале 5-го этапа имеется оптимальный план развития ТСОТБО с жестко согласованными материальными и энергетическими потоками между источниками отходов, транспортной схемой, местами переработки и конечного размещения отходов. При этом реализуемый в определенные виды ТЭР топливный потенциал отходов привязывается к конкретным потребителям. Данное решение отражается на электронной карте в ГИС- модели ТСОТБО.

На пятом этапе оптимальная структура ТСОТБО транслируется в динамическую модель ТСОТБО, в которой проводится серия вычислительных экспериментов. Полученная при балансовой оптимизации в рамках дискретных шагов планирования (трех-, пятилеток) структура ТСОТБО проверяется на модели с непрерывным временем с шагом дискретизации в 1 день. При этом в динамической модели снимаются упрощения балансовой модели и ТСОТБО исследуется как сложная система [20], с учетом обратных связей [21] с добавлением информационных потоков и экономических сценариев. Таким образом, оптимальная структура ТСОТБО проверяется на динамическую устойчивость, выявляются ее «узкие места» и системные особенности поведения. В результате имитационных экспериментов получаются данные о сложных материальных, энергетических и денежных потоках между системой и ее окружением и влияние системы на окружающую среду с учетом накопительного эффекта. Полученная информация с пространственной привязкой в ГИС-модели представляет собой исчерпывающую основу для разработки долгосрочной научно-обоснованной программы развития ТСОТБО с учетом использования их топливного потенциала.

При различных вычислительных экспериментах при варьировании внешних условий часть решений, как правило, будет повторяться вне зависимости от текущего сценария. Такие инвариантные решения на шестом этапе сводятся в блок первоочередных мероприятий программы развития ТСОТБО. Именно такой подход обеспечит возможность управления ТСОТБО по принципу опережения, как этого требуют современные экологические реалии [22].

Использование топливного потенциала бытовых отходов коломескгого района

На сегодняшний день, по согласованию с членами рабочей группы по координации обращения ТБО Министерства экологии и природопользования Правительства Московской области, Коломенский район выбран для апробации разработанных инструментальных средств и подготовки с их помощью пилотной «Научно-обоснованной районной программы обращения твердых бытовых отходов до 2020 года с учетом эффективного использования их энергетического потенциала». Далее приведены расчетные данные, полученные при работе с моделью обращения ТБО Коломенского района.

Рис. 2. Географическая основа

Модель коломенского района реализована в геоинформационной системе. Географическая карта и полученная после первых трех этапов алгоритма ГИС модель ТСОТБО приведены на рисунках 2 и 3.

После реализации 4-го этапа полученная при оптимизации схема анализируется с помощью сводных диаграмм материальных и энергетических балансов локально для узлов модели и агрегировано для всей территории в целом. Начинается анализ решения со сводной диаграммы дисконтированных денежных потоков ТСОТБО (Рис. 4. ). В ней для каждого периода планирования (для каждой трехлетки из всего 15-летнего промежутка с 2006 по 2020 г.) отражаются денежные потоки, разделенные по трем категориям: процессы, транспорт и запасы. На приведенной диаграмме видно, что основные расходы приходятся на развитие мощностей по энергетической переработке отходов (мусоросжигательный завод с утилизацией тепла), и на строительство нового санитарного полигона в районе существующего полигона «Воловичи» (узел 209), который исчерпывает в этот 2-ой период свои емкости.

Рис. 3. Выделенные источники ТБО, плотность образования и транспортная схема с узлами

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

После анализа поведения целевой функции, которая и представляет собой совокупные затраты, производится анализ материально-энергетических балансов. Сначала агрегированных (Рис. 6. ) а потом локальных по узлам (Рис. 5. ). Локальные узловые материально-энергетические балансы представляют собой системы показателей прихода и расхода в каждую статью баланса участвующими в ней процессами, транспортом или запасами. Проанализировать представленный блок массового баланса для узла 209 для 2-го периода планирования (2009-2011 гг.), можно проследив движение по стрелкам (Рис. 5. ).

В первом столбце справа указаны статьи баланса в сокращенных названиях («Отх-ТБОгор» - городские ТБО; «Отх-ТБОсел» - сельские ТБО; «Отх-ОФ-гТБО30» - остаточная фракция городских ТБО, после рециклинга на глубину 30%, т.е. извлечения сопутствующих товарных фракций, представленных балансом «Спт-Ф-гТБО30» и т.д.). В следующем столбце указаны знак и числовое значение правой части баланса [17]. Для отходов правые части представляют собой жесткие равенства («=0»), что означает, что в увязанном балансе оптимальной структуры ТСОТБО приход каждого вида отхода (ввоз, сбор, образование при переработке других видов отходов, извлечение из запасов предыдущего периода) должен равняться в узле его расходу (вывозу, переработке, запасу на следующий период). Для «полезных» балансов (например сопутствующая фракция черного метала, при рециклинге) правые части балансов представляют собой неравенства вида «>0», что означает, что приход в такой баланс может быть больше чем расход. Далее на левой части диаграммы следует столбец со значением дополнительных переменных линейной задачи [18], которые показывают значения избытка или дефицита баланса в зависимости от условия задачи (значения правой чисти баланса). На примере приведенного локального баланса (Рис. 5. ) такая дополнительная переменная для сопутствующей фракции утильных компонентов рециклинга городских ТБО («Сп-Ф-гТБО30») равна 42,24 тонны и означает массу этих отходов накопленную в рассматриваемом узле модели №209 за приведенный 2ой период планирования (2009-2011гг.).

Рис. 4. Дисконтированные денежные потоки территориальной системы обращения бытовых отходов, тыс. руб.

Рис. 5. Пример анализа материального баланса узла 206 второго периода планирования (2009-2011 гг.)

Следующий столбец - двойственная оценка баланса, который означает во сколько реально обходится 1 тонна для ТСОТБО в рамках текущей оптимальной структуры. Например, такая оценка для городских ТБО в балансе узел (см. «Отх- ТБОгор») равна «-1,94 руб./тонна». Это значит, что если значение баланса изменить на 1 тонну (т.е. присвоить правой части значение «=1», вместо «=0», и таким образом разрешить 1 тонне ТБО остаться в узле и не направлять на дальнейшую переработку), то совокупные дисконтированные издержки (т. е. значение целевой функции), уменьшится на 1,94 руб. Принципиальной особенностью этого показателя является то, что двойственная оценка (или «теневая цена» [17]), аккумулирует в себе затраты и прибыли по всему «жизненному циклу» данного вида отхода учитывая сбор за прием отхода, затраты на транспортировку до этого узла, затраты на необходимые мощности рециклинга, которые влечет за собой эта 1 тонна, дальнейшие затраты на установленную мощность мусоросжигательного завода, и затраты на захоронение шлака, после сжигания, а также прибыль после реализации сопутствующих фракций материальных продуктов и ТЭР.

Следующие три столбца относятся к подробному представлению составляющих приходных и расходных частей баланса. Для процессов эти три столбца содержат последовательно: класс процесса - тип процесса - режим работы процесса. Для транспорта: номер узла источника потока (или номер узла получателя) - класс транспорта - вид транспорта.

На приведенном балансе для 209 узла на период 2009-20011 гг. по стрелкам можно проследить последовательность анализа «жизненного цикла» в узле потока городских ТБО. В приходную часть баланса (стрелка 1) ввозятся 140,8 тыс. тонн из узла 41 грузовым автомобильным транспортом (см. отражение транспортных потоков в ГИС-модели Рис.7). Ввезенный поток отходов полностью расходуется комплексом рециклинга (стрелка 2), который относится к классу процессов «Механическая и биологическая переработка отходов» (на диаграмме «Отх-МБ-пер.»). Данный процесс параллельно участвует в приходных частях двух других материальных балансов. В баланс сопутствующей фракции рециклинга («Спт-Ф-гТБО30») он выдает 42,24 тыс. тонн (стрелка 3), а в баланс остаточной фракции рециклинга («Отх-ОФ-гТБО30»), т.е. оставшийся поток отходов после рециклинга 98,56 тыс. тонн (стрелка 4). Далее, так как остаточная фракция относится к отходам, и баланс представляет жесткое ограничение («=0»), то остаточная фракция перерабатывается на мусоросжигательном заводе (на диаграмме: «МСЗ с утил. тепла, нов. 200 тыс. т/год»), который относится к классу процессов «Использование топливного потенциала отходов» (на диаграмме: «Отх-в-ТЭР») и перерабатывая пришедшие 98,56 тыс.тонн остаточной фракции выдает в баланс шлака 18,88 тыс. тонн (стрелка 6), что вместе со шлаком остальных режимов МСЗ, составляет 22,45 тонн (стрелка 7), которые размещаются новом санитарном полигоне.

Таким образом, анализируются все материальные и энергетические балансы для каждого из 42-х узлов модели и для каждой из принятых пяти трехлеток планирования (с 2006-2020 гг.).

Рис. 6. Сводный массовый баланс системы за промежуток планирования (15-лет), тыс. тонн.

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

При этом частный материальный баланс для конкретного узла и конкретного периода выглядит так, как представлено на рисунке (Рис. 5.)

Рис. 7. Вариант распределения транспортных потоков по видам отходов для периода одного из сценариев

В результате оптимизационных расчетов использование топливного потенциала для выработки электрической и тепловой энергии в структуре ТСОТБО модели имеет место. Прежде всего это обусловлено малой емкостью существующего полигона «Воловичи» (500 тыс.м3), который полностью заполняется за два первых периода планирования (с 2006 по 2011 гг.), а также тем, что энергетическому использованию ТБО путем сжигания в топках котлов противопоставляется альтернатива - санитарный полигон, затраты на который (денежные, материальные и энергетические) взяты с учетом самых высоких экологических требований. В результате, в оптимальной структуре часть отходов сжигается с выработкой энергии. Однако предварительно поток отходов подвергается рециклингу на глубину 30%. Тем самым лишь оставшиеся по массе 70% отходов подлежат дальнейшей переработке, что пропорционально снижает затраты на установленные мощности таких капиталоемких объектов как МСЗ и санитарный полигон. Также рассмотренный в модели вариант приготовления специального топлива из отходов RDF (Residue Derived Fuel [23]) с высокой теплотой сгорания не нашел применения, в следствии того, что остаточная фракция после приготовления такого топлива была отнесена к отходам (а не к товарным фракциям, как в случае ре- циклинга), что потребовало дальнейших затрат на ее переработку.

Заключение

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

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

Пример использования методологии представлен для модели ТСОТБО Коломенского района. В процессе расчетов на 15-летнюю перспективу:

•  определена оптимальная структура системы с обоснованием транспортных потоков;
•  определены требуемые объемы финансирования по этапам реализации программы;
•  обоснованы временные рамки ввода в действие новых предприятий;
•  определен потенциал возможной экономии ТЭР при использовании топливного потенциала отходов в рамках обоснованной оптимальной структуры ТСОТБО.

М.А. Ушаков, Д.А. Николаев, Т.А.Степанова
Московский энергетический институт (технический университет),
Россия Институт проблем энергетической эффективности (ИПЭЭФ),
кафедра ЭВТ

Список литературы

1. Концепция перехода РФ на модель устойчивого развития. - Утв. Указ. Презид. РФ от 01.04.96 N 440
2. Проблема твердых бытовых отходов: комплексный подход/ Эколайн,ЕС0Ь001Л 1996 http://www.ecoline.ru/mc/books/tbo
3. ГОСТ 30772-2001 Ресурсосбережение. обращение с отходами. Термины и определения. - Введ. 01.07.2002 - М.,2002. - 20 с.
4. Информационная технология исследования развития энергетики/ Л. Д. Криворуцкий, Л.В. Массель. - Новосибирск: Наука. Сиб. издательская фирма РАН, 1995. - 160 с.
5. Оптимизация развития и функционирования автономных энергетических систем/ А.М. Клер, Н.П. Дека- нова, Б.Г. Санеев и др. - Новосибирск: Наука, 2001. - 144 с.
6. A Multiobjective Approach for Solid Waste Management/ R. Minciardi, M. Paolucci, M. Robba, R. Sacile
7. A Simulation System for Waste Management - From System Dynamics Modelling to Decision Support/ P.A. Wager and L.M. Hilty (www.iemss.org/iemss2002)
8. An Optimal Model and Its Application for the Management of Municipal Solid Waste from Regional Small Cities in China/ ISSN 1047-3289 J. Air & Waste Manage. Assoc. 54:191-199
9.  Степанова Т.А. Ушаков М. А. Методы эффективного вовлечения энергонесущих промышленных и бытовых отходов в топливно-энергетические балансы районов/ Вестник МЭИ. - 2005. - №1. - С. 32-39.
10. Степанова Т.А. Ушаков М.А. Экономикоматематическая модель для обоснования вовлечения энергонесущих промышленных и бытовых отходов в топливно-энергетические балансы предприятий и населенных пунктов/ Информационные технологии и математическое моделирование: Мат-лы 3-й всерос. науч.- приктич. конф. (11-12 дек.2004г.) - Томск: Изд-во Том. ун-та, 2004. - Ч.2 - С.156-159
11. Ушаков М. А. Степанова Т.А. Разработка методологии и инструментальных средств обоснования вовлечения энергонесущих промышленных и бытовых отходов в топливно-энергетические балансы районов и населенных пунктов/ Радиоэлектроника, электротехника и энергетика: Тез.докл. 11-ой Международ. науч.-тенич. конф. студ. и аспир. (Москва, 1-2 марта). - М.: ЗАО «Знак», 2005. - Т.2. - С. 374-375
12. Цветков В.Я. Геомаркетинг: прикладные задачи и методы. - М.: Финансы и статистика, 2002. - 240 с.
13. Имитационное моделирование производственных систем/ Под общ. ред. А. А. Вавилова. - М.: Машиностроение; Берлин: Техника, 1983. - 416 с.
14. Технология системного моделирования/ Е.Ф. Аврамчук, А.А.Вавилов, С.В.Емельянов и др.; Под общ. Ред. С.В. Емельянова и др. - М.: Машиностроение; Берлин: Техник, 1988. - 520 с.
15. Сборник удельных показателей образования отходов производства и потребления. - Утв. 07.03.1999. - М: НИЦПУРО при Минэкономике и Минприроде России, 1999 - 36 с.
16. Методические рекомендации по оценке объемов образования отходов производства и потребления. - М: ГУ НИЦПУРО, 2003
17. Таха, Хэмиди, А. Введение в исследование операций, 6-е издание.: Пер. с англ. - М.: Издательский дом Вильямс, 2001. - 912 с.
18. Вентцель Е.С. Исследование операций. Задачи, принципы, методология. Учеб.пособие для студ. вузов. - 2-е изд., стер. - М.: Высш.шк., 2001. - 208 с.
19. Об отходах производства и потребления в Московской области - Закон Московской области от 8 ноября 2001 г. № 171/2001-03
20. Жилин Д.М. Теория систем: опыт построения курса.- 2-е изд., испр. - М.: Едиториал УРСС, 2004. - 184 с.
21. Информационные технологии в управлении качеством среды обитания: Учеб. пособие для студ. высш. пед. учеб. заведений/ В. Е. Гершензон, Е. В. Смирнова,В.В. Элиас; Под ред. В.Е. Гершензона. - М.: Издательский центр «Академия», 2003. - 288 с.
22. Экологические проблемы Московской области и пути их решения/ Качан А.С. // Журнал «Устойчивое развитие. Наука и Практика». - 2003. - №3. - с. 25-28
23. Шершнев Е.С., Ларионов В.Г., Куркин П.Ю. Сжигание ТБО с получением электроэнергии и теплоты // Экол. и пром-сть России. - 1999. - Апрель. - С.44-47