Научная статья на тему 'МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ УЛЬТРАЗВУКОВОЙ ОБРАБОТКИ АСФАЛЬТОСМОЛОПАРАФИНОВЫХ ОТЛОЖЕНИЙ'

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ УЛЬТРАЗВУКОВОЙ ОБРАБОТКИ АСФАЛЬТОСМОЛОПАРАФИНОВЫХ ОТЛОЖЕНИЙ Текст научной статьи по специальности «Химические технологии»

CC BY
232
41
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ / УЛЬТРАЗВУК / УЛЬТРАЗВУКОВАЯ / УСТАНОВКА / НЕФТЬ / АСФАЛЬТОСМОЛИСТЫЕ ПАРАФИНОВЫЕ ОТЛОЖЕНИЯ / НАСОСНОЕ ОБОРУДОВАНИЕ / СКВАЖИНА / РЕЗЕРВУАР / НЕФТЕПРОВОД / ВТОРИЧНОЕ СЫРЬЕ

Аннотация научной статьи по химическим технологиям, автор научной работы — Хурамшина Регина Азатовна, Соколова Виктория Владимировна

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по химическим технологиям , автор научной работы — Хурамшина Регина Азатовна, Соколова Виктория Владимировна

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

MATHEMATICAL MODELLING OF ULTRASONIC TREATMENT OF ASPHALTENE DEPOSITS

Objective. The removal of asphaltene deposits at oil and gas facilities is one of the urgent and important problems and requires significant material and labor costs. It is possible to reduce costs by creating and implementing effective technical means, which requires an in-depth study of the processes of organic matter deposition at oil and gas facilities and their use as a secondary raw material. Methods. This paper discusses modern views on the state of the problem of asphaltene deposits in oil shipping and storage equipment and possible ways to solve it. The paper provides an overview of various ways to clean shipping and storage objects from asphaltene deposits: chemical (adding additives, solvents), thermal (heating by special devices or injection of superheated steam during exploitation), mechanical (using scrapers and pistons), and refers, among other things, to scientific works on the use of ultrasound to accelerate the removal of deposits. Results. The paper considers methods for removing deposits, as well as using the positive effect of the removed layer as a secondary energy source. A procedure for model calculation of the use of ultrasonic equipment to remove deposits has been developed. As a result, the deposit melting front velocity was determined depending on the duration of exposure. Conclusion. Taking into account the positive world experience, the level of development of the ultrasonic method for removing asphaltene deposits in the oil and gas industry and the use of asphaltene deposits as a secondary raw material, this area needs further development. The widespread implementation of equipment and, from the standpoint of rational use of natural resources, the use of deposits as a secondary raw material will increase cost efficiency and equipment efficiency, and reduce environmental impact.

Текст научной работы на тему «МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ УЛЬТРАЗВУКОВОЙ ОБРАБОТКИ АСФАЛЬТОСМОЛОПАРАФИНОВЫХ ОТЛОЖЕНИЙ»

Вестник Дагестанского государственного технического университета. Технические науки. Том 48, №2, 2021 Herald of Daghestan State Technical University. Technical Sciences. Vol.48, No.2, 2021 _http://vestnik.dgtu.ru/ISSN (Print) 2073-6185 ISSN (On-line) 2542-095Х_

Для цитирования: Р.А. Хурамшина, В.В. Соколова. Математическое моделирование ультразвуковой обработки асфальтосмолопарафиновых отложений. Вестник Дагестанского государственного технического университета. Технические науки. 2021; 48 (2): 60-72. DOI:10.21822/2073-6185-2021-48-2-60-72

For citation: R.A. Khuramshina, V. V. Sokolova. Mathematical modelling of ultrasonic treatment of asphaltene deposits. Herald of Daghestan State Technical University. Technical Sciences. 2021; 48 (2): 60-72. (In Russ.) D0I:10.21822/2073-6185-2021-48-2-60-72

ИНФОРМАТИКА, ВЫЧИСЛИТЕЛЬНАЯ ТЕХНИКА И УПРАВЛЕНИЕ COMPUTER SCIENCE, COMPUTER ENGINEERING AND MANAGEMENT

УДК 622.692

DOI: 10.21822/2073-6185-2021 -48-2-60-72

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ УЛЬТРАЗВУКОВОЙ ОБРАБОТКИ АСФАЛЬТОСМОЛОПАРАФИНОВЫХ ОТЛОЖЕНИЙ Р.А. Хурамшина, В.В. Соколова

