Способ обнаружения предвестников чрезвычайных ситуаций на линейной части подземного магистрального продуктопровода

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

 

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

Известен способ обнаружения изменения состояния участка трубопровода по появлению виброакустического сигнала в оболочке трубы, формирующегося при установке в нее врезки [Защита трубопроводов от несанкционированных врезок / А.А. Казаков // Системы безопасности. - 2008. - №5. - С.150-154]. Недостатком этого способа является запоздалое появление предупреждающего сигнала, исключающее возможность предотвратить нарушение целостности трубы злоумышленниками и сопутствующую этому чрезвычайную ситуацию.

Известен способ обнаружения «аварийно-опасного» участка трубопровода, основанный на возбуждении ударных виброакустических импульсов в оболочке трубы с помощью приваренных к ней звукопроводящих стержней с последующим определением отношения резонансной частоты диагностируемого трубопровода к эталонной [Пат. 2350833 РФ, МПК F17D 5/00. Способ контроля и диагностики состояния трубопровода [Текст] / Толстунов С.А., Мозер С.П., Толстунов А.С]. В основу способа положено известное соотношение зависимости резонансной частоты пластины ƒ₀ от ее толщины h:ƒ₀=c/2h, c - скорость распространения продольных волн в трубопроводе. Использование этой закономерности для выявления земляных работ в охраняемой зоне не представляется возможным.

Известна заявка №2006137406/28 от 23.10.2006 (дата публикации 27.04.2008) на способ и устройство дальнего обнаружения утечек в трубопроводе. Согласно заявке, в перекачиваемом продукте создаются периодические волны давления, которые регистрируются на другом конце контролируемого участка. По искажению регистрируемой волны судят о наличии утечки на этом участке. Способ не позволяет фиксировать изменения, происходящие за пределами оболочки трубопровода.

Известен способ обнаружения утечек на трубопроводном транспорте углеводородов, основанный на регистрации и анализе инфразвуковых сигналов в перекачиваемом продукте [Пат. US 666861982. Pattern matching for real time leak detection and location in pipelines (Распознавание образов для детектирования в реальном времени факта и локализации врезок в трубопроводы) и в различных вариантах исполнения описанный в http://acoustic-solution-intl.com/faq_index.htm; www.grouplb.com; http://torinsk.ru/publication/25-mpp2007.html и др. Недостаток способа - регистрируется факт нарушения целостности трубопровода, а не подготовительные работы по установке врезки или боеприпаса.

Из известных технических решений наиболее близким по совокупности существенных признаков к заявляемому является способ обнаружения изменений параметров среды в окружении заглубленного магистрального трубопровода [Пат. 2463590 РФ МПК G01N 29/04. Способ обнаружения изменений параметров среды в окружении заглубленного магистрального трубопровода/Епифанцев Б.Н., Федотов А.А.].

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

Целью изобретения является исключение процесса создания физических эталонов предвестников ЧС на заглубленных продуктопроводах при сохранении надежности обнаружения на необходимом уровне в любое время года.

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

$${\epsilon }_j\ge {\displaystyle \underset{\pm N_j\Delta C_j}{\overset{\pm \infty }{\int }}P\left(\sqrt{N_j{\sigma }_{\Delta E}}\right)d\Delta E}$$

по задаваемой вероятности ложных решений ε_j для j-го предвестника ЧС, расчетной оценке ожидаемого сигнала ΔC_j в точках регистрации, среднеквадратичному отклонению регистрируемых отсчетов указанных интегралов ΔE, т.е. σ_ΔE, в качестве эталонов используют уровни ±N_j ΔC_j, а решение о появлении j-го предвестника ЧС принимают при превышении накопленного за цикл результата установленных для него эталонных уровней.

Сущность изобретения поясняется нижеследующим описанием и прилагаемыми к нему чертежами. На фиг.1 представлен вариант структурной схемы, реализующей предлагаемый способ. Фиг.2 поясняет алгоритм определения числа накоплений интегралов разностей регистрируемых сигналов.

Обозначения на фигурах: 1 - изображение продуктопровода; 2 - окружающий продуктопровод грунт; 3, 4 - датчики виброакустических колебаний в оболочке трубы; 5 - преобразователь электрических сигналов в упругие колебания в заданной точке трубы; 6 - ключ, обеспечивающий пропускание входного сигнала на выход в течение действия управляющего сигнала по другому входу; 7 - схема вычитания; 8 - линия задержки; 9 - генератор «прозванивающих» импульсов; 10 - интегратор входного сигнала в течение существования «прозванивающего» импульса; 11 - вычислитель среднеквадратичного отклонения проинтегрированных разностей сигналов от датчиков 3, 4; 12 - микроконтроллер; 13, 14 - накопители входных сигналов; 15, 16, 17, 18 - решающие устройства; 19, 20 - перемножители; 21 - громкоговоритель; A, B - смежные участки продуктопровода, C - устройство принятия решений; ΔC_y, ΔC_вр - приращения регистрируемого сигнала при появлении на участке продуктопровода соответственно утечки или изменений, связанных с подготовкой врезки; ε_y, ε_вр - задаваемые (приемлемые) вероятности ложных решений, обеспечиваемые устройством, N_j, - необходимое число накоплений для обеспечения обнаружения соответствующего предвестника ЧС с вероятностью ложных решений ε_j, «Син» - синхронизация работы (сброс, пуск) устройства принятия решений C, P(…) - плотности вероятностей проинтегрированных разностей сигналов от датчиков ΔE(t_i) в моменты t_i.

Известно, что амплитуда сигнала на выходе датчиков 3, 4 содержит информацию о параметрах перекачиваемого продукта, характеристиках трубы и окружающей его среды [Буденков С.А. и др. Оценка возможностей метода акустической эмиссии при контроле магистральных трубопроводов // Дефектоскопия, 2000, №2, С.29-36]. Переизлучение упругих колебаний оболочки в жидкость определяется ее плотностью, турбулентная среда, которая представляет собой перекачиваемый продукт, порождает флуктуации амплитуды регистрируемых сигналов. Корни деревьев в ветреную погоду являются источниками упругих колебаний в почве, проникающими из окружающей среды в оболочку продуктопровода. Эти и ряд других источников неинформативных сигналов (шумов) определяют дисперсию амплитуд регистрируемых сигналов, поступающих на схему вычитания 7.

Импульсы от генератора 9 преобразуются в упругие колебания элементом 5, регистрируются датчиками 3, 4, вычитаются один из другого схемой 7, поступают через ключ 6 на интегратор 10, отсчеты получаемых интегралов разностей указанных импульсов поступают на схемы 11, 13, 14. Поскольку при распространении упругих колебаний по продуктопроводу происходит их задержка относительно возбуждаемых импульсов генератора 9, введена линия задержки 8 для управления ключом 6. Интегрирование получаемых разностей увеличивает отношение сигнал/шум, подавляются высокочастотные составляющие шумов.

На схему 11 возложена функция вычисления среднеквадратичного отклонения поступающих отсчетов σ_N:

$${\sigma }_N=\sqrt{\frac{1}{N}{\displaystyle \sum _{t=1}^N{({\displaystyle \underset{t_i}{\overset{t_i+\tau }{\int }}\Delta U(t)dt-m_N^{*}})}^2}},{\displaystyle \underset{t_i}{\overset{t_i+\tau }{\int }}\Delta U\left(t\right)dt=\Delta E\left(t_i\right)},$$

где N - число поступивших отсчетов в схему 11, ΔU(t) - сигнал на выходе схемы 7 в момент существования импульса на выходе схемы 8, τ - длительность импульса задающего генератора 9, $$m_N^{*}$$ - оценка математического ожидания ΔE(t_i), которое при отсутствии предвестников ЧС устанавливается равным 0. Число N определяется по задаваемой точности определения σ_N схемой 11 с учетом количества импульсов, поступающих на ее вход со схемы 8.

По завершению формирования значения σ_N включается в работу микроконтроллер 12. Его задача по заданным ε_y, ε_вр, ΔC_y, ΔC_вр и σ_N определить значения чисел накапливаемых сигналов ΔE(t_i) для обнаружения соответствующих предвестников ЧС: «утечки» или «подготовка к установке врезки». В общем случае предвестников может быть больше, устройство на фиг.1 отражает случаи обнаружения обозначенных отклонений от нормы на одном из контролируемых участков A, B.

Алгоритм работы микроконтроллера 12 можно пояснить с помощью изображенных на фиг.2 плотностей вероятностей сигналов ΔE(t_i). При отсутствии их накопления (N=1) плотность вероятностей имеет вид Р(ΔЕ(t_i), N=1). Шумы, порождаемые турбулентными пульсациями давления в перекачиваемом потоке, по своей природе имеют нормальный закон распределения вероятностей и при $$m_N^{*}\approx 0$$ , характеризуются на выходе схемы 11 только дисперсией $${\sigma }_N^2$$ . Накопление независимых отсчетов $${\displaystyle \sum _{i=1}^N\Delta E\left(t_i\right)}$$ увеличивает (согласно положениям теории вероятностей) среднеквадратическое отклонение обозначенной усредненной суммы в $$\sqrt{N}$$ раз (см. фиг.2). В то же время накопление сигнала ΔC_y·/(ΔC_вр) увеличивает итоговый сигнал в N раз. Поэтому увеличивая N, можно в принципе выявить сколь угодно малый полезный сигнал ΔC_j. Это известный вывод теории обнаружения сигналов. Применительно к рассматриваемой задаче не ясно, на каком количестве накоплений следует остановиться для ее решения.

С физической точки зрения формирование шурфа на смежном участке (А или В) и подготовка обнажившейся трубы для установки приспособления для создания врезки ведет к резкому сокращению переизлучения распространяющихся колебаний по продуктопроводу в этой зоне и появлению положительного приращения ΔC_вр в сечении регистрации (формулы для расчета этого приращения приводятся в [Меркулов Л.Г. Затухание нормальных волн в пластинах, находящихся в жидкости//Акустический журнал, 1964, т.Х, вып.2, С.206-211; Буденков Г.А. и др. Оценка возможностей метода акустической эмиссии при контроле магистральных трубопроводов//Дефектоскопия, 2000, №2, С.29-36]). При утечке окружающий продуктопровод грунт будет пропитываться вытекающим продуктом. В результате переизлучение распространяющейся энергии упругих колебаний увеличивается, сигнал на выходе системы регистрации уменьшается (-ΔC_y). Существующие расчетные методики позволяют оценить значения ΔC_вр и ΔC_y. Но есть проблема обнаружить эти отклонения в условиях, когда отношение сигнал/шум на входе датчиков оказывается существенно меньше единицы. Единственно видимый путь решить обозначенную задачу - провести накопление малых отклонений ΔC_вр, ΔC_y. По условиям задачи допускается проведение нескольких тысяч накоплений и, следовательно, можно увеличить исходное отношение сигнал/шум до 50-60 раз. Противоположное требование: чем больше времени будет затрачено на реализацию операции накопления, тем больший объем грунта будет загрязнен вытекающим продуктом, т.е. интервал накопления следует уменьшить.

Компромиссный вариант решения задачи просматривается при анализе фиг.2.

При решении задач обнаружения (распознавания) задаются приемлемыми для потребителя вероятностями ложных решений ε_вр (ε_y) и пропуска цели. При накоплении сигнал увеличивается по закону NΔC_вр(NΔC_y), среднеквадратическое отклонение шумов - по закону $$\sqrt{N}{\sigma }_N$$ . Вероятность ложной тревоги

$$$$ $${\epsilon }_j\ge {\displaystyle \underset{\pm N_j\Delta C_j}{\overset{\pm \infty }{\int }}P\left(\sqrt{N_j{\sigma }_N}\right)d\Delta E}$$ .

При знании оценки энергии шумов $${\sigma }_N^2$$ (определяется экспериментально при каждом цикле принятия решений схемы 11), заданной вероятности ложных решений ε_j (например, 0,01) и расчетной оценке сигнала ΔC_j, от предвестника ЧС при нормальном распределении $$P\left(\sqrt{N\cdot {\sigma }_N}\right)$$ последовательное задание числа N_j позволяет определить его величину при решении приведенного уравнения. Эта функция возложена на микроконтроллер 12.

Найденные значения N₁ и N₂ поступают на устройство принятия решений с одновременным запуском его в рабочее состояние (сигнал «Син»). Соответствующие накопители 13 и 14 производят суммирование N₁ и N₂ входных отсчетов, при превышении этими суммами порогов N₁ΔC_y и N₂ΔC_вр решающие устройства 15, 16, 17, 18 формируют сигналы «Решение».

Пороги принятия решений определяются числами N₁ и N₂, которые вычисляются по текущей дисперсии $${\sigma }_N^2$$ . Изменения режима перекачки, усиление ветра меняют показатель σ_N, а через него N_1, N₂. В результате поддерживается заданная вероятность ложных решений при обнаружении исходных изменений на продуктопроводе не ниже заданных при любых отклонениях параметров системы перекачки продуктов и погодных условий.

Способ обнаружения предвестников чрезвычайных ситуаций на линейной части подземных магистральных продуктопроводах, основанный на возбуждении в выбранном сечении оболочки трубы прозванивающих импульсов упругих колебаний, регистрации их на удалении от сечения возбуждения, накоплении определенного числа зарегистрированных импульсов в очередном цикле последовательно принимаемых решений, формировании эталонов предвестников чрезвычайных ситуаций и принятии решений по результатам сравнения результатов накопления с эталонами, отличающийся тем, что регистрацию импульсов проводят в двух сечениях оболочки продуктопровода, удаленных на одинаковые расстояния по обе стороны от сечения их возбуждения, результат накопления получают в виде суммы отсчетов интегралов от разностей этих импульсов, причем число накоплений N_j в цикле определяют в результате решения уравнения
,
где ε_j - вероятность ложных решений для j-го предвестника чрезвычайной ситуации, N_j - число накоплений для j-го предвестника чрезвычайной ситуаций; ΔC_j - расчетная оценка уровня ожидаемого сигнала в точках регистрации, - нормальная плотность вероятностей с нулевым математическим ожиданием, σ_ΔE - среднеквадратичное отклонение регистрируемых отсчетов указанных интегралов ΔE, в качестве эталонов используют уровни ±N_jΔC_j, а решение о появлении j-го предвестника чрезвычайной ситуации принимают при превышении накопленного за цикл результата установленных для него эталонных уровней, причем решение «Врезка» трансформируют в сигнал тревоги через установленный на каждом участке громкоговоритель, а сигналы всех принимаемых решений передаются на мнемосхему в службу безопасности по каналам телемеханики.