Neural network modeling for process of oil dehydration on the basis of incomplete experimental data

Abstract


The practical experience of synthesis of the model for the water-oil emulsion separation process is considered. In modern literature, the process does not have an accurate analytical description, and its parameters are nonlinearly depend on each other. The source of the data for the model is a series of laboratory experiments conducted using certified oil test methods. The resulting tables have gaps due to an irregular step in the values of the parameters. To simulate a nonlinear process, the Levenberg-Marquardt method was used to train the neural network. An analysis was made that determined the quality criteria for the model. To automate the modeling process, the MATLAB software complex was used, as well as the Neural Network Toolbox. The training of the neural network on the original data set showed that, due to missing parameter values and a small training sample size, it is not possible to obtain an adequate neural network model without preliminary data preparation. To fill in the missing data values, the method of interpolation with Hermite splines was used in the conditions of constraints by the minimum and maximum values of the parameters from the data set. With the help of this method, it was possible to expand the training data set for the neural network by 40 times. Further training of the neural network showed that the reproduction of the training sample is performed with high accuracy, and the quality criteria of the model show acceptable values. Based on the obtained four-dimensional model, it is possible to conduct various types of analysis, both expert and automated. The results of the analysis can be the basis for issuing recommendations on the organization and conducting of the technological process of oil treatment, as well as the synthesis of an adaptive control system with a reference model.