Уфимский государственный нефтяной технический университет, 450062, г. Уфа, ул. Космонавтов 1, Россия

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

Ключевые слова: математическое моделирование, ультразвук, ультразвуковая установка, нефть, асфальто-смолистые парафиновые отложения, насосное оборудование, скважина, резервуар, нефтепровод, вторичное сырье

MATHEMATICAL MODELLING OF ULTRASONIC TREATMENT OF ASPHALTENE DEPOSITS R.A. Khuramshina, V. V. Sokolova

Ufa State Petroleum Technological University, 1 Kosmonavtov St., Ufa 450064, Russia

Abstract. Objective. The removal of asphaltene deposits at oil and gas facilities is one of the urgent and important problems and requires significant material and labor costs. It is possible to reduce costs by creating and implementing effective technical means, which requires an in-depth study of the processes of organic matter deposition at oil and gas facilities and their use as a secondary raw

60

material. Methods. This paper discusses modern views on the state of the problem of asphaltene deposits in oil shipping and storage equipment and possible ways to solve it. The paper provides an overview of various ways to clean shipping and storage objects from asphaltene deposits: chemical (adding additives, solvents), thermal (heating by special devices or injection of superheated steam during exploitation), mechanical (using scrapers and pistons), and refers, among other things, to scientific works on the use of ultrasound to accelerate the removal of deposits. Results. The paper considers methods for removing deposits, as well as using the positive effect of the removed layer as a secondary energy source. A procedure for model calculation of the use of ultrasonic equipment to remove deposits has been developed. As a result, the deposit melting front velocity was determined depending on the duration of exposure. Conclusion. Taking into account the positive world experience, the level of development of the ultrasonic method for removing asphaltene deposits in the oil and gas industry and the use of asphaltene deposits as a secondary raw material, this area needs further development. The widespread implementation of equipment and, from the standpoint of rational use of natural resources, the use of deposits as a secondary raw material will increase cost efficiency and equipment efficiency, and reduce environmental impact.

Keywords: mathematical modeling, ultrasound, ultrasonic installation, oil, asphalt-resinous paraffin deposits, pumping equipment, well, reservoir, oil pipeline, secondary raw materials

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

Чаще всего особенностью данных месторождений является высокая вязкость нефтей и проблема образования асфальтосмолопарафиновых отложений (АСПО) в насосно-компрессорных трубах, резервуарах для хранения нефти и магистральных трубопроводах.

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

Рис.1. Пример асфальтосмолопарафиновых отложений (АСПО) Fig. 1. An example of asphalt-resin-paraffin deposits (ARPD)

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

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

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

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

Известно, что в настоящее время, как в России, так и за рубежом для контроля АСПО в нефтегазовой отрасли существуют и активно применяют несколько наиболее известных способов, подробно описанных в статье [2], которые позволяют с определенной эффективностью предупредить или устранить данную проблему (рис 2).

Рис 2. Депарафинизация АСПО Fig. 2. Dewaxing of ARPD

Методы исследования. Самым распространённым методом удаления АСПО в силу своей дешевизны и простоты реализации является «горячая» обработка. Однако из-за высокой степени плотности и высокой температуры плавления отложений асфальтенового типа эти методы теряют эффективность с ростом доли асфальтенов в составе отложений.

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

Применение того или иного греющего кабеля определяется способом добычи нефти. Для скважин, оснащенных ШГН, нагреть скважинную жидкость можно с помощью нагревательного кабеля, проложенного только снаружи НКТ (рис. 3), так как внутри НКТ находится штанга. Для скважин, оснащенных ЭЦН, а также фонтанных и газлифтных нагреть скважинную жидкость можно с помощью самонесущего нагревательного кабеля, опускаемого в НКТ (рис. 4б) через лубрикатор.

a) 6)

Рис. 3. Расположение нагревательных кабелей в скважине Fig. 3. Location of heating cables in the well а) скважина с ШГН, б) скважины с ЭЦН, фонтанные и газлифтные: 1 - насосно-компрессорные

труба; 2 - штанга насоса; 3 - кабель; 4 - обсадная колонна a) well with sucker-rod pumping pump, b) wells with ESP, flowing and gas-lift: 1 - tubing; 2 - pump rod;

3 - cable; 4 - casing

в) r)

