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

Полный текст

Одной из актуальных проблем промысловой подготовки нефти является то, что в определенные периоды времени не удается получить товарную нефть, соответствующую нормативным требованиям. К основным причинам можно отнести изменение состава и расхода исходной сырой нефти, меняющиеся погодные условия и различные нарушения в работе оборудования [1]. Все это приводит к необходимости постоянной корректировки режима работы установки подготовки нефти (УПН) [2]. Однако поиск необходимых значений параметров режима сопряжен с проведением большого количества сложных технологических расчетов [3–5], что затруднительно в условиях промысла.

Ранее в работах [6–7] было предложено одно из решений, способное оказать персоналу УПН поддержку в поиске новых значений параметров технологического режима. Данная статья является развитием этого решения. В ней рассматривается создание адаптивной системы поддержания качества (АСПК) подготовки нефти к транспортировке. АСПК представляет собой специализированную компьютерную программу, которая интегрируется в действующую АСУТП установки. Ее целью является прогнозирование значений давления насыщенных паров (ДНП) и содержания сероводорода в товарной нефти и выдача рекомендаций по корректировке режима работы УПН в случае их отклонения от нормативных значений.

Технологическая схема процесса подготовки представлена на рис. 1. Сырая нефть проходит три ступени сепарации С1, С2, С3, на которых последовательно понижается давление. В результате из нее удаляется попутный нефтяной газ, содержащий сероводород. Для доведения ДНП нефти до нормативного значения перед третьей ступенью сепарации С3 предусмотрен ее нагрев в теплообменниках Т1 и Т2. Теплообменник Т1 служит для предварительного подогрева за счет рекуперации тепла. Окончательный нагрев нефти в Т2 производится с помощью промежуточного теплоносителя, приходящего из печи П1. Повышение температуры также улучшает процесс отстаивания воды в отстойнике О1. Оставшийся сероводород нейтрализуется химическим реагентом, подаваемым насосом Н1 в емкость Е1.

Рис. 1. Схема технологического процесса подготовки нефти

Схема интеграции АСПК в систему сбора и подготовки нефти показана на рис. 2. Объектом управления является технологический процесс подготовки нефти. Его спецификой является то, что показатели качества целевой продукции не могут быть измерены непосредственно в ходе процесса, а определяются в химических лабораториях. Поэтому управление процессом по отклонению этих параметров от нормативных значений в настоящее время не представляется возможным.

В этой связи предлагаемая АСПК представляет собой систему управления по возмущающим воздействиям. Непосредственно измеряемыми параметрами возмущающих воздействий являются значения расходов нефти Q_q ($q=\overline{\mathrm{1,}{N}_q}$ – номер скважины, N_q – общее количество скважин), поступающие с групповых замерных установок (ГЗУ). Возмущающие воздействия со стороны окружающей среды и изменения эффективности работы оборудования могут быть косвенно отражены реальными значениями непосредственно измеряемых параметров технологического режима. Для сепаратора С1 такими параметрами являются температура T₃, давление Р₂ и расход G₂ выходящих потоков (индексы соответствуют номерам потоков на рис. 1). Аналогично для сепараторов С2 и С3 такими параметрами будут Т₅, Р₄, G₄ и Т₈, Р₁₄, G₁₄ соответственно. Для теплообменников T1 и T2 это температуры выходящих потоков T₆, T₇ и T₉, T₁₁ соответственно.

Химические составы пластовых нефтей ${\overline{x}}_q=\left(x_{q1},x_{q2}\mathrm{,...,}{x}_{qNj}\right)$ (N_j – количество веществ в химическом составе) очень медленно изменяются с течением времени, поэтому принимаются в АСПК постоянными. Однако при необходимости их значения могут быть скорректированы в процессе работы АСПК (например, при вводе в эксплуатацию новых скважин).

Рис. 2. Схема интеграции АСПК в систему сбора и подготовки нефти

Данные по показателям качества товарной нефти (давление насыщенных паров Р_нп и содержание сероводорода x_Н2S), периодически получаемые из химических лабораторий, вносятся в АСПК и используются для уточнения параметров математической модели.

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

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

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

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

