УДК 622.279, 620.193
НОВЫЙ МЕТОДИЧЕСКИЙ ПОДХОД К ПРОГНОЗИРОВАНИЮ ПОДВЕРЖЕННОСТИ ОБОРУДОВАНИЯ ВЫСОКОТЕМПЕРАТУРНЫХ ГАЗОКОНДЕНСАТНЫХ СКВАЖИН
УГЛЕКИСЛОТНОЙ КОРРОЗИИ
Пономарев Александр Иосифович1,
Иванов Николай Валерьевич2,
Юсупов Александр Дамирович1,
1 Уфимский государственный нефтяной технический университет, Россия, 450062, г. Уфа, ул. Космонавтов, 1.
2 Российский государственный университет нефти и газа (НИУ) им. И.М. Губкина, Россия, 119991, г. Москва, Ленинский пр., 65.
Актуальность. Перспективы увеличения добычи газового конденсата - ценного сырья для нефтепереработки и нефтехимии - в значительной мере связаны с дальнейшим освоением ресурсов газоконденсатных залежей в ачимовских отложениях севера Западной Сибири. Большие глубины, аномально высокое пластовое давление и высокая температура низкопродуктивных пластов наряду с присутствием в составе пластового газоконденсатного флюида диоксида углерода выдвигают, в первую очередь, повышенные требования к качеству заканчивания скважин и операций гидроразрыва пласта, надежности конструкций и материалов подземного оборудования скважин. Этими факторами обусловлены высокая стоимость реализации проектов и эксплуатационные риски добычи газа и конденсата, в связи с чем обоснование безопасных условий эксплуатации скважин в ачимовских отложениях является актуальной научно-технической задачей.
Цель: определить предельные значения термобарических параметров технологических режимов эксплуатации ачимовских газоконденсатных скважин с содержанием в добываемой продукции диоксида углерода, обеспечивающих физико-химические условия отсутствия электрохимической коррозии забойного оборудования - хвостовиков на длительный (20-летний) период. Объект: пластовый газоконденсатный флюид и забойное оборудование трех ачимовских газоконденсатных скважин с хвостовиками из углеродистой стали, нестойкой к углекислотной коррозии.
Метод: моделирование фазового поведения добываемого пластового флюида в скважинных условиях в среде программного обеспечения ГазКондНефть.
Результаты. Термодинамическими расчетами фазового поведения пластовой газоконденсатной смеси с учетом ее влаго-содержания показано, что снижение пластового давления в зоне расположения трех рассматриваемых скважин при разработке участка ачимовской залежи на протяжении 20 лет при проектных технологических режимах их эксплуатации со временем приводит к образованию двухфазной смеси «газ - нестабильный конденсат» в забойных термобарических условиях скважин. При этом водная жидкая фаза вследствие высокой температуры потока на забое скважин не образуется в течение всего расчетного периода. Гидродинамические расчеты параметров восходящего потока газожидкостной смеси показали, что высокие скорости потока скважинной продукции обеспечивают условия полного и непрерывного выноса нестабильного конденсата потоком газа с забоя на поверхность по каждой из рассматриваемых скважин в течение всего 20-летнего периода, тем самым предотвращаются физико-химические условия образования на поверхности хвостовиков электролита и протекания углекислотной коррозии.
Ключевые слова:
Газоконденсатная смесь, диоксид углерода, хвостовик скважины, параметры технологического режима эксплуатации скважины, симулятор, фазовые диаграммы, линии точек росы воды и углеводородов, режим течения многофазного потока, углекислотная коррозия, антикоррозионный режим эксплуатации скважины.
Введение
Вопросы углекислотной коррозии как подземного, так и наземного оборудования и промысловых трубопроводов были актуальны в 60-е гг. прошлого столетия при разработке газоконденсатных месторождений Краснодарского и Ставропольского краев [1]. Тогда отмечалась неравномерность коррозии насосно-компрессорных труб (НКТ) на некоторых скважинах. Чаще всего фиксировалась коррозия большей интенсивности в приустьевой зоне, чем в забойной, а на некоторых скважинах коррозия башмака и нижних НКТ не отмечалась вовсе. Наблюдаемые коррозионные поражения НКТ рядом исследователей [2, 3] объяснялись тем, что в стволе сква-
жин при снижении температуры и давления имеет место выделение из добываемого газа свободной водной жидкой фазы, присутствие которой инициирует процесс коррозии. При этом известно, что конденсационные воды сами по себе обладают небольшой коррозионной агрессивностью по отношению к металлу, но разрушительное действие этих вод многократно увеличивается при насыщении их углекислым газом. Ввиду отсутствия в тот период времени адекватных средств моделирования фазового поведения газоконденсатных смесей с учетом их влагосодержания, а также гидродинамических методов расчета многофазных вертикальных и горизонтальных потоков, исследования по определению гра-
DOI 10.18799/24131830/2021/06/3235
49
ничных термодинамических условии существования антикоррозионного режима эксплуатации скважинного и поверхностного оборудования не проводились.
Ачимовские отложения в настоящее время являются основным ресурсом прироста запасов и добычи углеводородного сырья (УВС) и, в первую очередь, газового конденсата в ЗападноИ Сибири [4]. Один из крупнейших - Уренгойский лицензионный участок - в свою очередь поделен на шесть эксплуатационных участков [5, 6].
Подземная конструкция большинства скважин на одном из эксплуатационных участков Уренгойского лицензионного участка (эксплуатационная колонна, хвостовик, насосно-компрессорные трубы с комплексом подземного оборудования) имеет коррозионно-стойкое исполнение. Такое проектное техническое решение обусловлено потенциальной коррозионной агрессивностью пластового флюида из-за содержания в его составе диоксида углерода. Также известно, что на Уренгойском лицензионном участке существует проблема углекислотной коррозии наземного оборудования. Так, элементы устьевой обвязки скважин и отдельные участки газосборных трубопроводов, выполненных в некоррозионно-стойком исполнении, за десятилетний период эксплуатации подверглись воздействию углекислотной коррозии [7-9].
Вместе с тем определение предельных значений термодинамических параметров технологических режимов эксплуатации ачимовских газоконденсатных скважин, обеспечивающих физико-химические условия отсутствия электрохимической коррозии забойного оборудования - хвостовиков, представляет большой практический интерес. При существовании таких режимов в течение длительного периода эксплуатации (несколько десятилетий) можно рассматривать вопрос о внесении изменений в проектную конструкцию новых скважин с применением хвостовиков в некоррозионно-стойком исполнении, что позволит частично снизить стоимость скважин.
Методика проведения исследования
Типовая конструкция скважин одного из эксплуатационных участков Уренгойского лицензионного участка [10] приведена на рис. 1 и включает:
• кондуктор (Б=324 мм), перекрывающий много-летнемерзлые породы и неустойчивые четвертичные отложения;
• техническую (промежуточную) колонну (0=245 мм), перекрывающую сеноманский горизонт с аномально низкими пластовыми давлениями;
• эксплуатационную колонну (Б=178 мм), перекрывающую пласты с нормальными и аномально низкими пластовыми давлениями с установкой башмака колонны на 20 м выше продуктивного пласта ачимовских отложений;
• хвостовик (0=114 мм), расположенный в интервале продуктивных пластов;
• лифтовую колонну (Б=89 мм), на которой установлен комплекс подземного оборудования (пакер, клапан-отсекатель, датчики давления и температуры и др.).
Рис. 1. Типовая конструкция ачимовских скважин Fig. 1. Typical well design for the achimov deposits
В 2020 г. в рамках программы исследовательских работ в эксплуатацию были введены три скважины, конструкция которых предусматривает хвостовики, выполненные в нестойком к коррозии исполнении -из углеродистой стали группы прочности Р110. При поверхностном анализе потенциальных рисков, связанных с процессами коррозии хвостовиков, напрашивался вывод о протекании более интенсивной уг-лекислотной коррозии металла хвостовиков в сравнении с коррозией металла поверхностного оборудования ввиду более высоких температур и парциальных давлений углекислого газа, наблюдаемых на забое.
Однако электрохимия углекислотной коррозии такова, что для протекания процесса коррозии металла необходимым условием является наличие на поверхности металла водной жидкой фазы. Считается, что пленка воды на поверхности металла толщиной 20-30 молекул является достаточной для протекания электрохимической коррозии [11], к которой относится углекислотная коррозия, наблюдаемая на наземных объектах добычи УВС ачимовских отложений. Поэтому присутствие воды на забое скважины только в паровой фазе исключает коррозию хвостовиков, несмотря на более высокие температуры и парциальное давление углекислого газа по сравнению с условиями устья скважины, где наблюдаются интенсивные коррозионные процессы.
Для достижения поставленной цели определены следующие задачи исследования: • для трех опытных скважин определить фазовое состояние воды на забое и устье при составах сква-жинной продукции, изменяющихся в течение разработки участка залежи по мере снижения пластового и забойного давлений на 20-летний период;
• предложить рекомендации по противокоррозионной защите и коррозионному мониторингу скважин, оборудованных хвостовиками из углеродистой стали.
Математическое моделирование фазового равновесия для систем природных углеводородов заключается в использовании различных уравнений состояния многокомпонентных систем, описывающих соотношения между параметрами углеводородной смеси в состоянии парожидкостного равновесия [12]. Уравнение состояния Ван-дер-Ваальса послужило основой для всех существующих уравнений и их модификаций, таких как: Редлиха-Квонга [13], Ли-Кеслера [14], Пенга-Робинсона [15, 16] и других. Выбор определенного уравнения для расчета фазового равновесия обусловлен диапазонами давления и температуры, а также составом природных газоконденсатных смесей, в который входят углеводороды различного строения (парафиновые, нафтеновые, ароматические) и неуглеводородные вещества (сероводород, диоксид углерода, азот и др.). Все современные программные комплексы, предназначенные для моделирования технологии добычи и подготовки углеводородного сырья, включают уравнение состояния Пенга-Робинсона в той или иной модификации. Адекватность уравнения состояния Пенга-Робинсона применительно к расчету процессов добычи и подготовки углеводородного сырья подтверждается работой [17]. Кроме того, в работах [18-22] дается рекомендация к использованию уравнения состояния Пенга-Робинсона для моделирования фазового поведения многокомпонентных углеводородных смесей, в т. ч. при пластовых условиях. По этой причине все выполненные технологические расчеты проводились с использованием уравнения состояния Пенга-Робинсона.
В качестве инструментария исследования требовалось определить программный комплекс технологического моделирования, адекватно описывающий фазовое поведение добываемой из ачимовских скважин
углеводородной смеси с учетом ее влагосодержания. Проведя анализ применения различных комплексов, описанных в работе [23], для расчетов была выбрана компьютерная программа ГазКондНефть.
В технологической модели «забой - устье скважины» производились расчеты фазового состояния пластового флюида, для чего использовались следующие исходные данные по исследуемым скважинам на 20-летний прогнозный период:
• компонентно-фракционные составы добываемой углеводородной смеси по состоянию на 2020 г.;
• проектные значения технологических режимов эксплуатации скважин: депрессии на пласт, деби-ты газоконденсатной смеси и воды.
Методика исследования заключается в последовательном решении следующих задач:
• определение термобарических условий выделения жидкой углеводородной фазы при изменении состава добываемого пластового флюида и параметров технологических режимов скважин в процессе истощения залежи;
• определение фазового состояния воды и возможности ее существования в забойных и устьевых условиях в водной жидкой фазе;
• определение гидродинамических условий полного и непрерывного выноса жидкости с забоя скважины на поверхность.
По результатам решения этих задач проводится анализ фазового состояния добываемого флюида на забое и устье каждой из скважин с точки зрения возможности формирования пленки или слоя водного раствора электролита на внутренней поверхности хвостовика и устьевой обвязки.
Полный состав газоконденсатной смеси для расчета фазового состояния добываемого пластового флюида и проектные параметры технологического режима скважины № 1 на 2020 г. приведены соответственно в табл. 1, 2.
Компоненты и фракции Components and fractions Содержание, % мас. Content, wt. % Фракции Fractions Содержание, % мас. Content, wt. % Фракции Fractions Содержание, % мас. Content, wt. %
Азот/Nitrogen 0,1590 F9(130-140 °С) 1,7182 F26(300-310 °С) 2,39-10-3
Диоксид углерода Carbon dioxide 1,3544 Fw(140-150 °С) 0,6482 F27(310-320 °С) 1,57^0-3
Метан/Methane 42,4308 F„(150-160 °С) 1,0133 F2S(320-330 °С) 1,03^0-3
Этан/Ethane 8,9696 F12(160-170 °С) 0,3451 F29(330-340 °С) б^ЫО-4
Пропан/Propane 8,6921 F13(170-180 °С) 0,3794 F30(340-350 °С) 3,22-Ш-4
Изобуган/Isobutane 3,3817 FM(180-190 °С) 0,2779 F31(350-360 °С) 1,63^10-4
Норм. буган/n-butane 5,1723 F15(190-200 °С) 0,1469 F32(360-370 °С) 8,68^10-5
Изопентан/i-pentane 2,6813 F16(200-210 °С) 0,1100 F33(370-380 °С) 3,8M0-5
Норм. пентан/n-pentane 3,1376 F„(210-220 °С) 0,0899 F34(380-390 °С) 1,91^10"5
F1(45-60 °С) 1,6847 F18(220-230 °С) 0,0443 F35(390-400 °С) 9,62-Ш-6
F2(60-70 °С) 2,7769 F19(230-240 °С) 0,0334 F36(400-410 °С) 5,15^10-6
F3(70-80 °С) 1,3186 F20(240-250 °С) 0,0277 F37(410-420 °С) 2,62^10-6
F4(80-90 °С) 2,2671 F21(250-260 °С) 0,0173 F3S(420-430 °С) 1,46^10-6
F5(90-100 °С) 3,0199 F22(260-270 °С) 0,0132 F39(430-440 °С) 6,52^10-7
F6(100-110 °С) 4,2434 F23(270-280 °С) 0,0095 F40(440-450 °С) 2,28 •ÎO-7
Fy(110-120 °С) 2,6816 F24(280-290 °С) 0,0063 F41(450-460 °С) 1,14^10-7
F8(120-130 °С) 1,1382 F25(290-300 °С) 0,0041 F42(460-470 °С) 4,65^10-s
Таблица 1. Компонентно-фракционный состав углеводородной смеси по скважине № 1 на 2020 г. Table 1. Component-fractional composition of the hydrocarbon fluid for well no. 1 in 2020
Таблица 2. Параметры проектного технологического режима скважины № 1 в 2020 г.
Table 2. Parameters of the design technological regime of well no. 1 in 2020
Параметр Parameter Значение Value
Расход углеводородной смеси, тыс. м3/ч (тыс. м3/сутки) Hydrocarbon fluid consumption, thousand m3/hour (thousand m3/day) 9,59 (230,16)
Расход пластовой воды, м3/сут Formation water consumption, m3/day 1,05
Забойное давление, МПа Bottomhole pressure, MPa 38,72
Температура на забое Bottomhole temperature, °С 108,00
Устьевое давление, МПа Wellhead pressure, MPa 28,35
Температура на устье Wellhead temperature, °С 45,76
Результаты расчетов
Данные табл. 1, 2 были экспортированы в симуля-тор для решения двух первых задач моделирования
фазовых превращений добываемого из скважины № 1 флюида. Изменение компонентно-фракционного состава добываемой газоконденсатной смеси с укрупненными фракциями по скважинам № 1-3 по годам разработки приведено в табл. 3.
По результатам расчетов также определены области двухфазного состояния воды и углеводородов в газоконденсатной смеси, ограниченные соответствующими линиями точек росы (нулевыми изоплерами), для изменяющихся в процессе разработки компонентно-фракционных составов добываемого пластового флюида. При этом следует отметить, что в расчетах на симуляторе использовались полные составы углеводородной смеси, аналогичные приведенным в табл. 1. Смещение границ двухфазной области для воды и углеводородов по состоянию на 2020, 2030 и 2040 гг. по скважинам № 1-3 представлены на диаграммах рис. 2-4, соответственно. На диаграммах также указаны термобарические условия на забое и устье в соответствии с проектными технологическими режимами эксплуатации скважин на вышеупомянутые годы разработки участка ачимовской залежи.
Таблица 3. Изменение компонентно-фракционного состава добываемой газоконденсатной смеси с укрупненными
фракциями по скважинам № 1-3 (2020, 2030, 2040 гг.) Table 3. Change in the component-fractional composition of the produced gas-condensate mixture with enlarged fractions for wells no. 1-3 (in 2020, 2030, 2040)
Компонент Component Содержание компонента/фракции, % масс. Component/fraction content, wt. % (excel кроме воды) 3t water)
Скважина № 1 Well no. 1 Скважина N° 2 Well no. 2 Скважина № 3 Well no. 3
2020 2030 2040 2020 2030 2040 2020 2030 2040
Метан/Methane 42,43 51,29 54,87 47,13 54,77 58,13 45,90 52,18 57,93
Этан/Ethane 8,97 10,28 10,67 8,57 9,75 10,32 8,38 9,35 10,28
Пропан/Propane 8,69 9,11 9,03 6,80 7,55 7,90 6,67 7,31 7,87
Изобуган/Isobutane 3,38 3,26 3,10 2,27 2,48 2,56 2,23 2,42 2,55
Норм. Буган/n-butane 5,17 4,73 4,39 3,18 3,42 3,49 3,13 3,35 3,49
Изопентан/i-pentane 2,68 2,18 1,92 1,42 1,48 1,48 1,40 1,47 1,48
Норм. пентан/n-pentane 3,14 2,46 2,13 1,59 1,65 1,63 1,57 1,65 1,63
Ф1 (фракция 45-100 °С) F1 (fraction 45-100 °C) 11,07 7,21 5,86 5,62 5,27 4,74 5,61 5,49 4,78
Ф2 (фракция 100-200 °С) F2 (fraction 100-200 °C) 12,59 7,45 5,93 12,25 9,41 6,94 12,41 10,79 7,11
Ф3 (фракция 200-300 °С) F3 (fraction 200-300 °С) 0,36 0,22 0,19 5,38 1,96 0,79 5,77 3,15 0,85
Ф4 (фракция 300-470 °С) F4 (fraction 300-470 °С) 0,01 0,00 0,00 4,19 0,42 0,09 5,38 1,08 0,10
Азот/Nitrogen 0,16 0,20 0,21 0,19 0,22 0,23 0,18 0,21 0,23
Углекислый газ/Carbon dioxide 1,35 1,60 1,69 1,40 1,61 1,71 1,37 1,54 1,71
Вода, г/м3 флюида Water, g/m3 of fluid 4,58 6,50 8,45 4,62 6,44 8,55 4,80 6,24 8,07
Для протекания углекислотной коррозии необходимо наличие на поверхностях хвостовика и некоррозионно-стойкого внутрискважинного оборудования, контактирующих с восходящим потоком добываемого флюида, электролита за счет образования слоя или пленки водной жидкой фазы, насыщаемой диоксидом углерода. Поэтому важно определить: а) в какой фазе (паровой или жидкой) находится вода в потоке газоко-нденсатной смеси; б) в какой фазе (газовой или жидкой) находятся углеводороды; в) если имеет место образование водной и углеводородной жидких фаз, то в каком количественном соотношении они находятся, какого типа эмульсии они образуют «вода в масле» или
«масло в воде»; г) какой режим восходящего многофазного потока обеспечивается в хвостовике, лифтовых трубах и на устье скважин при выполнении проектных технологических параметрах их эксплуатации.
Анализ результатов расчетов фазового поведения добываемого пластового флюида (газоконденсатной смеси с учетом влагосодержания) показал, что жидкая водная фаза из-за высокой температуры в пласте и скважинах не образуется ни в одной из исследуемых скважин в течение ближайшего 20-летнего периода их эксплуатации (рис. 2-4, табл. 4). Жидкая углеводородная фаза с незначительным содержанием растворенной в ней воды (0,461.. .0,754 % мольн.) появ-
ляется на забое скважин № 1 и № 3 через 20 лет - к 2040 г., а в скважине № 2 уже через 10 лет - к 2030 г. При таких низких концентрациях в жидкой углеводородной фазе вода находится в молекулярно-растворенном состоянии и не образует с диоксидом углерода электролит. При этом избыточная масса воды находится в паровой фазе в составе газа.
Теперь важно выяснить, каковы гидродинамические условия течения восходящего потока газокон-денсатной смеси в скважине при соблюдении проектного технологического режима эксплуатации: не происходит ли накопление жидкости на забое в хвостовике, в котором объемные скорости потока характе-
ризуются наименьшими значениями в скважине из-за наибольшего давления и большим по сравнению с лифтовой колонной диаметром. При скоплении углеводородной жидкости на забое пусть с небольшим содержанием воды, растворенной в жидкой углеводородной фазе, со временем может образоваться сначала эмульсия типа «вода в масле», частичный объем которой при дальнейшем накоплении может подвергнуться инверсии со сменой на тип «масло в воде». В этом случае неизбежно образование электролита с протеканием углекислой коррозии незащищенного металла скважинного оборудования, в рассматриваемом здесь случае - хвостовика.
100 120 Температура, °С
Рис. 2. Границы двухфазной области для воды и УВ и термобарические условия на забое и устье скважины № 1 по
состоянию на 2020, 2030 и 2040 гг. Fig. 2. Boundaries of the two-phase region for water and hydrocarbons (HC) and thermobaric conditions at the bottomhole and wellhead of well no. 1 as of2020, 2030 and 2040
170
Температура, °C
Рис. 3. Границы двухфазной области для воды и УВ и термобарические условия на забое и устье скважины № 2 по
состоянию на 2020, 2030 и 2040 гг. Fig. 3. Boundaries of the two-phase region for water and HC and thermobaric conditions at the bottomhole and wellhead of well no. 2 as of2020, 2030 and 2040
170 220
Температура, "С
Рис. 4. Границы двухфазной области для воды и УВ и термобарические условия на забое и устье скважины № 3 по
состоянию на 2020, 2030 и 2040 гг. Fig. 4. Boundaries of the two-phase area for water and HC and thermobaric conditions at the bottomhole and wellhead of well № 3 as of2020, 2030 and 2040
Таблица 4. Расчетные физико-химические свойства скважиной продукции (2020, 2030, 2040 гг.) Table 4. Estimated physical and chemical properties of the well production (in 2020, 2030, 2040)
Физические свойства Physical properties Скважина/Well
1 2 3
2020 2030 2040 2020 2030 2040 2020 2030 2040
Наличие жидкой углеводородной фазы на забое скважины Presence of liquid hydrocarbon phase at the bottomhole нет no нет no да yes нет no да yes да yes нет no нет no да yes
Наличие жидкой водной фазы на забое скважины Presence of liquid aqueous phase at the bottomhole нет no нет no нет no нет no нет no нет no нет no нет no нет no
Плотность смеси, кг/м3 Mixture density, kg/m3 340,25 142,65 51,74 325,23 155,30 57,74 345,20 228,87 112,59
Плотность газовой фазы, кг/м3 Gas phase density, kg/m3 340,25 142,65 51,30 325,23 151,44 55,95 345,20 228,87 111,43
Плотность жидкой фазы, кг/м3 Liquid phase density, kg/m3 - - 612,35 - 557,74 630,89 - - 583,14
Молекулярная масса смеси, г/моль Mixture molecular weight, g/mol 27,44 24,12 23,05 26,59 23,50 22,40 27,18 24,44 22,46
Молекулярная масса газовой фазы, г/моль Gas phase molecular weight, g/mol 27,44 24,12 22,89 26,59 22,96 21,82 27,18 24,44 22,26
Молекулярная масса жидкой фазы, г/моль Liquid phase molecular weight, g/mol - - 86,35 - 71,20 92,11 - - 77,57
Мольная доля растворенной воды во всем объеме жидкости Dissolved water mole fraction in the entire volume of liquid - - 4,63 10-3 - 6,96 10-3 4,61 ■ 10-3 - - 7,54-10-3
Для оценки режима течения газожидкостного потока из симулятора фазового поведения газоконден-сатной смеси были выгружены значения вязкостей газа и конденсата, поверхностного натяжения на границе газ-конденсат, коэффициента сверхсжимаемости газа, плотностей газовой и жидкой фаз для термобарических условий на забое и устье каждой из трех скважин для составов и проектных технологических режимов на 2020, 2030 и 2040 гг. Режим течения определялся по результату сравнения значения безразмерной скорости смеси [24]:
w,
= 4(fr + Q /Pi - Рд\°25/Рд\05 nd2 \ аа ' \р,/
да / \ри
со значением скорости нижней границы существования кольцевого режима:
wb = (0,82 - 0,0017^-°'б^ю^+и5^1-«),
wc г
где ql и qg - дебиты конденсата и газа, приведенные к скважинным термобарическим условиям, м3/сутки; Р1 и рг - плотности конденсата и газа в скважинных термобарических условиях, кг/м3; d - диаметр хвостовика или лифтовых труб, м; а - поверхностное
натяжение на границе газ-конденсат, Н/м; — = -
—
соотношение динамической вязкости газа и конденсата в скважинных термобарических условиях;
С
ß =-
С + q,
- расходное объемное газосодержание.
Результаты расчета скоростей w* и wь показали, что в хвостовике, а тем более в лифтовых трубах на забое каждой из скважин, газожидкостной поток в рассматриваемый период разработки (до 2040 г.) движется в дисперсно-кольцевом режиме с превышением скорости wb более чем на порядок (табл. 5).
Таблица 5. Значения расходного газосодержания, безразмерной скорости нижней границы существования кольцевого
режима и безразмерной скорости потока газоконденсатной смеси в хвостовике (2020, 2030, 2040 гг.) Table 5. Values of the flow rate gas content, dimensionless velocity of the lower boundary of the existence of the annular regime and dimensionless flow rate of the gas-condensate liquid in the tailpipe (in 2020, 2030, 2040)
Год разработки Development year Скважина/Well
1 2 3
ß wb w* ß Wb w* ß Wb w*
2020 1,0000 - - 1,0000 - - 1,0000 - -
2030 1,0000 - - 0,9905 1,239 8,848 1,0000 - -
2040 0,9992 0,834 15,635 0,9970 0,890 9,235 0,9975 0,892 29,568
Кроме того, для определения гидродинамических условий полного и непрерывного выноса жидкой фазы с забоя на поверхность размерные скорости потока газа в хвостовике и башмаке лифтовых труб
4 С
W =
nd2ß
сравнивались с критическои скоростью газа опрокидывания жидкой пленки, или скоростью реверса, при котороИ газ с дисперсной фазой жидкости в ядре потока движется вверх, а масса жидкости в кольцевой пленке, условно говоря, зависает, вычисляемой по формуле А.А. Точигина [24]
. 0,25
-0,5
Р- ■
Wrev = 3,3 М-)
Результаты расчетов показали, что для всех трех скважин в течение рассматриваемого периода разработки до 2040 г. выполняется неравенство
ю >
т. е. проектные режимы эксплуатации скважин обеспечивают полный и непрерывный вынос жидкой фазы и из хвостовика, и из башмака лифтовой колонны (табл. 6), т. к. скорость взвешивания частиц жидкости, диспергированной в ядре потока, всегда меньше скорости реверса пленки жидкости при дисперсно-кольцевом режиме течения вертикального потока.
Поскольку нестабильный конденсат не обладает электропроводимостью, пленка на внутренней поверхности хвостовика и лифтовых труб не является электролитом, и при непрерывной эксплуатации скважин в период до 2040 г. углекислотная коррозия хвостовиков исключается. При остановках скважин возможно стекание конденсата на забой, однако и в этом случае из-за незначительного содержания водной фазы в газоконденсатной смеси может образоваться эмульсия типа «вода в масле», что также не приведет к протеканию процесса коррозии.
В рамках данного исследования также был проведен анализ влияния количества воды конденсационного происхождения в составе добываемого флюида
на смещение границы перехода воды в двухфазное состояние. Моделирование фазового поведения газоконденсатной смеси осуществлялось для компонентно-фракционного состава, приведенного в табл. 1, для которого изменялось удельное содержание конденсационной воды. Результаты расчетов приведены на рис. 5, на котором также показаны термобарические параметры на забое и устье скважины № 1 в 2020 г.
Таблица 6. Сравнение скоростей потока газа в хвостовике со скоростью опрокидывания кольцевой пленки конденсата при двухфазном течении газоконденсатной смеси Table 6. Comparison of the gas flow rates in the tailpipe with the overturning rate of the annular condensate film in the two-phase flow of the gas-condensate fluid
Год разработки Development year Скважина/Well
1 2 3
Wrev W Wrev W Wrev W
2020 - - - - - -
2030 - - 0,72 1,17 - -
2040 1,32 3,75 1,29 2,17 0,86 4,66
Из рис. 5 следует, что увеличение удельного содержания воды смещает границу ее фазового состояния в сторону более высоких температур. При удельном содержании воды свыше 5,7 г/м3 уже в 2020 г. не исключается появление жидкой водной фазы с образованием эмульсий типа «вода в масле» и возникновение рисков углекислотной коррозии хвостовика из-за возможной инверсии эмульсии при дальнейшем росте содержания воды в продукции и проявлении ее электрической проводимости.
Поэтому мониторинг добычи воды раздельно по рассматриваемым скважинам составляет обязательную часть программы исследовательских работ, а в случае отклонения от проектных значений с резким ростом дебита воды потребуется техническая диагностика состояния хвостовика и проведение повторного моделирования фазового поведения газоконденсатной смеси с учетом актуализированных показателей.
Содержание воды в пластовом флюиде
^—4,6 г/гмЗ 5,7 г/мЗ - —6,8 г/мЗ 7,9 г/мЗ Забой
Г/ —1—1—1—
Устье III / //
S3
ЖИДКОСТЬ + ПАР ПАР
/ / / г / f/
/ // / / / / /
0 20 40 60 80 100 120 140
Температура, °С
Рис. 5. Смещение границ двухфазной области воды, содержащейся в газоконденсатной смеси скважины № 1 с компонентно-фракционном составом (табл. 1), при разных удельных содержаниях воды Fig. 5. Displacement of the two-phase water region boundaries contained in the gas-condensate fluid of well no. 1 with a component-fractional composition (table 1), at different specific water contents
Заключение
Для определения предельных значений термобарических параметров технологических режимов эксплуатации высокотемпературных ачимовских газоконденсатных скважин с содержанием в добываемой продукции диоксида углерода, обеспечивающих физико-химические условия отсутствия электрохимической коррозии забойного оборудования -хвостовиков из углеродистой стали на длительную перспективу, предложен комплексный методический подход, включающий моделирование фазового поведения влажной газо-конденсатной смеси на забое и устье скважин и гидродинамический расчет многофазного потока в скважине.
На основе предлагаемого подхода изучено влияние изменения термобарических условий эксплуатации трех опытных скважинах на проектных технологических режимах в процессе разработки газоконден-сатной залежи на фазовое поведение пластового флюида в интервале установки незащищенных от уг-лекислотной коррозии хвостовиков с учетом влагосо-держания добываемой продукции. Расчеты показали, что в двух из трех скважин из-за некоторого различия исходных составов газа и технологических режимов жидкая углеводородная фаза - нестабильный конденсат начинает выделяться только к двадцатому году эксплуатации, а по одной - к десятому году; при этом формирования жидкой водной фазы в забойных условиях не наблюдается ни по одной из скважин в течение всего прогнозного периода.
Учитывая, что вода не образует химических соединений с углеводородами, и при низкой концентрации в нестабильном конденсате (порядка 0,6 % мольн.) она находится в молекулярно-растворенном состоянии, сконденсированная жидкость не обладает электрической проводимостью. С другой стороны, продуктивные характеристики рассматриваемых трех скважин при проектных технологических режимах и диаметрах
хвостовиков и лифтовых колонн обеспечивают в течение всего прогнозного периода значения скорости восходящего потока пластового флюида, превосходящие критические для полного и непрерывного выноса жидкости (в случае ее образования) на поверхность, т. е. накопления ее на забое не происходит. Физико-химические свойства пластового флюида при высокой температуре пласта и гидродинамические условия на забое и устье каждой из скважин позволяют прийти к заключению о невозможности протекания электрохимической коррозии хвостовиков в течение 20-летнего периода во всех трех скважинах.
Однако в случае значительного превышения удельного содержания конденсационной или иной воды в продукции скважин необходимо предусмотреть проведение мероприятий по технической диагностики состояния хвостовиков и повторное моделирование с учетом отклонений от проектных технологических режимов эксплуатации скважин. Поэтому освоение скважин с хвостовиками из углеродистой стали после операции ГРП необходимо проводить в сжатые сроки, а в процессе ГРП использовать водные растворы с добавкой ингибиторов коррозии либо иные реагенты, не обладающие коррозионной агрессивностью.
Таким образом, предлагаемый подход позволяет:
• обосновывать антикоррозионные режимы на стадии проектирования скважин, когда в проектных решениях возможна замена части подземного оборудования из дорогих коррозионностойких сталей на более дешевые углеродистые стали без снижения уровня надежности и безопасности;
• прогнозировать возможность протекания процессов углекислотной коррозии не только подземного оборудования, но и трубопроводов обвязки скважин и системы сбора газа, выполненных из углеродистой стали.
СПИСОК ЛИТЕРАТУРЫ
1. Маркин А.Н., Низамов Р.Э. С02-коррозия нефтепромыслового оборудования. - М.: ОАО «ВНИИОЭНГ», 2003. - 187 с.
2. Кузнецов В.П. Основные факторы углекислотной коррозии газоконденсатных скважин и их изменение в процессе эксплуатации месторождений: дис. ... канд. техн. наук. - Краснодар, 1974. - 143 с.
3. Легезин Н.Е. Противокоррозионная защита систем добычи, сбора и транспорта природного газа с применением ингибиторов: дис. ... д-ра техн. наук. - М, 1988. - 284 с.
4. Освоение участков ачимовских отложений ООО «Газпром добыча Уренгой» / А.Ю. Корякин, И.В. Игнатов, А.Ю. Неудахин, М.Г. Жариков, С.А. Скрылев, А.Н. Нестеренко, В.П. Тюрин // Научный журнал российского газового общества. - 2017. -№ 3. - С. 21-28.
5. Курчиков А.Р., Бородкин В.Н., Хромцова А.В. Условия формирования и атлас текстур пород ачимовского комплекса севера Западной Сибири. - Новосибирск: Изд-во СО РАН, Зап.-Сиб. фил. ин-та нефтегазовой геологии и геофизики им. А.А. Трофимука СО РАН, 2010. - 130 с.
6. Корякин А.Ю., Кобычев В.Ф. Перспективные направления развития Уренгойского комплекса // Сборник научных трудов ООО «Газпром добыча Уренгой». - М.: ООО «ИД Недра», 2018. - С. 350-355.
7. Yusupov A. Mathematical simulation of the rate of carbon dioxide corrosion at the facilities of Gazpromdobycha Urengoy LLC // E3S Web of Conferences. - 2019. - V. 121. - P. 1-5.
8. Koriakin A. Carbon dioxide corrosion at the objects of the second district of Achimovsk deposits of Urengoy oil and gas bearing complex // International Journal of Mechanical Engineering and Technology. - 2018. - V. 9. - P. 1073-1080.
9. Опыт подбора ингибиторов коррозии для защиты от углекис-лотной коррозии объектов второго участка ачимовских отложений Уренгойского нефтегазоконденсатного месторождения / А.Ю. Корякин, Д.В. Дикамов, И.В. Колинченко, А.Д. Юсупов, Д.Н. Запевалов, Р.К Вагапов // Оборудование и технологии для нефтегазового комплекса. - 2018. - № 6. - С. 48-55.
10. Совершенствование конструкций скважин на ачимовские залежи Уренгойского нефтегазоконденсатного месторождения / А.Ю. Корякин, М.Г. Жариков, И.А. Яскин, С.А. Скрылев, Р.П. Гресько // Газовая промышленность. - 2017. - № 9. -С. 76-79.
11. Скорчеллетти В.В. Теоретическая электрохимия. - Л.: Госхи-миздат, 1963. - 609 с.
12. Modeling binary mixtures of water + light hydrocarbon using the perturbed-chain statistical associating fluid theory with induced association: improvement in describing all equilibrium phases /
P.T. Sugata, S.K. Jeffrey, D.V. Steven, M.L. Rosaly // ACS Earth and Space Chemistry. - 2019. - V. 3. - Iss. 11. - P. 2569-2581.
13. Soave G. Equilibrium constants from a modified Redlich-Kwong equation of state // Chemical Engineering Science. - 1972. -V. 27. - P. 1197-1203.
14. Li C., Jia W., Wu X. Application of Lee-Kesler equation of state to calculating compressibility factors of high pressure condensate gas // Energy Procedia. - 2012. - V. 14. - P. 115-120.
15. Peng D.Y., Robinson D.B. A new two-constant equation of state // Industrial & Engineering Chemistry Fundamentals. - 1976. -V. 15. - P. 59-64.
16. Robinson D.B., Peng D.Y., Chung S.Y. The development of the Peng-Robinson equation and its application to phase equilibrium in a system containing methanol // Fluid Phase Equilibria. -1985. - V. 24. - P. 25-41.
17. Модификация кубических уравнений Пенга-Робинсона и Брусиловского для описания поведения воды и метанола / Ю.В. Калиновский, А.В. Минеев, А.И. Пономарев // Нефтегазовое дело. - 2006. - Т. 4. - № 1. - С. 293-297.
18. Пономарев А.И., Калиновский Ю.В. Комплексирование газо-конденсатных и гидродинамических исследований скважин для определения характеристик многофазного потока // Наука и техника в газовой промышленности. - 2017. - № 2 (70). -С. 7-16.
19. High pressure phase behavior modeling of asymmetric alkane + alkane binary systems with the RKPR EOS / M.D. Cismondi, M.V. Galdo, M.J. Gomez, N.G. Tassin, M. Yanes // Fluid Phase Equilibria. - 2014. - V. 362. - P. 125-135.
20. Калашников О.В. Моделирование фазового поведения углеводородов: выбор уравнения состояния // Экотехнологии и ресурсосбережение. - 2003. - № 1.- С. 22-30.
21. Lopez-Echeverry J., Reif-Acherman S., Araujo-Lopez E. Peng-Robinson equation of state: 40 years through cubics // Fluid Phase Equilibria. - 2017. - V. 447. - P. 39-71.
22. Peng D.Y., Robinson D.B., Bishnoi P.R. The use of the Soave-Redlich-Kwong equation of state for predicting condensate fluid behavior // 9th World Petroleum Congress. - Tokyo, 1975. -P. 1-12.
23. Literature review on modeling and simulation of energy infrastructures from a resilience perspective / J. Wanga, W. Zuo, L. Barbarigos, X. Lu, J. Wang, Y. Lind // Reliability Engineering & System Safety. - 2019. - V. 183. - P. 360-373.
24. Движение газожидкостных смесей в трубах / В.А. Мамаев, Г.Э. Одишария, О.В. Клапчук, А.А. Точивин, Н.И. Семенов. -М.: Недра, 1978. - 270 с.
Поступила: 29.07.2020 г.
Информация об авторах
Пономарев А.И., доктор технических наук, профессор, заведующий кафедрой разработки и эксплуатации газовых и нефтегазоконденсатных месторождений Уфимского государственного нефтяного технического университета.
Иванов Н.В., аспирант кафедры машин и оборудования нефтяной и газовой промышленности Российского государственного университета нефти и газа (НИУ) им. И.М. Губкина.
Юсупов А.Д., аспирант кафедры разработки и эксплуатации газовых и нефтегазоконденсатных месторождений Уфимского государственного нефтяного технического университета.
UDC 622.279, 620.193
NEW METHODICAL APPROACH TO THE CO2-CORROSION ESTIMATION OF HIGH-TEMPERATURE GAS-CONDENSATE WELLS EQUIPMENT
Alexander I. Ponomarev1,
Nikolay V. Ivanov2,
Alexander D. Yusupov1,
1 Ufa State Petroleum Technological University, 1, Kosmonavtov street, Ufa, 450062, Russia.
2 Gubkin Russian State University of Oil and Gas, 65, Leninsky avenue, Moscow, 119991, Russia.
The relevance. The prospects for increasing gas condensate - valuable raw material for oil refining and petrochemicals, are largely associated with the further development of the achimov gas condensate deposits in the north of Western Siberia. Large depths, abnormally high reservoir pressure and high temperature of low-productivity formations, along with the presence of carbon dioxide in the gas-condensate formation fluid, put forward, first of all, increased demands on the quality of well completion and hydraulic fracturing operations, the reliability of structures and materials of underground well equipment. These factors are responsible for the high cost of project implementation and operational risks of gas and condensate production, and therefore justification of safe operating conditions for achimov is an urgent scientific and technical task.
The main aim of the research is the limiting thermobaric parameters determination of the technological operating modes of the achimov gas condensate wells, taking into account the production of CO2 in products.
Objects: reservoir gas-condensate fluid and bottomhole equipment of three achimov gas-condensate wells with carbon steel tailpipes, unstable to CO2-corrosion.
Methods: modeling the phase behavior of the produced reservoir fluid in bottomhole conditions by the GasCondNeft software. Results. Phase behavior thermodynamic calculations of the reservoir fluid, taking into account its moisture content, showed that a reservoir pressure decrease at the location of the three considered wells during the development of the achimov deposits for 20 years under project operation modes leads to the formation of a two-phase fluid (gas-unstable condensate) in bottomhole thermobaric conditions. Moreover, the aqueous liquid phase due to the high temperature of the flow at the well bottomhole is not formed during the entire calculation period. Hydrodynamic parameters calculations of the gas-liquid upward flow showed that high flow rates of well products provide conditions for the complete and continuous removal of unstable condensate by the gas flow from the bottom to the surface of all considered wells for the entire 20-year period, thereby preventing the physical and chemical conditions of electrolyte formation on tailpipes surface and carbon dioxide corrosion.
Key words:
Gas condensate fluid, carbon dioxide, tailpipe, well operation parameters, simulator, phase diagrams,
water and hydrocarbon dew point lines, multiphase flow regime, carbon dioxide corrosion, anticorrosive well operation.
REFERENCES
1. Markin A.N., Nizamov R.E. CÛ2-korroziya neftepromyslovogo oborudovaniya [CO2-corrosion of oilfield equipment]. Moscow, VNIIOENG Publ., 2003. 187 p.
2. Kuznetsov V.P. Ûsnovnye factory uglekislotnoy korrozii gazokon-densatnykh skvazhin i ilh izmenenie v protsesse ekspluatatsii mes-torozhdeniy Dis. Kand. nauk [Main factors of carbon dioxide corrosion of gas condensate wells and their changes during field operation. Cand. Diss.]. Krasnodar, 1974. 143 p.
3. Legezin N.E. Protivokorrozionnaya zashchita system dobychi, sbora i transporta prirodnogo gaza s primeneniem ingibitorov. Dis. Dokt. nauk [Corrosion protection of natural gas production, gathering and transportation systems using inhibitors. Dr. Diss.]. Moscow, 1988. 284 p.
4. Koryakin A.Yu., Ignatov I.V., Neudakhin A.Yu., Zharikov M.G., Skrylev S.A., Nesterenko A.N., Tyurin V.P. Development of blocks of achimov deposits Gazprom dobycha Urengoy LLC. Nauchny zhurnal rossiyskogo gazovogo obshchestva, 2017, no. 3, pp. 21-28. In Rus.
5. Kurchikov A.R., Borodkin V.N., Khromctsova A.V. Usloviya formirovaniya i atlas tekstur porod achimovskogo kompleksa severa Zapadnoy Sibiri [Formation conditions and atlas of rocks of
the Achimov complex textures in the north of Western Siberia]. Novosibirsk, SO RAN Publ., 2010. 130 p.
6. Koryakin A.Yu., Kobychev V.F. Perspektivnye napravleniya razvitiya Urengoyskogo kompleksa [Prospective directions of development of the Urengoy oil and gas complex]. Sbornik nauch-nylh trudov OOO «Gazprom dobycha Urengoy» [Collection of research papers of Gazprom dobycha Urengoy LLC]. Moscow, Nedra Publ., 2018. pp. 350-355.
7. Yusupov A. Mathematical simulation of the rate of carbon dioxide corrosion at the facilities of Gazprom dobycha Urengoy LLC. E3S Web of Conferences, 2019, vol. 121, pp. 1-5.
8. Koriakin A. Carbon dioxide corrosion at the objects of the second district of Achimovsk deposits of Urengoy oil and gas bearing complex. International Journal of Mechanical Engineering and Technology, 2018, vol. 9, pp. 1073-1080.
9. Koryakin A.Yu., Dikamov D.V., Kolinchenko I.V., Yusupov A.D., Zapevalov D.N., Vagapov R.K Experience of corrosion inhibitors selection to protect the objects of the second site of achimovsky deposits of Urengoy oil and gas-condensate field from carbon acid corrosion. Equipment and technologies for oil and gas complex, 2018, no 6,pp. 48-55. In Rus.
10. Koryakin A.Yu., Zharikov M.G., Yaskin I.A., Skrylev S.A., Gres-ko R.P. Improving the well design of the Achimovsk deposit of the
Urengoy oil and gas condensate field. Gas Industry, 2017, no. 9, pp. 76-79. In Rus.
11. Skorchelletti V.V. Teoreticheskaya elektrokhimiya [Theoretical electrochemistry]. Leningrad, Goskhimizdat Publ., 1963. 609 p.
12. Sugata P.T., Jeffrey S.K., Steven D.V., Rosaly M.L. Modeling binary mixtures of water + light hydrocarbon using the perturbed-chain statistical associating fluid theory with induced association: improvement in describing all equilibrium phases. ACS Earth and Space Chemistry, 2019, vol. 3, pp. 2569-2581.
13. Soave G. Equilibrium constants from a modified Redlich-Kwong equation of state. Chemical Engineering Science, 1972, vol. 9, pp. 1197-1203.
14. Li C., Jia W., Wu X. Application of Lee-Kesler equation of state to calculating compressibility factors of high-pressure condensate gas. Energy Procedia, 2012, vol. 14, pp. 115-120.
15. Peng D.Y., Robinson D.B. A new two-constant equation of state.
Industrial & Engineering Chemistry Fundamentals, 1976, vol. 15, pp. 59-64.
16. Robinson D.B., Peng D.Y., Chung S.Y. The development of the Peng-Robinson equation and its application to phase equilibrium in a system containing methanol. Fluid Phase Equilibria, 1985, vol. 24, pp. 25-41.
17. Kalinovskiy Yu.V., Mineev A.V., Ponomarev A.I. Modification of the Peng-Robinson and Brusilovsky cubic equations to describe the behavior of water and methanol. Oil and gas business, 2006, vol. 1, pp. 293-297. In Rus.
18. Ponomarev A.I., Kalinovskiy Yu.V. Integration of gas condensate and hydrodynamic studies of wells to determine the characteristics of multiphase flow. Science and Technology in the Gas Industry, 2017, no. 2, pp. 7-16. In Rus.
19. Cismondi M.D., Galdo M.V., Gomez M.J., Tassin N.G., Yanes M. High pressure phase behavior modeling of asymmetric alkane + alkane binary systems with the RKPR EOS. Fluid Phase Equilibria, 2014, vol. 362, pp. 125-135.
20. Kalashnikov O.V. Modeling the phase behavior of hydrocarbons: the choice of the equation of state. Ecotechnology and resource saving, 2003, vol. 1, pp. 22-30. In Rus.
21. Lopez-Echeverry J., Reif-Acherman S., Araujo-Lopez E. Peng-Robinson equation of state: 40 years through cubics. Fluid Phase Equilibria, 2017, vol. 447, pp. 39-71.
22. Peng D.Y., Robinson D.B., Bishnoi P.R. The use of the Soave-Redlich-Kwong equation of state for predicting condensate fluid behavior. IX World Petroleum Congress. Tokyo, 1975. pp. 1-12.
23. Wanga J., Zuo W., Barbarigos L., Lu X., Wang J., Lind Y. Literature review on modeling and simulation of energy infrastructures from a resilience perspective. Reliability Engineering & System Safety, 2019, vol. 183, pp. 360-373.
24. Mamaev V.A., Odishariya G.E., Klapchuk O.V., Tochivin A.A., Semenov N.I. Dvizhenie gazozhidkostnylh smesey v trubakh [The movement of gas-liquid mixtures in pipes]. Moscow, Nedra Publ., 1978. 270 p.
Received 29 July 2020.
Information about the authors
Alexandr I. Ponomarev, Dr. Sc., professor, Ufa State Petroleum Technological University. Nikolay V. Ivanov, postgraduate student, Gubkin Russian State University of Oil and Gas. AlexandrD. Yusupov, postgraduate student, Ufa State Petroleum Technological University.