Измерение и расчет фазовых соотношении в условиях равновесия для многокомпонентных смесей

ГЛАВА VI ИЗМЕРЕНИЕ И РАСЧЕТ ФАЗОВЫХ СООТНОШЕНИИ В УСЛОВИЯХ РАВНОВЕСИЯ ДЛЯ МНОГОКОМПОНЕНТНЫХ СМЕСЕЙ

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

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

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

§ 1. АНАЛИЗ ГАЗА

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

1.    Низкотемпературную ректификацию.

2.    Аналитическую масс-спектрометрию.

3. Инфракрасную абсорбционную спектрометрию (ИКС).

4.    Газовую хроматографию.

Эти методы определения состава углеводородного газа или жидкости будут кратко описаны в данной главе.

Общепринятые методы анализа газа путем абсорбции компонента в растворе описаны в учебниках по физической химии [V. 58] и в стандартах ASTM (например, метод D1136 [V. 6]). Плотность газа можно установить взвешиванием определенного объема [V 6,    1070-52].

Метод измерения плотности газа является вполне удовлетворительным анализом двухкомпонентных смесей, но при этом нельзя не учитывать отклонения от законов идеального состояния [IV’. 83]. Для анализа потоков стали применять элементы, измеряющие удельную теплопроводность, а также аппаратуру для измерения других характеристик [VI. 9].

НИЗКОТЕМПЕРАТУРНАЯ РЕКТИФИКАЦИЯ

До конца двадцатых годов не существовало систематических методов лабораторного определения состава природных газов. После сожжения пробы, при котором измерялись уменьшение объема вследствие сгорания компонентов газа и содержание двуокиси углерода в продуктах сгорания, компоненты природного газа выражались в метановом и этаповом эквиваленте. В. Дж. Подбельняк (W. J. Podbielniak) сконструировал и применил низкотемпературную фракционирующую колонку, в которой в качестве охлаждающего агента, создающего орошение, использовался жидкий воздух или азот (VI. 41).

На рис. VI. 1 показана одна из первых моделей аппарата для низкотемпературной ректификации. Для обеспечения орошения жидкого метана в. колонках последние помещались в посеребренные вакуумные рубашки.

Температуру кипения измеряли термопарой, установленной на верху колонны. Объем отогнанного газа определяли по повышению давления в вакуумирован-ном приемнике, измеряемому при помощи манометра. По очереди отгоняли метан, этан, пропан, изобутан, н-бутан и т. п. Когда точка кипения достигнет комнатной температуры (в обычных условиях изопентан кипит при 28° С), становится необходимым создать в колонке вакуум для предотвращения конденсации отходящих паров в соединительных линиях. Несколько изомеров гексана имеют относительно близкие температуры кипения и поэтому их объединяли в одну гекса-новую фракцию. Помимо гексанов, может испариться и остаток, состоящий из гептана и высших углеводородов, если предел выкипания этого остатка относительно низок. Чаще, однако, гептан и высшие углеводороды оставались на дне колонки в виде жидкости, плотность которой определяют путем взвешивания известного объема в пикнометре, а молекулярный вес — по понижению температуры замерзания бензола. Если имеется жидкий остаток, объем всех компонентов газа и объем жидкости выражается числом молей, откуда легко уже определить мольные доли каждого компонента. Из подаваемого в колонку газа обычно удаляют (при помощи твердой каустической соды) двуокись углерода и сероводород, так как эти соединения образуют твердые отложения во фракционной колонке. Содержание этих газов определяют на отдельном приборе путем абсорбции раствором едкого кали, для чего нужна вто-

Рис. VI. 1. Одна из первых моделей для низкотемпературной фракциояировки, применявшаяся для анализа природных газов (Л) и смесей летучих жидкостей (В и Л).

/ — к вакууму; 2—1 л\ 3 — 2 л; 4 — манометры; 5 — запорный клапан; 6 — крошка пемзы; 7 — стальной контейнер с образцом»' 8 — к милливольтметру; 9— уравнительный сосуд.

рая проба газа. Содержание сероводорода можно установить и при помощи раствора йода [VI. 2].

С момента появления аппаратуры для низкотемпературной ректификации в нее были внесены значительные усовершенствования [VI. 44]:    улучшена насадка,

введены автоматическая регистрация объема и температуры отогнанного газа, автоматическое регулирование охлаждения, увеличена скорость перегонки (рис. VI. 2 и VI. 3). Серьезным шагом вперед явилось измерение теплопроводности отходящих паров в качестве средства распознавания границ кипения фракций [VI. 43]. Американская Ассоциация газового бензина (NSAA) производит испытания, целью которых является определение точности анализа методом низкотемпературной ректификации, ибо этот метод часто является основой контрактов на поставку [VI. 35]. Объем образца газа, отбираемого для низкотемпературного анализа, обычно исчисляемого литрами, а продолжительность анализа часами.

Состав сложных жидких смесей, содержащих компоненты от метана до углеводородов с температурой кипения 121° С и выще, можно определить двухступенчатой низкотемпературной ректификацией. Г азы вплоть до пентанов анализируются подобно природному газу, а остаток образца ректифицируют во фракционирующей колонне с конденсацией отходящего продукта. Для жидкой части можно построить кривую разгонки, состав жидкой части можно также характеризовать последовательным рядом фракций, каждая из которых выкипает в узком температурном интервале. Фракции эти различаются по молекулярному весу и плотности. Используемые температурные пределы выкипания фракций иногда соответствуют температурам кипения нормальных предельных углеводородов, и тогда эти фракции именуют октановой, нонановой, декановой и т. п. В табл. VI. 1 показан состав пластовой жидкости, определенный комбинированным методом анализа [VI. 27].

Таблица VI. 1 Состав углеводородной смеси, содержащей газообразные и высококипящие компоненты [VI. 27] Состав пластовой нефти и газа, находящихся в равновесии при абсолютном давлении, равном 269 ат, и температуре 99° С

* Интервал выкипания фракции соответствует интервалу температур кипения нормальных предельных углеводородов с числом атомов, равным номеру фракции.

Рис. VI. 2. Современная аппаратура для низкотемпературной ректификации.

о — внешний вид* 6 — схема потока. I *— гидравлическо-электронный прибор для регулирования а контроля расхода дистиллята; 2 — парообразный дистиллят; «? — конгуры для измерения теплопроводности; 4 — датчик; 5 — регистрирующие манометр, кабельная связь и механизм, перемещающий картограмму; 6 — манометр колонны; 7—приборы, регулирующие давление в колонне и охлаждение флегмы; «5 — автоматическая подача образцов газа; 5 —двигатель; 10 теплопроводность; // — температура; 12 — вакуум; 13 — удаление паров; 14 — воздух; 15 — выпускное отверстие; 16 —. термостатическая воздушная баня.

Рис. VI. 3. Запись на картограмме, полученная в результате низкотемпературной ректификации.

Образец —. иефтезаводский крекинг-газ. /-—разница в разрывах характерной кривой разгонки объясняется удержанием жидкости

во флегме: 2 — кривая теплопроводности дистиллята.

МАСС-СПЕКТРОМЕТР

Масс-спектрометр ионизирует молекулы газа в условиях низкого давления и классифицирует образующиеся молекулярные и осколочные ионы в соответствии с их массами [VI. 12, VI. 1, VI. 6, VI. 63]. Схема действия масс-спектрометра показана на рис. VI. 4. Подлежащий анализу газ в количестве около 0,2 см³ при давлении

'.:2

11 у

16

1 ат помещают в баллон напуска. Масс-спектрометр может иметь также устройства, обеспечивающие испарение жидкостей с целью последующего анализа их в парообразном состоянии при низком давлении. Из баллона напуска погок газа через диафрагму (10—15 мк) поступает в ионизационную камеру, в которой происходит ионизация и диссоциация молекул и образование осколочных и молекулярных ионов. Ионизированные частицы под действием высокого напряжения приобре. тают ускорение и проходят через щель в магнитное поле. Эти частицы отклоняются магнитным полем в соответствии с их массой и фокусируются на выходную щель и коллектор. Ток в цепи коллектора усиливается и регистрируется. Изменяя ускоряющее напряжение ионов, можно получить спектр заряженных частиц с массой от 1 до 250 и более.

Фотография аналитического масс-спектрометра фирмы «Консолидейтед инжиниринг корпорейшн» показана на рис. VI. 5.

Регистрирующий прибор фиксирует высоту пики (рис. VI. 6) в зависимости от различных соотношений между массой и зарядом. Соотношение между высотой пики, выраженной в единицах шкалы, и давлением образца, выраженным в микронах, в первоначально поступившем объеме является постоянным. Э^т0 соотношение называется чувствительностью Совокупность относительных высот пик по сравнению с высотой какого-либо пика, взятого в качестве эталона, является масс-спектром, являющимся функцией структуры молекулы. Различные молекулы одинакового молекулярного веса различаются распределением интенсивностей ионов в масс-спектре. Аналитический масс-спектрометр

Характерный образец: газ, жидкость Образец

i