Введение При добыче нефти важной задачей является отделение воды с растворенными в ней солями с целью предотвращения коррозии технологического оборудования и снижения расхода энергии на транспортировку балластной жидкости. Для обезвоживания водонефтяной эмульсии к ней подмешивают и равномерно распределяют по всему объему химический реагент - деэмульгатор, который представляет собой поверхностно-активное вещество. Процесс протекает за счет уменьшения силы поверхностного натяжения капель воды, находящихся в нефти, благодаря чему капли воды объединяются и оседают на дно, образуя сплошной слой. За время процесса не происходит химических реакций, то есть не происходит изменения молекулярной структуры веществ. Химический состав нефти и пластовой воды, а также наличие в них механических примесей обусловлены геологическими факторами различных географических регионов. Реологические свойства (вязкость, текучесть) и эмульгирующие способности различаются в зависимости от места добычи. Современные системы управления процессом обезвоживания не учитывают весь комплекс факторов, влияющих на скорость отделения воды от нефти [1-3]. Создание автоматизированной системы, предсказывающей состояние процесса на основании эталонной модели, могло бы повысить качество управления таким многопараметрическим процессом, как разделение водонефтяной эмульсии [4, 5]. 1. Постановка задачи Необходимо построить модель процесса разделения водонефтяной эмульсии, идентифицировать управляющие параметры и определить управляемую величину. Модель должна быть пригодна для работы в широком диапазоне параметров с достаточной для выполнения анализа и прогнозирования точностью [6, 7]. Также необходимо определить критерии, по которым будет оцениваться качество модели. Таблица 1 Результаты теплохимического обезвоживания нефти (обводненность 19,8% об., t = 40 °С) Наименование реагента Расход реагента Содержание воды в нефти , % об., после отстаивания за время , мин Отделилось воды , %, за время отстаивания(Т), мин г/м3 жидк. г/т нефти 15 30 60 90 120 15 30 60 90 120 Без реагента 0 0 0 0 0 0 0 19,8 19,8 19,8 19,8 19,8 ДЭ-1 15 21 0 2,5 5,0 7,6 10,1 19,8 19,4 19,0 18,5 18,1 30 42 0 5,2 10,4 15,6 20,8 19,8 19,0 18,1 17,4 16,3 60 84 2,5 12,6 27,8 35,4 40,4 19,4 17,7 15,1 13,8 12,8 ДЭ-2 15 21 10,4 23,4 39,1 44,3 49,5 18,1 15,9 13,1 12,1 11,1 30 42 36,8 52,6 65,8 71,7 73,7 13,5 10,5 7,8 6,7 6,1 60 84 67,7 88,5 88,5 9,1 91,1 7,4 4,0 2,8 2,3 2,3 Для построения модели процесса термохимического обезвоживания взяты отчеты лабораторных исследований нефтей различных месторождений [8]. Источником данных послужил отчет об исследования качества разделения водонефтяной эмульсии при воздействии деэмульгаторов ДЭ-1 и ДЭ-2. Эксперимент проводился с нерегулярным шагом по входным параметрам. Данные представлены в табличном виде с группировкой относительно различных значений действующих параметров. Данные для температуры нагрева 40 °С представлены в табл. 1. Значение параметра E рассчитывается по формуле , (1) где - исходная обводненность; - обводненность после эксперимента. Полнота таблицы экспериментальных данных является определяющим фактором при проведении анализа. Пропуск или недостоверность значений параметров могут привести к искажению результата процесса моделирования [9]. Наибольшую эффективность обезвоживания в представленной серии экспериментов показал деэмульгатор под обозначением ДЭ-2. Его показатели качества обезвоживания и рассмотрены в качестве источника данных. 2. Управляющие параметры Действующие параметры - температура подогрева эмульсии, количество введенного деэмульгатора в пересчете на массовый эквивалент смеси (расход в граммах на тонну), а также время отстаивания. Моделирование производится при условии, что эмульсия имеет вид «вода в нефти» (эмульсия второго рода) и процесс разделения водонефтяной эмульсии происходит однообразно с одинаковой скоростью при любом процентном содержании воды - от 0% до точки инверсии фаз (50-60 %об.), когда эмульсия меняет вид на «нефть в воде» (эмульсия первого рода). Для построения модели введен относительный параметр - эффективность разделения эмульсии , который рассчитывается по формуле (1) и принимает значения от 0 до 1. Параметр описывается нелинейной функцией вида , (2) где - расход реагента; -температура нагрева; -время отстаивания. Исходная таблица экспериментальных данных для моделирования имеет размерность 110×4, то есть 110 сочетаний значений параметров , , , . При анализе экспериментальных данных прослеживается прямая корреляция между расходом деэмульгатора и долей отделившейся воды, а также между температурой нагрева и долей отделившейся воды.Таким образом, модель процесса должна быть монотонно возрастающей в пределах рабочих значений параметров.Следовательно, качество модели определяется следующими параметрами: 1) минимальное значение дисперсии непрерывной функции; 2) неотрицательные значения частных производных на всем диапазоне моделирования ; 3) дополнительные критерии определения качества модели. 3. Построение модели экспериментальных данных В качестве инструмента создания модели экспериментальных данных технологического процесса предложено использовать искусственную нейронную сеть (ИНС). В отличие от других методов моделирования, данный способ позволяет найти зависимости в больших массивах данных, предоставляя возможности для дальнейшего подробного анализа полученной модели. Обучение нейронной сети [10-14] производилось с помощью модуля NeuralNetworkToolbox программного комплекса MATLAB. Как дополнительный критерий качества при обучении рассчитывался параметр Performance, представляющий собой интегральный показатель среднеквадратичной ошибки (функция MSE) по следующей формуле: , (3) где - значение обучающего параметра; - значение параметра на выходе нейросети; - количество обучающих примеров. Для расчета стандартного отклонения функции двух переменных используется расчет отдельно по каждой из них по формуле , (4) где - массив данных; - количество обучающих примеров. В табл. 2 представлены значения критериев качества для нейросетевых моделей с различным количеством нейронов в скрытом слое. Таблица 2 Показатели качества модели Кол-во нейронов Performance Дисперсияпо Q Дисперсия по t Знак первой производной 1 0.0177 9.2437 16.1609 Положительный 2-5 ~0.0001 9.1857 23.2762 Неотрицательный 10-20 ~10-12 10.041 20.1697 Неотрицательный 21-100 ~10-13 11.491 27.3278 Переменный При нейросетевом моделировании возникают следующие проблемы: 1) при 1 нейроне в скрытом слое моделирование производится неудовлетворительно (рис. 1,а), так как показатель MSE = 0.0177, что сравнимо с порядком величины ; 2) при небольшом количестве нейронов в скрытом слое (2-5) моделирование также производится неудовлетворительно (рис. 1,б). Показатель MSE ≈ 0.0001. В тоже время показатель отклонения по увеличивает свое значение; 3) при среднем количестве нейронов в скрытом слое (10-20) моделирование в точках известных значений производится удовлетворительно, MSE ≈ 10-12. В то же время значения отклонений высокие. На графике в местах отсутствия обучающей выборки наблюдаются провалы и пики (рис. 1,в); 4) при большом количестве нейронов (21-100) показатель Performance (2) принимает наименьшее значение ~10-13. Функция перестает быть монотонно возрастающей, так как меняется знак первой производной. На графике в областях отсутствия обучающей выборки пики и провалы еще больше (рис. 1,г). Рис. 1. Графики с выходов нейросетей с различным количеством нейронов в скрытом слое: а) нейросеть с 1 нейроном; б) нейросеть с 5 нейронами; в) нейросеть с 20 нейронами; г) нейросеть со 100 нейронами Таким образом, с прямой задачей заполнения значений между известными наборами данных нейросеть не справляется. Необходимо либо использовать нейросеть с малым количеством нейронов, либо изменить архитектуру нейросети, либо увеличить обучающую выборку. Для заполнения отсутствующих данных в таблице существуют различные методы [15]: 1) заполнение пропущенных значений усредненными данными; 2) метод ближайших соседей; 3) регрессионный метод; 4) интерполяционный метод. Для расширения обучающего набора данных был применен метод поэлементной интерполяции функции кубическими полиномами Эрмита. Отличие данного метода заключается в том, что интерполируемая функция задается не только набором данных, но и их первыми производными, что усложняет вычисления, но дает преимущество в точности. Функция имеет следующий вид: (5) где переменная представляет собой замену параметра из набора данных по формуле , (6) переменная является центральной разностной производной: . (7) При фиксировании двух параметров - температуры и времени отстаивания - получается функция одного параметра - зависимость эффективности обезвоживания от дозировки реагента-деэмульгатора. В программном комплексе MATLAB для интерполяции сплайнами Эрмита используется функция pchipinterp из пакета CurveFittingToolbox, которая лучше вcего справляется с интерполяцией данных в границе экспериментальных значений. При необходимости на небольшом доверительном интервале с помощью pchipinterp можно проводить и экстраполяцию. Рис. 2. График после расширения данных по оси расхода деэмульгатора: температура 40°С, время отстаивания 120 мин Набор данных после расширения по оси представлен на рис. 2, 3. Рис. 3. Экспериментальные данные после расширения по расходу деэмульгатора при различном времени отстаивания (0, 15, 30, 60, 90, 120 мин) Аналогичным образом проведено расширение для значений второго параметра - температуры нагрева эмульсии. При этом фиксируются значения времени и расхода реагента-деэмульгатора. Финальный этап расширения обучающего набора данных - это фиксирование значений расхода реагента-деэмульгатора и температуры нагрева эмульсии для интерполяции по времени отстаивания. 4. Обучение нейросети на расширенных данных Изначальный массив экспериментальных данных имел размер 110×4. После расширения данных полученный массив имеет в своем составе 4275×4, то есть 17100 элементов, среди которых 4275 целевых выходных значений, на основе которых производится обучение нейронной сети. Для этого используется модуль NeuralNetworkToolbox из состава программного комплекса Matlab. Подбор оптимальной структуры нейронной сети произведен экспериментально. Обучение проведено методом Левенберга - Марквардта. Параметр Performance = 3.18·10-6. Поверхность, которую генерирует обученная нейронная сеть, представлена на рис. 4. Показатели качества обученной нейросети, обученной на расширенном массиве данных, представлены в табл. 3. Рис. 4. Выходные данные обученной на расширенном наборе данных нейросети Таблица 3 Показатели качества итоговой нейросети Кол-во нейронов Performance Дисперсияпо Q Дисперсия по t Знак первой производной 100 5.23*10-7 8.9224 15.3233 Неотрицательный В результате расширения данных обученная нейросеть корректно воспроизводит экспериментальные данные в рабочем диапазоне параметров. Показатели качества имеют удовлетворительные значения, функция монотонно возрастает, показатели отклонений имеют наименьшее значение по сравнению с исходным набором данных. Заключение С помощью методов интерполяции и возможностей нейронных сетей по моделированию графиков удалось создать модель процесса в виде MISO-системы на основе экспериментальных данных, количество которых изначально не позволяло сразу приступить к обучению нейронной сети и получить представительную модель. Итоговая точность модели процесса в рамках предполагаемых способов использования считается приемлемой. Полученная нами нейросетевая модель процесса обезвоживания нефти может использоваться в качестве эталонной модели в адаптивной системе управления на реальном производственном объекте.