Рис. 4. Поперечное сечение кабелей Fig. 4. Cross-section of cables а) и б) плоские нагревательные кабели; в) комбинированный нагревательно-силовой кабель; г) самонесущий нагревательный кабель: 1 - нагревательная жила, 2 - двухслойная изоляция из полимерного материала с обмоткой, 3 - подушка под броню, 4 - броня из стальной оцинкованной профилированной ленты; 5 - изоляция нагревательная жила из полимерного материала; 6 - три силовые жилы; 7 - изоляция силовых жил; 8 - металлический теплоотводящий экран; 9 - оболочка, 10, 11 - двухслойная броня из стальных оцинкованных проволок a) and b) flat heating cables; c) combined heating and power cable; d) self-supporting heating cable: 1 - heating conductor, 2 - two-layer insulation made of polymer material with a winding, 3 - cushion under the armor, 4 - armor made of galvanized steel profiled tape; 5 - insulation heating core made of polymeric material; 6 - three power cores; 7 - insulation of power conductors; 8 - metal heat sink; 9 - sheath,

10, 11 - double-layer armor made of galvanized steel wires

На рис. 5 дан сравнительный анализ двух вариантов. Видно, что при мощности 100 Вт/м температура нефти в НКТ составит 47°С, в то время как при нагреве самонесущим кабелем, расположенным в НКТ, - 43°С при мощности 24 Вт/м. Следовательно, нагрев кабелем, расположенным внутри НКТ, требует в несколько раз меньшей мощности, чем нагрев кабелем, расположенным снаружи НКТ.

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

" SB

а б

Рис. 5. Распределение температуры в поперечном сечении скважины Fig. 5. Temperature distribution in the cross-section of the well а) при мощности 100Вт/м для кабеля, расположенного с внешней стороны НКТ; б) при мощности 24Вт/м

для кабеля, расположенного внутри НКТ a) at a power of 100 W / m for a cable located on the outside of the tubing; b) at a power of 24W / m for a cable

located inside the tubing

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

Химические методы широко используются для удаления АСПО, но довольно дорогостоящие и оказывают большое экологическое воздействие на окружающую среду. Помимо этого агрессивные реагенты сокращают срок службы металлических частей скважин, причиняя прямые убытки добывающим компаниям. Химические методы позволяют увеличить нефтеотдачу на 25...30 % и имеют экономическую эффективность - 648 р./т [0].

Тепловые - весьма эффективные методы. Они основаны на способности парафина плавиться при температурах выше 50°С и стекать с нагретой поверхности. Осаждение парафина сильно зависит от температуры, поэтому термические методы могут быть очень эффективны как для предотвращения, так и для устранения АСПО [6]. Недостатки термических методов заключаются в высокой энергоемкости, электро- и пожароопасности, ненадежности и низкой эффективности.

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

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

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

ударных волн разрушаются агрегаты твердой фазы. УЗ диспергирование проходит успешно при условии, что напряжение, действующее на частицу при захлопывании кавитационной полости, выше реальной прочности измельчаемых материалов. Таким образом, активно воздействуя на кинетику химических реакций и обеспечивая стимуляцию тепло- и массообменных процессов, УЗ способствует увеличению производительности различных технических систем и снижению их энергоемкости и, за небольшим исключением позволяет достичь высокой дисперсности практически любых твёрдых материалов [7]. Интересным вопросом является возможность применения ультразвука совместно с другими методами.

Модельный расчет применения ультразвукового оборудования для удаления АСПО. Предлагается, что в резервуар с нефтью и отложениями помещается ультразвуковой излучатель, соответственно он занимает положение на границе фаз - «нефть-отложение». Излучатель можно помещать в резервуар, который выведен из эксплуатации, так и в резервуар во время эксплуатации. Далее излучатель включается, и вокруг него нефть с отложениями нагревается, происходит расплавление АСПО, и радиус границы фаз увеличивается.

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

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

C^ + ^f^^ (1)

где T - температура нефти/АСПО в исследуемом объеме; r - радиус (расстояние до начала координат); cv - объёмная теплоёмкость нефти/АСПО; X - коэффициент теплопроводности нефти/АСПО.

В рамках данной модели предполагается зависимость температуры только от радиуса и времени T(r, t) (сферическая симметрия). Параметр a - коэффициент затухания ультразвуковых волн в нефти/АСПО; хо - коэффициент, отражающий свойства среды и параметры ультразвуковых волн. Для расчёта запишем уравнение в виде

