Современные сооружения биохимической очистки (БХО) сточных вод крупных промышленных предприятий и населенных пунктов представляют собой сложный комплекс взаимосвязанных объектов основного и вспомогательного назначения. Можно выделить следующие функциональные группы объектов: производственные, транспортно-складские, подсобные объекты инженерного обеспечения и ремонтной службы, объекты административно-хозяйственного и бытового назначения. Выбор площадок для строительства сооружений БХО сточных вод, планировка и застройка их территорий регламентированы требованиями СНиП [1-3] и рекомендациями, приведенными в работе [4]. Выбор оптимальных объемно-планировочных решений по размещению станции БХО сточных вод невозможен традиционными ручными методами. Повышение качества проектных работ с одновременным сокращением сроков проектирования возможно только на основе широкого использования в процессе поиска оптимальных проектных решений современной вычислительной техники, что в свою очередь невозможно без разработки математических моделей, методов и алгоритмов решения соответствующих задач.

Формирование генерального плана станции по чистки сточных вод с использованием средств вычислительной техники представляет собой сложную задачу размещений зданий и сооружений различных функциональных групп с учетом укрупненных технологических, инженерных и транспортных коммуникаций и природно-климатических особенностей территориального района [5, 6]. Эффективность решения такой задачи зависит от большого числа факторов, значительная часть которых трудно формализуема. Качество решения задачи, с точки зрения проектировщика, зависит от числа размещаемых объектов и их габаритов, количества минимальных и максимальных разрывов между объектами, подлежащих обязательному соблюдению, заполненности территории объектами других производств, конфигурации существующих на площадке магистралей и т.п.

Если решение задачи осуществляется для площадки с уже размещенными объектами, то актуальным становится вопрос о делении или блокировке объектов. В этом случае может возникнуть необходимость в разрыве технологической схемы и, как следствие этого, изменении его аппаратурного оформления [7, 8].

Осуществим формализацию задачи размещения объектов станции на генплане. Для этого введем ряд допущений и обозначений.

1.     

2.      Размещение объектов осуществляется в двумерном пространстве с метрикой

где (xo_i, yo_i), (xo_j, yo_j) - координаты центра объектов с номерами i, j;

 N - число размещаемых объектов.

3.      Объекты задаются в форме прямоугольников или окружностей.

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

Определение 1. Под элементарной (одиночной) трассой будем понимать отдельную коммуникацию.

Определение 2. Под групповой трассой будем понимать совокупность отдельных коммуникаций, расположенных в одном канале или на одной эстакаде.

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

Элементарные трассы либо сами связывают объекты, либо частично или полностью вливаются в групповые трассы.

Элементарную трассу представим в виде:

Для r-ой групповой коммуникации имеем:

где: N_k1, N_k2 - соответственно число групповых и элементарных трасс;

 n¢ _r, n² _j - соответственно число фрагментов r-ой групповой и j - ой элементарной трасс;

lk¢ _r, lk² _j - соответственно ширина r-ой групповой и j-ой элементарной трасс;

 x¢ _r, y¢ _r, h¢ _r; x² _j, y² _j, h² _j - координаты осевой линии трасс.

где S¢ _r - число элементарных коммуникаций, образующих групповую трассу;

 d_S, d_j - диаметр элементарных коммуникаций с номерами s, j;

 d¢ _S - соответственно величина санитарно-технического и противопожарного разрыва между соседними коммуникациями в канале (эстакаде);

 Z¢ _or, Z''_oj - соответственно ширина зоны обслуживания для r-ой групповой трассы и j-ой элементарной трасс;

 hk¢ _r, hk¢ '_j - высота канала для трасс с номерами r, j.

5.      Территория разбита на кварталы (модули) с заданными размерами. В процессе размещения объектов допускается изменение размеров кварталов.

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

7.      Подъездные пути проходят (преимущественно) по границам модулей.

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