Разработанная АСПК выдает рекомендации для ЛПР по температуре теплоносителя на выходе из печи $T_{печь}^{рек}$ и расходу реагента для нейтрализации сероводорода $G_{реаг}^{рек}$. Значения давлений на ступенях сепарации поддерживаются соответствующими подсистемами стабилизации имеющейся АСУТП на уровне регламентных значений и не меняются с течением времени. Процесс сепарации является адиабатическим, поэтому каждому составу поступающей порции нефти будет соответствовать вполне определенное значение температуры теплоносителя и расхода реагента. Они могут быть определены с помощью математической модели, описывающей квазистационарный процесс подготовки нефти данного состава.

В ходе оценки временных параметров АСПК, во-первых, прогнозирует составы поступающей нефти ${\overline{x}}_{сум}^{расч}$ и соответствующие показатели качества готовой продукции $P_{нп}^{расч}$ и ${х}_{H2S}^{расч}$ с упреждением на некоторое время. Для этого АСПК по модели квазистационарного процесса выполняет ряд расчетов на ожидающиеся в ближайшее время составы нефти. Исходя из этих прогнозов ЛПР имеет возможность по своему запросу получить от АСПК рекомендуемые значения параметров технологического режима $T_{печь}^{рек}$ и $G_{реаг}^{рек}$ для любой ожидающейся порции нефти для времени t_выбр.

Во-вторых, АСПК определяет время внесения этих изменений ( $t_{печь}^{рек}$ и $t_{реаг}^{рек}$) в технологический режим. В силу того, что [8]:

1) изменение расходов нефти, поступающей на установку, как правило, происходит не чаще, чем раз в сутки;

2) влияние погодных условий характеризуется преимущественно суточной периодичностью (день/ночь);

3) диапазон изменений температур теплоносителя, как правило, находится в пределах 10 °С, что позволяет достаточно быстро менять режим,

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

Исходя из вышесказанного алгоритм для ЛПР строим следующим образом.

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

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

Схема АСПК промысловой подготовки нефти показана на рис. 3. Данные о расходах Q_q с ГЗУ поступают в модуль смешения нефтей идентифицируемой математической модели (ИММ). На их основании рассчитывается мольный расход $G_{сум}^{расч}$ и ${\overline{x}}_{сум}^{расч}$ состав  поступающей на подготовку суммарной нефти, которые передаются в расчетный модуль технологического процесса ИММ. В него же поступают данные о текущих значениях параметров технологического режима (P₂, T₃, G₂, P₄, T₅, G₄, P₁₄, T₈, G₁₄, T₆, T₇, T₉, T₁₁). На их основе прогнозируются показатели качества готовой продукции: давление насыщенных паров $P_{нп}^{расч}$ и содержание сероводорода ${х}_{H2S}^{расч}$.

Эти данные поступают в режимный модуль ИММ для сравнения с нормативными значениями $P_{нп}^{норм}$ и ${х}_{H2S}^{норм}$. В случае их отклонения определяются новые значения управляющих параметров: температуры теплоносителя на выходе из печи $T_{печь}^{рек}$ и расхода реагента для нейтрализации сероводорода $G_{реаг}^{рек}$. В этом же модуле определяются времена внесения изменений в режим работы установки $t_{печь}^{рек}$ и $t_{реаг}^{рек}$. В конечном итоге вся информация передается персоналу УПН, который принимает решение и вносит необходимые изменения в технологический процесс с помощью действующей системы АСУТП.

При периодическом появлении новых результатов лабораторных анализов качества товарной нефти P_нп и x_H2S и их отклонения от расчетных значений модули идентификации ${\gamma }_{C_1-C_4}$ и ${\gamma }_{H_2{S}_{}}$ проводят корректировку поправочных коэффициентов (см. ниже).

Описанный алгоритм работы АСПК представлен на рис. 4.

Особенностью модуля смешения нефтей ИММ является то, он отражает различное удаление нефтедобывающих скважин от УПН и смешение в трубопроводах системы сбора потоков добытой жидкости в различных сочетаниях. В этой связи для прогнозирования состава и расхода суммарной сырой нефти необходимо учитывать задержки t_сбор, связанные с движением жидкости по трубопроводам [9].

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

Рис. 3. Схема АСПК промысловой подготовки нефти:

t_тек – текущий момент времени; t_выбр – выбранный ЛПР момент времени;

МИС1, МИС2, МИС3 – модули идентификации параметров сепараторов С1, С2, С3;

МИТ1, МИТ2 – модули идентификации параметров теплообменников Т1, Т2

Рис. 4. Алгоритм работы АСПК

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

