CC BY
71
УДК 620.193
С. А. Долгих, Ф. Ш. Шакиров, С. С. Виноградова, Б. Л. Журавлев, И. Н. Галиев
ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ РАСЧЕТНОЙ МОДЕЛИ КАТОДНОЙ ЗАЩИТЫ
ОБСАДНЫХ КОЛОНН ПО РЕЗУЛЬТАТАМ ИССЛЕДОВАНИЙ В ПРОМЫСЛОВЫХ УСЛОВИЯХ
Ключевые слова: обсадная колонна, катодная защита, ток защиты, смещение потенциала.
Рассмотрена стандартизованная методика определения параметров катодной зашиты обсадных колонн. Разработано программное обеспечение для реализации методики. Проведены расчеты параметров катодной защиты для двух обсадных колонн скважин НГДУ «Бавлынефть». Показана связь срока эксплуатации скважины с параметрами катодной защиты.
Keywords: well casing, cathodic protection, protection current, potential displacement.
Standardized method of determining the parameters of cathodic protection of an oil casing are consider. Software has been developed to implement the methodology. The calculations of the parameters of cathodic protection for two of the well casing NGDU "Bavlyneft" were made. Displayed communication parameters of cathodic protection on the life of the well.
Ресурс обсадных колонн повышают, осуществляя ряд комплексных мероприятий, включающих использование катодной защиты [1,2].
Основным параметром катодной защиты является величина защитного тока. Существует ряд методов, позволяющих определить величину требуемого тока катодной защиты [3,4]: метод поляризационных кривых (электрического каротажа); метод определения профиля падения напряжения на эксплуатационной колонне; метод расчета сдвига потенциала в стволе скважины и величины сопротивления в системе скважина/грунт; метод моделирования катодной защиты скважины.
В данной работе рассматривается метод расчета сдвига потенциала в стволе скважины и сопротивления в системе скважина/грунт
Цель работы заключалась в разработке программного обеспечения для расчетно-экспериментальной методики определения параметров катодной защиты обсадных колонн с учетом конструкции скважин; проведении расчетов параметров модели и сопоставлении результатов расчетов и результатов экспериментальных исследований в промысловых условиях.
Метод расчета сдвига потенциала в стволе скважины и сопротивления в системе скважина/грунт, позволяет, зная величину защитного тока и смещение потенциала на устье скважины, рассчитать сдвиг потенциала на забое скважины [4] или, зная требуемую величину сдвига потенциала на забое скважины, рассчитать ток защиты. В соответствии с этим методом алгоритм расчета параметров катодной защиты состоит из четырех этапов.
Первый этап расчета включает определение геометрических и электрических параметров конструкции в соответствии с типом и характеристиками скважины. Для обсадной колонны, состоящей из трех участков, поперечное сечение металлической конструкции на разных уровнях определяется согласно схеме, представленной на рис. 1 [4].
Рис. 1 - Схема определения поперечного сечения обсадной колонны, состоящей из трех участков (слева направо): эксплуатационной колонны; эксплуатационной колонны и кондуктора; эксплуатационной колонны, кондуктора и направления
Формула определения общей площади поперечного сечения участков имеет вид:
= ■ d ■ е, м2 где d - наружный диаметр трубы; е - толщина стенки трубы.
Продольное омическое сопротивление обсадной колонны определяется по формуле:
rk =
Р-1 Sk
Ом
м
-1
(единица) - длина участка.
В систему электрических параметров при расчетах входят коэффициент затухания ак и характеристическое сопротивление ук,
определяемые по соответствующим формулам:
[м-1]
У к =
co -гк,Ом
ж - d,
где гк - продольное омическое сопротивление; гсо -удельное сопротивление покрытия; dk - наружный диаметр участка обсадной колонны.
Второй этап включает в себя испытания с подачей поляризующего тока (1п), при котором измеряют сдвиг потенциала на устье скважины (ип) и рассчитывают соотношение ип/1п, .
Третий этап методики предусматривает определение методом последовательных
а.. = r.
k
r
co
приближений величины удельного сопротивления покрытия (rco).
Расчеты проводят, начиная с нижнего участка, принимая произвольное смещение потенциала на забое U0 (например, U0 = 1).
Для колонны, состоящей из трех участков, используют три системы уравнений [4]. Для первого участка (забой):
fU1 = U0 • cosh (b1)
А \ I1 = ^ • sinh(b1)
I Yi
Для второго участка:
{U2 = U1 • cosh (b2) + Y2 • I1 • sinh (b2)
U
I2 = —- • sinh( b2) +11 • cosh( b2 )
Y2
Для третьего участка:
С С U3 = U2- cosh (Ьз) + Y3 • I2 • sinh (Ьз)
С -i U2
I I3 =—2• sinh(b3) +12-cosh(b3) L Y3 В этих уравнениях
bk = ak -Lk
где ak - коэффициент затухания; Lk - длина участка.
В данной методике предполагается, что величина rco, выражающая удельное сопротивление покрытия, неизменна по всей поверхности скважины для всех участков.
Для решения системы уравнений выражения для Ui и I1 системы А подставляют в систему В и соответственно выражения для U2 и I2 , полученные в системе В, подставляют в систему С. В результате подстановок получают систему С:
С
' <
72 cosh(b1) • cosh(b2 ) • cosh(b3 )+ • sinh(b1) • sinh(b2 )x
n
U3 = U0 • 73 xcosh(b3 )+ •cosh( b1 )sinh(b2 )• sinh(b3 )+
72
73 + sinh( b1 )'Cosh(b2 )• sinh(b3 )
n
1 Y2
-cosh(b1 )cosh(b2 )sinh(b3 )+-sinh(b1 )sinh(b2 )x
Y3 Y1Y3
I3 = Uq• xsinhb )+-cosh(b1 )• sinh(b2 )cosh(b3 )+-sinh(b1 )x
Y2 Y1
xcosh(b2 )• cosh(b3 )
В рассматриваемой модели принято, что значение соотношения (и3/!3 - отношение сдвига потенциала на устье скважины к току защиты) не зависит от величины и0, а зависит только от значений параметров Ь-|, Ь2, Ь3, у-ь у2, у3, которые являются функциями гсо. Подставляя параметр гсо в выражение для и3 и !3 и приравнивая соотношение и3/!3 к полученному экспериментально значению этого соотношения, методом последовательных
приближений рассчитывают значение гсо.
Четвертый этап методики включает подстановку значения «гсо» во все выражения для Ь-^
Ь2, Ь3, у1, у2 и Y3 в системе С/. Полученные уравнения связывают значения смещения потенциала на устье скважины (и3) со смещением потенциала на забое (и0), и величину смещения потенциала на забое (Ы0) с величиной тока защиты (|3):
из=К||*ио !з=К2*ио
ы (13 {
где К1 - коэффициент пропорциональности между сдвигами потенциала на устье и на забое; К2 -коэффициент пропорциональности между сдвигом потенциала на забое и током защиты
Используя данную систему уравнений, рассчитывают смещение потенциала на устье и забое скважины при любой заданной силе тока.
Для реализации рассмотренного алгоритма расчета параметров катодной защиты разработано программное обеспечение с использованием объектно-ориентированного языка программирования Delphi 7.
Диалоговое окно программы условно можно разделить на три области (рис.2).
диаметр, м длина, м
ИСХОДНЫЕ ДАННЫЕ
1. Конструкция скважины
число колонн N 1. нап равл яющая колонна кондуктор
промежуточная колонна обсадная колонна-хвостовик эксплуатационная колонна
2. Результаты испытаний
потенциал на устье скважины до подачи тока, В сипа тока А
потенциал на устье скважины после отключения тока, В РАСЧЕТНЫЕ КОЭФФИЦИЕНТЫ ПРОПОРЦИОНАЛЬНОСТИ Коэффициент пропорциональности КТ(В/В) между сдвигами потенциала на забое и на устье (из=К1'иа] I
Коэффициент пропорциональности К2(А/В) между сдвигам потенциала на забое и током защиты (13=4\£*иС) РАСЧЕТ ПАРАМЕТРОВ КАТОДНОЙ ЗАЩИТЫ ° Вариант!
Расчет сдвига потенцигла за забое иО;В;=
при заданном токе защиты 13(А)= счаетд |
Г Вариант2 в
Расчет тока защиты 13[А|=
при заданном сдвиге потенциала на заЗое иО (В)= В адить ]
Рис. 2 - Диалоговое окно программы для определения параметров катодной защиты обсадной колонны, состоящей из пяти участков
Область ввода исходных данных:
- ввод параметров колонны (число участков от одного до пяти);
- ввод параметров каждого участка колонны (м) -диаметр, длина и толщина стенки;
- ввод значения потенциала на устье скважины до подачи тока (В);
- ввод значения силы тока (А);
- ввод значения потенциала на устье скважины после отключения тока (В);
Область результатов расчета коэффициентов пропорциональности:
- коэффициент пропорциональности К1 между сдвигами потенциала на устье и на забое скважины;
- коэффициент пропорциональности К2 между сдвигом потенциала на забое и током защиты.
Область расчета параметров катодной защиты:
- Вариант 1 - Расчет сдвига потенциала на забое при заданном токе защиты;
- Вариант 2 - Расчет тока защиты при заданном сдвиге потенциала на забое.
Для выбора варианта расчета выделяют "Вариант 1" или "Вариант 2". В случае выделения "Вариант 1" задают ток защиты !з(А) и определяют сдвиг потенциала на забое и0 (В), в случае
выделения "Вариант 2" задают сдвиг потенциала на забое U0 (В) и определяют ток защиты !з(А).
Сопоставление результатов расчетов и результатов экспериментальных исследований в промысловых условиях провели для обсадной колонны скважины № 7263г НГДУ «Бавлынефть» (2012 г. пуск в эксплуатацию), состоящей из трех участков: направление, кондуктор, эксплуатационная колонна (рис. 3).
Направление
Кондуктор
Эксплуатационная колонна
ЕГ
та
Рис. 3 - Трехступенчатая обсадная колонна
На первом этапе исследования были проведены испытания с подачей тока поляризации (!п=6А), при этом потенциал на устье скважины до подачи тока (и0, В) составлял -0,757 В (МСЭ), а после отключения тока ,В) составил 2,9 В. Рассчитанное значение соотношение (Л, В/А) при этом составило 0,357.
Результаты определения соотношение (Д, В/А) и удельного сопротивления покрытий ( гсо 2
Ом . м ) при значениях силы тока защиты от 3 до 10 А приведены в табл.1. В рассматриваемой методике предполагается, что значения соотношения между смещением потенциала на устье скважины и силой тока защиты, а так же значения удельного сопротивления покрытия не изменяются при изменение силы тока защиты. Однако испытания в промысловых условиях и результаты расчета показали, что значения этих параметров модели зависят от величины силы тока защиты при этом коэффициент вариации составляет 0,11.
После подстановки значения «гсо» во все
выражения для Ь1, Ь2, Ь3, у1, у2 и Y3 были получены уравнения, связывающие значения смещения потенциала на устье скважины (и3), со смещением потенциала на забое (и0), и величину смещения потенциала на забое (и0) с величиной тока защиты (!3) при рассматриваемых значениях силы тока:
= Ир ■ 1.25
,21
Для силы тока 3 А: Для силы тока 4 А: Для силы тока 5А: Для силы тока 6А: Для силы тока 7А:
Ш, = VE
U =
tó = UP ■ 1, Ц = Up * í № = Ц r ц = Ц * s;
ft, =4e.-1,
Ц = Up * i № = ц -r Ц = ÜP * i!
17 17 141? 19 33
Для силы тока 8А: _
си, =
Для силы тока 9А: ]» _
ГО. =
Для силы тока 10А: 1т —
Ц 1,19 Ге * 3,57 = С в-JJ
UP ■ 1,12 UL + 4U1
Результаты расчета сдвига потенциала на забое (Uo) с использованием значения силы тока защиты U., =f(l3) и сдвига потенциала на устье скважины líc =f(U3) и усредненных значений
ю,фф_ — ■
табл.1.
Таблица 1 - Результаты определения сдвига потенциала на устье скважины, измеренные и рассчитанные соотношения между сдвигом потенциала и током защиты
Сила Сдвиг А,
тока потенц В/А rco , =f(Ua), =f(Ia),
(Ja, иала Ом ■ В В
А) (U3 В) м2
3 0,783 0,261 373 0,65 0,80
4 1,243 0,311 450 1,04 1,06
5 1,843 0,369 541 1,54 1,33
6 2,143 0,357 523 1,88 1,59
7 2,363 0,338 492 1,97 1,86
8 2,663 0,333 485 2,22 2,12
9 2,843 0,316 458 2,37 2,39
10 2,963 0,296 428 2,47 2,65
Данные табл. 2 показали, что при использовании в расчетах средних значений коэффициентов К1 и К2, результаты расчетов сдвига потенциала на забое и0=1(из) близки к результатам расчетов сдвига потенциала на забое и0=Т(!з).
В качестве примера на рис. 4 показано диалоговое окно программы с вариантом расчета сдвига потенциала на забое при токе защиты 6А.
Оценку степени влияния срока эксплуатации колонны на значения параметров модели провели на основании расчетов для обсадной колонны скважины № 22505 НГДУ «Бавлынефть» (1980 г. ввод в эксплуатацию), состоящей из трех участков: направление, кондуктор, эксплуатационная колонна (рис. 3).
Сдвиг потенциала, измеренное
соотношение(д В/А) и удельное сопротивления покрытия для для обсадной колонны скважины № 22505 составили: сила тока (^=6 А), потенциал после отключения тока, (и1= -0,97В), сдвиг потенциала (из =0,26В), измеренное соотношение (Д=0 043 В/А), удельное сопротивление покрытия (гсо =30 Ом* м2).
Уравнения для расчета смещения потенциала на устье и забое скважины при любой заданной силе тока имеют вид: Ц, = иЕ -2М (э = Ч * №№
D
R
A
Результаты расчета сдвига потенциала на устье (из,В) и на забое (и0,В) с использованием этих уравнений представлены в табл.2.
PjnCTA Ркчгт C^riiTiTh Ji р Ы ' i
ИСХОДНЫЕ ДАННЫЕ 1, Конструкция скважины
число колонн ZI направляющая колонна кондуктор
диаметр, м
|o.«Ä
толщина, м
Jö.oon
эксплуатационная колонна К1*5 2. Результаты испытаний потенциал на устье скважины до подачи тока. В сила тока А
потенциал на устье скважины после отключения тока. В РАСЧЕТНЫЕ КОЭФФИЦИЕНТЫ ПРОПОРЦИОНАЛЬНОСТИ Коэффициент пропорциональности К1(В/В) между сдвигами потенциала на забое и на устье (иЗ=К1*иО) г"
Коэффициент пропорциональности К2(А/Б) между сдвигам потенциала на забое и током защиты (13=К?Ц0) I3-19 РАСЧЕТ ПАРАМЕТРОВ КАТОДНОЙ ЗАЩИТЫ
Вариант!
Расчет сдвига потенциала за забое У0(В)=1! н
при заданном токе защиты 13(А)= о«*«
Вариант? I
Расчет тока защиты 13(А)= I
при заданном сдвиге потенциала на забое иО (В)= В
Рис. 4 - Результаты расчетов «Вариант 1» скважине № 7263г НГДУ «Бавлынефть»
на
Таблица 2 - Результаты определения потенциала на устье и на забое скважины
сдвига
Сила Сдвиг Сдвиг
тока потенциала на потенциала на
защиты устье забое
А) (из В) (И В)
3 -0,128 0,048
4 -0,171 0,064
5 -0,213 0,080
6 -0,256 0,096
7 -0,299 0,112
8 -0,341 0,128
9 -0,384 0,144
10 -0,426 0,160
Сопоставление данных табл. 1 и 2 показывает существенное различие в сдвигах потенциалов, как на
устье так и на забое скважин при одной и той же силе тока защиты для скважины № 7263г НГДУ «Бавлынефть» (2012 г. ввод в эксплуатацию) и для скважины № 22505 НГДУ «Бавлынефть» (1980 г. ввод в эксплуатацию). Это различие в основном связано, с тем что удельное сопротивление покрытия для обсадной колонны скважины № 7263г
НГДУ «Бавлынефть» составляет 469 Омм , а для обсадной колонны скважины № 22505 НГДУ
«Бавлынефть» - 30 Ом м .
Выводы
1. Разработано программное обеспечения для расчетно-экспериментальной методики определения параметров катодной защиты обсадных колонн с учетом конструкции скважин.
2. Установлено, что соотношение между смещением потенциала на устье скважины и силой тока защиты, а так же рассчитанное значение удельного сопротивления покрытия зависят от величины силы тока защиты.
3. Показано, что рассматриваемая методика позволяет эффективно оценивать удельное сопротивление покрытия колонны, которая в свою очередь зависит от возраста скважины и от условий эксплуатации.
Литература
1. Долгих, С.А. Расчет распределения плотности тока защиты по глубине обсадной колонны нефтяной скважины / С. А. Долгих, Р. А. Кайдриков, Б. Л. Журавлев, В.Э. Ткачева // Вестник Казан. технол. ун-та.
- 2012. - № 20. - С. 191-193.
2. Долгих, С.А. определение тока защиты и расчет смещения потенциала на забое обсадной колонны нефтяной скважины / С.А. Долгих, Б.Л. Журавлев, Р.А. Кайдриков, В. Э. Ткачева // Вестник Казан. технол. ун-та.
- 2012. - № 22. - С. 241-244
3. EN 12954, Cathodic protection of buried or immersed metallic structures. General principles and application for pipelines, 2001.
4. EN 15112, Äußerer kathodischer Korrosionsschutz von Bohrlochverrohrungen; Deutsche Fassung EN 15112:2006
© С. А. Долгих - асп. каф. технологии электрохимических производств КНИТУ, с.н.с. ТатНИПИнефть г.Бугульма, dolg@tatnipi.ru, Ф. Ш. Шакиров - зав. лаб. защиты от коррозии института ТатНИПИнефть г.Бугульма, fsharirov@tatnipi.ru, С. С. Виноградова - канд. техн. наук, доц. каф. технологии электрохимических производств КНИТУ, vsvet@kstu.ru; Б. Л. Журавлев - д-р хим. наук, проф. той же кафедры, bgyr@kstu.ru; И. Н. Галиев - асп. той же кафедры.