Пусть U - заданная область, в которой необходимо разместить N соединенных определенным образом объектов, с размерами xp_i и yp_i, (xd_i - для объектов, проекция которых на плоскость задана в виде окружности), i=1,:,N. Качество размещения будем оценивать величиной критерия F, характеризующего стоимость коммуникаций и занятой ими и размещаемыми объектами площади. Значения координат размещаемых объектов xo_j, yo_i и трасс коммуникаций T_j, T° _r, r=1,:,N_k1, j=1,:,N_k2, удовлетворяющие условиям M(xo_i, yo_i, T_j, T° _r)=0 и доставляющие минимум критерию оптимизации F, будем считать решением задачи. Условия M(xo_i, yo_i, T_j, T° _r)=0, i=1,:,N, j=1,:,N_k2, r=1,:,N_k1, характеризующие правила размещения объектов и коммуникаций, являются моделью проектного решения генерального плана станции БХО сточных вод.

Представим критерий F в следующем виде:

где S₁, S₂ - соответственно стоимость территории, занятой под объекты и коммуникации; S₃ - стоимость коммуникаций; S₄ - удорожание (удешевление) стоимости объектов от деления их на части (блокировки с другими объектами).

Составляющие критерия определяются по приведенным ниже вычислениям.

где su¢ _i - стоимость единицы площади территории, занятой под i-й объект.

где M¢ _j, M² _t - число кварталов, по территории которых проходят групповые и одиночные трассы; su² _jm, su² ¢ _tl - соответственно стоимость единицы площади, занятой под групповые и одиночные трассы.

где N_k - число всех коммуникаций;

lz_j, sk_j - соответственно длина и стоимость единицы длины j-ой коммуникации;

 sl_j - затраты на эксплуатацию и техническую реализацию одного поворота для j-ой коммуникации;

 n1_j - число поворотов j-ой коммуникации.

где S0_i - стоимость i- го объекта;

 e _i - коэффициент, учитывающий блокировку (деление объекта на части);

e _i<0 - i- ый объект блокируется с другими объектами;

e =0 - i-ый объект размещается отдельно стоящим с учетом всех санитарно-технических и противопожарных ограничений;

e _i>0 - i-ый объект делится на части.

Введем следующие обозначения:

1.     

Территория (область U), на которой осуществляется размещение объектов задается в форме односвязного многоугольника с вершинами:

где xu_l, yu_l - координаты l-ой вершины, l=1,:,KU; KU - количество вершин многоугольника.

Область U делится на кварталы (модули), с учетом уже сформировавшейся части района строительства, ограниченные прямыми, параллельными осям координат:

где G - стороны квартала (модуля); N_x - количество модулей по оси 0X; N_y - количество модулей по оси 0Y; YU_max, YU_min, (XU_max, XU_min) - соответственно максимальная и минимальная ординаты (абсциссы) вершин многоугольника.

Стоимость единицы площади для каждого модуля представим в виде матрицы:

где SUa b - стоимость единицы площади для a ,b - го модуля.

Преобладающее направление векторов в области U зададим вектором a U, образующим с осью 0X угол a . За линию главного фасада станции примем прямую, параллельную оси 0X.

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

где k_g - номер внешней коммуникации;

 xk_g, yk_g - координаты точек подвода внешних коммуникаций;

N¢ _k - число внешних коммуникаций.

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

Параметры "истоков" коммуникаций представим в виде матрицы:

где k¢ _h - номер объекта для h-ой коммуникации;

 N² _k - число коммуникаций, связывающих объекты;

 xk¢ _h, yk¢ _h - координаты центра подвода h-ой коммуникации к k¢ _h - му объекту.

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

Функциональные параметры коммуникаций зададим следующим образом:

где fk_h - функциональный признак h-ой коммуникации (водопровод, линия связи и т.п.);

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

4.     