Периодически отВираемыи образец расширяется до газообразного состояния при ~~ 0,05мм pm cm Просашбание газа

(Молекулы)+(е) ^\(аон)~+{ион)**{мткуло$* L    ? (атом)+(е) ~ J

сн_?+е-—сн}+сн?+сн$ + СН⁺+С⁺+Н⁺

+10006 -10006 _и

I    Ион*\ Напряжение

I    ускорения

<5    ионоб

Г-5

10

Г-7

13

IV-

17

Рис. VI. 4. Схема потока в масс-спектрометре, /—сосуд емкостью 3 л; 2— диафрагма; 3 — ионизирующая камера; 4 — поле электродвижущей силы;

5 — входная щель: 6 — ускоритель частиц; 7 — выходная щель; 8 — коллектор и усилитель; 9 — регистрирующий прибор; 10— ускоренные частицы; //— тяжелые ионы, 12 — легкие ноиы; J3 — коллектор;

14 — сигнал в усилитель; 15 — ионы в фокусе; 16 — сопротивление.

обеспечивает постоянство масс-спектров и чувствительности для каждого вида молекул независимо от времени и парциального давления компонентов в смеси. Американский нефтяной институт (API) составил справочные таблицы масс-спектров и чувствительностей для большого числа углеводородов [VI. 1]. Для выполнения аналитической работы масс-спектрометр калибруют по эта-

Рис. VI. 5. Приборный сгол масс-спектрометра фирмы «Консолидейтед» (Consolidated), модель 21-103В.

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

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

определить состав сложных смесей, состоящих из двадцати и более компонентов. Более того, в благоприятных обстоятельствах прибор обнаруживает следы веществ, концентрация которых составляет несколько частей на 1 млн. В табл. VI. 2 приведены результаты анализов искусственной смеси, которые проводились периодически в течение 6 месяцев, и анализов образцов, взятых из газового пробоотборника и отправленных в лабораторию наряду с другими обычными образцами без указания, что они собой представляют.

Таблица VI. 2 Параллельные анализы газовых смесей масс-с пектрометром

А, Анализы искусственной смеси газа, проводившиеся периодически в течение 6 месяцев

Мольная доля

28

24

25

26    27

Масса /заряд

Продолжение табл. VI. 2 Б. Анализы образцов, последовательно взятых из газового пробоотборника

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

ИНФРАКРАСНАЯ СПЕКТРОСКОПИЯ

В инфракрасном спектрометре получают спектр поглощения инфракрасных лучей образцом газа или жидкости [VI. 22]. Различные молекулы поглощают инфракрасные лучи различных частот, что позволяет определить содержание каждой из них в образце. Оптическая схема инфракрасного спектрометра приведена на рис. VI. 7.

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

Инфракрасным спектрометром можно анализировать газы, жидкости, прозрачные твердые вещества или ве-

Рис. VI. 7. Схема оптической системы инфракрасного спектрометра [VI. 22].

/ — деталь гребневидной заслонки; 2 — приемное зеркало;

3 — болометр; 4 — зеркало, установленное под углом 45°;

5— выходная щель; 6 — призма с углом 60°; 7 — зеркало длины волн; 8 — прицельное зеркало; 9 — входная щель; 10 — вогнутое зеркало; // — сосуд с образцом; 12 — источник инфракрасных лучей; 13—гребневидная заслонка; 14 — вращающаяся отражательная заслонка; 15 — вращающийся столик; 16 — синхронный двигатель; 17 — сосуд с образцом, принятым за эталон.

щесТва в растворителях. При анализе жидкостей берут образцы объемом 1—2 см³ или меньше, а при анализе газа — образцы объемом 200 см³ или меньше.

Поглощение лучей различной длины волны регистрируют и сопоставляют со стандартными величинами поглощения, указанными на шкале спектрометра или взятыми из опубликованных работ [VI. 1].

ГАЗОВАЯ ХРОМАТОГРАФИЯ

Газовая хроматография — сравнительно новый метод анализа. Своим названием она обязана изменению окраски адсорбентов при поглощении определенных соединений. Во время анализа смеси разделяются на отдельные компоненты, при этом компоненты поглощаются на поверхности твердого адсорбента и каждый из них десорбируется в потоке газа--носителя [VI. 42,

VI. 18, VI. 43, VI. 8, VI. 17]. По теплопроводности выходящего потока газа можно судить о содержании каждого из компонентов.

При хроматографии пар адсорбируется или абсорбируется на поверхности пористого твердого вещества. Газ-носитель переносит индивидуальные соединения вперед по слою, где они распределяются по степени их адсорбируемости. К тому в^меии, когда соединения достигнут конца колонки, они оказываются уже раз-

Рис. VI. 8. Схема газового хроматографа.

/—•змеевик с твердым адсорбентом или же с твердым материалом, покрытым нелетучим жидким абсорбентом;    2 —¦ термостатируемая

баня:    3— впуск образца:

4 —элемент для измерения теплопроводности; 5 — регистрирующий вольтметр; 6 —¦ расходомеры:    7    —    поток    ге

лия.

деленными и следуют один за другим в устройство, в котором измеряется их теплопроводность. Основное оборудование, необходимое для разделения и анализа, показано на рис. VI. 8. Проявляющей колонной может служить трубка диаметром 6,4 мм и длиной 1,2—1,8 м. Она заполняется гранулированным твердым веществом, например кирпичной крошкой, поверхность которой покрыта нелетучим абсорбентом, например триизобутиле-ном, диоктилфталатом или силиконом. На устройстве, измеряющем теплопроводность, показывается концентрация соединений в выходящем потоке газа-носителя. Переменными факторами при хроматографии являются диаметр насадки колонки, размеры колонки, температура, состав газа-носителя, скорость продвижения газа-носителя и объем образца. Для каждого данного вида анализа подбирают оптимальное сочетание параметров и тарируют аппарат по известным образцам.

Процедура анализа заключается в следующем.

1.    Потоком газ а-носителя, например гелия, стабилизируется колонна при постоянной температуре.

2.    Добавляют во входящий поток гелия порцию газа.

3.    Регистрируют изменение теплопроводности во времени.

Хроматограмма показана на рис. VI. 9. Продолжительность анализа зависит от длины колонки и ее содержимого, скорости движения газа-носителя и объема

Рис. VI. 9. Хроматограмма деэтанизиро-ванного газового бензина.

Колонна высотой 1,8 м; поглотитель — трикрезил-фосфат на целине; температура 50° С; газ-носи-тель — гелий. / — «-гептан; 2—н-тексан; 3— изо-гексан; 4— «-пентач; 5 — изопентан; 6—пропаи;

7 — изобутаи; 8 — «-бутан; 9 — начало

образца. Соединения идентифицируют на основании времени, которое затрачивается на прохождение через слой.

Эта техника анализа непрерывно усовершенствуется, особенно в части методики выбора надлежащей адсорбционной колонны. Типовая аппаратура, применяемая в настоящее время, показана на рис. VI. 10.

ОПРЕДЕЛЕНИЕ СОДЕРЖАНИЯ ГЕЛИЯ

Часто при изыскании источников получения гелия необходимо определить содержание его в природном газе. Хорошие результаты при этом дает использование масс-спектрометра. Горное бюро [VI. 3] разработало метод анализа на гелий, принятый в качестве стандартного. Этот метод состоит в адсорбции всех прочих газов на активированном угле при температурах жидкого азота и в измерении давления, создаваемого гелием.

ОПРЕДЕЛЕНИЕ СОДЕРЖАНИЯ СЕРОВОДОРОДА

