CC BY
141
УДК 621.644.8
Е. В. Врагова 1, Л. А. Скляров 2
1 Институт почвоведения и агрохимии СО РАН ул. Советская, 18, Новосибирск, 630099, Россия E-mail: vragovae@rambler.ru
2 Новосибирский филиал Государственной академии профессиональной подготовки и повышения квалификации руководящих работников и специалистов инвестиционной сферы Комсомольский пр., 4, Новосибирск, 630004, Россия
E-mail: courses@nsunet.ru
ОБНАРУЖЕНИЯ УТЕЧЕК ГАЗА ИЗ МАГИСТРАЛЬНЫХ ГАЗОПРОВОДОВ В ТЕПЛОВОМ ПОЛЕ ИЗЛУЧЕНИЯ ЗЕМНОЙ ПОВЕРХНОСТИ
Представлен анализ новых методов расчета проявления утечек газа из трубопровода в тепловом поле излучения земной поверхности. Показано, что даже малые трещины в трубах приводят к появлению на поверхности грунта тепловых аномалий, позволяющих локализовать место утечки беспилотным летательным аппаратом с установленным на борту ИК радиометром.
Ключевые слова: моделирование, аварии, трубопроводный транспорт, температурные поля, тепловизор.
Постановка проблемы
Территория России характеризуется разветвленной газотранспортной сетью (рис. 1). Эти магистрали являются источниками повышенной опасности как с точки зрения экологии, так и возможного создания чрезвычайной ситуации при возникновении аварий на газопроводах. Предвестниками таких опасностей могут являться утечки газа, возникающие из-за механических повреждений трубопровода или его коррозии. Раннее обнаружение утечек является актуальной задачей. Возникает острая необходимость в разработке информационной системы, способствующей скорейшему обнаружению и устранению утечек газа. Прежде чем представить такую систему, предлагаем вашему вниманию настоящую работу, посвященную разработке инструментального метода обнаружения утечек газа из трубопроводов, который бы решал одну из основных проблем диагностики их состояния.
Анализ последних исследований и публикаций
Основным методом обнаружения мест утечек в настоящее время является патрулирование газопроводов с применением переносных газоанализаторов. Этот метод связан с большими трудозатратами и отличается невысокой оперативностью. Более оперативными и перспективными представляются методы обнаружения утечек, основанные на анализе тепловых изображений земной поверхности в районе залегания трубопроводов [1; 2], тепловизорами, установленными на беспилотных летательных аппаратах. Принцип тепловизионной диагностики состоит в сравнении эталонного и анализируемого полей температуры. Аномалии температуры служат индикаторами дефектов, а величина температурных сигналов и их поведение во времени лежат в основе оценок параметров объектов. Над диагностируемой трассой запускают беспилотный летательный аппарат, оборудованный системой, состоящей из видеокамеры, тепловизора и системы спутниковой навигации. Беспилотный летательный аппа-
ISSN 1818-7900. Вестник НГУ. Серия: Информационные технологии. 2009. Том 7, выпуск 4 © Е. В. Врагова, Л. А. Скляров, 2009
рат, перемещаясь по заданному маршруту, осуществляет запись и передачу на Землю цифровой диагностической информации в виде последовательности привязанных к карте изображений. Стационарный комплекс на Земле в режиме реального времени осуществляет сшивку изображений в ленту, а также обработку и анализ полученной цифровой информации с целью выявления мест разрывов и утечек.
Рис. 1. Разветвленная система транспортировки газа от добычи до потребителя
Постановка задачи и ее решение
Метод основан на обнаружении тепловых аномалий на поверхности почвы, вызываемых резким перепадом давления газа при выходе из канала утечки и появлением локальных тепловых градиентов в слое почвы над трубопроводом. Рассмотрим задачу оценки влияния этого процесса на тепловое поле в грунте.
Для систем теплоснабжения, паропроводов и газопроводов трубопроводы, чаще всего, заложены в грунте - массиве. Объем неограниченного массива имеет коэффициент теплопроводности X и бесконечные размеры по всем направлениям. Трубопроводы можно рассматривать как действующие сосредоточенные линейные, положительные и отрицательные источники теплоты длиною Ь (рис. 2).
Рис. 2. Расположение источника и стока теплоты: (+Q) - положительные источники теплоты; (-0 - отрицательные «стоки» теплоты; К - место совмещения изотерм
По сравнению с протяженностью Ь неограниченного массива источник (+О) и «сток» (-О) расположены близко друг к другу, и поэтому путь между ними является линией наименьшего термического сопротивления. В связи с этим вся теплота, выделяемая источником, будет полностью поглощаться «стоком». Если бы источник и «сток» не воздействовали друг на друга, то каждый из них создавал бы в теле массива температурное поле в виде концентрических изотерм. Тогда
пф/-тъ),{ -т;)
(+О) =
1 2к'
; (-О) =
1п-
2к
1п %
г
1п %
откуда Т3 = Т + О-—.
31 2пЬк
Если бы температурное поле в массиве формировалось только источником или только «стоком», то в точке К была бы температура Г3 или же Т3 . При одновременном действии
источников и «стоков» результирующее температурное поле получается путем сложения температурных полей, возбуждаемых в массиве отдельными источниками и «стоками» в предположении, что они не мешают друг другу. Если обозначить температуру источника теплоты через Г0 , а «стока» - Т0, то действительная температура в точке К определится из выражения
Т = Т + Т - Т = Т + Т -
1К 13^13 10 11 Т11
О
2пЬк
1п
Г г" ^
V Г" г3у
Для расчетов условно полагают, что теплота от цилиндрического газопровода (+О), заложенного в грунт, передается не в окружающую среду, а забирается отрицательным источником теплоты (-О) (рис. 3).
Рис. 3. Расположение газопровода: Ь - длина цилиндра с радиусом Я; Н - глубина заложения; X - коэффициент теплопроводности массива; Т0 - температура поверхности слоя
В этом случае «сток» (-О) размещен симметрично источнику (+О) и окружен точно таким же массивом, а реальная картина температурного поля в грунте при этом не нарушается. Учитывая, что
(+О) = пЬ(ТК - Т0) ; (_д) = пЬ(Т00 - Т'к ) ,
— 1п А ' —1п к '
2Х г 2Х "
и используя принцип совмещения температурных полей, можно выразить температуры в любой точке К грунта:
О к ,, ,, О г" ТК = 1п —; Тк = 1п—.
2п1Х г 2пЬ'к к
Суммируя температуры ТК и Тк , получим:
О г'
Т = Т +1п —
1к ^о
2кЬХ
или
Т(х,у) =То +-О-
2кЬХ х2 + (к - у)2
Таким образом, температурное поле в грунте становится определенным, если замерена температура на поверхности массива То и известен поток теплоты О. С другой стороны, тепловой поток цилиндром может быть рассчитан по температуре поверхности массива и еще по одной, любой, температуре в грунте:
О = кЬ[Т(х,у)-То] пЬ[Тт - То]
1 Ух2 + (к + у)^±1пГ2 к - IV 2^1П х2 + (к - у)2 Я )
Естественно ожидать, что максимальная температура грунта в точке п (х = 0, у = к + Я) равна
Тт = Т0 + 1пГ 1 + 2 к ] . т 0 2лЬЛ, ^ Я)
Если в массиве будет заложено два теплопровода, то соответственно им появится и два стока. Температурное поле в грунте получится уже как результат совмещения четырех температурных полей.
Расчет температурного поля методом релаксации. Релаксация - процесс установления термодинамического равновесия или восстановления неустойчивого положения в физической системе. Метод релаксации используется для решения задач стационарной теплопроводности в телах сложной конфигурации, когда при расчете температурного поля дифференциальное уравнение теплопроводности не поддается аналитическому решению.
Расчет температурного поля методом релаксации удобно иллюстрировать на примере, когда теплота распространяется в двух измерениях. Сечение тела (рис. 4, 5) обычно разбивается релаксационной решеткой на ячейки квадратной формы (Ах = Ду). Тело имеет глубину Ь и коэффициент теплопроводности X.
В дальнейшем допускается следующее:
1) процесс теплопроводности концентрируется в стержнях релаксационной решетки, и чем меньше будут размеры ячейки, тем выше точность вычислений, но количество расчетов при этом увеличится;
2) по каждому стержню релаксационной решетки передается в точности такое же количество теплоты, которое в действительности передается через элемент с размерами Ах и Ду;
3) в качестве расчетного соотношения может быть использована формула расчета теплого потока через плоскую стенку:
г
к
Ок =7 Р (Т - Ч),
о
где 1 - отмечает наибольшее значение температуры; к - наименьшее значение температуры.
Рис. 4. Ячейка релаксационной решетки
Рис. 5. Проекция ячейки релаксации
Приняв эти допущения, можно рассчитать количество теплоты, которое протекает по каждому стержню:
• для горизонтального стержня 5 = Ах, Р = АуЬ
к
Ок =—АуЬ(Тг - 1к);
Ах
• для вертикального стержня 5 = Ау, Р = АхЬ
к
Ок = -к АхЬ(Т - Ч); Ау
• для любого стержня
Ок =кЬ(Тг - ). (1)
При этом могут иметь место различные схемы подвода и прохождения теплоты через ячейку релаксации.
1. Теплота идет от точки 1 к узловой точке 0, а от точки 0 расходится в направлении точек 2, 3, 4. Уравнение баланса теплоты имеет вид О1^0 = О0^2 + О0^3 + О0^4. Согласно (1)
кьТ -Т)=кь(Т0 -г2)+кь(Т0 -Т^)+кьТ -Т4),
откуда T0 =1 (T + T2 + T3 + T4 ).
2. Теплота идет от точек 1 и 2 к точке о, а от точки о расходится к точкам 3 и 4. Уравнение баланса имеет вид О1^0 +О +О о^4, или
ХЬ (Т - То)+и (Т - То) = ХЬ (То - Тз) + ХЬ (То - Т4),
откуда T0 =1 (T + T2 + T3 + T4 ).
3. Теплота идет от точек 1, 2, 3, 4 к точке 0 и там взаимно уничтожается. Уравнение баланса тепла имеет вид Q^0 + Q2^0 + Q3^0 + Q4^0 = 0, или
XL(T^-T0H^T -T0) + XL(T3 -T0) + ХЬ(Т4 -T0) = 0, откуда T0 =1 (T + T2 + T3 + T4 ).
Следовательно, какова бы ни была схема прохождения теплоты, температура в узловой точке квадратной релаксационной ячейки всегда будет равна среднему арифметическому значению из температур, окружающих эту точку.
Для каждой узловой точки релаксационной решетки существует закон релаксации:
T + T + T + T
AP = T + 7 + 7 + 7 -T0 = 0. (2)
4 0
Принцип релаксации заключается в следующем.
1. Тело сложной конфигурации разбивается на релаксационную решетку. Исходя из предварительного объема знаний приближенно задают значения температур в узлах решетки.
2. Эти приближения проверяются от точки к точке в соответствии с требованием закона релаксации, и устанавливают точки, в которых наблюдается наибольшее отклонение от закона релаксации (2).
3. Начиная с мест (точек) наибольшего отклонения вносятся поправки А, чтобы удовлетворить закону релаксации.
4. Эти исправления в свою очередь вызывают новые отклонения в соседних точках решетки, и возникает необходимость в повторной коррекции.
5. Повторная коррекция вносится каждый раз последовательно в порядке убывающих отклонений и производится до тех пор, пока численные значения температур во всех точках сетки не придут во взаимное соответствие, т. е. везде будет соблюдаться закон релаксации.
Записав уравнение (2) для каждого из узлов тепловой сетки, получим систему, состоящую из числа линейных уравнений, равного числу узлов сетки. Для решения этой системы уравнений применяют различные численные методы, и в частности метод релаксации. Название метода происходит от латинского relaxation 'ослабление', означающего постепенный переход системы в равновесное состояние. Например, если температура в каком-либо узле сетки, зависящая от четырех соседних значений температуры, находится в равновесии с ними, то выполняется уравнение (2). Если она не находится в равновесии с соседними значениями температуры, то правая часть этого уравнения не будет равна нулю, т. е.
T + T2 + T3 + T4 - 4T0 = AT, (3)
где AT - остаток.
Для сведения к нулю правой части (остатка) каждого из уравнений системы, т. е. для приведения системы в равновесное состояние, и применяется этот метод.
Как уже сказано, глубинный профиль температуры в грунте, связанный процессами солнечного нагрева почвы и теплообмена поверхности с окружающей средой, может быть рассчитан в соответствии с законом Фурье [4].
л
Qtp = — FT (ti -13), (4)
о,-,-
где X - коэффициент теплопроводности грунта; 5- - расстояние между точками (¡^ и ;
FT - площадь теплообмена; ti, tj - температуры грунта в точках (¡^ и ^ — соответственно.
Наличие в деятельном слое грунта трубопровода приводит к нарушениям процесса теплопередачи [5] и вызывает искажения глубинного профиля этих изотерм. Этот процесс усиливается при отличии температуры газа в трубе от температуры грунта, что характерно для реальных условий прокачки газа при его транспортировке.
Для определения поля температур использован метод релаксаций. Для этого расчетная область с учетом симметрии относительно вертикальной оси, проходящей через центр трубы, разбивается сеткой на элементарные объемы. Для выбранной точки каждого объема составляется уравнение баланса тепла, позволяющее определить значение температуры.
Количество тепла, передаваемого от газа, движущегося в трубе, к грунту определялось по
формуле
где к =
О = kF (t -1 ),
(5)
1
аг
- коэффициент теплоотдачи; а - коэффициент теплоотдачи от газа к
+ Я.
стенкам трубопровода; ЯЕ - суммарное термическое сопротивление стенки трубопровода и участка грунта между внешней поверхностью трубы и расчетной точкой (¡^ в грунте; Fg -теплоотдающая поверхность трубопровода; tg , ti - температуры движущегося по трубопроводу газа и точки {¡^ грунта соответственно.
При возникновении утечки из газопровода массовая скорость истечения газа определяется отношением атмосферного давления и давления в трубе.
Так, если перепад давлений атмосферного воздуха Ратм (Па) и давления в газопроводе Ргаза (Па) удовлетворяет неравенству:
Р,
Р.
<
' 2 ^
у +1
т
У-1
то в начальный момент времени истечение газа происходит со звуковой скоростью (такой режим наблюдается в случае перекачки газа под давлением более 0,18 МПа). Здесь у = ср / су - отношение теплоемкостей.
Массовая скорость звукового истечения газа определяется как
^ =
(
2
у+1
^у-Т Р
У + 1
кг/с,
где
Ргаза - давление в трубопроводе (Па); Т - температура в газопроводе (К);
Я = 8314 - универсальная газовая постоянная (кг-м2/К-кмоль-с2); М - молярная масса (кг/кмоль).
Чем больше перепад давления, тем больше скорость и массовый расход газа при истечении, и скорость истечения становится равной так называемой критической скорости, равной местной скорости звука
= , к-,
Я V м
где к - величина показателя адиабаты газа.
Для газа в модельном случае возьмем метан с М = 16, примем к = 1, и тогда wg = 22,8 кг/с.
1
У
Расход G = w S ,
^r^ g трещ '
где S - площадь трещины в газопроводе (м2).
Количество тепла, отводимого выходящим через отверстие газом, может быть определено, согласно [3], по формуле
Qi = Ggcgt, (6)
где cg - удельная теплоемкость газа Дж/(кг*К).; t - температура газа (К).
Утечка в расчетах имитировалась стоком тепла, помещенным в точку Qi , расположенную
в нижней части трубы. Интенсивность стока из отверстия размерами (1х0,1) мм2 в соответствии с [3] оценена как Qi = -33,15 Вт для температуры газа +24 °С и давления в трубе 70 атм.
Коэффициент теплопроводности грунта в расчетах принимался равным 1,49 Вт*м-1 "град-1. Температура воздуха в расчетах принималась равной te = 20 °С. Расположение стока тепла в верхней части увеличит «тепловой провал», соответствующий месту утечки.
Пространственное распределение теплового поля на поверхности может быть получено путем вращения рисунка результирующего распределения температур вокруг оси, проходящей через центр трубопровода.
Программный комплекс разработан на языке Java с помощью NetBeans (c).
Выводы
Оценки вероятности обнаружения, проведенные по методике, представленной в [5] для беспилотных летательных аппаратов с ИК аппаратурой (например, GasFindIR), показывают возможность правильной идентификации утечек с высот, характерных для патрульных полетов 200-5оо м. С помощью разработанного комплекса программ были рассчитаны типичные варианты возможных повреждений газопроводов. Полученные выводы на основе анализа численных расчетов и атлас тепловых полей могут быть использованы при проведении реальных тепловизионных съемок трубопроводов.
Список литературы
1. Викторов В. А., Лункин Б. В., Совлуков А. С. Радиоволновые измерения параметров технологических процессов. М.: Энергоатомиздат, 1989. 208 с.
2. Неразрушающий контроль. Россия, 1900-2000 гг.: Справочник / В. В. Клюев, Ф. Р. Со-снин, С. В. Румянцев и др.; Под ред. В. В. Клюева. М.: Машиностроение, 2001. 616 с.
3. ВулисЛ. А. Термодинамика газовых потоков. М.; Л.: ГЭИ, 1950. 304 с.
4. Исаченко В. П., Осипова В. А., Сукомел А. С. Теплопередача. М.; Л.: Энергия,1965. 424 с.
5. Павлов Н. И., Эльц Е. Э. Обнаружение температурных аномалий, обусловленных заглубленными в грунт инородными объектами // Изв. вузов. Радиоэлектроника. 2007. Т. 5о, № 3. С. 12-2о.
Материал поступил в редколлегию 26.10.2009
E. V. Vragova, L. A. Sklayrov
DETECTION OF LEAKS GAZA FROM THE MAIN GAS PIPELINES IN THE THERMAL FIELD OF RADIATION OF THE TERRESTRIAL SURFACE
Analysis of existing methods and a new methodic of calculation of leaks of gas from the pipeline in a thermal field of radiation of a terrestrial surface is considered. It is shown, what even small cracks in pipes lead to occurrence on a surface of a ground of the thermal anomalies, allowing localizing a leak place a pilotless flying machine with established onboard infra-red radio metre.
Keywords: modeling, accident, pipeline transport, temperature fields, teplovizor.