При решении задачи необходимо знать расположение источников энергии (ЛЭП, магистральные трубопроводы, природные водоемы - приемники очищенных сточных вод), расположенные за пределами области U. Зададим эти источники в виде:

где L - количество источников, kwl - количество вершин для l - го источника;

 xwu l , ywu l - координаты u - ой вершины для l - го источника.

5.     

Транспортные коммуникации, проходящие в области U, определим в форме:

где xpj _t, ypj _t - координаты точек, через которые проходит осевая линия t-ой коммуникации;

N_p - число транспортных коммуникаций;

fp_t - функциональный признак t-ой коммуникации (автодорога, коллектор и т.п.);

hp_t - ширина t-ой коммуникации.

6.      На территории строительства существуют зоны, где строительство недопустимо вообще (наличие водоемов, действующих объектов и т.п.). Запретные зоны определим в виде многоугольников:

где kz_g - количество вершин g-ой зоны запрета; Q - количество зон запрета; xz_sg, yz_sg - координаты s-ой вершины для g-ой зоны запрета.

7.     

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

где kz^'_q - количество вершин q-ой зоны для целесообразного размещения отдельных объектов;

 Q^' - количество зон;

 xz^'_pq, yz^'_pq - координаты p-ой вершины для q-ой зоны;

fp^'_q - функциональный признак q-ой зоны.

8.     

Матрицы L¢ ₁ и L² ₁ задают минимальные допустимые расстояния между объектами и границей области U:

где l1¢ _ic - минимальное допустимое расстояние между i-ым и c-ым объектами;

 l1² _i - минимальное допустимое расстояние от i-го объекта до границ области U.

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

где l1² ¢ _1i - номер объекта "истока" i-ой коммуникации;

 l1² ¢ _2i - номер объекта "стока" i-ой коммуникации;

 l1² ¢ _3i - максимальная длина i-ой коммуникации; N_md - число коммуникаций, для которых вводится данное ограничение;

 d - номер производственного комплекса;

 D - число производственных комплексов.

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

где l2¢ _ij - минимальное допустимое расстояние между i-ой и j-ой коммуникациями.

где l2² _ij - минимальное допустимое расстояние между i-м объектом и j-ой коммуникацией.

Пункты 1-10 характеризуют исходную информацию для решения рассматриваемой задачи.

Рассмотрим правила размещения объектов и коммуникаций:

Условие 1. Необходимость соблюдения санитарно-технических и противопожарных разрывов:

а) между объектами: " i,c Î [1,N]

Это условие можно представить неравенством:

б) между коммуникациями и объектами, не являющимися точками "истока" и "стока" для

соответствующих коммуникаций:

где xc_j, yc_j, hc_j - координаты точки

_jÎ T_j или c_jÎ T° _j.

Аналогично условию 1а, условие 1б запишем в форме:

 в) между коммуникациями:

или

Условие 2. Образование застойных зон с повышенной концентрацией выбрасываемых в атмосферу веществ (хлора и др.) исключает размещение объектов в секторе с углом a 0_i, образуемым линией фасада i-го объекта (осью 0X относительной системы координат для объекта) с направлением преобладающих ветров:

где Na - число объектов, для которых вводится данное ограничение.

Условие 3. Прокладка коммуникаций осуществляется на заданных уровнях:

где E_j - число возможных уровней прохождения j-ой коммуникации;

h° _jt - t-ый уровень прокладки j-ой коммуникации, знак (+) или (-) определяет расположение трассы относительно поверхности земли.

Условие 4. Соблюдение максимальных расстояний между объектами отдельного производственного комплекса обеспечивается выполнением неравенства:

или, что аналогично:

Условия (20)-(25), характеризующие правила размещения объектов и коммуникаций на генплане предприятия, и необходимая для них исходная информация (4)-(19) представляют собой математическую модель проектного решения:

Задача формирования генерального плана предприятия относится к классу комбинаторных задач. Исходя из особенностей математической модели формирования генерального плана сооружений БХО сточных вод предлагается алгоритм решения задачи, базирующейся на совмещении процесса размещения объектов и трассировки коммуникаций. Укрупненная блок-схема алгоритма представлена на рис. 1.

1.      Из множества размещаемых объектов выделяются группы, относящиеся к отдельным производственным комплексам (блокам). Если площадь, занимаемая блоком с номером d, d=1,:,D, отнесенная к площади квартала, значительно меньше принятой для подотрасли минимальной плотности застройки кварталов, то делается попытка объединения d - го комплекса с другими производственными комплексами или объектами других функциональных групп с учетом правил размещения объектов (20) - (25), плотность застройки которых, при их отдельном размещении в квартале минимальна. Объединение производится исходя из связей между блоками. Очередность размещения блоков определяется из следующих соображений: блок, занимающий большую площадь с учетом величин санитарно-технических разрывов между объектами (15) - (17), имеет более высокий приоритет при размещении. В случае, если два или несколько комплексов при плотном размещении занимают одинаковую площадь, более высокий приоритет при размещении имеет комплекс с большей удельной стоимостью коммуникаций с внешними источниками энергии или транспортными коммуникациями (точками подвода коммуникаций к границе области размещения) (7).

Для выбранного d-го производственного комплекса определяется квартал (несколько кварталов, если комплекс занимает большую площадь) на территории которого он будет размещен. Критерием оценки выбора квартала является стоимость оптимальных инженерных и транспортных трасс коммуникаций от границы квартала до внешних источников энергии и транспортных коммуникаций, то есть внешних коммуникаций для рассматриваемого комплекса, исходя из допущений, принятых при описании модели (условие 1). Оптимальными считаются трассы, проведенные внутри прямоугольника, охватывающего "истоки" и "стоки" коммуникаций. Поиск оптимального квартала осуществляется методом координатной релаксации. В качестве начального приближения выбирается квартал, до границ которого коммуникация с максимальной удельной стоимостью имеет минимальную длину. Выбор квартала для размещения следующего производственного комплекса определяется исходя из минимума стоимости коммуникаций с внешними источниками и объектами уже размещенного (размещенных) комплекса (комплексов) (3). При наличии зоны целесообразного размещения объектов определенного назначения (14) процесс компоновки начинается с данных объектов и квартала, включающего соответствующую зону.

2.      Так как полученное размещение комплексов из-за ряда допущений может быть не оптимальным, произведем оптимизацию начального размещения комплексов с помощью итерационного алгоритма, основанного на парных перестановках комплексов. Исходя из начального размещения (или результата предыдущей итерации) выбираются d-ый и j-ый комплексы (d¹ j, d, j=1,:,D) и меняются местами; рассчитывается значение критерия оптимизации (суммарной стоимости коммуникаций, для которых изменяются координаты "истоков" ("стоков")); если оно меньше предыдущего, то производится обмен. Если для d-го комплекса отсутствует улучшение критерия при всех перестановках с другими комплексами осуществляется переход к отысканию "наилучшего" места для следующего комплекса и т.д. В связи с тем, что для отдельных комплексов (блоков) существуют типовые проекты (например, для блока двух коридорных аэротенков с размерами коридора м и вторичных горизонтальных отстойников - проект 902-2-442.87 [4].), в которых приводятся рекомендуемые варианты размещения сооружений внутри блока с трассировкой отдельных внутренних коммуникаций. В этом случае можно задав вариант размещения внутри блока и обеспечив его ориентацию в области U, нужно выполнить трассировку внешних коммуникаций (см. пункт 5).

