УДК658.012.011.56:658.512
СТОХАСТИЧЕСКИЕ МОДЕЛИ И МЕТОДЫ В СИСТЕМАХ ЭЛЕКТРОХИМИЧЕСКОЙ ЗАЩИТЫ ПОДЗЕМНЫХ ТРУБОПРОВОДОВ
ТЕВЯШЕВ А.Д., СТРИЖАК Л.В.,ТКАЧЕНКО А.А.
Приводятся стохастические модели и методы решения задачи оперативного планирования режимов работы систем электрохимической защиты подземных трубопроводов в условиях риска и неопределённости. Ээффек-тивность моделей и методов иллюстрируется на примере решении задачи оперативного планирования режимов работы систем электрохимической защиты участка газораспределительных сетей высокого давления г. Харькова.
Введение
В настоящее время в состав ГТС Украины входят: газопроводов высокого давления протяженностью свыше 37,5 тыс.км и газораспределительных трубопроводов низкого давления свыше 255 тыс. км [1].
Трубопроводы, уложенные в грунт, подвергаются воздействию различных физических и грунтоклиматических факторов, в той или иной степени влияющих на их состояние. Нарушения изоляции, воздействие давления грунта, влияние блуждающих токов в земле, значительные перепады температур и внутреннего давления ведут к существенному сокращению срока службы трубопроводов. Поэтому одним из решающих факторов, определяющих надежность трубопроводных систем на практике, является состояние их антикоррозионной защиты. На протяжении последних нескольких лет почвенная коррозия становится причиной около 30% отказов на трубопроводах. В результате аварии только стоимость капитального ремонта 1 км магистрального газопровода диаметром 1000 мм составляет около $1,5 млн грн, не считая прочих негативных последствий, имеющих зачастую намного больший финансовый эффект.
Как правило, на практике не удается добиться полной сплошности изоляционного покрытия: различные виды покрытия имеют различную диффузионную проницаемость и поэтому обеспечивают различную изоляцию трубы от окружающей среды. В процессе строительства и эксплуатации в изоляционном покрытии возникают трещины, задиры, вмятины и другие дефекты. Наиболее опасными являются сквозные повреждения защитного покрытия, где, практически, и протекает грунтовая коррозия. В ДК «Укртрансгаз» более половины трубопроводов находятся в эксплуатации более 20 лет. Это означает, что вероятность возникновения коррозионного разрушения газопроводов таких участков возрастает ввиду тенденции существующих на сегодняшний день изоляционных материалов к ухудшению диэлектрических свойств с течением времени.
Активный метод защиты от коррозии путем катодной поляризации основан на снижении скорости растворения металла по мере смещения его потенциала коррозии в область более отрицательных значений, чем естественный потенциал. Прилагая между поверхностью металла трубы и грунтом электрический ток, необходимо достигнуть снижения потенциала в дефектных местах изоляции трубы до значения ниже критерия защитного потенциала, равного 0,85 В. В результате этого скорость коррозии снижается до 10 мкм в год.
Для повышения оперативности сбора, обработки и анализа информации о степени защищенности объектов магистральных газопроводов от коррозии и, прежде всего, для обеспечения дистанционного мониторинга и регулирования защитных параметров установок катодной защиты используются системы анализа и оперативного управления режимами работы станциями электрохимической защиты (ЭХЗ) подземных трубопроводов.
В качестве основных элементов системы ЭХЗ подземных трубопроводов используются станции катодной защиты (СКЗ), которые обеспечивают необходимое значение защитного потенциала на всей протяженности защищаемого подземного трубопровода. В свою очередь величина защитного потенциала «труба-земля» является случайной величиной, зависящей от многочисленных факторов. Основными из них являются: фактические режимы работы СКЗ; фактическое состояние изоляционного покрытия трубопровода в любой точке рассматриваемого участка в каждый момент времени и соответствующее ему переходное сопротивление изоляционного покрытия; фактическое значение сопротивления грунта вдоль трассы трубопровода. Это приводит к необходимости рассматривать систему ЭХЗ подземных трубопроводов как стохастическую систему, функционирующую в стохастической среде. Известно [2], что наиболее эффективным способом управления такими системами являются двухэтапные процедуры, включающие в себя решение двух задач:
- оперативное планирование режимов работы СКХ на основании прогнозов (условных математических ожиданий) основных возмущающих факторов системы;
- стабилизация величины защитного потенциала после наблюдения фактических реализаций возмущающих факторов.
В данной работе приведены стохастические модели и методы решения задачи оперативного планирования режимов работы системы СКЗ, обеспечивающей ЭХЗ заданного участка подземных трубопроводов.
Цель исследования - повышение эффективности работы системы электрохимической защиты подземных трубопроводов путем решения задачи оперативного планирования режимов работы СКЗ трубопроводов.
РИ, 2008, № 2
61
1. Математическая модель распределения защитного потенциала «труба-земля»
Детерминированная модель
Для участка подземного трубопровода с подключенными к нему СКЗ разность потенциалов «труба-земля» в каждой точке на всей протяженности трубопровода для каждой СКЗ складывается из двух составляющих: положительного потенциала грунта ЦТр(х), создаваемого электрическим полем анодных зазем-лителей, и отрицательного потенциала поверхности трубы UT(x), возникающего за счет протекания катодного тока вдоль трубопровода.
Изолированный подземный трубопровод, защищаемый станциями катодной защиты, можно представить в виде протяженного электрического проводника постоянного тока с утечкой (рис.1) [3].
Величина Rпр [Ом/м] - электрическое сопротивление единицы длины протяженного проводника с утечкой (удельное продольное сопротивление, погонное сопротивление).
Величина Rиз [Омм] - переходное сопротивление утечке, характеризует электрическое сопротивление между отрезком протяженного проводника и соприкасающейся с ним средой.
При наличии какой угодно приложенной к протяженному проводнику нагрузки в любой точке х проводника имеем, потенциал UT(x) и протекающий по проводнику ток i, зависящие от координаты х.
Рассмотрим приведенную на рис.1 схему замещения протяженного проводника с утечкой. Пусть эта схема представлена в виде последовательно соединенных
продольных элементов сопротивлений R прАх , к кон-
цам которых параллельно подключены элементы пе-
R
из
реходных сопротивлений . Пусть, кроме этого,
на участке ОС имеется некоторый промежуток АВ, внутри которого к протяженному проводнику не приложено никакой нагрузки. Ввиду того, что промежуток АВ не загружен, ток i на всем этом промежутке не имеет разрыва, т.е. изменяется непрерывно.
Ах
■Ax-
Ax
• 2Л. і + Аі
RM
с R Ax b
А і
t
2
Rsl
Ax
М
ж
2
х
а
2
Рис.2. Отрезок ab схемы замещения протяженного проводника с утечкой
Выделим на промежутке АВ достаточно малый отрезок ab длиной Дх. Пусть начальная точка а отрезка ab имеет координату х. Следовательно, координата точки b есть х + Дх. К средней точке с отрезка ab подключен один из элементов переходных сопротив-
R
RU3
лений ——, в связи с чем между точками а и с и между точками с и b включены по продольному
элементу сопротивления, имеющие величину
R пр Ах 2
На рис. 2, для удобства показан отрезок ab в более крупном масштабе, где
ab = Дх, (1)
ac = cb = ■
Дх
(2)
Пусть в точке а имеются потенциал UT и ток i, протекающий от точки а к точке с, а в точке b -потенциал UT + Ди^ и ток i + Ді, протекающий от с к b, т.е. приращению Дх координаты х соответствуют приращения AUt и Ді потенциала Ut и тока i. Если ток направлен слева направо, т.е. если он имеет положительный знак, то Ut > Ut + AUt ,
вследствие чего AUt < 0 , т.е. приращение потенциала имеет отрицательный знак, т.е. знаки величин AUt и i противоположны. Как видно из рис. 2, это падение напряжения равно:
-Д UT =
. RпрАх i 2
+ (i + Д i)
R прАх 2
, (3)
A/a
1—1 . 1-І гЧ І—і
j Rn»Ax
-і I 2 J 2 J
Ax
ШШШШШШШШШЖШШШШШШШШШШШШШЖШШШЖіШШ; М Рис.1. Схема замещения протяженного проводника с М
утечкой
где Ді представляет собой ток на участке
ab = Дх между отрезком проводника и эквипотенциальной поверхностью ММ, т.е. ток
Ді есть приращение тока на участке Дх, которое будет положительным в случае его
входа в отрезок Дх (ток подтекания), и отрицательным в случае его выхода из отрезка Дх (ток утечки).
Величина тока д. определяется потенциалом Ut c точки с и переходным сопротивле-
62
РИ, 2008, № 2
R„
нием
, причем при положительном значении по-
Дх
тенциала UTc приращение Ді будет отрицательным, а при отрицательном потенциале VTc приращение положительно:
где коэффициент затухания а определяется выраже-
_ lR пр
нием а — a R , а постоянные А и В определяются
V RU3
для каждого участка непрерывного изменения тока из условий на границах того же участка.
ИТc — ИТ _ i‘
Ді = -
И
Тс
Ит — -
R пр Ах 2 ’
R пр Лх
R
из
Ах
Ru
■■Ах
(4) Значения положительного потенциала грунта И 1гр (х),
создаваемого электрическим полем анодного заземления і-й СКЗ, определяется выражением [4]:
(5)
Игр (х)
1 Ргр
2л(Ьі — ai)
Таким о бразом, при известных значениях потенциала и тока в точке подключения СКЗ можно последовательно получить значения потенциала Ит и тока і в каждой точке рассматриваемого участка.
Представим уравнения (3) и (5) следующим образом:
диТ
Дх
= -Rпрі -
R
пр
2
Ді,
(6)
(х - аі) + У(х ~ аі)2 + yf
(х - Ьі) +7(х - Ьі)2 + yA _ ’
(12)
аі, Ьі, yA - координаты расположения протяженного анодного заземления i-й СКЗ.
Функция распределения потенциала «труба-земля» Ит_з(х) для n-СКЗ находится следующим образом:
Ді _ Ит ^ Rпрі д
Дх R из 2RU3
n n
Ит-з (х) = Е иТ(х) -Е игр(х). (із)
і=1 і=1
Если отрезок Дх неограниченно уменьшать, устремляя его к значению дифференциала dx координаты x , то приращения ДИт и Ді будут соответственно стремиться к дифференциалам <іИт и d. Вторые слагаемые правых частей уравнений (6) и (7) превратятся в величины высшего порядка малости по сравнению с первыми слагаемыми, а сами уравнения из конечноразностных превратятся в дифференциальные.
Функция ИТ(х) - распределения отрицательного потенциала поверхности трубы каждой і-й СКЗ является решением системы дифференциальных уравнений:
dUTr)+R пр і =0, (8)
di ! иТ(х) _ 0 ^х RU3 (9)
і - функция распределения значений силы тока.
Решение системы дифференциальных уравнений в общем виде:
ИТ(х) = Ae
(10)
ЛіІ^пр~~^и
r(Ae
—х-
-Be
),
(11)
Стохастическая модель
Основными случайными величинами, влияющими на значение защитного потенциала «труба-земля», являются: фактическое состояние изоляционного покрытия трубопровода в точке x в момент времени t и соответствующее ему переходное сопротивление изоляционного покрытия, а также фактическое значение сопротивления грунта вдоль трассы трубопровода, зависящее от влажности почвы, ее состава, плотности, температуры и т.п.
Пусть P) - вероятностное пространство, где П - пространство элементарных исходов; F - а -алгебра подмножеств П; P - вероятностная мера на F.
Определим параметры RU3 и Ргр как случайные величины на пространстве (fi,F,P): Ргр = Ргр (WA),
RU3 = RU3 (<Д t) , где ш ЄО .
Как уже рассматривалось ранее, потенциал «труба-земля» в точке х представляет собой сумму потенциалов каждой СКЗ и является разницей отрицательного потенциала трубы и положительного потенциала грунта. Отрицательный потенциал трубы зависит от случайной величины сопротивления изоляции ИТ (х) = ИТ (х, RU3 (ш, t)). Положительный потенциал грунта зависит от случайной величины удельного сопротивления грунта Ир (х) = И\р (х, ргр (w,t)). Таким образом, потенциал «труба-земля» также представляет собой случайную величину:
ИТ-3(х, ^w) = Е ИТ(х, RU3(^,t)) - Е игр (х, Ргррt)).
і=1 і=1
РИ, 2008, № 2
63
2. Стохастические показатели оценки качества и эффективности функционирования системы ЭХЗ подземных трубопроводов
Повысить качество и эффективность функционирования систем ЭХЗ подземных трубопроводов - это значит обеспечить наиболее эффективную защиту подземных трубопроводов от электрохимической коррозии при одновременном уменьшении затратэлек-троэнергии на СКЗ.
Система ЭХЗ полностью выполняет свое функциональное назначение на интервале управления [0,T] в том и только в том случае, если [5]:
Vt Є [0, T] Л Vx Є L АУш ЄП,
P(W-3 < Ut(t, x, ш) < U+_3) = Po. (14)
Для реальных систем ЭХЗ условие (14) оказывается чрезвычайно жестким, и для достаточно больших интервалов управления [0,Т] ему не удовлетворяет практически ни одна реальная система ЭХЗ. Это означает, что для реальных систем ЭХЗ вероятность выхода значения разности потенциалов «труба-земля» за допустимые границы отлична от нуля.
Качество функционирования систем ЭХЗ на интервале времени [0,T] будем определять на основе локальных критериев, количественно выражающих целевое назначение системы, режим р аботы и параметры СКЗ. Оценки качества функционирования систем ЭХЗ на основе таких критериев должны в явном виде включать требования вида (14) и рассматриваться как стохастический функционал, определенный на пространстве случайных функций переходного сопротивления изоляционного покрытия подземного трубопровода и сопротивления грунта вдоль трассы трубопровода. В качестве конкретных локальных показателей качества функционирования систем ЭХЗ могут быть использованы необходимые статистические характеристики этого функционала, такие как: математическое ожидание и дисперсия.
Рассмотрим однониточный участок подземного трубопровода [0,L], на котором установлено n СКЗ с координатами xj, i = 1,2,...n. Назовем i-м участок
трубопровода между [xj_i, xj J точками подключения (i — 1 )-й и i-й СКЗ.
Функция распределения потенциала Ut _з (t, x, ш) на интервале |xj_і, xjJ между i-1-й и i-й СКЗ представ-
ляет собой выпуклую по x є |xi, xi+iJ функцию, имеющую два минимума в точках xj_i и xj и один максимум в некоторой внутренней точке xj интервала [xi-b xi ].
Локальной «диктующей точкой» системы ЭХЗ в момент времени t є [0,T] на интервале [xj_i, xjJ называется такая точка x*(t, ш) є |xi_i, xi J, в которой функция Ut_з (t, x, ш) достигает максимума.
Локальной «диктующей точке» соответствует точка максимума случайной величины Ut_з (t, x, w). Практически это означает, что локальная «диктующая точка» также является случайной величиной и для Уш Є ПЛ Vt Є [0,T] может принимать значения из некоторого открытого интервала Xj = (x^, xj ^), являющегося подмножеством множества [xj_i, xj J, т.е. X* c[xb xi+i] .
Подмножество Xj c [xj+i, xj] называется множеством «диктующих точек» системы ЭХЗ на интервале
[xi-l, xi] .
Глобальной «диктующей точкой» систем ЭХЗ на интервале |xj_i, xj J называется такая точка x* є X* , для которой выполняется условие
Vt Є [0, T] Л Vw Є Q :
*
xj = argmaxUT _з (t, x, ш) (,5)
xeX* . (15)
Знание местоположения глобальной « диктующей точки^) на интервале [xi_i, xiJ позволяет получить оценки качества функционирования системы ЭХЗ на этом интервале на основании оценок вида (14), вычисленных при x = x* .
Так как функция распределения потенциала Ut-3(t,x,w) на интервале |xi_i, xiJ достигает минимального значения только в точках xj_i и xj, оценки качества функционирования системы ЭХЗ будут иметь вид:
G : Vt Є [0,T] Л Vw єП,і = 1,n ;
P(Ut-3(xj, t, ш) < U^_3) > a , (16)
P(UT-3(xi, 1, w) > UT-3) ^ P . (17)
Критерием оценки качества функционирования системы ЭХЗ является ограничение вероятности выхода значения потенциала «труба-земля» за пределы верхней границы U^_3 = —0.85 [B] с вероятностью не ниже а = 0.98 ; нижней границы U^_3 = —1.15 [B] с вероятностью ниже, чем р = 0.95.
3. Математическая постановка задачи оперативного планирования режимов работы СКЗ подземных трубопроводов
Для рассматриваемого однониточного участка подземного трубопровода энергетические затраты на работу i-й СКЗ в момент времени t є [0, T] определяются мощностью по постоянному току на выходе преобразователя СКЗ Сэксп, [кВт ч]:
Ciэксп = Pi(t) ■ t , (18)
64
РИ, 2008, № 2
где мощность i-й СКЗ Pj(t), [кВт] определяется выражением:
pi(t) = ui«P (t)ii
(19)
здесь Uinp
Unp (t) = Ii
Ri л + Riпз + Zi
напряжение на
Ц л 1 л пз
выходе преобразователя i-й СКЗ [В] в момент времени t; Ii - величина силы тока на выходе i-й СКЗ [А] в момент времени t; Riл - сопротивление соединительных линий i-й СКЗ [Ом]; Ria3 - сопротивление анодного заземления i-й СКЗ [ Ом];
zi = yjRjпр ■ RiH3 ■ (лБ)-1 - характеристическое сопротивление трубопровода [Ом]; Ri пр - продольное сопротивление подземного трубопровода, [Ом/м] на
участке [xi_i, xiJ; Ri из - среднее значение переходного сопротивления изоляционного покрытия, [Омм2] на участке [xi_i, xi J.
Переходное сопротивление изоляции является случайной величиной Ri из = Ri из (t, ш), следовательно, и характеристическое сопротивление трубопровода представляет собой также случайную величину
Zi = Zi (t, ш) .
ного потенциала «труба-земля», которая должна удовлетворять критериям, определяющим область G.
Задача оперативного планирования режимов работы СКЗ в стохастической постановке будет иметь вид:
П
М Е Ji(T,ra) ^ mm (22)
ш i=i IeG v '
На основании показателей оценки качества функционирования системы ЭХЗ область допустимых решений G определяется системой неравенств
G : Vt Є [0, T] Л Vw Є П,і = 1, n ,
>|c P(Ut-3(xj, t, w) < U^_3) > a , (23)
P(UT-3(xi,tiw) > UT-З) ^ P , (24)
0 < Ii < Iimax, (25)
где I - вектор значений силы тока СКЗ I = [Ii, I2,_, 1П]; U+_3 = —0.85 - максимальное зна-
чение защитного потенциала [В]; Up_3 =—1.15 -минимальное значение защитного потенциала [В]; Imax - максимальные значения силы тока СКЗ [А].
4. Детерминированный эквивалент стохастической задачи оперативного планирования режимов работы СКЗ подземных трубопроводов
Как следует из (19), мощность i-й СКЗ Pi является в этом случае функцией от случайно величины Zj(t, ш) и, в свою очередь, также будет случайной величиной, т. е.
Pi(t,W) = I2[Riл + Ri пз + zi(t,w)]. (20)
В этом случае в качестве интегрального показателя оптимизации работы i-й СКЗ на интервалах времени [0, T] будем использовать выражение
T
Ji(T, ш) = J Pi(t, w)tdt =
0
T
= /I2[Riл + Ri„ + zi(t w)] tdt . (21)
0
При построении детерминированного эквивалента задачи стохастического программирования (22)-(25) будем предполагать, что случайные величины сопротивления изоляции R^ (w) и сопротивления грунта Ргр Е) распределены по нормальному закону распределения с известными параметрами, т.е.
R«s (w) = N(R«5, ^R„ ) , Ргр (w) = ^(Ргр, арр ) , где
Rm , Ргр - математические ожидания случайных
величин сопротивления изоляции и удельного сопро-
22
тивления грунта соответственно, cr , ар - дисперсия случайных величин сопротивления изоляции и удельного сопротивления грунта соответственно. На заданном интервале управления [0, T] параметры распределения случайных величин R из (w) и ргр Е) будем считать постоянными.
Численное значение Ji(T, ш) характеризует суммарные затраты мощности i-й СКЗ на интервале времени [0, T] и, как следует из (21), является случайной величиной.
Так как Ji (T, ш) - случайная величина, то для получения операциональной системы интегральных показателей эффективности функционирования i-й СКЗ необходимо использовать ее математическое ожидание.
Получили задачу стохастического программирования М-типа с вероятностными ограничениями. Целевая функция данной задачи является суммой энергетических затрат по эксплуатации каждой станции катодной защиты на промежутке времени [0, T]. Регулирование мощности СКЗ определяет величину защит-
Целевая функция задачи имеет вид:
П
M Е Ji(T,ra) ^ min
w i=1 IgG
После подстановки выражения для вычисления затрат мощности i-й СКЗ и выражения характеристического сопротивления участка трубопровода целевая функция задачи будет иметь вид:
n [T „
Е M j/I2[Ri л + Ri„ +
i=1 ш [ 0
^ ^Rinp ’ Riu? (^) ’ (^D) 1 ]tdt| ^ min . (26)
J IgG
После вычисления определенного интеграла выражение (26) можно представить в виде:
РИ, 2008, № 2
65
E M-flifR;л + Ri +
1=1 ы
+ 4Ri«p ■ RiH3 (w) ■ (^D)-1 ]-yJ^ • (27)
Для получения математического ожидания случайной величины воспользуемся неравенством Йенсена.
Для нелинейной функции f(x, гі(ш)) и случайной величины г)(ш) имеет место неравенство Йенсена:
Mf(x,^(w)) > f(Mf(x, ^(ш)))
если функция f(x, ^(w)) выпукла, и
Mf(x , ц(ш)) < f(Mf(x, ^(w)))
если функция f(x, ^(w)) вогнута.
алов, создаваемую каждой СКЗ, и является разницей отрицательного потенциала трубы и положительного потенциала грунта:
П
ит_з(x, t, ш) = Е UT(x,R „з ц t)) -
1=1
П
Е Егр (x, Ргр Eit))
1=1
Найдем распределение случайной величины отрицательного потенциала трубы UT(x* , R„3 (w, t)) i-й СКЗ в «диктующей точке» xJk при известном распределении случайной величины сопротивления изоляции R„3 (w) = N(R„3, <rR ). С помощью критерия Пирсона проверки статистической гипотезы о распределении случайной величины принимаем, что случайная величина UT(x*,R„3(w,t)) имеет нормальное распределе-
Построим график зависимости функции затрат мощности i-й СКЗ (21) от значения случайной величины
RU(Е (рис. 3).
ние:
UT(x; ,R„3(w,t)) = N(UT(x; ), ^U^(x;)) • (29)
Найдем распределение случайной величины положи-
тельного
потенциала грунта Up (x, ргр (w,t)) i-й СКЗ
«диктующей точке» xk при известном распределении
удельного сопротивления грунта Ргр (w) = N(P^ Как видно из выражения (12):
) •
игр (x, Ргр (w,t))
1 Ргр (w,t)
* ln
2:r(b; — a;) (x - a;) + 4(x - ai)2 + yA
(x - b;) + 4(x “ bi)2 + УіА
Рис. 3. График зависимости функции затрат мощности i-й СКЗ от значения случайной величины RiH3 (w)
Случайная величина (x, ргр (ш, t)) является линей-
ной функцией от нормально распределенного аргумента удельного сопротивления грунта т.е. является нормально распределенной случайной величиной:
в
Из графика зависимости видно, что функция, определяющая затраты при эксплуатации СКЗ, вогнута, т.е. имеет место неравенство:
Е M jl2 [Riл + Ri„ +4RnP ■ Ri„3 И ■ КТ1 ] T^J <
2 [Ri л + Ri„ +fi-R~(^T ] ■
<E
i=1
. (28)
Таким образом, оценка затрат при эксплуатации СКЗ, полученная в результате решения задачи стохастического программирования, будет больше реальных затрат.
Рассмотрим неравенства (23)-(25), определяющие область G задачи стохастического программирования. Для построения детерминированного эквивалента задачи необходимо нахождение функции распределения разности потенциалов «труба-земля» Ut-З(x, t, ш) в «диктующих точках».
Как уже рассматривалось ранее, потенциал «труба-земля» в точке х представляет собой сумму потенци-
Ер(xbРгр(w,t)) = ВДр(xk)!CTU;(x*k)). (30)
Таким образом, потенциал «труба-земля» в «диктующих точках» представляет собой сумму нормально распределенных случайных величин. При этом случайные величины сопротивления изоляции RU3 (ш) и удельного сопротивления грунта Ргр Е) являются независимыми случайными величинами:
иТ-З (x; 1ш) = N(UT-З (x; )i °Ut 3 (x*)) .
При известных параметрах математического ожидания и дисперсии случайной величины потенциала «труба-земля» в «диктующих точках» выражение (23) можно представить в виде:
>|с I
P(UT-3(xi Л » < uр_З) =
= р(ит-3(xi,t!w) - Ur-3 ^ 0 = (31)
Ut _3 (x*) - U+_3
a
UT-3 (xi )
> a .
66
РИ, 2008, № 2
Неравенство (31) представим в виде:
ф~1(а) ■ V3(x*) -“(UT(x*) - ^-3), (32) где Ф-1(а) - обратная функция Лапласа. Выражение (24) можно представить в виде:
P(ut-3(хіЛ w) > Uт-3) =
= p(ut-3(xi,t!ш) - ut-3 > 0) = (33)
= Ф
UT-3(xi) - Ur-3
> /3
Ut _3 (x-)
неравенство (33) представим в виде:
ф-V) ■
< ut_3 (xi) - U7_3 . (34)
Ut-3 (xi)
Таким образом, детерминированный эквивалент задачи (22)-(25) принимает вид:
Е
i=1
I? (t) [Riл + Ri„ Rinp ■ RU ■ (xD) 1 ] '
“ ss •<]---->+
T
сто пеш-иі
ITS у
I
t
Э V
Ї
min
IeG ,
РДул. ПЕрелслсеоап 11В
(35)
ОРПК,№рж.р-Н
t
№ ESS sy0 ЧП'БС-Сери
5t 1T2 1T1 1TD T 169 Ul65
Г Г “ЖЬ”Г
иОО ■АСС'ллЭоАоч.л Z33S І I I I
26П У' V у! '/ 1ST 498 у'
А
Т^Лсботш
X 377
Ч----IX
1
26І
К • "s'-
05 ^ и
СКЗ I 22ІІ I
ПРП 4-Залютио \/ Vn
РП ЗСсртіроБИ
і
ДсфЮДЕТ^КТЕДЛіА
1ТТ
“4>У-
_ рдп.Лвдное
-к--------[>
-<1 \* скз
325 у
4 881 ЛГНКС З.Леанм ЕЕ1
ПРП S , Н. Е-ЗЕіфІЯ
РД О АО 'Л ала-Xqi ыкс' 116+
G : i = 1, n : Ф 1(а) ■ а.
UT-3 (xi )
<
<-l Ur_3(xi) -Щ_3
Ф G' aUr 3(x) ^
п. Н.Высокий ZS7 255
*<]----->'4-
п. НщчьДВ67 257
*<]----->>
385
к. Отдьк
Z3rmi »<Л ’>..•>
J Т18
0 < Ii < Iim
33? Пс+итилдаа З
9TT CTicprttoa ІАіалллсг Э7Т
-к----------1>*
317
*<]------
Втсріерілт Беглолсек
ІШі
*<]—•"-• "В.ГВДгф" 1Ш43
(36)
(37)
(38)
Задача (35)-(38) относится к классу задач нелинейного математического программирования с ограничениями в виде нелинейных неравенств и двусторонней ограниченностью переменных. Для ее решения использован метод множителей Лагранжа.
5. Результаты решения задачи оперативного планирования режимов работы СКЗ подземных трубопроводов
Результаты решения задачи иллюстрируются на примере участка газораспределительных сетей высокого давления г. Харькова (рис. 4): длина участка 41510 м; количество СКЗ - 5; координаты расположения СКЗ от начала участка: 2000 м, 13690 м, 20900 м, 30063 м, 39513 м на участке имеется трубопровод различного диаметра: [0, 19950] D = 0,529 м;
[19950, 30050] D = 1,020 м; лгнкс2т.гач1»ч- [30050, 41510] D = 0,720 м;
А
^ изоляционное покрытие - заводское напы-
ленное, полимерное, усиленного типа, среднее значение переходного сопротивления изоляции RU3 = 5 ■ 103 [Ом м2]; продольное
сопротивление
Rпр = 6.6 ■ 10
-6
Рис. 4. Участок газораспределительных сетей высокого давления г. Харькова
.и ^
трубопровода [Ом/м]; сопротивление анодного заземления СКЗ R^ = 1.37 [Ом]; сопротивление соединительных линий СКЗ
67
n
для оптимальных значений силы тока СКЗ
Рис. 6. Распределение потенциала «труба-земля» для оптимальных значений силы тока СКЗ при различных реализациях случайных величин
R л = 0.23 [Ом]; среднее удельное сопротивление грунта по трассе составляет ргр = 40 [Ом.м]; анодное заземление - несмещенное, т.е. aj = xf —100 [м], bj — Xf А 100 [м] j = 1..5 ; длина цепочки анодного заземления LA = 200 [м] j = 1..5 ; расстояние до анодного заземления yA = 220 [м] j = 1..5 ; максимальная сила тока на выходе СКЗ 12.5 А; Umjn = —0.85 - минимальный защитный потенциал [В]; Umax = —1.15 - максимальный защитный потенциал [В].
Результатом решения задачи является оптимальное значение силы тока СКЗ на заданный период времени. Для описанных исходных данных получили такие значения силы тока СКЗ: І8к^[1]=7.032А,
^[2]=6.790А, Iskz[3]=7.154 А, ^[4]=6.650А, М^[5]=6.287А. На рис. 5 показан график распределения защитного потенциала «труба-земля» для полученных значений силы тока СКЗ трубопровода:
На рис. 6 изображены графики распределения потенциала «труба-земля» для оптимальных значений силы тока СКЗ при различных реализациях случайных величин сопротивления изоляции и удельного сопротивления грунта.
Выводы
В результате решения задачи оперативного планирования режима работы системы ЭХЗ были получены оптимальные значения силы тока СКЗ в условиях неполной информации о сопротивлении изоляционного покрытия трубопровода и удельного сопротивления грунта, показана эффективность решения данной задачи.
Научная новизна: на примере участка газораспределительных сетей высокого давления г. Харькова впервые предложены стохастические модели и методы решения задачи оперативного планирования работы системы ЭХЗ подземных трубопроводов.
Практическая значимость: результаты решения задачи будут использованы в системе анализа и оперативного управления режимами работы системой ЭХЗ подземных трубопроводов, разработка которой направлена на повышение эффективности работы и сокращение затрат на эксплуатацию систем ЭХЗ подземных трубопроводов.
Литература: 1. Стратегія розвитку нафтогазового комплексу України на період до 2030р. та дальшу перспективу. 2. Юдин Д.Б. Математические методы управления в условиях неполной информации. М.: Сов.радио, 1974. 3. Лорткипа-нидзе Б.Г. Цепи постоянного тока с утечкой и электрическая зашита от коррозии. Изд-во АН Груз. ССР, Тбилиси, 1962. 4. Тевяшев А.Д., Ткаченко В.Ф., Попов А.В., Стрижак Л.В. Стохастический подход к постановке и решению задачи оперативного планирования режима работы системы электрохимической защиты подземных трубопроводов // Восточно-Европейский журнал передовых технологий. 2005. № 3/2(15). С. 94-98. 5. Тевяшев А.Д., Ткаченко В.Ф., Попов А.В. Система интегральных показателей качества и эффективности функционирования систем электрохимзащиты подземных трубопроводов // Восточно-Европейский журнал передовых технологий. 2006. № 1/2(19). С. 21-26.
Поступила в редколлегию 11.06.2008
Рецензент: д-р физ.-мат. наук, проф. Дорошенко В. А.
Тевяшев Андрей Дмитриевич, д-р техн. наук, профессор, заведующий кафедрой ПМ ХНУРЭ. Научные интересы: проблемы моделирования и оптимизации технологических процессов в системах энергетики. Адрес: Украина, 61166, Харьков, пр. Ленина, 14, тел. (057) 702-14-36, е-mail: tevjashev@kture.kharkov.ua.
Стрижак Людмила Викторовна, инженер кафедры прикладной математики ХНУРЭ. Научные интересы: проблемы электрохимической защиты трубопроводов. Адрес: Украина, 61166, Харьков, пр. Ленина, 14, тел. (057) 702-16-59, е-mail: lvst@mail.ru.
Ткаченко Александр Анатольевич, начальник отдела организации продаж услуг ОАО «Укртелеком», Харьковский филиал. Научные интересы: управление в технических системах. Адрес: Украина, 61005, Харьков, пл. Свободы, Госпром, 3 подъезд, 3 этаж.
68
РИ, 2008, № 2