«V. (l+X^e-J1^ = (U^2%)) (2)

Будем использовать классические разностные схемы для численного решения параболического уравнения общего вида

^'«■^(^■«Э+Я'''«' (3)

при проведении основного расчёта

c(r,t) = cva (l + Xol ■e-ar)~1r2, k(r,t)=Kr2, f = 0. (4)

Введём на отрезке [r0,rmax] равномерную сетку с шагом h = (rmax — г0)/п и узлами xi = ro + ih (i = 0,..., n).

Аналогично для j-го шага по времени вводим tj = jx, j = 1,..., m, шаг по времени равен т = Т/т.

Разностная схема будет по очереди вычислять значение температуры в узлах:

ui = T(xi,tj) (5)

Основная часть консервативной разностной схемы второго порядка [17] может быть записана в следующем виде:

Pi.j

J+1 J u< -u

L = aAJ'uJ'+1 + (1- а)АЫ) + pUj,

(6)

где оператор Л аппроксимирует приток тепла в узел } из соседних узлов

tiuk

(7)

и имеет вид

коэффициенты р^у, а^^, соответствуют средним значениям с, к, / на отрезках и могут быть вычислены различными способами [17].

В данной работе используется простейший стандартный вариант:

_ c{xi-^tj+2)+c{xi+!2,tj+2) Pi, j = 2 '

<i'j = 2

aiJ =

k(xi-h,tj+^)+c(xi,tj+^)

(8) (9) (10)

Уравнения (6) дополняются начальными и граничными условиями для конкретной задачи. При переходе на каждый следующий шаг по времени решается система линейных алгебраических уравнений относительно и0 — иП при известных и0 ...иП.

Вид такой системы определяется из (6). Основная разностная схема с весом о (обычно о = 0.5) имеет вид:

I ^оТ ( , Л ¡+1 от ¡+1 от ¡+1 Ь (Ри +Т2( аи + ))и1 -Т2 а1+1и{+1 — Т2 = ьи,

h2K~ i'] , - M'jjj^ h2

(l-a)Tf^ l \\ „ J (1-(T)T

hi = {pij - (aij + ai+u)) ui

h2

h2

( ai+iuj+i + aiu]i-1) + T<pij. (11)

Так как данное уравнение содержит только неизвестные с номерами / — 1,1,1 + 1, матрица системы является трёхдиагональной и для решения применяется стандартный метод прогонки [17]. Граничные условия первого рода (для случая известной температуры с возможной зависимостью от времени) принимают вид

uo = TL(tJ+i)' uh = TR(tj+1). Граничное условие второго рода

- =QL(t)

x=Ro

(12) (13)

может быть с точностью до h2 аппроксимировано выражением [17]:

g 1 h0 +(1 - a) = -Qlj --Wj+-Po,

h

I+1 I u0 - u

2

2

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

°'QL,j = QL(tj+l)' (14)

что приводит к первому уравнению линейной системы вида

■ и0з+1 — I аии1+1 = Яи + ^ а1Аи1 — и0з) + 2к]

Аналогичное условие на правой границе интервала (—к—)

У дх) Х=К

h

2Т'

= QR(t),

(15)

(16)

может быть записано в виде

--an,juJn+\ + (han,j + ~hpn,j) uii+1 = QR,j--~a1,j(ui - uCi) + hfnj + ~xp0,nuL (17)

Условие фиксированного потока тепла через сферу

Q

du

r = Ro:TU

XF

x=R0

с учетом к = Ar2 F = 4nR0

(18)

(19)

(20)

соответствует подстановке

2

т

Qw = (21>

Одним из часто применяемых методов при решении задач с фазовыми переходами является метод с переменным шагом по времени. В каждый момент tj положение фронта плавления считается совпадающим с одним из узлов сетки: s(tj) = xk и шаг по времени выбирается таким образом, чтобы в момент tj+1 фронт оказался в соседнем узле (с некоторой погрешностью аппроксимации): s(tj) = xk+1 (при тепловом потоке направо вдоль оси). Положение фронта задаётся функцией s (t), уравнение теплопроводности (1) решается в двух областях (до и после фронта плавления):

Rmin <r < s(t): X = Ян, cv = cVh, (22)

S(t) < Т < Rmax: Я Яа, ^v ^va• Один из путей определения желаемого шага по времени - использование разностной аппроксимации в условия Стефана

( Ан дн дк) =-Lpa —, (23)

\ н дг a дг J r=„ft) Иа dt' v '