Для определения содержания сероводорода в природном газе широко используется метод Тутвейлера (Tutweiler) [VI. 2]. Сероводород растворяют в воде и воздействуют на него стандартным раствором йода, используя в качестве индикатора крахмал. Обычно содержание H2S выражается в гранах серы на 100 куб. фут. газа (г/100 л³); 590 гран серы на 100 м³ (или 1360 г/100 м³ соответствуют 1% от объема (при 15,6° С и абсолютном давлении, ранном 1,03 ат).

§ 2. ИЗМЕРЕНИЕ РАВНОВЕСИЯ ПАР —ЖИДКОСТЬ

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

Аппаратура, примененная Сиднеем Юнгом (Sidney Jonng) [VI. 68] для определения упругости пара и плотности насыщенных фаз тридцати чистых веществ, показана на рис. VI. 11. Вещество «запирается» ртутью в стеклянных трубках, окруженных термостатирующими паровыми рубашками; достижение температуры начала кипения и точки росы устанавливается визуально. Кей (Kay) [III. 30] применял аналогичную аппаратуру при своих исследованиях двойных систем. Он вставлял небольшой стальной шарик в капиллярную стальную трубку и перемешивал жидкость при помощи соленоида, находящегося снаружи обогревательной рубашки.

Для приведения в равновесие многокомпонентных смесей при давлениях 350 ат и выше использовалась аппаратура различных типов. Толстостенные стальные сосуды, вмещающие 500—1000 см³, могут быть снабжены внутренней мешалкой [III. 26]. как показано на рис. VI. 12, и прикреплены к цапфе для качания [111. 68]. Для достижения в таких сосудах равновесия можно также пропускать пар через жидкость [VI. 11]. Оборудование этого типа позволяет получать равновесные образцы пара и жидкости для анализа.

Сейдж и Лес и (Sage and Lacey) и их сотрудники разработали аппаратуру для комбинированного определения объемных свойств и коэффициентов распределения пара и жидкости [VI. 55].

Общая схема этой аппаратуры показана на рис. VI. 13. Равновесный сосуд 1 погружают в масляную ванну с мешалкой, температуру которой регулируют с помощью обычных устройств, принятых для этой цели. Ртуть может быть введена в сосуд 1 из резервуара 2 через вентиль 3,

Рис. VI. 11. Аппарат Сиднея Юнга для определения свойств газа и жидкости.

На верхнем конце вертикального стержня 4, входящего в сосуд 1 снизу, находится указатель уровня ртути 5 с электрическим контактом. Электроконтакт и идущий от него по стержню электропровод изолированы от стержня и сосуда. Стержень 4 выведен из заполненной ртутью части сосуда через сальник, имеющийся в нижней стенке камеры 6. Этот сальник позволяет стержню двигаться вертикально лишь с очень ничтожной утечкой ртути из равновесного сосуда. Под сальником стержень присоединен к ходовому винту 7, который приводится в движение электродвигателем через червяк, находящийся в сцеплении с муфтой ходового винта. Счетчик числа оборотов 8, находящийся на валу червячной передачи, определяет после надлежащей калибровки положение точки контакта уровнемера 5 с точностью до 0,05 мм.

Во избежание изменений общего объема системы сосуда при подъеме или спуске стержня 4 к нижнему концу ходового винта присоединяют стержень такого же диаметра, что и стержень 4. Стержень 10 проходит через сальник в заполненную ртутью камеру 11, соединяющуюся с равновесным сосудом. Это компенсационное устройство не только устраняет какое бы то ни было влияние движения стержня на общий объем заполненной жидкостью системы, но служит также для предотвращения больших перемещений уровнемера 5, которые могли иметь место при изменении давления в сосуде. Благодаря компенсатору 10 при повышении давления наблюдаются лишь такие деформации, как упругое сжатие в винте 7 и упругое удлинение в окружающем его корпусе.

В стенку равновесного сосуда вделаны вентили 12 И 13. Через них сосуд присоединен к вакуум-насосу и К

Рис. VI. 12. Установка Катца и Гахмуса (Katz and Hachmuth) для изучения сщесей в равновесных условиях [111.26].

1— резервуар с ртутью; 2— насос для подачи ртутн;

3 — ртуть; 4 — сжатый азот; 5 — ртутный резервуар высокого давления; 6 манометр; 7 — к приборному щиту; 8 — термометр; 9 — манометр с надлежащим интервалом измерения; 10 — газовая фаза;

И—жидкая фаза; 12 — трубка для взятия образца газа; 13 — трубка для взятия образца жидкости;

14—мешалка; 15 — заземление; 16 — двигатель мешалки термостата; 17 — масляный или водяной термостат;

18 — двигатель мешалки сосуда; 19 — сосуд; 20 — нагреватель на 1000 вт\ 21 — к анализатору газа; 22 — к анализатору жидкости; 23 — вентиль понижения давления;

24 — вентиль для загрузки и очистки.

аппарату для ввода или извлечения образцов (через вентиль 13, а также вентиль 15 или 16). Вентиль 15 позволяет откачивать содержимое нижнего манифольда независимо от равновесного сосуда.

Установление равновесия в сосуде ускоряется при помощи спиральной мешалки 17, которая сконструирована так, что свободное поперечное сечение внутри сосуда остается неизменным в любой точке, за исключением точки, расположенной близ самого днища и самой крышки. Эта спиральная мешалка приводится в действие трубчатым валом, к которому внутри корпуса прикреплен якорь из мягкого железа. Корпус изготовлен из нержавеющей стали. Снаружи корпуса находится вращающийся электромагнит, вращающий спиральную мешалку со скоростью около 100 об/мин.

Детали уровнемера показаны на рис. VI. 14, а равновесного сосуда — на рис. VI. 15. Для определения давления применяют грузовой манометр. Предельное давление равно 700 кГ/см², а интервал температур, который обеспечивает термостат, — от —1,1° С до +238,8" С с точностью +0,06° С.

Объемные и фазовые измерения производятся на взвешенном количестве материала. Газ поступает из

резервуара 19 через вентиль 20, а жидкость — из сосуда 21 через вентиль 16 (см. рис. VI. 13). Общий объем углеводородов в равновесном сосуде определяется посредством указателя ртути 5, а поверхность раздела углеводородного газа и жидкости — при помощи теплового индикатора. Объем системы, заполненный углеводородами, регулируют подачей и откачкой ртути через вентиль 3.

Результаты, полученные на описанном выше оборудовании для смесей, состав которых дан в табл. VI. 3, приведены на рис. VI. 16 и VI. 17.

Методы определения равновесий в системе газ — жидкость при низких температурах рассмотрены в работе Блумера и Перента (Bloomer and Parent) [III. 7].

В двухфазном потоке, движущемся по трубам к сепаратору, работающему при постоянных температуре и давлении, фазы обычно находятся в состоянии, близком к равновесному. Результаты анализа образцов двух фаз, взятых из промысловых сепараторов, приведены в табл. VI. 4.

‘Средний молекулярный вес 137,1.

Состав образцов, взятых из промысловых сепараторов

Таблица VI. 4

§ 3. КОЭФФИЦИЕНТЫ РАСПРЕДЕЛЕНИЯ (КОНСТАНТЫ РАВНОВЕСИЯ) СЛОЖНЫХ СМЕСЕЙ

Первые диаграммы для определения констант равновесия были основаны на экстраполяции фугетивно-стей и их соотношений [III. 67, III. 41, VI. 30]. Такие диаграммы назывались «диаграммами идеальных констант равновесия», поскольку они строились для углеводородных смесей, подчиняющихся законам идеальных растворов. На рис. VI. 18—VI. 20 приведены идеальные константы равновесия, вычисленные Брауиом и Холкомбом (Brown and Holcomb) [VI. 14] для пропана, изобутана и я-бутана.

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

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

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

Для смесей оклахомской нефти и природного газа эта работа была проведена Катцем и Хахмусом (Katz and Hachmuth) [111.26]. Диаграммы рис. VI. 21—V. 29, а применимы для многих нефтей среднеконтинентальной платформы. Диаграммы для изобутана, я-бутана, изо-пентана и я-пентана были составлены в 1953 г. Хиндсом и Берроусом (Hinds and Burrows). Константа равновесия для гептана, умноженная на коэффициент 0,15, дает величину, достаточно точно выражающую константу равновесия для суммы гептанов и более тяжелых углеводородов оклахомской нефти. При использовании диаграмм для систем иефть — природный газ для вычисления констант равновесия метана, а также гептана и более тяжелых углеводородов может потребоваться

Рис. VI. 13. Аппарат конструкции Сейджа и Леси для определения объемных свойств и фазового поведения [VI. 55].

Рис. VI. 14. Деталь уровнемера [VI. 55].

Рис, VI. 16. Объемное поведение смеси природного газа и яефти [VI. 55].

/ — точка начала кипения; 2 — две фазы при 37.8° С; 3 — жидкость при 37.8° С.

Рис. VI. 17. Объемное поведение углеводородной смеси, обнаруживающее ретроградную точку .росы [VI. 55].

/ — газовая и жидкая фазы; 2 — жидкая фаза; 3 — газовая фаза; 4 — ретроградная точка росы,

К-у/.х WO г

80 60 50 40 30

20

10

8 6

5

41 3

2

1

0,8

0,6

0,5

0,4

0,3

0,2

0,1

W 80    120    160    200    240    280

Температура, °C

-40    0

Рис. VI. 18. Идеальные константы равновесия для пропана [IV. 14].

W 80    120    160    200 2kO 28<

0,06

0,05

0,04

0,03

К⁼УР>

Температура,°С

Рис. VI. 19. Идеальные константы равновесия для изобутана [IV. 14].

Температура^С

Рис. VI. 20. Идеальные константы равновесия для «-бутана [IV. 14].

100 п во -

30    40    50    60    70    80    90

Температура, °С

K-?

1000т 800

600

-30    -20    -10

по т

to

<Б

-30 -20    -10    0    10    20    30    kO    50    60    70    80    90    100    110    120    130 Температура, °С

Рис. VI. 22 Коистаиты равновесия для этана в оклахомской нефти [111.26].

К-?

^к~х

Температурке

„ к=У/х

N0

N0

СЛ

н-лснтана

изопентана

-20    -10    0    10    20    30    Ц-0    50    60    70    80    90

Температура, °С

Рис. VI. 28. Константы равновесия для гексанов в оклахомской нефти [III. 26].

к=у/х

30    -20    -10 О W го 30 ч-о 50    60    70    80    90

Температура, °С

Таблица VI. 5

Экспериментальные данные по равновесиям фаз сложных углеводородных смесей

некоторая модификация этих диаграмм [VI. 53]. В работе Стендинга и Катца [III. 68] показано, что состав нефти имеет большое значение при высоких давлениях. Константы, полученные для нефти типа оклахомской, можно с достаточной достоверностью применить при расчетах, связанных с абсорбцией маслами. Данные Роланда, Смита и Кавелера, выполненные в виде графиков, могут быть использованы для расчета сепарации газоконденсатной смеси [VI. 49].

Как видно из табл. VI. 5, были исследованы и условия равновесия неуглеводородных компонентов — двуокиси углерода, азота и сероводорода. Диаграмма, составленная Поттменом (Poettmann) для двуокиси углерода в нефти, приведена на рис. VI. 30 [VI. 47]. Поттмен производил измерения при содержании двуокиси углерода в газовой фазе в 6—11 % мол. В сырой нефти из месторождения Биллингс растворялось больше двуокиси углерода, чем в конденсате из Эратского месторождения [VI. 45]. Якоби и Рзаса (Jacoby and Rzasa) [VI. 29] определили составы фаз системы нефть—природный газ (молекулярный вес С₇ + равен 227, удельный вес 0,894), в которой содержалось азота, двуокиси углерода и сероводорода по 5% мол. Величины К для этих компонентов сопоставлены с величинами К для метана и этана на рис. VI. 31—VI. 33. Диаграммы для азота и сероводорода в этой системе приведены на рис. VI. 34 и VI. 35. Вагтборг (Vagtborg) [VI. 62] изучал естественную систему нефть — природный газ, содержавшую значительное количество сероводорода и небольшое — двуокиси углерода и азота. Он осуществлял дифференциальное испарение нефти, насыщенной при абсолютном давлении 175 ат и 67,8° С. Осредненные данные по дифференциальным составам фаз приведены в табл. VI. 7, а константы равновесия (все при 67,8° С)—в табл. VI. 8.

Таблица VI. 6

Опубликованные работы, в которых выведены корреляции констант равновесия жидкой и парообразной

фаз

К=у/х

-10    О    10    20    30    40    50    60    70    80    90    100    110

Т емпература,°С

Рис, VI. 30. Константы равновесия для двуокиси углерода в нефти Биллингса [VI. 47].

10 100 Абсолютное давление,ат

10

Рис. VI. 31. Константы равновесия для компонентов нефти при 37,8° С [VI. 29].

/ — сероводород; 2 — этан; 3— двуокись углерода; 4 — метан; 5 — азот*

СЭ

СХ

СЗ

5

сз

5

со

з:

о

си

:э

Рис. VI. 32. Константы равновесия для компонентов нефти при 65,50° С [VI. 29].

/ — сероводород; 2 — этан*, 3 — двуокись: углерода; 4 — метан; 5 =— азот.

сз

са.

сэ

е

а:

сз

е

Со

а;

о

5:

100

давление,ат

Лбсолютное

Рис. VI. 33. Константы равновесия для компонентов нефти при 93,3° С.

/ — сероводород; 2 — этан; 3 — двуокись углерода; 4 — метан; 5 азот.

ДВсалютнив даВление,от

го за т so 60 ю so зо wo по по по Температура/С-

Рис. VI. 34. Константы равновесия азота в нефти (смесь С-6) [VI. 29].

20 30 40 50 60 70 80 90 100 НО 120 130 Температура, °С

Рис. VI. 35. Константы равновесия сероводорода в нефти (смесь С-6) [VI. 29],

0,01

1,0    2,0 3,04,0 6,08,010    20 30 40 6080100    200300    600

moo

800

600

CO n

о * 3

m

300

2 oo

юо

80

60

w

30

20

I 2

В ,

las

0,6

0,4

0,3

0,2

0,1

0,08

0,06

0,04

0,03

0,02

Абсолютное давление,am

Рис. VI. 36, Константы равновесия при 48,9° С.

/ — двойные системы Ci—С_х — пС₄; С¹!— nC_t; С_г—пС4, С₈ — лС,; 2 природный газ — нефть (Катц и Гахмус); «?-— природный газ— конденсат (Роуленд, Смит и Кейвлер);

4 — природный газ — нефть (Стендинг и Катц); 5 — метан; б—этан; 7 — пропан; 8 — бу-таны; 9 — пентаны; 10 — гексаны; — гептаны.

Таблица VI. 7

Состав фаз при выделении растворенного газа из Пластовой жидкости при 67,8° С [VI. 62]

Абсолютное давление, ат

* Состав жидкой фазы при атмосферном давлении определен экспериментально. Состав пара выражен «сглажен

ными» данными.

** Состав жидкости при температуре начала кипения дни.

(абсолютное давление 175). Определен путем экстраполя-

Таблица VI. 8 Константы равновесия для пластовой жидкости при 67,8° С [VI. 62]

КОРРЕЛЯЦИЯ КОНСТАНТ РАВНОВЕСИЯ

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

Хенсон и Браун (Hanson and Brown), [VI. 26] сопоставили различные системы при температуре 48,9° С (рис. VI.36). Они показали, что при одинаковом давлении схождения данный компонент даже при различном составе смеси будет иметь одинаковые константы равновесия. Первоначально применявшийся метод определения констант равновесия для систем с летучестью, промежуточной между исследованными системами, заключался в интерполяции диаграмм, подобных показанной на рис. VI. 36. Уайт и Браун (Whitte and Brown) [VI. 65] предложили простую методику отыскания кон-ртант равноиесия при известных кажущихся величинах давления схождения. Кривая идеальной константы равновесия по этому методу строится в координатах log К—log р. По таким диаграммам получаем величины минимума констант равновесия и давления, необходимого для минимума К, как функцию давления схождения и упругости паров данного компонента. Схема такого построения кривой константы равновесия показана на рис. VI. 37.

щи

Рис. VI. 37. Определение величин К по давлению схождения.

/ — идеальная К; 2 — минимум К; 3 — по закону Рауля; 4 — давление минимума К',

5 — давление схождения.

Давления схождения или геометрические места критических точек двойных систем показаны на рис. III. 14. В общем можно констатировать, что в многокомпонентных системах повышение давления схождения вызывается либо добавлением к системе высоколетучих веществ, например азота и водорода, либо уменьшением содержания промежуточных компонентов. Высокое давление схождения приводит к замедленному повышению величин К для менее летучих компонентов и к ускоренному снижению величин К для более летучих компонентов при давлениях, при которых они не следуют больше законам идеальных растворов. Хеддн (Hadden) [VI. 24] составил корреляцию констант равновесия на основе давления схождения и метода приближения к давлению схождения [VI. 25]. Уинн (Winn) [VI. 67] разработал номограмму (рис. VI. 38) для нахождения констант равновесия при давлении схождения 350 ат.

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

"1 ^f л

emtpou**

®- I

6ем%ьпв

Абсолютное давление,а

Абсолютное давление схождения Рси.ат

tr.

C\j    С\|    ГЪ    frj    t»J

(    С\1    f>.'    Q1*    ж«ч*    г-^*

10    0.2 Ц25 0.3 0,35 0*50,5 0,6 0,70,д 1,0    1,5    2,0    2,5 3.0 3,5 40 50 6 7

Абсолютное даВление 6 системе р, ат

PhCj VI. 39. График для определения сеточного давления при давлении схождения,

отличном от 350 кГ/см² абс. [VI. 67].

Рис. VI. 40. Показатель степени летучести в уравнении (VI. 1) [VI. 67].

1 •— чистые углеводороды; 2 — нефтяные фракции (погон уд. веса 0,78)* на рис. VI. 39. Методика определения К в этом случае заключается в следующем.

120 160 200 2Ь0 280 320 360 Температура нипеная при атмосферном давлении,

1.    На рис. VI. 38 находим точку давления—температуры системы. Проводим прямую линию через эту точку и К = 1,0.

2.    Находим точку пересечения этой линии с кривой давления, соответствующего давлению на сетке рис. VI. 39.

3.    Проведя через эту точку пересечения и точку компонента прямую линию, определяем К.

Для компонентов с температурой кипения выше 98,9° С (гептан С₇) константы определяем по следующему уравнению:

^(V,'¹⁾

где К_п—константа равновесия нефтяной фракции или смеси с высокой температурой кипения; Кг — константа равновесия этана при температуре, давлении и давлении схождения в системе; Кч — константа равновесия гептана при температуре, давлении и давлении схождения в системе; Ь — показатель степени летучести (рис. VI. 40).

Давления схождения Уинн (Winn) определяет по критическим кривым двойных систем, приведенным на рис. VI. 41. Для тройных и сложных систем Уинн предлагает следующее.

-40    0    40    80    120    160    200    240    280    320    360    400    4-4-0    480

Температура,°С

Рте. VI. 41. Давление схождения двойных углеводородных систем [VI. 67].

W00 в ООО 5000 4000

0    4-0    80    120    1Б0    200    74Q    280

Температура, °С

По рис. VI. 41 при температуре системы нужно найти давления на двух кривых, являющихся геометрическим местом критических точек двойных смесей’ самого легкого компонента с промежуточным и тяжелым ^компонентами. Если температура системы выше критической температуры промежуточного компонента, берем критическое давление этого компонента. Умножая эти две величины давления на соответствующие весовые доли промежуточного и тяжелого компонентов в жидкой фазе, складываем результаты и делим их на сумму весовых долей промежуточного и тяжелого компонентов в жидкой фазе.

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

(MlW'i + (Рс, vh^W2 + (Pc. vh У₃ + • • •

^Рс- ⁰ ~    W₂    +    W₃ + . . .

где р _Cf0 — абсолютное давление схождения системы в ат- {р _{c v}) 1—давление на кривой критических точек Самого легкого компонента и компонента, обозначенного, данным подстрочным значком при температуре системы; W1 — весовая доля/компонента, обозначенного дан ным подстрочным значком, в жидкой фазе. Индексы 1,

2, 3 ... обозначают 1 — самый легкий компонент, 2 — следующий, более тяжелый компонент и т. п.

Рзаса, Гласс и Опфелль (Rzasa, Glass and Opfell) [VI. 51] разработали корреляцию величин К для сложных нефтяных систем, по которой давление схождения определяется из температуры и произведения молекулярного веса фракции гептана и более тяжелых углеводородов на ее удельный вес (рис. VI. 42).

Рис. VI. 42. Определение давления схождения по температуре и свойствам фракции гептанов и более тяжелых углеводородов [VI; 51 ].

Ч В 8 10    20    М 6080100 200 400600

Абсолютное давление,ат

Рис. VI. 43. Константы равновесия для различных величин давления схождения при 126,7° С [VI. 51],

/ — метан; 2 — этан; 3— пропан; 4 — изобутан; 5 — «-бутан; б —I изопентан; 7— н-пентан; 8 — н-гексан; 9 — н-гептан; 10 — н-октан; // — н-нонан; 12 — «-декан.

Константы равновесия при различных давлениях схождения и температуре 126,7° С приведены на рис. VI. 43.

Для определения равновесий в системе природный газ—жидкость широкс применяют диаграммы фугетив-ности, составленные фирмой «Келлогг» [VI. 13]. Константу равновесия определяют по уравнению

(VI. 2)

где Цх — фугетивность компонента, деленная на его мольную долю в жидкой фазе; f/y — фугетивность компонента, деленная на его мольную долю в парообразной фазе.

Соотношение между фугетивностью и мольной долей зависит от средней мольной температуры кипения фазы и выражается рядом диаграмм, каждая из которых дает это соотношение для данного компонента при данном давлении. Применение диаграмм (примером может являться диаграмма для метана, приведенная на рис. VI. 44) осуществляется по методу последовательных приближений. Для определения К исходят из четырех переменных: температуры, давления, мольных средних температур кипения жидкости и пара.

Другие методы корреляции и работы, в которых они приведены, перечислены в табл. VI. 6.

КОРРЕЛЯЦИЯ ОРГАНИКА—БРАУНА (ORGANICK—BROWN), ДИАГРАММЫ NGAA

(Американской ассоциации газового бензина)

Органик и Браун [VI. 37] разработали корреляцию констант равновесия, в которой кажущееся давление схождения вычисляется по средней мольной температуре кипения паровой фазы и средневзвешенному эквиваленту молекулярного веса жидкости. На диаграммах К выражены как функции температуры, давления и давления схождения.

Давление схождения определяют по двум диаграммам (рис. VI. 45 и VI. 46). Для данных температуры и давления равновесия находят среднемолекулярную температуру кипения пара и средневзвешенный молекулярный вес жидкости М_в. Задаются давлением схождения или корреляционным давлением р_к и определяют соотношение р/р_к. По рис. VI. 45 находят показатель степени в члене (ЮО/Мв)".

По рис. VI. 46 определяют ординату р_к (ЮО/Мв)", сопоставляя ее с величиной, вычисленной по заданному значению рк и М_ъ. Этот метод последовательных приближений в сочетании с повторными вычислениями фазовых равновесий очень трудоемок. Однако, пользуясь этим методом, можно получить точные величины К при любом давлении схождения.

Можно построить график константы равновесия как функции нормальной температуры кипения компонента (NBP) при постоянной температуре и давлении. На рис. VI. 47—VI. 49 приведены такие графики, построенные для температуры 93° С и нескольких давлений, при трех давлениях схождения с использованием констант Органика и Брауна.

Если расчет ведется не для нормальных парафиновых углеводородов, а для других соединений, Органик и Браун рекомендуют применять эквивалентный молекулярный вес (рис. VI. 50). Данные о влиянии растворителя для ароматических и нафтеновых абсорбционных масел приведены в работе Киркбрайда и Бертетти (Kirkbride and Bertetti) [VI. 31].

Вычислять давления схождения следует лишь в тех случаях, когда другие расчеты неприменимы. Американская ассоциация газового бензина (NQAA) выпустила сборник диаграмм констант равновесия и кривых давления схождения наиболее часто встречающихся смесей (рис. VI. 51—VI. 53). Диаграммы К при абсолютном давлении схождения, равном 140 ат, для некоторых компонентов природного газа приведены на рис. VI. 54—VI. 66. Имеются также диаграммы абсолютных давлений схождения, равных 42, 56, 70, 210, 280, 350, 700 и 1400 [V. 1].

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

Роланда (Roland) [VI. 50] показывают, что по мере приближения к точке росы константы могут и не сойтись и не иметь тенденцию к схождению (рис. VI. 67). Это явление можно отнести за счет присутствия в нефти асфальтенов. Точка росы в их присутствии становится нечеткой, ибо при этом остается пластическая твердая асфальтовая фаза. Отсюда следует, что для сложных систем природного газа и нефти вместо термина «давление схождения» правильнее применять «кажущееся давление схождения», или «корреляционное давление».

jjx, am

350

“-70-60    -40    -20    0    20    40    60    80    100    120    W    160    180    200

Температура, °С

Рис. VI. 44. График Келлога: фугетивность для метана при абсолютном давлении, равном 70 [VI, 13].

КОМПОНЕНТЫ С ВЫСОКОЙ ТЕМПЕРАТУРОЙ КИПЕНИЯ

В пластах высокого давления высококипящие компоненты при температуре пласта являются летучими. Гоффман, Крамп и Хокотт (Hoffman, Crump and Hocott) [VI. 27] приводят некоторые характеристики этих высоко-кипящих компонентов в пластовых условиях. Летучесть абсорбционного масла при некоторых условиях характеризуется в работе Веббера (Webber) [VI. 64].

Поттмен и Мейленд (Poettmann and Mayland) вывели корреляцию констант равновесия высококипящих компонентов [VI. 46] на основе диаграммы фугетивности, построенной Киркбрайдом.

На рис. VI. 68 приведена диаграмма для компонентов с нормальной температурой кипения, равной 260°С. Компоненты рассматриваются раздельно, поскольку характер компонента влияет на давление схождения. На рис. VI. 69—VI. 71 представлены данные Гоффмана, Крампа и Хокотта [VI. 27], совмещенные с диаграммами Поттмена и Мейленда.

Показатель п

0,2

0,1----------1

О 0,1    0,2    0,3    0,4-    0,5    0,6 OP 0,8    0,9    1,0р/р_к

Рис. VI. 45. График 1 для определения давления схождения [VI. 37].

Равновесное давление ^к Корреляционное давление На кривых приведена среднемольная температура кипения паров (МАВР).

0    Q2    0,4-    0,5    0,8    10

Рис. VI. 46. График 2 для определения давления схождения.

ю

СО

-100 -80 -4-0    0    b-0    80    120 tBO 200 24-0

нормальная температура кипения₁⁰ С

Рис. VI. 47. Константы равновесия как функция нормальной температуры кипения компонента при абсолютном давлении схождения, равном 280 кГ/сж², и температуре 93° С.

0,0001

-120 -80 -4-0    0    40    80    100    160    200    24-0

Пормольноя температура кипения, °С

опою

0008

0,0006

0,000b

0,0003

00002

Рис. VI. 48. Константы равновесия как функция нормальной температуры кипения компонента при абсолютном давлении схождения, равном 420 кГ/см², и температуре 93°С.

ooow

0,0008

0,0006

0,0004-

0,0003

0,0002

0.0001

-120 -во -40    0    4О 80    120

iso гоо 2ьо °С

нормальная температура кипения,

Рис. VI. 49. Константы равновесия как функция нормальной температуры кипения компонента при абсолютном давлении схождения, равном 700 кГ/см², и температуре 93° С.

200    400    600    800    1000

Эквивалентный молекилярнь/й Вес

Рис. VI. 50. Эквивалентные молекулярные веса углеводородов [VI. 37].

80 120 160 Температура,°С

Рис. IV. 51. Давления схождения для типовых смесей в стабилизационной колонне [VI, 40].

Абсолютное давление схождения,am

ЮОО 900 800 7 00 600

500

Ш

3 00

Абсолютное давление схождения,am

160    200    2W

Температура, °С

Рис. VI. 52. Давления схождения для типовых смесей в абсорбционных колоннах [VI. 40].

Абсолютное давление схождения, am

Рис. VI. 53. Давления схождения для типовых смесей в десорбционной колонне [VI. 40].

Абсолютное довление, ат

30    50    70    100    300    500    1000    3000    5000    10000

Абсолютное давление, am

fOOcrrт elTb-

6 5 U

3 2

W

8

6

5

U

3

2

1,0

8

6

5

2

о;

8

б

5

U

3

0,1

8

6

5

4

3

2

0,001    __________________________________________U

0,7 /    2    3    4 5 6 8 /0    20    30 4-0 60 80 (00    200 300 500 700

Абсолютное давление, ат

0.7 I 2    3    С    5    6    7    8    10    20    30 UO 60 80 WO

Абсолютное давление, am

200 300    500    700

2 J 4 5 6 7 8 10    20    30 UO 60 80100    200    300    500    700

Абсолютное давление, am

Абсолютное давление, am

0,0001

0,001

0,0008

0,0006

0,0005

0,0004

0,0003

0,0002

1,0    2,0    3,0    4,05,0 7,0 10,0    20    30 4050 70 1П0    200    300    400    600

Абсолютное давление, ат

1.0    2.0    3,0    if,0 5,0 7,0 10,0    20 30 4050 70 /00

Абсолютное давление,am

200 300 т 600

0,0008

0,0006

0,0005

0,0004

0,7 t 2    3    U    5 6 78910    20    30 UO 50 70 100    200 300 500

Абсолютное давление, am

700

1    2    3    4    5    6    7    8    910    20    30    4050 60 80100    200    300    4    00

Абсолютное давление, ат

Абсолютное давление, am

Рис. VI. 67. Константы равновесия для природного газа и нефти из месторождения

Биллингс.

/ — постоянный состав; 2—различные смеси; «У —данные Катца и Хахмуса; 4 — равновесия, полученные при абсолютной упругости паров, равной 656,2 кГ!смЪ

°С

537 °С 482° О

Ц-26 °С

0,0001

0,6 1,0

3,0    6,0    10,0    Kt

0,03    0,06    0,1

0,0003 0,0006 0,001    0,003    0,006    0,01

Рис. VI. 68, График констант равновесия компонентов с нормальной температурой кипения 260° С. Характеристический фактор: / — 10,5; 2—11,5; 3 —12,5.

Рис. VI. 69. Константы равновесия компонента с нормальной температурой кипения

N9° С.

Рис. VI. 70. Константы равновесия компонента с нормальной температурой кипения

260° С,

0,00010,0002 0,00050,001 0,002 0,005 0,01 0,02    0,05    0,1    0,2    0,5    1,0    2    5    10    К-.

Рис. VI. 71. Константы равновесия компонента с нормальной температурой кипения

204° С.

§ 4. ВЫЧИСЛЕНИЕ РАВНОВЕСИЙ СИСТЕМ ПАР — ЖИДКОСТЬ

Для смесей, состав которых известен, температуру и давление точек начала кипения и точек росы, а также состав фаз смеси в данных условиях можно вычислить по константам равновесия. Количество и состав каждой фазы вычисляют по методу последовательных приближений [VI. 30, 111.27].

Точке начала кипения (или, как ее иначе называют, точка появления первого пузырька) соответствуют температура и давление, при которых сумма членов Кх равна единице:

ЪКх = Ху = 1,    (VI.    3)

где К — константа равновесия; х — мольная доля компонента в жидкой фазе; у — мольная доля в паровой фазе.

Таблица VI. 9 Вычисление давления точки начала кипения смеси при 49° С (измеренное абсолютное значение равно 85,75 ат)

Примечание. Путем экстраполирования получаем точку начала кипения, при температуре 49° С, равную 83,5 ат.

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

Точка росы — это температура, при которой образуется первая капля жидкости, и ей соответствуют температура и давление, при которых сумма членов у/К равна единице:

1у/К=1х=1.    (VI.    4)

Пример расчета точки росы приведен в табл. VI. 10.

Таблица VI. 10

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

Примечание. Путем интерполяции находим, что температура точки-росы при 9,1 ат равна 62,6° С.

Равновесное или однократное испарение можно рассчитать для тех смесей, которые при определенных температуре и давлении разделяются на две фазы (рис. VI. 72). При этом весь пар находится в равновесии со всей жидкостью. Установление равновесия может иметь характер непрерывного или периодического про-

Подставляя значение х из (VI. 5) в (VI. 7), получим zF=-^L + yV.    (VI.    S)

При F = 1

Если числитель и знаменатель умножим на К и разделим их на V, то

Кг

zF

У —

L/K+V- V K + L/V

Кг

У—'

(*-1)1/+1 Для всех компонентов

F Кг V' K + L/V

Аналогично, подставляя значение (VI. 7), получим

F г    г

х

1у= 1 = Х-

I ' К +L/V- (/С — 1) V + 1 •

Кг

V-X _{K + LjV} -^_{K + LlV} . (VI. 14)

Зная состав исходной смеси и константы равновесия для избранных температуры и давления, задаемся величиной V. Чтобы определить, правильно ли взята эта величина₁ вычисляем

Mz/(K+L/V) или lz/(K+L/V).

Диаграмма принятых и вычисленных величин V приведена на рис. VI. 73. Следует отметить, что кривые приближаются к начальной точке, даже если правильное решение находится где-то вдоль линии, проходящей под углом 45°. Пример расчета приведен в табл. VI. 11. Иной путь решения также по методу последователь-(VI. 7) ных приближений получаем при другом сочетании урав-

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

Обозначим F — число молей исходной смеси; V — число молей пара; L — число молей жидкости, находящейся в равновесии; г — мольная доля компонента в исходной смеси; х— мольная доля компонента в жидкости; у — мольная доля компонента в паре; i — любой компонент.

По равновесию

— =Kl.    CVI- 5)

^Х1

По общему материальному балансу

F=L+V.    (VI. 6)

По материальному балансу данного компонента

Рис. VI. 72. Равновесное испаре-    Рис.    VI.    73.    Решение    уравнения

ние при постоян-    равновесного    испарения,

ных температура и дазлении.

1 —исходная смесь [г]; 2 — пар; 3 — жидкость.

zF — xL + yV.

(VI. 9)

(VI. 10)

(VI. 11) у в уравнение

(VI. 12)

Х*= 1 = 2-у- д._|_ _L/v .    (VI.    13)

Пусть F— 1 моль, тогда по уравнениям (VI. 11) и

(VI. 13)

Таблица VI. 11

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

равном 112 am и 48,9° С

Примечание. Путем экстраполяции получим С, « 196; удельный вес С, + = 0,850.

V г= 0,887 мол, пара, Z, = 113 мол. жидкости. Молекулярный вес

Пример вычисления равновесного испарения методом повторных попыток (метод Ньютона)

V₂ = 0,8 +

0,22616

1,99409 V₂ = 0,8 + 0,1134 V₂ = 0,913

l/„ = 0,913 .

0,12251

5,01925 V_a = 0,913 — 0,0244 V_a = 0,8886

Продолжение табл. VI. 12

V = 0,8886 —

0.01466

3,84803

V    = 0,8886 — 0,00381

V    = 0,88479

V    = 0.8848

0

(VI. 16)

1,0

dV

нений (VI. 5), (VI. 6) и (VI. 7). Если вычтем уравнение (VI. 12) из (VI. 11), то получим

z(K— ])

. (VI.lt,

Если V определена верно, правый член уравнения [VI. 15] равен нулю.

Выбор величины V основан на определении угла наклона кривой при первоначально заданной величине V (рис. VI. 74). Угол наклона кривой при V получаем, дифференцируя по V выражение г(К- 1)

г (К— I)²

Величина z(K—l)/[V(K—1) + 1], полученная при первой попытке, и угол наклона кривой позволяют вычислить Л V, которую приплюсовывают к величине, заданной при первой попытке:

[V (K — 1) + l]² '

•V(K- 1)+ 1

= X

Рис. VI. 74. Решение уравнения равновесного испарения (VI. 15) по методу Ньютона.

Профиль кривой может отклоняться от прямой линии между первой и второй точками (рис. VI. 74) и для получения точного ответа может понадобиться третья попытка. Величину V₃ получаем из V2 точно так же, как V₂ была получена из V|. Пример подобного решения приведен в табл. VI. 12,

.    z(*-l)

' V(K- 1)+1

.    z(/C-l)²    ¦

' IV/ (iC— 1) + l]²

(VI. 17)

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

§ 5. ДИФФЕРЕНЦИАЛЬНОЕ ИСПАРЕНИЕ

Интегрируя это выражение, получаем !п (x_f/x₀)

(VI. 21)

К—I

K-l

(VI. 22)

(VI. 23)

(VI. 24) (VI. 25)

1

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

ydW = Wx — (W — dW) (x— dx),    (VI. 18)

где W — число молей жидкости; * — мольная доля данного компонента в жидкости; у — мольная доля данного компонента в паре.

Раскрывая скобки, опуская dWdx и подставляя у = = Кх. получаем

xdW + Wdx = KxdW.    (VI. 19)

После преобразования и интегрирования при допущении, что К — const, получим

I

CdW № ~К

(VI. 20)

Индекс 0 — первоначальное состояние и индекс f — окончательное состояние.

1-я попытка

X = 0,9885

2-я попытка

Q    0,1113    0,8120    27,07    6,19    5,19    0,3395    0,0378    8,13    0,307

С₃    0,888?    0,8120    27,07    0,617    —0,3831    1,0828    0,962    0,720    0,693

268    2    = 0,9998    ?=1,000

Wj_

W о

h

wv

In

Xo

W,\K-

-¹ =1 X()    1^0/

В конце вычисления у нас будет

Х.*^ = 1,

lK_fxi = 1.

Пример

Какое количество пара нужно дифференциально удалить при постоянной температуре (37,8° С) для того, чтобы получить жидкость с абсолютным давлением точки начала кипения, равным 21,1 ат, из жидкости, состоящей из 11,13% мол. метана и 88,87% мол. пропана. Точке начала кипения жидкости соответствует абсолютное давление, равное 33 ат при 37,8° С.

Решение

Применяем величины К из рис. III. 31.

Поэтому нужнс дифференциально испарить 18,8% жидкости.

Аналогичные уравнения можно вывести, заменив К отношением K_n/Ki- Таким образом, для приведенного

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

Пример

Определить конечное давление, необходимое для дифференциального испарения 25% мол. смеси, состоящей при температуре 37,8°С из 32% мол. метана и 68% мол. пропана. Точке начала кипения соответствует абсолютное давление, равное 70,3 ат.

Решение

Применяем величины К из рис. III. 31

Первая попытка Допускаем, что абсолютное давление р_те— 63 ат.

X = 1,0035

Вторая попытка Принимаем 'абсолютное давление равным р_те =

= 62 ат.

X = 1,000

Определяем конечное давление путем вычисления точки начала кипения конечной жидкости.

Первая попытка Предположим, что абсолютное давление p_f- = 56 ат.

S = 0,9675

Вторая попытка Предположим, что абсолютное давление р₍ =54 ат.

Pf — конечное абсолютное давление, равное 54 ат.

§ 6. ПРИМЕНЕНИЕ ЭЛЕКТРОННЫХ ВЫЧИСЛИТЕЛЬНЫХ МАШИН ДЛЯ РАСЧЕТОВ ИСПАРЕНИЯ

Расчеты испарения можно вести при помощи логарифмической линейки, однако этот метод очень трудоемок, а результаты получаются неточные. Настольные счетные машины повышают точность результатов, но все же расчеты остаются весьма трудоемкими. Поэтому для расчетов испарения, связанных с многочисленными повторениями большого числа сложных вычислений, все шире применяются цифровые электронные вычислительные машины [VI. 39, VI. 48, VI. 36, VI. 16].

Цифровая электронная вычислительная машина [VI. 7, VI. 28, VI. 66, VI. 15, VI. 5, VI. 4] производит арифметические действия с очень большой скоростью, превышающей 100 000 действий в секунду. В машину вводят исходные данные и команды с перфокарты или перфорированной ленты, которые она запоминает, затем она может выполнять команды по заданной программе. По существу машина может только складывать и сравнивать число с нулем. Все вычисления могут быть сведены к арифметическому действию сложения. Результаты вычислений, хранимые в ячейках памяти, печатаются в виде таблиц или же выдаются машиной в виде перфокарт или ленты.

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

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

Уже из этих кратких предварительных замечаний видно, что многие проблемы нецелесообразно решать при помощи электронно-вычислительной машины, при использовании которой все переменные должны быть связаны посредством математических уравнений или сведены в таблицы. Если подобные уравнения (или их приемлемые приближения) трудно получить, следует рассмотреть вопрос о применении моделирующих устройств. Но, несмотря на эти ограничения, вычислительные машины все более широко применяются при инженерных расчетах. Ниже рассматриваются принципы работы, относящиеся ко всем электронно-вычислительным машинам и лишь в отдельных случаях только к машинам модели «Дататрон» [VI. 4].

ПРОГРАММИРОВАНИЕ

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

Научиться программированию нетрудно, ибо всякий, кто решал ту или иную задачу, фактически уже программировал ее. Программирование сводится по существу к написанию инструкции по решению данной задачи для расчетчика, не знающего ни математики, ни инженерного дела, не. '-меющего пользоваться настольной вычислительной машиной. Поскольку наблюдается тенденция к использованию одноадресных электронно-вычислительных машин с хранимой в памяти программой, в данном разделе будут рассматриваться машины этого класса. Изменить программу для двухадресной машины нетрудно, по существу двухадресная машина может использоваться как одноадресная. Программист может рассматривать вычислительную машину просто как так называемый «черный ящик» и в основном не обращать внимания на детали механизма, выполняющего вычисления. Устройство электронно-вычислительных машин описывается здесь лишь в объеме, необходимом для выполнения программистом его работы.

ЗАПОМИНАЮЩЕЕ УСТРОЙСТВО ИЛИ ПАМЯТЬ

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

РЕГИСТРЫ И КОД

Программист должен также знать, что арифметическое устройство вычислительной машины имеет ряд регистров, где числа, над которыми производятся действия, появляются на различных этапах вычислений. Эти регистры аналогичны имеющимся у настольной счетной машины. Действия, производимые над числами этими регистрами, почти сходны с соответствующими действиями настольных счетных машин, но с одним важным исключением: ответ всегда появляется в одном и том же регистре. Этот регистр имеет в 2 раза больше разрядов, чем любая ячейка памяти, плюс знаковый разряд. При сложении числа ответ (или сумма плюс знаковый разряд) появляется на левой половине регистра. Когда число в этом регистре умножается на другое число, взятое из памяти, произведение плюс знаковый разряд занимают весь регистр. С другой стороны, когда число в этом регистре делится на число, взятое из памяти, частное плюс знаковый разряд появляются в левой половине регистра, а остаток в правой. Этот регистр для удобства разделен на две части. Левую часть, или ответную, называют регистром А. Правую часть, где после деления появляется остаток, называют регистром R. Число, прежде чем над ним может производиться какое-либо действие, должно находиться в регистре А. После выполнения действия ответ автоматически появляется в некоторой форме в регистрах А и R и над ним можно производить последующие действия. Если сумма двух чисел имеет большое число цифр, не помещающихся в регистре А, регистр отбрасывает разряды, находящиеся справа. Это отбрасывание опять-таки аналогично действию настольного счетного устройства с тем исключением, что электронно-вычислительная машина в данном случае автоматически остановится и вспыхнет световой сигнал, означающий, что произошло «переполнение». Поэтому данные следует располагать в каждой ячейке памяти или в регистре А так, чтобы предотвратить возникновение «переполнения». Некоторые вычислительные машины имеют другие регистры, которые могут значительно облегчить задачу программиста, но рассмотрение их выходит за границы этого элементарного изложения, Регистр А и каждая ячейка памяти содержат несколько цифровых разрядов плюс знаковый разряд. Эти цифровые разряды в совокупности со знаковым разрядом называются кодом. Код для большинства вычислительных машин состоит из десяти цифровых разрядов и знакового разряда. Длина кода является фактором, определяющим число цифр, над которыми производятся действия.

Рис, VI. 75. Блок-схема вычислений на цифровой вычислительной машине.

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

(VI. 28) (VI. 29) (VI. 30)

БЛОК-СХЕМА

F'

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

Уравнение (VI. 15) для расчета испарения имеет форму, удобную для решения его на электронно-вычислительной машине:

^y-'Zx = 'Lyf^-₎^l)_+l =F(V) = 0.

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

Рэчфорд и Райс (Rachford and Rice) [VI. 48] предложили использовать метод 20-ти итераций (попыток). Значения функции F(V) определяются для ряда значений V, выбранных таким образом, чтобы сузить область изменения переменной V, в которой находится решение уравнения [VI. 15]. Для определения корня функции F(V) возьмем сначала величину V, равную 0,5. Из рис. VI. 74 можно видеть, что если Ёг (К — 1)/[V(К — 1) + 1] — величина положительная, то значение V слишком мало, чтобы удовлетворить уравнению (VI. 15), следовательно, величина V должна быть больше 0,5. Если же значение функции отрицательное, то величина V слишком велика, чтобы удовлетворить уравнению (VI. 15), т. е. V должно быть меньше 0,5. Первая величина V сузила интервал изменения переменной V, в котором лежит решение уравнения (VI. 15), до половины предыдущего интервала изменения от 0 до 1. Теперь, если для нахождения значения функции F( V) взять другую величину V, находящуюся в середине интервала величин, в котором находится решение уравнения (VI. 15), то знак функции F(V) покажет, в какой из половин оставшегося интервала находится решение уравнения (VI. 15). Повторяя эту процедуру 20 раз, можно определить величину V, которая дает значение функции F(V), меньшее величины 10^-8. Определив величину V, являющуюся решением уравнения (VI. 15), используют эту величину в уравнениях (VI. 10) и (VI. 12). Хотя этот метод был бы слишком трудоемким для настольного счетного устройства, электронно-вычислительная машина может при его использовании получить решение через 1 мин после ввода в нее исходных данных и программы.

Команда

Очистить, сложить

Сложить

Другим подходом к задаче определения V является так называемый метод Ньютона. В данной точке определяется тангенс угла наклона функции F(V) к V. Функция F( V), деленная на тангенс угла наклона, служит поправкой, которую нужно добавить к первому взятому значению V. Новое значение V будет гораздо ближе к величине, полученной при решении уравнения (VI. 15). Этот метод состоит в следующем:

Очистить, вычесть

F (У) = 0,    (VI.    26)

F (V) = F (V;) + (У_/+1 - Vjt F‘ (V,) + . . . , (VI. 27)

где /—номер итерации. Опуская члены второго порядка, получаем

F (V,)

^Vj+\ — ^V j = — F’ (VJ "

F(V,)

^1/ж

^<V/)—    1)    П²

Так, из уравнения (VI. 29) можно найти ряд последующих приближений из величины V, пока функция F(V) не станет меньше какого-то заданного числа. Для данной задачи по методу Ньютона получают ответ после нескольких итераций. Этот метод дает сходимость быстрее, чем метод 20 итераций. Однако в большинстве случаев удовлетворительными являются оба метода. Лишь при более сложной задаче, когда уравнение приходится решать много раз, лучше переходить с программы, составленной для метода 20 итераций, на программу, составленную по методу Ньютона. Блок-схема этапов вычислений по методу Ньютона показана на рис. VI. 75. Цель этой блок-схемы — дать программисту возможность увидеть процесс вычисления в целом и установить, какие промежуточные результаты следует сохранить для последующего использования (например, результаты в блоке С, которые потребуются в блоках / и N). Блок-схема позволяет также легко перегруппировать этапы вычислений, сводя к минимуму число операций, и тем самым получить оптимальную программу.

КОМАНДЫ

Этапы блок-схемы после ее составления переводятся в команды для вычислительной машины. Команда состоит из трех частей: 1) адреса ячейки памяти, в которой будет храниться команда; 2) кода операции; 3) адреса числа, над которым должно быть произведено действие. Некоторые из возможных команд перечислены в табл. VI. 13.