Iliya V Artyushkin

Joint Stock CompanyInstitute for Design and Research in Oil Industry «Giprovostokneft»

93, Krasnoarmeyskaya st.,Samara,443041, Russian Federation Postgraduate Student

  1. Позднышев Г.Н. Стабилизация и разрушение нефтяных эмульсий. - М: Недра, 1982.
  2. Тронов В.П. Промысловая подготовка нефти. - Казань, 2000.
  3. Путохин В.С. Математическое моделирование технологического процесса обезвоживания нефти на промыслах // Нефть и газ. - М.:МИНХ и ГП, 1977. - С. 37-42.
  4. Веревкин А.П.,Ельцов И.Д., Зозуля Ю.И., Кирюшин О.В.Оперативное управление технологическими процессами под готовки нефти по технико-экономическим показателям // Автоматизация, телемеханизация и связь в нефтяной промышленности. - 2006. - № 3. -ОАО «ВНИИОЭНГ». - С. 48-53.
  5. Андреев Е.Б.,Ключников А.И., Кротов А.В., Попадько В.Е., Шарова И.Я.Автоматизация технологических процессов добычи и подготовки нефти и газа: Учеб.пособие для вузов. - М.: ООО «Недра-Бизнесцентр», 2008. - 399 с: ил.
  6. Артюшкин И.В. Возможность создания комплексной экспертной автоматизированной системы управления процессом промысловой подготовки нефти // Нефтяное хозяйство. - 2016. - Вып. 6. - С. 29-31.
  7. Артюшкин И.В.,Максимов А.Е.Разработка автоматической системы управления процессом термохимического обезвоживания нефтяных эмульсий на основе искусственной нейронной сети // Вестник Самарского государственного технического университета. Сер. Технические науки. - 2017. - Вып. 1 (53). - С. 7-15.
  8. Бортников А.Е.,Кордик К.Е., Савиных А.В., Ницин А.С. Некоторые результаты лабораторных экспериментов по разрушению водонефтяной эмульсии под воздействием равномерного электрического поля// Геология, геофизика и разработка нефтяных и газовых месторождений. - 2013. - Вып. 9. - С. 48-56.
  9. Jerzy W. Grzymala-Busse, Teresa Mroczek. Definabilityin Mining In complete Data // Procedia Computer Science, Volume 96, 2016. - pp. 179-186.
  10. Аксенов С.В.,Новосельцев В.Б.Организация и использование нейронных сетей (методы и технологии). - Томск, 2006.
  11. Каплан Р. Основные концепции нейронных сетей. - М., 2001.
  12. Neural Network Software, About NeuroSolutions. URL: http:// www.neuroproject.ru/aboutproduct.php (дата обращения 02.09.2018).
  13. ХайкинС. Нейронныесети: полныйкурс = NeuralNetworks: AComprehensiveFoundation:2-е изд. - М.: Вильямс, 2006.
  14. Рутковская Д.,Пилиньский М., Рутковский Л.Нейронные сети, генетические алгоритмы и нечеткие системы. - М., 2006.
  15. Сташкова О.В.,Шестопал О.В.Использование искусственных нейронных сетей для восстановления пропусков в массиве исходных данных // Известия вузов. Северо-Кавказский регион. Сер. Техническиенауки. - Ростов-н/Д, 2017. - Вып. 1. -С. 37-42.

Views

Abstract - 13

PDF (Russian) - 4

Cited-By


PlumX

Dimensions

Refbacks

  • There are currently no refbacks.

Copyright (c) 2018 Samara State Technical University

Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies