Коксование углей

Размещено на http://www.allbest.ru/

ВВЕДЕНИЕ

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

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

Кожухотрубчатые теплообменники - наиболее распространенная конструкция теплообменной аппаратуры.

По ГОСТ 9929 стальные кожухотрубчатые теплообменные аппараты изготовляют следующих типов (рисунок 1):

- ТН - с неподвижными трубными решетками;

- ТК - с температурным компенсатором на кожухе;

- ТП - с плавающей головкой;

- ТУ - с U-образными трубами;

- ТПК - с плавающей головкой и компенсатором на ней.

В зависимости от назначения кожухотрубчатые аппараты могут быть теплообменниками, холодильниками, конденсаторами и испарителями; их изготовляют одно- и многоходовыми.



Рисунок 1 Основные типы кожухотрубчатых теплообменных аппаратов: а - с неподвижными решетками (ТН) или с компенсаторами на корпусе; б - с плавающей головкой; в - с U-образными трубками

Кожухотрубчатый аппарат с неподвижной трубной решеткой типа ТН (рисунок 2) имеет цилиндрический кожух 1, в котором расположен трубный пучок 2; трубные решетки 3 с развальцованными трубками крепятся к корпусу аппарата. С обоих концов теплообменный аппарат закрыт крышками 4. Аппарат оборудован штуцерами 5 для теплообменивающихся сред; одна среда идет по трубкам, другая проходит через межтрубное пространство.



Рисунок 2 Одноходовой теплообменник типа ТН вертикального исполнения

Теплообменники этой группы изготовляют на условное давление 0,6…4,0 МПа, диаметром 159…1200 мм, с поверхностью теплообмена до 960 м²; длина их до 10 м, масса до 20 т. Теплообменники этого типа применяют до температуры 350 °С.

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

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

Достоинством аппаратов этого типа является простота конструкции и, следовательно, меньшая стоимость.

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

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

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

Вертикальный кожухотрубчатый теплообменник типа ТК (рисунок 3) отличается от теплообменника типа ТН наличием вваренного между двумя частями кожуха 1 линзового компенсатора 2 и обтекателя 3.



Рисунок 3 Вертикальный кожухотрубчатый теплообменник типа ТК

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

При установке линзового компенсатора на горизонтальных аппаратах в нижней части каждой линзы сверлят дренажные отверстия с заглушками для слива воды после гидравлических испытаний аппарата.

Теплообменники с U-образными трубками типа ТУ имеют одну трубную решетку, в которую завальцованы оба конца U-образных трубок, что обеспечивает свободное удлинение трубок при изменении их температуры.

Недостатком таких аппаратов является трудность чистки внутренней поверхности труб, вследствие которой они используются преимущественно для чистых продуктов.

Такие аппараты (рисунок 4) состоят из кожуха 2 и трубного пучка, имеющего одну трубную решетку 3 и U-образные трубы 1. Трубная решетка вместе с распределительной камерой 4 крепится к кожуху аппарата на фланце.

Рисунок 4 Теплообменник у U-образными трубами

Теплообменники этого типа могут быть в горизонтальном и вертикальном исполнении. Их изготовляют диаметром 325…1400 мм с трубами длиной 6...9 м, на условное давление до 6,4 МПа и для рабочих температур до 450 °С. Масса теплообменников до 30 т.

Для обеспечения раздельного ввода и вывода теплоносителя в распределительной камере предусмотрена перегородка 5. Теплообменники типа ТУ являются двухходовыми по трубному пространству и одно- или двухходовыми по межтрубному пространству. В последнем случае в аппарате установлена продольная перегородка, извлекаемая из кожуха вместе с трубным пучком.

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

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

Поскольку механическая очистка внутренней поверхности труб в аппаратах типа ТУ практически невозможна, в трубное пространство таких аппаратов следует направлять среду, не образующую отложений, которые требуют механической очистки.

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

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



Рисунок 5 Горизонтальный двухходовой теплообменник с плавающей головкой