дг ' а дг,

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

Для вывода расчётной формулы интегрируем уравнение (3) сначала по переменной г (предполагаем зависимости с, к только от г) по двум отрезкам [И0,5(0] и [з(^, Rmаx]:

r=s(t)-0 ' дг r=R0 JRo

r=s(t)+0

д£Тс(гтг + ^(сТ)\ =к(гл)д- -к(гл)д- „ +

д: 3(:) м *г=5(0+0 дг г=Ктах дг г=5(0+0 3(:)

(24)

Просуммируем эти равенства друг с другом с учётом непрерывности профиля температур и известного разрыва теплового потока (условия Стефана, к = Лг2)

дт

дт

r=s(t)-0 дг

d

k(r,t_ „-k(r,t)- =-Lpar(25)

r=s(t)+0

d

В итоге получим соотношение

£ ¡RR™ c(r)Tdr =Ql-Qr- Lpas21 + ¡^ fdr. (26)

Интегрируя его по отрезку [0, tj+1\, запишем в случае известных постоянных потоков (граничных условий второго рода)

¡RR™ c(r)T(r, t)dr - ¡RR™ c(r)T(r, 0)dr = = (Ql - Q«)tj+i - Lp,,s2(s(t) - s(0)} +

+'j0'"Cfdr- ' (27)

откуда легко выразить перейти к приближённому соотношению для вычисления tJ+1 .

При движении фронта слева s(tj) = Xj и отсутствии источников тепла (f = 0) получим

tj+i = (Ql - Qr)-1 (h Ih Pi,ju|+1 - h $*=i piju? + paL(xj+i)2(xj+i - s0)) (28) и шаг по времени

Ti+i = ti+i-tj. (29)

Если тепловой поток на границах не постоянный, например, вычисляется приближенно по условию фиксированной температуры, то вместо полного интегрирования Ql - QR по времени можно ограничиться отрезком интегрирования [tj, tj+1\ и аналогично записать

Ъ+i = (Ql - Qr)-1 (h HiPijU^1 - h ZUpiX + PaL(xj+i)2h). (30)

Процесс вычислений для перемещения фронта плавления на один шаг сетки заключается в составлении и решении системы описанной выше системы уравнений с трёхдиагональной матрицей с дополнительной фиксацией температуры Tc в узле Xj+1 для некоторого шага ,

(1)

вычисления нового шага по времени т-+1 по формуле (8) или (9) (по полученным ранее значениям температур в узлах сетки u), и повторении процесса до достижения установившегося шага по времени:

ЛI) _ (1-1) Lj+1 Lj+1

< £Tj+1 )

(31)

Обсуждение результатов. Для численного решения задачи Стефана разработана программа на языке Python, основанная на использовании классической консервативной разностной схемы второго порядка из [17]. Проведём численный расчёт для следующих параметров: Q=2000 Вт; г0=0,2 м; L = 170 кДж/кг; To = 20 °С; ТКр = 30 °С; Хн=Ха = 1,4 Вт/м°С; cVH = cm = 2 кДж/м3°С. (рис. 6)

Рис. 6. Программа на языке Python Fig. 6. Python program

Пример профиля температуры приведен на рис. 7.

Рис. 7. Пример профиля температуры Fig. 7. Example of a temperature profile

Завышенная температура объясняется достаточно произвольно заданным значением теплового потока. Перемещение фронта плавления изображено на рис. 8:

Рис. 8. Перемещение фронта плавления Fig. 8. Moving the melting front

Изменим граничные условия, зафиксировав про г = R0 температуру T н = 45 °С (тепловой поток становится переменным, для согласования шага по времени используется формула (9)).

Графики перемещения фронта плавления приведены на Рис. 8. Пунктиром изображено более быстрое перемещение изотермы T = Ткр = 30°С без учёта теплоты плавления вещества (в предположении L = 0)

Рис. 9. Перемещение фронта плавления для граничного условия c температурой TL =45°С Fig. 9. Moving the melting front for the boundary condition with temperature TL = 45 °C

69

Пример профиля температур для данного граничного условия приведён на рис. 10.

Рис. 10. Пример профиля температуры Т(г) для граничного условия Ti =45°С Fig. 10. Example of a temperature profile T (r) for the boundary condition T_L = 45 °C

Вывод. Разработана математическая модель применения ультразвукового оборудования для удаления АСПО. Приведены результаты математического моделирования распространения тепла.

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

Библиографический список:

1. Галиакбарова, Э.В, Бахтизин Р.Н, Надыршин, А.С., Галиакбаров, В.Ф. Безопасное энергетически эффективное исключение осадконакопления при хранении нефти в резервуарных емкостях // Электронный научный журнал «Нефтегазовое дело». 2015, №4. С. 142-148.

2. Иванова, Л.В, Буров Е.А., Кошелев В.Н. Асфальтосмолопарафиновые отложения в процессах добычи, транспорта и хранения // Электронный научный журнал «Нефтегазовое дело» 2011, №1. С. 268-284.

3. Hamida T. Effects of Ultrasonic Waves on Immiscible and Miscible Displacement in Porous Media/T. Hamida, T. Babadagli//2005 SPE Technical Conference and Exhibition. Dallas, Texas, 2005.

4. Goswant, Т. К. How to Get Rid of Paraffin Deposition in Pipe Line / Т. К. Goswant // Chem. Age India. - 1969. Vol. 20. No 4. рр. 15-16.

5. Malykh, N. On Sonocapillary Effect / N. Malykh, V. Petrov // XIII Session of the Russian Acoustical Society. Moscow, 2003.

6. Злобин А.А. К вопросу о механизме действия ингибиторов для защиты от асфальтосмолопарафиновых отложений (АСПО)/А.А. Злобин, И. Р. Юшков // Вестник ПНИПУ. Геология. Нефтегазовое и горное дело. 2011. С. 78-83.

7. Бергман Л. Ультразвук и его применение в науке и технике / Пер. с немец. М.: Иностр. лит., 1957. 726 с.

8. Abramov V.O., Ultrasonic technology for enhanced oil recovery from failing oil wells and the equipment for itsimplemention / Abramov V.O., Mullakaev M.S., Abramova A.V., Esipov I.B., Mason T.J. // Ultrasonics Sono-chemistry. 2013. Vol. 20. No. 5. рр. 1289-1295.

9. Abramov O.V, Extraction of bitumen, crude oil and its products from tar sand and contaminated sandy soil under effect of ultrasound / Abramov O.V., Abramov V.O., Myasnikov S.K., Mullakaev M.S. // Ultrasonics Sonochem-istry. 2009. Vol. 16. No. 3.pp. 408-416.

10. Хофштаттер Г., Павлов М.В., Мастобаев Б.Н. Применение ультразвука для разрушения асфальтомолопа-рафиновых отложений в трубопроводе при транспортировке нефти // Транспорт и хранение нефтепродуктов и углеводородного сырья, 2014. № 3. С. 6-9.

11. Павлов М.В., Применение ультразвука для очистки от асфальтосмолистых и парафиновых отложений на объектах транспорта и хранения нефти: дис. канд. техн. наук: 25.00.19 / Павлов Михаил Валентинович, М., 2019. 133 с.

12. Асылбаев Д.Ф., Муллакаев М.С., Абрамов В.О., Прачкин В.Г. Исследование комбинированного воздействия ультразвука и химических реагентов на реологические свойства нефти Боровского месторождения // Оборудование и технологии для нефтегазового комплекса. 2013. № 5. С. 34-36.

13. Beresnev I.A., Johnson P. A., Elastic-wave stimulation of oil production: A review of methods and results, Geophysics, vol. 59, no. 6, pp. 1000-1017, 1994.

14. F. Van der Bas and e. al., "Radial Near Wellbore Stimulation by Acoustic Waves," in SPE International Symposium and Exhibition on Formation Damage Control, Lafayette, 2004

15. Mason T., Advances in Sonochemistry // JAI Press. 1 - 1989, 323 p.; 2 - 1991. - 345 p.; 3 - 1993. - Р. 328; 4 -1996. - 351 p.

16. Price G. Current Trends in Sonochemistry. L.: Royal Soc. of Chemistry, 1990. 121 p.

17. Самарский А. А. Введение в теорию разностных схем [Текст]. М.: Наука, 1971. 552 с.

18. Zerroukat, Mohamed & Chatwin, Chris. (1997). Computational Moving Boundary Problem. Journal of Fluid Mechanics. 343. pp.407-407.

19. Gupta R. S., Kumar D. Variable time step methods for one-dimensional Stefan problem with mixed boundary condition //International Journal of Heat and Mass Transfer. 1981. Vol. 24. No. 2. pp. 251-259.

20. Кот В.А. Высокоточные полиномиальные решения классической задачи Огефана // Технические науки, 2017. №6. С. 112-122.

21. Муллакаев, М.С. Ультразвуковая интенсификация добычи и переработки нефти // М.: ОАО "ВНИИО-ЭНГ", 2014. 168. с.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

22. Чесноков А.В. Что ждет сервисные компании в нефтегазовой отрасли России?//Нефть и Газ Евразия. 2008. № 10, с. 10-15.

References:

1. Galiakbarova, E.V, Bakhtizin R.N, Nadyrshin, A.S., Galiakbarov, V.F. Bezopasnoye energeticheski effektivnoye isklyucheniye osadkonakopleniya pri khranenii nefti v rezervuarnykh yemkostyakh // Elektronnyy nauchnyy zhurnal «Neftegazovoye delo» 2015, №4. S. 142-148. [Galiakbarova, E. V., Bakhtizin, R. N., Nadyrshin, A. S., Galiakbarov, V. F. Safe energy-efficient exclusion of sedimentation during oil storage in reservoir tanks / / Electronic scientific journal "Oil and Gas Business" 2015, No. 4. pp. 142-148. (In Russ)]

2. Ivanova, L.V, Burov Ye.A., Koshelev V.N. Asfal'tosmoloparafinovyye otlozhe niya v protsessakh dobychi, transporta i khraneniya // Elektronnyy nauchnyy zhurnal «Neftegazovoye delo» 2011, №1. S. 268-284. [Ivanova, L. V., Burov, E. A., Koshelev, V. N. Asphalto-resin-paraffin deposits in mining, transport and storage processes // Electronic scientific journal "Oil and Gas Business" 2011, No. 1. pp. 268-284. (In Russ)]

3. Hamida, T. Effects of Ultrasonic Waves on Immiscible and Miscible Displacement in Porous Media / T. Hamida, T. Babadagli // 2005 SPE Technical Conference and Exhibition. Dallas, Texas, 2005.

4. Goswant, Т. К. How to Get Rid of Paraffin Deposition in Pipe Line / Т. К. Goswant // Chem. Age India. - 1969. Vol. 20. No 4. рр. 15-16.

5. Malykh, N. On Sonocapillary Effect / N. Malykh, V. Petrov // XIII Session of the Russian Acoustical Society. Moscow, 2003.

6. Zlobin, A. A. K voprosu o mekhanizme deystviya ingibitorov dlya zashchity ot asfal'tosmoloparafinovykh otlozheniy (ASPO) / A. A. Zlobin, I. R. Yushkov // Vestnik PNIPU. Geologiya. Neftegazovoye i gornoye delo. 2011. S. 78-83. [Zlobin, A. A. On the question of the mechanism of action of inhibitors for protection from asphalt-resin-paraffin deposits (ASPO) / A. A. Zlobin, I. R. Yushkov // Vestnik PNRPU. Geology. Oil and gas and mining. 2011. pp. 78-83. (In Russ)]

7. Bergman L. Ultrasound and its application in science and technology / Per. s. nemets. M.: Inostr. lit., 1957. 726 p.

8. Abramov V.O., Ultrasonic technology for enhanced oil recovery from failing oil wells and the equipment for itsimplemention / Abramov V.O., Mullakaev M.S., Abramova A.V., Esipov I.B., Mason T.J. // Ultrasonics Sonochemistry. 2013. Vol. 20. No. 5. pp. 1289-1295.

9. Abramov O.V, Extraction of bitumen, crude oil and its products from tar sand and contaminated sandy soil under effect of ultrasound / Abramov O.V., Abramov V.O., Myasnikov S.K., Mullakaev M.S. // Ultrasonics Sonochemistry. 2009. Vol. 16. No. 3.C. 408-416.

10. Khofshtatter G., Pavlov M.V., Mastobayev B.N. Primeneniye ul'trazvuka dlya razrusheniya asfal'tomoloparafinovykh otlozheniy v truboprovode pri transportirovke nefti // Transport i khraneniye nefteproduktov i uglevodorodnogo syr'ya, 2014. № 3. S. 6-9. [Hofstatter G., Pavlov M. V., Mastobaev B. N. Application of ultrasound for destruction of asphalt-paraffin deposits in the pipeline during oil transporta-tion//Transport and storage of petroleum products and hydrocarbon raw materials, 2014. No. 3. pp. 6-9. (In Russ)]

11. Pavlov M.V., Primeneniye ul'trazvuka dlya ochistki ot asfal'tosmolistykh i parafinovykh otlozheniy na ob"yektakh transporta i khraneniya nefti: dis. kand. tekhn. nauk: 25.00.19 / Pavlov Mikhail Valentinovich, M., 2019. 133 s. [Pavlov M. V., Application of ultrasound for cleaning from asphalt resinous and paraffin deposits at oil transport and storage facilities: dis. cand. tech. nauk: 25.00.19 / Pavlov Mikhail Valentinovich, M., 2019. 133 p. (In Russ)]

12. Asylbayev D.F., Mullakayev M.S., Abramov V.O., Prachkin V.G. Issledovaniye kombinirovannogo vozdeystviya ul'trazvuka i khimicheskikh reagentov na reologicheskiye svoystva nefti Borovskogo mestorozhdeniya // Oborudovaniye i tekhnologii dlya neftegazovogo kompleksa. 2013. № 5. S. 34-36. [Asylbaev D. F., Mullakaev M. S., Abramov V. O., Prachkin V. G. Investigation of the combined effect of ultrasound and chemical reagents

on the rheological properties of Borovskoye oil // Equipment and technologies for the oil and gas industry. 2013. No. 5(In Russ)]

13. Beresnev I.A., Johnson P. A., Elastic-wave stimulation of oil production: A review of methods and results, Geophysics, Vol. 59, No. 6, pp. 1000-1017, 1994.

14. F. Van der Bas and e. al., "Radial Near Wellbore Stimulation by Acoustic Waves," in SPE International Symposium and Exhibition on Formation Damage Control, Lafayette, 2004

15. Mason T., Advances in Sonochemistry // JAI Press. - 1 - 1989, 323 p.; 2 - 1991. - 345 p.; 3 - 1993. - Р. 328; 4 -1996. - 351 p.

16. Price G. Current Trends in Sonochemistry. L.: Royal Soc. of Chemistry, 1990. 121 p.

17. Samarskiy A. A. Vvedeniye v teoriyu raznostnykh skhem [Tekst]. - M.: Nauka, 1971. - 552 s. [Samarsky A. A. Introduction to the theory of difference schemes [Text]. M.: Nauka, 1971. 552 p. (In Russ)]

18. Zerroukat, Mohamed & Chatwin, Chris. (1997). Computational Moving Boundary Problem. Journal of Fluid Mechanics. 343. pp. 407-407.

19. Gupta R. S., Kumar D. Variable time step methods for one-dimensional Stefan problem with mixed boundary condition //International Journal of Heat and Mass Transfer. 1981. Т. 24. No. 2. pp. 251-259.

20. Kot V.A. Vysokotochnyye polinomial'nyye resheniya klassicheskoy zadachi Ctefana // Tekhnicheskiye nauki, 2017. №6. S. 112-122.

21. Mullakayev, M.S. Ul'trazvukovaya intensifikatsiya dobychi i pererabotki nefti // M.: OAO "VNIIOENG", 2014. 168. s. [Mullakaev, M. S. Ultrasonic intensification of oil production and processing / / Moscow: JSC "VNIIOENG", 2014. 168p. (In Russ)]

22. Chesnokov A.V. Chto zhdet servisnyye kompanii v neftegazovoy otrasli Rossii? // Neft' i Gaz Yevraziya. - 2008. № 10. S. 10-15 [Chesnokov A.V. What is waiting for service companies in the oil and gas industry of Russia? // Oil and Gas Eurasia. 2008. № 10. pp. 10-15. (In Russ)]

Сведения об авторах:

Хурамшина Регина Азатовна, аспирант, [email protected]

Соколова Виктория Владимировна, кандидат филологических наук, доцент, [email protected] Information about the authors:

Regina A. Khuramshina, Postgraduate, [email protected]

Viktoriya V. Sokolova, Cand. Sci. (Philoligical), Assoc. Prof., [email protected]

Конфликт интересов. Conflict of interest.

Авторы заявляют об отсутствии конфликта интересов. The authors declare no conflict of interest.

Поступила в редакцию 14.05.2021. Received 14.05.2021.

Принята в печать 01.06.2021. Accepted for publication 01.06.2021.

i Надоели баннеры? Вы всегда можете отключить рекламу.