На основе анализа аварийности, технологических нарушений при транспортировке газа и аналитического обзора вопросов технического диагностирования в работе [3] поставлена задача повышения эффективности диагностирования газораспределительных сетей на предмет возникновения неучтенного газа (утечек) при его транспортировке. На основании данного аспекта выявлена необходимость определения значимости факторов [2], влияющих на возникновение утечек природного газа.
Постановка задачи. Цель работы – организация и проведение экспертного опроса [4] для выявления значимости факторов возникновения утечек и неучтенного газа, а также оценка меры согласованности мнений экспертов. Интерпретацию результатов требуется представить в виде теоретико-множественной и ситуационно-фреймовой моделях.
На основании природы возникновения, в результате системного анализа процесса транспортировки газа и изучения измерений расхода, количества и учета природного газа [3], все факторы были условно разделены на три группы: утечки газа из газопроводов и сооружений на них; утечки, зависящие от физико-технических свойств газа и окружающей среды, неучтенный газ, зависящий от состояния учета [2].
Организация и проведение экспертного опроса
Экспертная группа составлена из 5 экспертов, принципиальное значение при выборе экспертов имели следующие факторы: опыт работы, компетентность, объективность. Каждый из экспертов имеет стаж работы в организациях, связанных с транспортировкой газа, более 5 лет, является высококвалифицированным специалистом оперативно-диспетчерских и аварийных служб линейно-промысловых управлений магистральных газопроводов.
В качестве формы опроса использован опросный лист. С целью обеспечения независимости оценок исключено взаимовлияние экспертов и снижено воздействие посторонних факторов. Работа организована таким образом, что эксперт находился в привычной для него обстановке, при отсутствии влияния на оценку других лиц [5].
На основании имеющейся информации синтезирована диаграмма «Количество неучтенного газа», представленная в [3] и объединяющая факторы, возникающие в каждой из групп. С помощью диаграммы экспертам предлагалось оценить степень возможности возникновения каждого из факторов в своей группе с помощью метода ранжирования. Данный метод является наиболее простым для проведения процедуры упорядочения объектов и не требует трудоемкого обучения экспертов. Метод идеально подходит для упорядочения объектов, число которых не превышает 15 [6]. Среди объектов нет одинаковых по сравниваемым показателям, т.е. нет эквивалентных объектов. В этом случае между ними существует только отношение строгого порядка. Пример результатов ранжирования для удобства обработки представлен в табл. 1.
После того как данные опросных листов собраны, проводится обработка полученных оценок. При ранжировании объектов используется мера согласованности мнений группы экспертов – дисперсионный коэффициент конкордации. При оценке объектов исследования эксперты расходятся во мнениях по решаемой проблеме. В связи с этим возникает необходимость количественной оценки степени согласия экспертов. Получение количественной меры согласованности позволяет обоснованно интерпретировать причины расхождения мнений.
Оценка согласованности мнений экспертов
Результаты опросов экспертов представлены в виде матрицы:
Для табл. 1 c представленными в ней результатами опроса экспертов в группе утечек, возникающих на газопроводах и сооружениях на них, вычислили ris для фактора возникновения утечки на стыковых соединениях
Рассчитаны суммы рангов для каждого фактора табл. 1:
ris2 = 25; ris3 = 24; ris4 = 8; ris5 = 19; ris6 = 8.
По результатам опроса экспертов в группе утечек, зависящих от физико-технических свойств газа и окружающей среды, получены суммы рангов для каждого из факторов: ris1 = 10; ris2 = 9; ris3 = 21 ris4 = 20. По результатам опроса экспертов в группе неучтенного газа, зависящего от состояния учета, получены суммы рангов для каждого из факторов: ris1 = 15; ris2 = 25; ris3 = 10. Проведена оценка математического ожидания для каждой группы факторов:
Рассчитан коэффициент конкордации W, значение которого должно быть больше заданного Wз (W > Wз). Обычно Wз = 0,5, т.е. при W > 0,5 выводы экспертов согласованы в большей мере, чем не согласованы. Таким образом, коэффициент конкордации для каждой из групп: W1 = 0,670857; W2 = 0,976; W3 = 0,881333 [1].
Таблица 1
Результаты ранжирования факторов
Утечки газа из газопроводов |
Э № 1 |
Э № 2 |
Э № 3 |
Э № 4 |
Э № 5 |
На стыковых соединениях |
4 |
3 |
3 |
5 |
6 |
Через раковины в полиэтиленовых и стенки асбестоцементных труб |
5 |
5 |
4 |
6 |
5 |
Коррозия труб |
6 |
6 |
5 |
4 |
3 |
Повреждения землеройными механизмами |
2 |
1 |
2 |
1 |
2 |
Утечки в арматуре |
3 |
4 |
6 |
2 |
4 |
Эксплуатационные работы |
1 |
2 |
1 |
3 |
1 |
Теоретико-множественная модель
Интерпретацией результатов будет служить теоретико-множественная модель [7–8]. Опишем подробно группы подмножеств, представленных в формуле (1) множества неучтенного газа и утечек.
Подмножество G – совокупность факторов возникновения утечек природного газа:
G = {g1, g2, g3, g4, g5, g6}, (1)
где g1 – фактор возникновения утечки на стыковых соединениях; g2 – фактор возникновения утечки через раковины в полиэтиленовых и стенки асбестоцементных труб; g3 – фактор возникновения утечки в результате коррозии труб; g4 – фактор возникновения утечки в результате повреждения газопровода землеройными механизмами; g5 – фактор возникновения в результате утечки в арматуре; g6 – фактор возникновения утечки в результате проведения эксплуатационных работ.
Подмножество F представляет собой совокупность факторов, влияющих на возникновение утечек в результате изменения физико-технических свойств газа и окружающей среды:
F = {f1, f2, f3, f4}, (2)
где f1 – фактор возникновения утечки в результате изменения температуры газа; f2 – фактор возникновения утечки в результате изменения давления газа; f3 – фактор возникновения утечки в результате изменения влажности газа; f4 – фактор возникновения утечки в результате изменения плотности газа.
Подмножество W представляет собой совокупность факторов влияющих на возникновение неучтенного газа из-за состояния учета:
W = {w1, w2, w3, w4, w5, w6, w7, w8, w9}, (3)
где w1 – фактор возникновения неучтенного газа в результате неточного или неправильного измерения температуры и давления газа и окружающей среды; w2 – фактор возникновения неучтенного газа в результате усреднения данных показаний измерительных приборов за слишком продолжительный промежуток времени, в течение которого расход сильно колеблется; w3 – фактор возникновения неучтенного газа в результате неточного определения или пренебрежения поправками, которые необходимо вводить в расчетные формулы; w4 – фактор возникновения неучтенного газа в результате неправильного или неточного определения фактической плотности учитываемого газа; w5 – фактор возникновения неучтенного газа в результате несоответствия учтенного счетчиками газа фактическому расходу, из-за чего при малом расходе счетчик дает заниженные показания или вообще не работает; w6 – фактор возникновения неучтенного газа в результате неправильного подбора диафрагм или дифманометров при измерении расходомерами, из-за чего перо «зашкаливает» или запись ведется в нижней части картограммы; w7 – фактор возникновения неучтенного газа в результате несвоевременного обнаружения счетчиков, дающих неправильные показания; w8 – фактор возникновения неучтенного газа в результате неправильной записи расхода газа по показаниям счетчиков; w9 – фактор возникновения неучтенного газа в результате технических и арифметических ошибок при обработке картограммы и при подсчетах.
Используя формулу , рассчитаем веса каждого фактора в группе. Произведя вычисления, построим диаграмму весов факторов возникновения утечек и неучтенного газа.
Нечеткая ситуационно-фреймовая модель
Для формализации полученных знаний применим нечеткую ситуационно-фреймовую модель [9]. Модель представлена в виде формализованного образования факторов возникновения утечек и неучтенного газа.
Модель состоит из фреймов. Фрейм имеет определенную внутреннюю структуру, состоящую из множества элементов, называемых слотами, которым также присваиваются имена. Именами слотов будут служить факторы из подмножеств G, W, F. Каждый слот представляется определенной структурой данных с конкретным названием фактора возникновения утечки или неучтенного газа. За слотами следуют текущие значения слотов. В значении слота представляется конкретная информация, относящаяся к объекту, описываемому этим фреймом. Значения слотов будут выглядеть как качественные оценки факторов возникновения утечек и неучтенного газа, т.е. «маленькая величина», «средняя величина», «большая величина».
Удобство применения данной модели заключается в том, что значением слота может быть практически что угодно: числа, формулы, тексты на естественном языке или программы, правила вывода или ссылки на другие слоты данного фрейма или других фреймов. В качестве значения слота может выступать набор слотов более низкого уровня, что позволяет реализовывать во фреймовых представлениях «принцип матрешки».
Диаграмма весов факторов модели возникновения утечек и неучтенного газа
Таблица 2
Типовые ситуации возникновения утечек и неучтенного газа на ЛЧ МГ
g1 |
средняя величина |
0,2 |
f3 |
большая величина |
0,35 |
… |
… |
… |
w7 |
маленькая величина |
0,08 |
Далее следуют демоны. Демоном называется процедура, автоматически запускаемая при выполнении некоторого условия. Демоны автоматически запускаются при обращении к соответствующему слоту. Демон является разновидностью связанной процедуры. Демонами фреймов будут являться весовые показатели факторов, найденные ранее и отраженные в диаграмме на рисунке. Пример структуры фрейма представлен в табл. 2.
Основным преимуществом фреймов как модели представления знаний является то, что она отражает концептуальную основу организации памяти человека, а также ее гибкость и наглядность. Значения слотов представляются в системе в единственном экземпляре, поскольку включаются только в один фрейм, описывающий наиболее понятия из всех тех, которые содержат слот с данным именем. Такое свойство систем фрейма обеспечивает экономное размещение базы знаний в памяти компьютера.
Так как значения слота имеют лингвистическую модификацию, то необходимо описать их при помощи нечеткого множества. Факторы возникновения утечек и неучтенного газа являются лингвистическими переменными. Представим их в виде терм-множеств. Например:
g1 = {маленькая величина, средняя иеличина, большая величина}.
Главной задачей при оперировании с нечеткими множествами является определение функции принадлежности. Введем переменные. Пусть K – универсальное множество и x – элемент K, описывающий лингвистическую переменную фактора возникновения утечки.
(4)
В качестве функции принадлежности используем симметричную гауссовскую функцию принадлежности, которая имеет вид
(5)
Все факторы возникновения утечек и неучтенного газа в физическом выражении обозначают величину объема V, измеряемую в м3. Из этого следует, что единицей измерения универсального множества и его элементов также будет м3.
Таблица 3
Границы значений универсального множества с указанием значений элементов
x, фактор возникновения утечки или неучтенного газа |
K, м3/ч |
x1, м3/ч |
x2, м3/ч |
x3, м3/ч |
g1 |
0...350 |
50 |
150 |
330 |
g2 |
0...500 |
100 |
200 |
450 |
g3 |
0...450 |
70 |
180 |
420 |
g4 |
0...100 |
10 |
50 |
90 |
g5 |
0...300 |
30 |
140 |
280 |
g6 |
0...150 |
20 |
70 |
140 |
f1 |
0...10 |
1 |
5 |
8 |
f2 |
0...5 |
0.5 |
2 |
4.5 |
f3 |
0...25 |
3 |
13 |
24 |
f4 |
0...20 |
2 |
10 |
18 |
w1 |
0...15 |
1...5 |
7 |
14 |
w2 |
0...40 |
5 |
20 |
35 |
w3 |
0..10 |
1 |
5 |
9 |
w4 |
0...60 |
6 |
30 |
55 |
w5 |
0...10 |
2 |
5.5 |
8.5 |
w6 |
0...50 |
4 |
25 |
47 |
w7 |
0...30 |
3 |
15 |
27 |
w8 |
0...80 |
9 |
40 |
75 |
w9 |
0...70 |
8 |
35 |
67 |
Расчет возможных потерь будет точным, если мы будем знать участок исследования, а именно – состав газопровода, его диаметр, давление, наличие сооружений на газопроводе, наличие и перечень измерительных приборов, с указанием их технических характеристик и т.д. Зададим условные границы значений, характеризующие потери за единицу времени м3/ч. В табл. 3 указаны возможные значения величины x.
Заключение
Полученные результаты вычислений коэффициента конкордации результатов опроса экспертов по группам свидетельствуют о согласованности суждений экспертов в полной мере, т.к значения коэффициента во всех случаях приближены к 1. Построена теоретико-множественная модель и рассчитаны веса факторов модели множества возникновения утечек и неучтенного газа. Для формализации полученных знаний применена нечеткая ситуационно-фреймовая модель.
Рецензенты:
Попов Г.А., д.т.н., профессор, заведующий кафедрой «Информационная безопасность», ФГБОУ ВПО «Астраханский государственный технический университет», г. Астрахань;
Ханова А.А., д.т.н., доцент кафедры «Прикладная информатика в экономике», ФГБОУ ВПО «Астраханский государственный технический университет», г. Астрахань.