Таблица VI. 13 Перечень команд для вычислительной машины

Объяснение

Арифметические команды

Очистить регистр А и прибавить содержимое ячейки, адрес которой указан в команде

Прибавить содержимое ячейки, адрес которой указан в команде, к содержимому регистра А

Очистить регистр А и вычесть содержимое ячейки, адрес которой указан в команде, из содержимого регистра А

Команда

Запись

ПРОГРАММИРОВАНИЕ

Прежде чем программировать задачу однократного испарения, приведем в качестве упрощенного примера программу для решения следующего уравнения:

x = (D + B)C,

где ?>= 110,2345, В = 22,445566 и С = 33,557799. Программа показана в табл. VI. 14 в одной из применяемых форм.

Таблица VI. 14 Программа для х = (D +В)С

Оста нов.

Переполнение «включено», знаки различны

Переход при числе, отличном от нуля

X О * ’К Н Я о> к zS з* x cq

о ^а cl и * о

и К

Команды нлн даннь

команда

Считывание, начало вычисления

0001

Переход при нуле

Условный

переход

Безусловный

переход

Считывание

Считать, начать вычисления

Начать

числения

Принятие решения и разветвление

Остановить работу машины.

Если знак регистра А отличается от знака ячейки, адрес которой указан в команде, переполнение показывает «включено»