В теплообменных аппаратах подобного типа трубные пучки сравнительно легко могут быть удалены из корпуса, что облегчает их ремонт, чистку или замену. Горизонтальный двухходовой конденсатор типа ТП состоит из кожуха 10 и трубного пучка. Левая трубная решетка 1 соединена фланцевым соединением с кожухом и распределительной камерой 2, снабженной перегородкой 4. Камера закрыта плоской крышкой 3. Правая, подвижная, трубная решетка установлена внутри кожуха свободно и образует вместе с присоединенной к ней крышкой 8 «плавающую головку». Со стороны плавающей головки аппарат закрыт крышкой 7. При нагревании и удлинении трубок плавающая головка перемещается внутри кожуха. Для обеспечения свободного перемещения трубного пучка внутри кожуха в аппаратах диаметром 800 мм и более трубный пучок снабжают опорной платформой 6. Верхний штуцер 9 предназначен для ввода пара и поэтому имеет большое проходное сечение; нижний штуцер 5 предназначен для вывода конденсата и имеет меньшие размеры.

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



РАСЧЕТНАЯ ЧАСТЬ

Исходные данные:

Исходная смесь - сероуглерод-тетрахлорид углерода

G_F = 5,2 т/ч

x_F = 42 % мол.

tисх. см. = 11°С

Примем температуру исходной смеси в резервуаре 11°С. Температура конденсации водяного пара при давлении 0,7 МПа составляет t_конд^п = 165°С [6, C. 549, таблица LV11].

Условимся все величины, связанные с исходной смесью отмечать индексом «ИС», а связанные с паром - «П», со стенкой - «ст».

Напишем температурную схему:

Насыщенный водяной пар 165°С < 165°С

Исходная смесь 11°С > 46,3°С (температура кипения сероуглерода)

Большая и меньшая разность температур между исходной жидкостью и греющим паром:

Дt_б = 165-11 = 154°С,

Дt_м = 165-46,3 = 118,7°С.

Средняя разность температур Дt_ср, входящая в уравнение теплопередачи, определяется как [6, C. 169, формула (4.78)]:

(1)

Вычислим по формуле (1) Дt_ср:

°С

Средняя температура исходной смеси [6, C. 82, формула (4.82)]:

°С

Расход исходной смеси:

Плотность исходной смеси при 28,7°С рассчитаем по формуле (2) [6, C. 512, таблица 1V]:

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

(3)

где V - объёмный расход жидкости, м³/с;

- скорость движения жидкости, м/с;

- площадь круга, м².

Объёмный расход исходной смеси по формуле (4):

Рассчитаем теплоёмкость исходной смеси по формуле (, [6, C. 562, рисунок X1])

Расход теплоты на нагрев исходной смеси [5, C. 45, формула (2.2)]:

(5)

где - удельное количество тепла, Дж/с;

- количество подаваемой смеси, кг/с;

- теплоёмкость тела, Дж/(кг•К);

- температурный напор, К.

Расход теплоты для рассчитываемого процесса по формуле (5):

Расход греющего пара с учётом 7 % потерь теплоты , - энтальпии пара и конденсата при давлении 0,7 МПа и температуре 165°С [6, C. 549, таблица LX11]:

Уравнение теплопередачи [6, C. 168, формула (4.72)]:

(6)

где Q - количество передаваемого тепла, Вт;

- коэффициент теплопередачи, Вт/(м²*К);

F - площадь теплопередающей поверхности, м²;

- температурный напор, К.

Ориентировочно определяем максимальную величину площади поверхности теплообмена. Минимальное значение коэффициента теплопередачи для случая теплообмена от конденсирующегося водяного пара к органическим жидкостям [6, C. 164, таблица 4.8]. При этом по формуле (6):

По формуле (67) найдём вязкость исходной смеси при температуре 28,7°С, [6, С. 556, рисунок V]:

Критерий Рейнольдса находится по формуле [4, C. 64, формула (1-64)]:

(7)

где Re - критерий Рейнольдса;

- скорость движения жидкости, м/с;

- эквивалентный диаметр трубопровода, м;

- вязкость, Па•с.

Для обеспечения турбулентного движения исходной смеси при Re>10000 скорость в трубах должна быть больше (выбираем трубы для теплообменника 25Ч2, как самые часто используемые в промышленности) по формуле (7), выражая, скорость:

Число труб 25Ч2 мм, обеспечивающих объёмный расход гептана при Re = 10000 [4, C. 228]:

Условию n<23,0 и F<2,4 м² удовлетворяет [4, C. 215, таблица 4.12] одноходовый теплообменник наружным диаметром 150 мм с числом труб n = 13, длина труб 2 м.

Уточняем значения критерия Re:

Критерий Прандтля [6, C. 151, формула (4.12)]:

(8)

где - критерий Прандтля, безразмерная величина;

- теплоёмкость, Дж/(кгК);

- вязкость, Пас;

- теплопроводность, Вт/(мК).

Значения теплопроводности сероуглерода и тетрахлорида углерода близки, поэтому для расчёта теплопроводности исходной смеси воспользуемся формулой ( при 28,7°С:

Рассчитаем критерий Прандтля для исходной смеси при температуре 28,7°С по формуле (8):

Теплоотдача при развитом турбулентном течении в прямых трубах и каналах (Re > 10000) [6, C. 152, формула (4.17)]:

(9)

где - критерий Нуссельта, безразмерная величина;

- критерий Рейнольдса, безразмерная величина;

- критерий Прандтля, безразмерная величина;

- критерий Прандтля, рассчитанный при температуре теплопередающей стенки, безразмерная величина;

- поправочный коэффициент, учитывающий влияние на коэффициент теплоотдачи отношения длины трубы L к её диаметру d.

Найдём коэффициент . Для этого найдём отношение

[6, C. 153, таблица 4.3]. Отношение примем равным 1,0 (с последующей проверкой). Рассчитаем по формуле (9):

Критерий Нуссельта [4, C. 305, формула (2-38)]:

(10)

где - критерий Нуссельта, безразмерная величина;

- коэффициент теплоотдачи, Вт/(м²•К);

- эквивалентный диаметр трубопровода, м;

- теплопроводность, Вт/(м•К).

Из формулы (10) выразим и найдём :

Необходимо найти коэффициент теплоотдачи при конденсации водяного пара на пучке горизонтальных труб. Расчёт осуществляется приближённо по формуле (11) [6, C. 229]:

(11)

где - коэффициент теплоотдачи со стороны греющего пара, Вт/(м²•К);

- коэффициент, зависящий от числа труб в теплообменнике;

- коэффициент, зависящий от температуры греющего пара;

- число труб в теплообменнике;

- длина труб в теплообменнике, м;

- расход греющего пара, кг/с.

, при отсутствии воздуха в паре, [6, С. 162,рисунок 4.7] при шахматном расположение в пучке и при числе труб по вертикали n_p = 5 для выбранного нами теплообменника [6, С. 215, таблица 4.12] [6, C. 162, таблица 4.6]. Таким образом, по формуле (11) находим:

Принимаем тепловую проводимость загрязнений со стороны греющего пара 1/r_загрИС = 1/5800 1/(Вт/(м²•К)), со стороны исходной смеси 1/r_загрП = 1/5800 1/(Вт/(м²•К)) [6, C. 531, таблица XXX1]. Коэффициент теплопроводности стали [6, С. 529, таблица XXV11]. Можем рассчитать [6, C. 229]:

Коэффициент теплопередачи [6, C. 168, формула (4.74)]:

(12)

Рассчитаем коэффициент теплопередачи по формуле (12):

Поверхностная плотность теплового потока:

Проверим принятое значение [6, C. 229]. Определяем:

Находим температуру стенки со стороны исходной смеси:

При температуре 140°С находим теплоёмкость (, [6, C. 562, рисунок X1]), вязкость, теплопроводность исходной смеси:

Ввиду близости расположения точек на номограмме [6, C. 556, рисунок V] принимаем значение вязкости . Ввиду параллельного расположения зависимостей теплопроводности гептана и октана и необходимости экстраполяции до температуры 140°С, принимаем. Таким образом, по формуле (8):

Проверяем принятое соотношение :

Было принято . Посчитаем разницу в процентах:

Разница в 1,7 процентов между принятым значением и рассчитанным укладывается в инженерную пятипроцентную погрешность. Расчёт коэффициента теплопередачи К закончен. Коэффициент .

Расчётная площадь поверхности теплообмена по формуле (6):

Коэффициент теплоотдачи , следовательно, определяющий диаметр труб - внутренний, то есть .

Теплообменник с длиной труб 2 м имеет площадь поверхности теплообмена:

Запас площади поверхности теплообмена:

Запас площади теплообмена достаточен.

Принимаем одноходовой кожухотрубчатый типа ТН (ГОСТ 15118-79, ГОСТ 15120-79, ГОСТ 151122-79) теплообменник с внутренним диаметром кожуха 150 мм, числом труб 13 и длиной труб 2 м. К рабочему теплообменнику устанавливается резервный теплообменник того же типа.

Расчёт гидравлических сопротивлений трубопровода и выбор насоса

Исходная схема трубопровода представлена на рисунке 6.



Рисунок 6 Расчётная схема: 1 - резервуар с исходной жидкостью; 2, 5, 6, 8, 9 - вентили; 3 - центробежный насос; 4 - мембрана; 7 - теплообменник; 10 - тарельчатая ректификационная колонна

кожухотрубчатый теплообменник трубопровод гидравлический

Расчёт гидравлических сопротивлений

Принимаем, что исходная смесь хранится при температуре 11°С. Плотность исходной смеси при температуре 11°С рассчитаем по формуле (2), [6, C. 512, таблица 1V]:

Исходная смесь движется в теплообменнике. Средняя температура смеси в теплообменнике 28,7°С (см. п. 3.2), при этой температуре .

Плотность исходной смеси после теплообменника (температура 46,3°С) рассчитываем по формуле (2), [6, C. 512, таблица 1V]:

Таким же образом определим вязкость смеси на различных интервалах движения исходной жидкости.

Вязкость исходной смеси при t=11°С, [6, С. 556, рисунок V]:

Вязкость исходной смеси при t = 28,7°С, .

Вязкость исходной смеси при t = 46,3°С, [6, С. 556, рисунок V]:

Рассчитаем объёмные расходы исходной смеси на различных участках (см. рисунок 3).

Объёмный расход до теплообменника:

Объёмный расход в теплообменнике: .

После теплообменника:

Определим примерный диаметр трубопровода. Пусть скорость до теплообменника и после него (в теплообменнике ) будет щ = 2 м/с, тогда рассчитаем примерный диаметр трубопровода по формуле (3):

до теплообменника

после теплообменника

Из стандартного ряда труб принимаем трубу 76Ч3,5.

Внутренний диаметр такой трубы [5, C. 17]:

Уточним скорость движения в принятой трубе по формуле (3), выразив скорость:

до теплообменника

после теплообменника

Определим режим движения жидкости по формуле (7):

до теплообменника

в теплообменнике (см. п. 3.2):

после теплообменника

Найдём коэффициенты шероховатости труб [5, C. 14]. Будем использовать новые трубы, так как ректификационная колонна - это дорогостоящее сооружение, которое должно работать непрерывно, окупая затраты. Поэтому недопустимо прекращение работы ректификационной колонны из-за проблем в трубопроводе.

Для стальных новых труб принимаем коэффициент шероховатости . Рассчитаем вспомогательный коэффициент е [5, C. 14]:

Проверим критериальные условия:

Таким образом, расчётная формула для расчёта коэффициента трения л [5, C. 14, формула (1.6)]:

(13)

где - коэффициента трения, безразмерная величина;

- вспомогательный коэффициент, безразмерная величина;

- критерий Рейнольдса, безразмерная величина.

Рассчитаем коэффициент л по формуле (13):

до теплообменника

в теплообменнике

после теплообменника

Найдём местные сопротивления для каждого участка (рисунок 10):

до теплообменника [3, C. 14-15, а также 4. С. 520, таблица X111]

в теплообменнике

после теплообменника

Формула для расчёта гидравлического сопротивления [6, C. 25, формулы (1.49) и (1.50)]:

(14)

где - полное гидравлическое сопротивление, Па;

- скорость движения жидкости, м/с;

- коэффициента трения, безразмерная величина;

- коэффициент местного сопротивления, безразмерная величина;

- длина участка трубопровода с определённым , м;

- диаметр трубопровода, м;

- плотность перекачиваемой жидкости, кг/м3;

, - давления в приёмном и исходном резервуарах, Па;

- высота поднятия жидкости, м.

С учётом того, что в исходном и приёмном резервуарах давление одинаковое, то расчётная формула выглядит следующим образом:

Из рисунка 3 находим сумму расстояний на различных участках трубопровода и задаёмся их значениями:

до теплообменника

в теплообменнике (длина труб теплообменника - 2 м, а также по 0,5 м на входное и выходное пространство)

после теплообменника

Таким образом, по формуле (14) получаем:

до теплообменника

в теплообменнике

после теплообменника

Общее гидравлическое сопротивление трубопровода:

Из соотношения ДРсети = сgh определим

h_сети = ДРсети/ (с•g) = 173328/(1437,5•9,81) = 12,5 м.

Будем считать, что характеристика сети представляет собой правильную параболу, выходящую из точки с координатами V = 0 м³/ч; h на которой известна точка с координатами (для бесперебойной работы подставляем наибольший объёмный расход - после теплообменника) и h_сети = 12,5 м. Найдем коэффициент параболы.

Общее уравнение параболы у = а*х² + b. Подставив значения, имеем 12,5 = а*3,71² + 12,4. Тогда, а = 0,0073. Уравнение у = 0,0073*х² + 12,4.

Возьмем несколько значений объемной производительности и определим напор h_сети (таблица 1).

Таблица 1

Зависимость напора сети от производительности насоса

Производительность, м3 / ч	Напор сети, м	Напор насоса, м
0	12,4	16
5	12,6	15,6
10	13,1	14,6
15	14,0	13,6
20	15,3	12,6
25	16,9	11,6
30	18,9	10,6

По данным таблицы 1 построим зависимости сети и насоса (рисунок 7).

Наиболее близкий по параметрам к расчетному насосу в каталоге является центробежный насос К50-32-125а производительностью 11,2 м³/ч и напором 16 м. Характеристика такого насоса с диаметром ротора 165 мм также приведена на рисунке 7. Расположение рабочей точки (2) относительно расчетной точки (1), показывает, что насос может преодолеть гидравлические сопротивления сети и подавать в неё заданную смесь на заданную высоту.



Рисунок 7 Совмещение характеристик сети и насоса

К установке предлагается центробежный насос К50-32-125а с диаметром ротора 165 мм, который обеспечит перекачивание 3,71 м³/ч смеси сероуглерода с тетрахлорид углеродом по трубопроводу из хранилища в теплообменник и до питающей тарелки колонны. При этом его к.п.д. составит 58%, мощность на валу насоса 1,1 кВт.



ЗАКЛЮЧЕНИЕ

Для нагрева 5,2 т/ч бинарной смеси сероуглерода с тетрахлорид углеродом с 11єС до 46,3єС требуется кожухотрубный теплообменник типа ТН (ГОСТ 15118-79, ГОСТ 15120-79, ГОСТ 151122-79) с диаметром кожуха D = 150 мм, числом труб n = 13, длинной труб l = 2 м, площадью теплопередачи F = 1,7 мс запасом поверхности 59,7%. Кроме рабочего устанавливается также резервный теплообменник того же типа.

Для подачи смеси из хранилища необходимо установить центробежный насос типа К50-32-125а с производительностью V₂ = 11,2 м/ч, напором Н = 16 м, мощностью N = 1,1 кВт, с КПД = 58 %.



СПИСОК ЛИТЕРАТУРЫ

1. Ахметов С. А. и др. Технология и оборудование процессов переработки нефти и газа. СПб., Недра, 2006, 868 с.

2. Филиппов В. В. Теплообмен в химической технологии. Теория. Основы проектирования: учеб. пособие. Самара, 2014, 197 с.

3. Справочник нефтехимика. В двух томах. Т. 1/ С. К. Огородникова. Л., Химия, 1978. 496 с.

4. Касаткин А. Г. Основные процессы и аппараты химической технологии. М; 1961.

5. Основные процессы и аппараты химической технологии: Пособие по проектированию /Под ред. Ю. И. Дытнерского, 2-е изд., М.: Химия, 1991. 496 с.

6. Павлов К. Ф., Романков П. Г., Носков А. А. Примеры и задачи по курсу процессов и аппаратов химической технологии. Учебное пособие для вузов/10-е изд., перераб. и доп. Л.: Химия, 1987. 576 с., ил.

Размещено на Allbest.ru

...