3.      В целях получения наиболее качественного размещения объектов внутри кварталов производится предварительная трассировка внешних коммуникаций с целью получения точек подвода коммуникаций к границам (вершинам) квартала. Прокладка трасс внешних коммуникаций на основе современных методов организации строительства и комплектования технологических процессов осуществляется по границам кварталов. Общая ширина зон, в пределах которых производится прокладка трасс одиночных T_j (j=1,:,N_k2) и групповых T° _r (r=1,:,N_k1) коммуникаций, определяется принятыми для каждой подотрасли промышленности минимальными предельными расстояниями от границ кварталов до размещаемых объектов. Осевые линии основных дорог проходят строго по границам кварталов. Предварительный анализ трасс коммуникаций производится исходя из предположения, что все трассы будут расположены оптимально. Это обеспечивается принятыми размерами каналов для прокладки коммуникаций. Алгоритм, по которому проводятся трассы коммуникаций, будет рассмотрен ниже. Из всех точек возможного подвода коммуникаций для каждого производственного комплекса выбираются те, для которых стоимость коммуникаций минимальна.

4.      В пределах выбранного квартала (кварталов) производится последовательное размещение объектов, относящихся к d-ому производственному комплексу (d=1,:,D). Последовательность размещения определяется исходя из размеров (площади) объектов, начиная с объекта, имеющего максимальную площадь.

Оптимизация размещения объектов осуществляется методом координатной релаксации. Идея метода заключается в следующем: для каждого i-го объекта производится циклическая релаксация по координатам. При перемещении, которое осуществляется с заданным шагом (обычно 6 м) выполняется проверка, не пересекается ли перемещаемый объект с остальными (20) и вычисляется стоимость его связей. Если это значение меньше, чем на предыдущем шаге, то процесс продолжается; если больше предыдущего, то движение по этой координате прекращается. После размещения i-го объекта переходим к i+1-му объекту и так выполняется весь цикл. Процесс поиска заканчивается, когда изменение местоположения всех объектов не приводит к улучшению критерия.

5.      При формировании генерального плана особый интерес представляет трассировка всех инженерных и транспортных коммуникаций и в первую очередь построение связывающих сетей однородных коммуникаций (деревьев Штейнера) [5].

Как показали экспериментальные исследования, приведенные при использовании метода, основанного на предварительном формировании сетей магистралей, порядок трассировки отдельных коммуникаций или их сетей практически не влияет на качество конечного результата. Однако для облегчения проведения трасс с использованием выше названного метода установим следующий порядок трассировки коммуникаций. В первую очередь трассировке подлежат отдельные коммуникации, имеющие наименьшие санитарно-технические разрывы с объектами (18) и число уровней прокладки (24), так как коммуникации, имеющие большее число уровней легче трассировать в случае пересечения с другими трассами на основном уровне. При равных условиях порядок трассировки определяется исходя из общего эвристического правила Айкерса: коммуникации трассируются в порядке приоритетных номеров, который равен числу "истоков" или "стоков" в прямоугольнике, в рамках которого проходят оптимальные трассы. Обычно короткие фрагменты коммуникаций трассируются первыми, далее трассируются окружающие их коммуникации. После прокладки одиночных коммуникаций T_j (j=1,..,N_k2) осуществляется трассировка коммуникаций, для которых возможна совместная прокладка в одном канале (эстакаде) T° _r (r=1,:,N_k1). Под сетью будем понимать некоторое число коммуникаций, имеющих одинаковый уровень проведения трассы, соизмеримые санитарно-технические разрывы с объектами (15) и сравнительно малые разрывы между собой (19). Последовательность укладки коммуникаций в канале (эстакаде) определяется из размеров охватывающего прямоугольника. Ширина каждого ортогонального фрагмента сети определяется после проведения всех коммуникаций с учетом зон обслуживания (1)-(2).

Рассмотрим приближенный метод построения сетей коммуникаций (деревьев Штейнера) в ортогональной метрике. Перед началом трассировки отдельной сети формируется сеть магистралей. Построение магистралей осуществляется следующим образом: в области U проводятся линии, отстоящие для проведения внешних коммуникаций на расстоянии, определяемом санитарно-техническими нормами для данной сети коммуникаций от линии размещения объектов в кварталах. Для внутренних коммуникаций это расстояние определяется санитарно-техническими нормами для объектов и коммуникаций (17), а так же расстоянием от границы области U (16) и уже проложенных трасс (18).