Проверить содержимое (не знак) регистра А:

а) если содержимое регистра А равно нулю, установить знак регистра А на нуль и продолжать счет в должной последовательности; б) если содержимое регистра А не равно нулю, передать управление ячейке, адрес которой указан в команде Проверить содержимое (не знак) регистра А:

а) если содержимое регистра А не равно нулю, продолжать счет в должной последовательности; б) если содержимое регистра А равно нулю, передать управление ячейке, адрес которой указан в команде Если переполнение показывает «включено», передать управление ячейке, адрес которой указан в команде, и вновь установить переполнение.

Если переполнение показывает «выключено», продолжать счет в должной последов атель ности Передать управление ячейке, адрес которой указан в команде Считываются несколько команд и исходных данных, число которых указано в скобках, и передаются в последовательно расположенные ячейки памяти, начиная с ячейки, адрес которой указан в команде Записать содержимое нескольких последовательно расположенных ячеек памяти, число которых указано в скобках, начиная с той ячейки, адрес которой указан в команде

Читать непрерывно исходные данные или команды и хранить их в последовательно расположенных ячейках памяти, начиная с ячейки, адрес которой указан в команде, пока не подойдет команда начать вычисления; передать управление ячейке, адрес которой указан в этой команде Начать вычисления, передав управление ячейке, адрес которой указан в команде