$Q_{тр}\left(t\right)=Q_{тр1}\left(t\right)+Q_{тр2}\left(t\right)$.(1)

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

$V_{тр}=\underset{t-{\tau }_{сбор}}{\overset{t}{\int }}\left[Q_{тр1}\left(t\right)+Q_{тр2}\left(t\right)\right]dt$. (2)

Так как расходы входящих потоков нефти Q_тр1(t) и Q_тр2(t) являются положительными величинами, то уравнение (2) всегда имеет единственное решение.

Рис. 5. Расчетная схема бинарного дерева системы сбора нефти

Рис. 6. Расчетный элемент системы сбора нефти

Вычисленное из (2) время t_сбор задержки (для каждого трубопровода оно свое) позволяет определить состав потока, выходящего из трубопровода в момент времени t:

${\overline{x}}_{тр}\left(t\right)=\frac{Q_{тр1}\left(t-{\tau }_{сбор}\right){\overline{x}}_{тр1}\left(t-{\tau }_{сбор}\right)+Q_{тр2}\left(t-{\tau }_{сбор}\right){\overline{x}}_{тр2}\left(t-{\tau }_{сбор}\right)}{Q_{тр1}\left(t-{\tau }_{сбор}\right)+Q_{тр2}\left(t-{\tau }_{сбор}\right)}$.(3)

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

Для трубопроводов, в которых не происходит смешения, а поступает нефть со скважины, для расчета принимается Q_тр2(t) = 0. Суммарные расход $G_{сум}^{расч}$ и состав ${\overline{x}}_{сум}^{расч}$ сырой нефти будут равны расходу и составу потока, выходящего из корневого трубопровода.

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

$P_{нп}^{расч}=P_{нп}^{расч}\left(G_{сум}^{расч},{\overline{x}}_{сум}^{расч},P_2,T_3,G_2,P_4,T_5,G_4,P_{14},T_8,G_{14},T_6,T_7,T_9,T_{11}\right)$,(4)

$x_{H2S}^{расч}=x_{H2S}^{расч}\left(G_{сум}^{расч},{\overline{x}}_{сум}^{расч},P_2,T_3,G_2,P_4,T_5,G_4,P_{14},T_8,G_{14},T_6,T_7,T_9,T_{11}\right)$.(5)

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

Математическая модель сепаратора определяет расходы G_m, G_n и составы ${\overline{x}}_m,$ ${\overline{x}}_n$  выходящих из него газа и жидкости в зависимости от условий в нем:

$G_k^{}=G_m^{}+G_n^{}$;(6)

$G_k^{}\cdot {х}_{kj}=G_m^{}\cdot {х}_{mj}+G_n^{}\cdot {х}_{nj}$;(7)

$G_k^{}\cdot H_k=G_m^{}\cdot H_m\left({Т}_m,P_m\right)+G_n^{}\cdot H_n\left({Т}_n,P_n\right)$;(8)

$G_m^{}=e_i\cdot G_k$;(9)

$G_n^{}=(1-e_i)\cdot G_k$;(10)

${\phi }_j=\frac{{х}_{mj}}{{х}_{nj}}={\gamma }_j{\eta }_i\cdot \frac{f_{nj}^{}\left(T_n,P_n\right)}{f_{mj}^{}\left(T_m,P_m\right)}$;(11)

$\sum _j\frac{x_j}{1+e_i\cdot \left({\phi }_j-1\right)}=1$,(12)

где H_k, H_m, H_n – энтальпии потоков;

e_i – мольная доля газа;

f_j – константы фазового равновесия веществ;

f _nj, f_mj – фугитивности веществ;

h_i – испарительный КПД;

g_j – поправочные коэффициенты;

i – номер аппарата;

j – номер вещества;

k, m, n – номера потоков (см. рис. 1).

Для сепаратора С1: i=1, k=1, m=2, n=3, для сепаратора С2: i=2, k=3, m=4, n=5, для сепаратора С3: i=3, k=12, m=14, n=8.

Математическая модель теплообменника [10] связывает температуры T_m, T_u выходных потоков с тепловой нагрузкой E аппарата:

$G_k^{}=G_m^{}$;(13)

$G_k^{}\cdot x_{kj}=G_m^{}\cdot x_{mj}$;(14)

$G_k^{}\cdot {Н}_k\left(T_k,P_k\right)+E_i=G_m^{}\cdot H_m\left(T_m,P_m\right)$;(15)

$G_n^{}=G_u^{}$;(16)