Из всех "истоков" и "стоков" сети распространяются лучи, параллельные осям координат, до ближайшего препятствия. Препятствием является отрезок ранее построенной сети того же направления, что и рассматриваемый луч, граница квартала или граница области U. Разметка магистралей осуществляется следующим образом. В начале всем магистралям, пересекающимся с исходным множеством лучей присваивается вес 1 и они помечаются номером соответствующего луча. На каждом последующем k - ом шаге все магистрали, еще не помеченные и пересекающиеся с магистралями (k-1) - го уровня, получают вес k и соответствующую отметку исходного луча. Процесс разметки магистралей на k - ом шаге осуществляется попеременно для каждой группы магистралей (k-1) - го уровня, соответствующих "истокам" ("стокам") сети. Если на некотором шаге ни одна из вновь отмеченных магистралей не имеет номера одного из исходных лучей, это является признаком того, что соединение соответствующего "истока" ("стока") невозможно. В этом случае рассматривается другой уровень прокладки рассматриваемой коммуникации (24), отличный от основного. При выполнении условия (22) разметка магистрали продолжается. Если в процессе разметки магистралей обнаружено пересечение магистралей, имеющих разные отметки, это означает, что может быть построено соединение соответствующих выводов сети. В этом случае строится фрагмент сети, после чего процесс разметки магистралей может быть продолжен. При построении фрагментов при равных возможностях осуществляется минимизация длины сети с учетом числа поворотов (3). Указанный процесс разметки и построения фрагментов продолжается до тех пор, пока не будут подсоединены все выводы сети, соединение с которыми возможно. Успешное построение каждой сети обеспечивает достаточно большая ширина каналов. Перед построением следующей сети осуществляется преобразование сети магистралей. Оно состоит в исключении частей фрагментов построенной сети и добавлении новых магистралей, параллельных этим фрагментам и отстоящих от них на расстоянии, определенном исходя из санитарно-технических норм для рассматриваемой коммуникации (18). Пример, иллюстрирующий приведенный алгоритм трассировки приведен на рис.2.

Рис. 2 Иллюстрация алгоритма трассировки.
1- граница квартала; 2 – территория для размещения объектов; 3 – магистрали возможных трасс между объектами А и В; 4 – оптимальная трасса.

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

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

1. СНиП 2.04.03-85. Канализация. Наружные сети и сооружения. / Издание официальное. - М.: Госстрой СССР, 1986.

2. СНиП 2.04.02-84. Водоснабжение. Наружные сети и сооружения. / Издание официальное. - М.: Госстрой СССР, 1985.

3. СНиП II-89-80 Генеральные планы промышленных предприятий. / Издание официальное. - М.: Госстрой СССР, 1980.

4. Каталог типовой проектной документации для строительства зданий и сооружений канализации. - М., 1991.

5. Зайцев И.Д. Теория и методы автоматизированного проектирования химических производств. - Киев: Наукова думка, 1981.

6. Минаков И.П., Рафалович И.И., Тимощук В.С. Использование ЭВМ при проектировании генеральных планов и объемно-планировочных решений. Л.: Наука, 1982.

7. Биологическая очистка производственных сточных вод: Процессы, аппараты и сооружения./ Яковлев С.В., Скирдов И.В., Швецов В.Н. и др.; Под ред. С.В. Яковлева. М.: Стройиздат. 1985.

8. Автоматизированный выбор оборудования системы очистки сточных вод. / Немтинов В.А., Мокрозуб В.Г., Субочев С.А. и др. // Труды Тамбовск. гос. техн. универ-та. Тамбов, 1999. вып. 4. С. 51.