Ввести программу с перфокарт в память вычислительной машины, начиная с ячейки памяти 0001, пока не дойдет очередь кода, в состав которого будет входить команда начать вычисления; хранить D в ячейке 0001 Хранить D в ячейке 0001 Хранить В в ячейке 0002 Хранить С в ячейке 0003 Очистить регистр и прибавить D Прибавить В (в регистре А будет тогда D + В) Умножить на С

Объяснен не

Объяснение

На первом этапе программы электронно-вычислительная машина вводит сообщаемую информацию в последовательно расположенные ячейки памяти (начиная с ячейки 0001). Когда ввод дойдет до последнего кода, на ленте или перфокарте, ввод прекратится и управление будет передано ячейке 004. В данном случае вычисления начинают по команде очистить регистр А и прибавить содержимое ячейки 0001 или D к регистру А. После того как вычислительная машина завершила это действие, выполняется следующая команда и содержимое ячейки 0002, т. е. В, прибавляется к регистру А, который будет тогда содержать сумму D + В. После этого вычислительная машина производит еще одно действие, выполняя команду умножить эту сумму на содержимое ячейки 0003. Теперь в регистрах А и R будет находиться произведение (D + В)С.

ПЕРЕСТАНОВКА ДЕСЯТИЧНОЙ ЗАПЯТОЙ

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