$G_n^{}\cdot {Н}_n\left(T_n,P_n\right)-E_i=G_u^{}\cdot H_u\left(T_u,P_u\right)$;(17)

$G_n^{}\cdot x_{nj}=G_u^{}\cdot x_{uj}$;(18)

$E_i^{}=K_i\left(r_i\right)\cdot F_i\cdot \Delta T_i$,(19)

где K – коэффициент теплопередачи;

r – термическое сопротивление загрязнений;

F – поверхность теплообмена;

DT – средняя разность температур;

i – номер аппарата (см. ниже);

k, m, n, u – номера потоков (см. ниже и рис. 1).

Для теплообменника Т1: i=1, k=6, m=5, n=7, u=8, для теплообменника Т2: i=2, k=6, m=9, n=10, u=11.

Математическая модель отстойника определяет расходы и составы выходящих из него жидкостей:

$G_n^{}=G_k^{}\cdot x_{k\beta }^{}$;(20)

$G_m^{}=G_k^{}-G_n^{}$;(21)

$x_{nj}^{}=\left\{\begin{array}{c}\mathrm{1,}приj=\beta ;\\ \mathrm{0,}приj\ne \beta ;\end{array}\right.$(22)

$x_{mj}^{}=\left\{\begin{array}{c}\frac{x_{kj}\cdot G_k}{G_m},приj\ne \beta ;\\ \mathrm{0,}приj=\beta ;\end{array}\right.$(23)

$G_m^{}\cdot H_m\left(T_m,P_m\right)=G_k^{}\cdot H_k\left(T_k,P_k\right)-G_n^{}\cdot H_n\left(T_n,P_n\right)$;(24)

где b – номер вещества H₂O. Для отстойника О1: k=9, m=12, n=13 (см. рис. 1).

Математическая модель точки смешения газов сепарации определяет расход и состав суммарного газа:

$G_k^{}=G_m^{}+G_n^{}$;(25)

$G_k^{}\cdot x_{kj}=G_m^{}\cdot x_{mj}+G_n^{}\cdot x_{nj}$;(26)

$G_k^{}\cdot H_k\left(T_k,P_k\right)=G_m^{}\cdot H_m\left(T_m,P_m\right)+G_n^{}\cdot H_n\left(T_n,P_n\right)$.(27)

Для точки смешения М: k=15, m=4, n=14 (см. рис. 1).

Математическая модель печи связывает параметры входного и выходного потоков теплоносителя и тепловую мощность:

$G_k^{}=G_m^{}$;(28)

$G_k^{}\cdot x_{kj}=G_m^{}\cdot x_{mj}$;(29)

$x_{mj}^{}=\left\{\begin{array}{c}\mathrm{1,}приj=\beta ;\\ \mathrm{0,}приj\ne \beta ;\end{array}\right.$(30)

$E_i^{}=G_k\cdot H_k\left(T_k,P_k\right)-G_m\cdot H_m\left(T_m,P_m\right)$.(31)

Для печи П1: i=1, k=10, m=11 (см. рис. 1).

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

$G_n^{}={\psi }_i{G}_k^{}\cdot x_{k\xi }^{}$,(32)

где $\psi$ – удельный расход реагента;

$\xi$ – номер вещества H₂S.

Для емкости нейтрализации Е1: i=1, k=7, n=16 (см. рис. 1).

В основу расчета фазовых равновесий и теплофизических свойств материальных потоков положено кубическое уравнение состояния Пенга – Робинсона [11–12]. Согласно этому уравнению зависимости фугитивностей и энтальпий от температуры определяются выражениями [10]:

$\ln f_j\left(T,P\right)=\frac{B_j}{B}\left(Z-1\right)-\ln \left(Z-B\right)+\frac{A}{\mathrm{4,828}B}\left[\frac{B_j}{B}-\frac{2{\displaystyle \sum _{k=1}^N{x}_j{\left(a{\alpha }_{jk}\right)}^{\mathrm{0,5}}}}{a\alpha }\right]\ln \frac{Z+\mathrm{2,414}B}{Z-\mathrm{0,414}B}$, (33)

$H\left(T,P\right)=H_{ид}\left(T\right)+RT\left[1-Z+\frac{A}{\mathrm{2,828}B}\left(1+\frac{D}{a\alpha }\right)\ln \frac{Z+\mathrm{2,414}B}{Z-\mathrm{0,414}B}\right]$. (34)

В уравнениях (33)–(34) параметры a, a, A, B, D определяются через температуры, давления и критические свойства веществ по известным формулам [11]. Сжимаемость фаз Z рассчитывается путем решения кубического уравнения [11]

$Z^3-\left(1-B\right)Z^2+\left(A-2B-3B^2\right)Z-\left(AB-B^2-B^3\right)=0$.(35)

Зависимость энтальпии идеального газа H_ид(T) от температуры представляет собой аппроксимацию экспериментальных данных для химических веществ [13].

В силу сложности моделируемого процесса и многокомпонентности фаз аналитическое решение уравнений (6)–(35) является крайне затруднительной задачей. Для численного решения используется специализированный программный продукт «МиР ПиА Процесс» [14]. Он позволяет проводить численное моделирование технологических процессов и аппаратов в различных отраслях промышленности [15–16]. Модель УПН, реализованная с помощью программного продукта «МиР ПиА Процесс», представлена на рис. 7.

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

Для идентификации моделей каждого из сепараторов C1, C2, C3 (модули МИС1 – МИС3 на рис. 3) используются измеренные данные о температуре T₃, T₅, T₈, давлении P₂, P₄, P₁₄ и расходе газа на выходе G₂, G₄, G₁₄. При этом значения температур T₃, T₅, T₈ и давлений P₂, P₄, P₁₄ напрямую вносятся в математическую

Рис. 7. Модель УПН, реализованная с помощью программного продукта «МиР ПиА Процесс».

 

модель (6)–(35). Затем определяется значение испарительного КПД ступени сепарации h_i, при котором расчетный расход газа $G_m^{расч}$ равен реальному G_m. Это равенство обеспечивается за счет уменьшения констант фазового равновесия f_j на величину испарительного КПД:

${\phi }_j={\gamma }_j{\eta }_i\frac{f_j^{ж}}{f_j^{г}}$.(35)

Типовая зависимость расчетного расхода газа от величины КПД приведена на рис. 8а. Монотонный характер этой зависимости позволяет найти единственное решение задачи идентификации модели сепаратора.

Схожим образом проводится идентификация модели теплообменников (модули МИТ1 – МИТ2, рис. 3). При этом для каждого из них используются измеренные значения температур выходных потоков. Здесь варьируемой величиной выступает коэффициент термических сопротивлений загрязнений r, влияющий на общий коэффициент теплопередачи аппарата. Величина коэффициента r находится из условия равенства расчетной выходной температуры потока реальному значению.

Типовая зависимость выходной температуры от коэффициента r также носит монотонный характер (рис. 8б), что позволяет найти единственное решение задачи.

Результатом решения уравнений (6)–(35) является определение значений показателей качества товарных потоков. Эти данные поступают в блок выбора режима для сравнения с нормативными значениями (см. рис. 3). В случае отклонения найденных значений от нормативных режимный модуль ИММ определяет температуру теплоносителя на выходе из печи $T_{печь}^{рек}$ и расход реагента для нейтрализации сероводорода $G_{реаг}^{рек}$. Эта задача имеет единственное решение и не требует применения многомерных методов поиска. Это возможно, так как $P_{нп}^{расч}$ зависит только от $T_{печь}^{рек}$, $G_{реаг}^{рек}$ а  рассчитывается исходя из остаточного содержания сероводорода при температуре $T_{печь}^{рек}$.

Поэтому в первую очередь определяется температура T_печь теплоносителя на выходе из печи, обеспечивающая требуемое значение давление насыщенных паров $P_{нп}^{норм}$. Решение этой обратной задачи осуществляется в режимном модуле ИММ методом многократного решения прямой задачи.

Рис. 8. Типовые зависимости, получаемые и используемые в ходе идентификации моделей сепаратора (а) и теплообменника (б)

На рис. 9 показан типовой график зависимости давления насыщенных паров от температуры теплоносителя. Монотонный характер представленной зависимости позволяет легко определить решение T_печь, например методом половинного деления.

Рис. 9. Зависимость давления насыщенных паров от температуры теплоносителя

После определения температуры теплоносителя $T_{печь}^{рек}$ рассчитывается количество реагента, необходимого для нейтрализации сероводорода $G_{реаг}^{рек}$, по формуле (31).

Рис. 10. Временные диаграммы изменения потоков УПН и транспортных задержек управляющих воздействий

Представленная выше математическая модель (6)–(35) описывает стационарный режим работы УПН, а корректировки технологического режима должны вноситься в разных точках технологической схемы. Время их внесения рассчитывается с учетом временной задержки на транспортировку. Это время определяется на основе модели идеального вытеснения и индивидуально для каждой конкретной УПН. На рис. 10 показаны типовые временные диаграммы изменения потоков УПН и транспортных задержек управляющих воздействий. Эти временные диаграммы относятся к простейшему случаю скачкообразного изменения состава ${\overline{x}}_{сум}^{расч}$ исходной нефти (график а). При этом изменение давления насыщенных паров товарной нефти $P_{нп}^{}$ и содержания сероводорода ${х}_{H2S}^{}$ в ней происходит синхронно в конце цикла подготовки (графики б и в). График г отображает время внесения корректировки расхода реагента $t_{реаг}^{рек}$. Изменение температуры нагрева нефти ${Т}_{печь}^{}$ осуществляется с учетом инерции контура нагрева (график д).

Если, как отмечалось ранее, инерционность подачи реагента в модели (6)–(35) не учитывается, то изменение расхода реагента происходит в момент поступления потока новой нефти в емкость Е1. Однако в силу определенной инерционности контура нагрева теплоносителя регулировка печи должна происходить заблаговременно. Тогда для определения времени упреждения применяется эмпирическая формула, описывающая изменение во времени температуры потока 9, выходящего из теплообменника Т2 (см. рис. 11), как апериодического звена с запаздыванием

$T_9\left(t\right)=T_9^{кон}-\left(T_9^{кон}-T_9^{нач}\right)\exp \left(-\frac{t-t_{печь}^{рек}}{{\tau }_{цирк}}\right)$,(36)

где $T_9^{нач}$, $T_9^{кон}$  – начальное и конечное значения температуры, соответствующие старому и новому технологическому режиму;

${\tau }_{цирк}$ – постоянная времени контура циркуляции теплоносителя, определяемая экспериментально для конкретной УПН.

Рис. 11. График переходного процесса в теплообменнике T2

Тогда время упреждения для печи будет равно (при допустимом отклонении температуры на 0,1 °С)

${\tau }_{печь}={\tau }_{цирк}\ln \frac{T_9^{кон}-T_9^{нач}}{\mathrm{0,1}}$.(37)

Как отмечалось выше, внедрение АСПК не отменяет необходимости периодического проведения лабораторных анализов. Их результаты используются для уточнения параметров g_j математической модели (6)–(35). Необходимость проведения корректировки обусловлена тем, что в ходе эксплуатации меняются параметры технологического оборудования, что влияет на интенсивность испарения отдельных веществ [17]. В случае отклонения расчетного давления $P_{нп}^{расч}$ от действительного $P_{нп}^{}$ проводится корректировка значения ${\gamma }_{C_1-C_4}$ для легких компонентов от метана (C₁) до бутана (C₄) включительно. В случае отклонения расчетного содержания сероводорода ${х}_{H2S}^{расч}$ от действительного ${х}_{H2S}^{}$ проводится корректировка значения g_H2S  для сероводорода. Решение этих обратных задач проводится аналогично решению обратной задачи по определению температуры теплоносителя на выходе из печи.

Наибольшую практическую пользу внедрение предложенной АСПК может принести для УПН, получающих сырую нефть с небольшого числа скважин. Эти установки являются наиболее чувствительными, так как изменение объемов добычи или остановка на ремонт даже одной скважины может оказать существенное влияние на состав исходной нефти. Однако и для больших УПН внедрение АСПК является актуальным, так как позволит избежать получения некондиционной нефти при авариях в системе сбора нефти, а также при резкой смене параметров окружающей среды.

Об авторах

Сергей Борисович Коныгин

Самарский государственный технический университет

Автор, ответственный за переписку.
Email: info@eco-vector.com

д.т.н., доцент, заведующий

Россия, Россия, 443100, г. Самара, ул. Молодогвардейская, 244

Денис Владимирович Коноваленко

Самарский государственный технический университет

Email: info@eco-vector.com

старший преподаватель

Россия, Россия, 443100, г. Самара, ул. Молодогвардейская, 244

Дмитрий Александрович Крючков

Самарский государственный технический университет

Email: info@eco-vector.com

к.т.н., доцент

Россия, Россия, 443100, г. Самара, ул. Молодогвардейская, 244

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

Дополнительные файлы