Для этого служат команды сдвига влево и вправо. Команда, находящаяся в ячейке 0007 (см. табл. VI. 14), перемещает произведение, находящееся в регистре, таким образом, что десятичная запятая оказывается на том же месте, что и в первоначально введенных в память данных. По следующей команде к последней цифре, находящейся в регистре А, прибавляется 1, если первая цифра в регистре R равна или больше 5. По следующим двум командам выдается (печатается) ответ (D + В)С. Команда, находящаяся в ячейке ООП, завершает процесс вычислений.

ПРОГРАММА ДЛЯ РАСЧЕТОВ ОДНОКРАТНОГО (КОНТАКТНОГО) ИСПАРЕНИЯ

В табл. VI. 15 показана программа для расчетов мгновенного испарения по методу Ньютона. Для упрощения берем систему, состоящую лишь из Десяти компонентов. Программа составлена так, что после ее ввода в вычислительную машину можно получить ряд решений, непосредственно считывая для каждого из них ^zn, К_п и первую предположительно взятую величину V. В ячейках памяти 0001—0010 находятся данные о составе исходной смеси z„, а в ячейках ООП—0020 — константы равновесия К_п для каждого компонента, расположенные в той же последовательности, что и данные

о составе каждого компонента в исходной смеси. Ячейки 0021—0090 оставлены для хранения промежуточных результатов. Величины, хранимые в ячейках 0091—0099, будут рассмотрены по мере их использования в программе. Данные о составе исходной смеси располагаются в ячейках памяти так, что десятичная запятая находится между первым и вторым разрядом слева. Десятичная запятая в константах равновесия ставится между третьим и четвертым разрядом слева, а в величинах V ставится перед первым разрядом слева.

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

Команды, находящиеся в ячейках 0100—0113, ясно показывают этапы вычислений. Прежде всего очищаются регистры и прибавляется величина, находящаяся в ячейке ООП. Поскольку в этой ячейке находится величина Ль на регистре А появляется Ки как указано в колонке объяснений. Результаты каждого этапа вычислений вплоть до ячейки 0113 показаны в колонке объяснений. Команды с 0100 до 0113 можно дополнительно повторять еще 9 раз, по разу для каждого компонента. При этом необходимо лишь изменить адреса команд в ячейках 0100, 0102, 0103, 0107, ОНО, 0111 и 0113. При каждом повторении этапа к адресам команд указанных ячеек нужно прибавить 1. Так, для Кг команда в ячейке 0100 будет «очистить, прибавить 0012», и она будет храниться в ячейке 0114.

Таблица VI. 15 Программа для расчета однократного испарения по методу Ньютона