Разработка эффективного маршрута газопровода с учётом социально-культурных факторов Текст научной статьи по специальности «Экономика и бизнес»

УДК 330.15

DOI: 10.18384/2310-6646-2018-3-94-103

РАЗРАБОТКА ЭФФЕКТИВНОГО МАРШРУТА ГАЗОПРОВОДА С УЧЁТОМ СОЦИАЛЬНО-КУЛЬТУРНЫХ ФАКТОРОВ

Новоселов АЛ.1, Новоселова ИЮ.1, Желтенков А.В.2

1Российский экономический университет имени Г.В. Плеханова 115054, г. Москва, Стремянный пер., д. 36, Российская Федерация 2 Московский государственный областной университет 141014, Московская область, г. Мытищи, ул. Веры Волошиной, д. 24

Аннотация. В статье рассматривается проблема построения эффективного маршрута с учётом капитальных и эксплуатационных затрат, историко-культурной экспертизы и загрязнения окружающей среды. Предлагается формула оценки суммарных затрат на построение и эксплуатацию отдельных фрагментов газопровода, включающая компенсационные выплаты населению. При этом используется нечёткий метод оценки затрат, позволяющий представить ожидаемые суммарные затраты в виде треугольных чисел. Оптимальный вариант маршрута газопровода предложено определять с помощью метода Шимбола-Оттермана в нечёткой форме. Для выбора окончательного решения предложено воспользоваться риск-функцией оценки затрат. В статье представлены алгоритмы предлагаемых модификаций указанных методов и дан численный пример.

Ключевые слова: негативное влияние, окружающая среда, ущерб, компенсационные выплаты, алгоритм Шимбола-Оттермана, риск-функция, нечёткие числа, треугольные числа, оптимальный маршрут.

OPTIMAL DESIGN OF A GAS TRANSMISION NETWORK REGARDING SOCIO-CULTURAL FACTORS

A. Novoselov1, I. Novoselova1, А. Zheltenkov2

1Plekhanov Russian University of Economics

36, Stremyanny ln, Moscow, 115054, Russian Federation

2Moscow Region State University

141014, Moscow region, Mytischi, ul. Vera Voloshina, house 24

Abstract. The article deals with the problem of designing an effective route taking into account capital and operating costs, historical and cultural expertise and environmental pollution. A formula is proposed for estimating the total costs for the construction and operation of gas pipeline sections including compensation payments to the population. At the same time, a fuzzy cost estimating method is used, which allows to present the expected total costs in the form of triangular numbers. It is proposed to calculate the optimal variant of a gas transmission network by using the Schimbel-Otterman method in a fuzzy form. To choose the final solution, it is suggested to use a risk-function of cost estimation. The algorithms of the proposed modifications of these methods are presented and a numerical example is given.

© CC BY Новоселов А.Л., Новоселова И.Ю., Желтенков А.В., 2018.

Keywords: negative impact, environment, damage, compensation payments, Schimbel-Otter-man algorithm, risk-function, fuzzy numbers, triangular numbers, optimal route.

Постановка проблемы

Россия обладает четвертью общемировых разведанных запасов природного газа и крупнейшей газотранспортной системой - её протяжённость составляет 168,3 тыс. км. Добываемый в России природный газ поступает в магистральные газопроводы, объединённые в Единую систему газоснабжения России (ЕСГ). Сеть газопроводов растёт в связи с расширением эксплуатируемых месторождений, расширением рынков сбыта природного газа. При прокладке маршрута нового газопровода учитываются геофизические факторы, протяжённость, влияние на экологическое состояние территории [1], памятники истории и культуры, а также влияние на среду обитания населения региона. Два последних аспекта представляют важные факторы, которые следует учитывать при разработке маршрута газопроводов, но строгих математических методов решения данной задачи до сих пор не было разработано.

В основе такой оценки лежит историко-культурная экспертиза, проводимая этнологами, археологами и архитекторами в местах предполагаемых транспортных коридоров [2]. Решения по выбору оптимального маршрута газопровода в местах проживания и ведения хозяйственной деятельности населения должны приниматься с учётом результатов этнологических экспертиз и минимизации вмешательства в исконную среду обитания, традиционное природопользование и образ жизни. Необходимы информирование местного населения и привлечение его представителей к процессу консультаций [3]. При выборе маршрута газопровода должна учитываться необходимость минимизации рисков негативного влияния на традиционную экономику и культурную жизнь населения. Возможности развития туризма также могут быть сокращены в результате прокладки газопровода. Совокупность приведённых выше факторов должна быть учтена при разработке оптимального маршрута газопровода наряду с капитальными затратами на его строительство и текущими затратами на его эксплуатацию.

Разработанная модель и алгоритм решения задачи

В основе выбора оптимального маршрута газопровода лежит ориентированный граф G, в котором отсутствуют обратные дуги. Каждая дуга (i, j) соответствует фрагменту газопровода, после которого начинается очередная альтернатива. Поскольку любому фрагменту газопровода (i, j) соответствуют затраты на строительство, эксплуатацию, компенсационные выплаты населению и вероятностную оценку штрафов и экстренных выплат населению в случае аварий, эти затраты следует рассматривать в динамике за период, охватывающий строительство и эксплуатацию газопровода. Для этого суммарные затраты следует просуммировать с учётом фактора времени, т. е. с использованием дисконтного множителя. Поскольку суммарные затраты относятся к фрагменту маршрута газопровода (i, j), то их обозначим V(. j. Значения каждого из составляющих суммарных затрат, относящихся к компенсационным выплатам населению и оцен-

ке ущерба окружающей среды, определяются экспертным методом. Поэтому в результате экспертной оценки находятся минимальное (min), ожидаемое (av) и максимальное (max) значение оцениваемого показателя. Тогда в год t величины оценки штрафов и экстренных выплат населению в случае аварий техногенного или природного характера - R, = Rj Rj)) ,акомпенсационные выплаты населению, связанные со строительством и эксплуатацией газопровода, -K,Kj^. Такое представление этих показателей соответствует нечётким (треугольным) числам A, которые задаются тремя значениями, которым соответствует оценка уверенности или функция принадлежности ц (Л), причём ц (Amin) = 0; ц (Aav) = 1; ц (Amax) = 0. Нечёткие числа удобно отображать графически, например, для треугольного числа А=(5; 8; 9) получаем графическое отображение, представленное на рис. 1.

1

г о =t о,в

а s"

1 ь о 0,6

s т 1 1 ГС 0,4

пп X 0,2

О ' *

4 5 6 7 В 9 10

Значение показателя

Рис. 1. Графическое отображение треугольного числа

Арифметические операции с нечёткими числами представлены, например, в [4]. Суммарные затраты для дуги (у) определяются по формуле как треугольное число:

f К,, - К,,)(! + г Г -f Zllf, (1 + г)'-' + f if,,, (1 + г )'

jr _ \ 'I п .л I/ t \L 1 X ' г-ж i-1 \1 1 X ' TT- Zt Ч1 t

V(ij) - .

t-T t-1 t-1 (1)

где 8^^ - эксплуатационные затраты на фрагменте: газопровода (у) е О в год 1;

Я.^ - нечёткая оценка штрафов и экстренных выплат населению в случае аварии техногенного или природного характера при эксплуатации фрагмента газопровода (у) е О в год £

Z,..• - капитальные затраты на строительство фрагмента газопровода (у) е О

1{у) В год Р;

к{^ - нечёткая оценка компенсационных выплат, связанных со строительством и эксплуатацией фрагмента газопровода (у), е О в год £

^ - текущий год;

(1, 2...Т) - период строительства фрагмента газопровода (у);

(Т + 1,... Т1) - период эксплуатации фрагмента газопровода (у).

ISSN 2072-8549

2

2018/№ 3

Для поиска варианта строительства газопровода ь = {1,....и}, соединяющего месторождение (начальное событие графа альтернатив строительства газопровода, равное единице) с конечным пунктом - магистральным газопроводом или терминалом отгрузки природного газа (конечное событие графа альтернатив строительства газопровода, равное п), с минимальными суммарными затратами, следует воспользоваться одним из следующих методов:

• методом Флойда-Уоршелла [5];

• методом Беллмана-Форда [6];

• методом Дейкстры [5];

• методом Джонсона [7];

• алгоритмом Левита [8];

• методом Шимбела [9];

• методом Оттермана [10].

Проведённый анализ методов поиска кратчайшего расстояния в орграфе позволяет утверждать, что наиболее целесообразным является использование комбинированного метода Шимбела-Оттермана. Особенность матричных операций в методе Шимбела заключается в реализации операции суммирования элементов а и Ь (с = а ® Ь), которые в общем виде представлены ниже:

где a и b - числа больше нуля.

С учётом того, что в поставленной задаче необходим поиск маршрута L = {1,....и} в условиях нечётко заданных суммарных затрат для каждого фрагмента маршрута, в данный алгоритм внесён ряд корректировок.

В предлагаемом алгоритме проводится поиск оптимального варианта маршрута газопровода для разных уровней функции принадлежности - от нуля до единицы отдельно для левой и правой границ треугольного числа. Ниже приводится расчёт для заданных значений функции принадлежности ц* и варианта границы треугольных чисел (левой или правой), характеризующих суммарные затраты V^ для отдельных фрагментов маршрута газопровода (ij) е G.

Шаг 1. Расчёт проводится для определённых суммарных затрат при заданных значениях функции принадлежности ц* и границы треугольных чисел (min, max), задающих суммарные затраты для фрагментов маршрута газопровода V^.

Шаг 2. Формирование матрицы смежности A, в которой при наличии дуги указываются суммарные затраты V/!, которые должны быть предварительно найдены по формуле (1):

с = <

0, если a > 0 л b = 0 0, если b > 0 л a = 0 0, если a = 0 л b = 0

(2)

a + b, если a > 0 л b > 0

0, в противном случае

V(ij)е G

Шаг 3. Присвоение B = A; p=1

Шаг 4. Расчёт матрицы Cp суммы весов между вершинами для всех путей, состоящих из двух дуг в соответствии с правилами (2). Элементы матрицы C определяются по формуле:

С" =1k

min К© h\akj® К >

' еСЛи Ьik > 0

(4)

0, если все akj © bk = 0

Шаг 5. Формирование матрицы трассировки Т?, показывающей предшествующую вершину для каждой достигаемой вершины, представленной в столбцах этой матрицы по правилу:

_ |к* исходяоиз _ kпш^{ац ©ЬЛ 1 а^ ©ЬЛ > 0),если За^ ©ЬЛ > 0 (5) I ° если все ак. © Ькк _ 0

Шаг 6. Проверка: есть ли в матрице Ср ненулевые элементы? Если есть, то переход к шагу 6; в противном случае - к шагу 7.

Шаг 7. Присв оение В = Л; р = р+1; переход к шагу 3.

Шаг 8. Определение минимального веса для лучшего варианта строительства газопровода:

=, mn, К}

(6)

Шаг 9. Формирование оптимального маршрута Ь ={1,....и} на основе первых строк матриц трассировки С', I = 1,2,..р*. Конец расчёта для заданной функции принадлежности и указанной границы треугольного числа, характеризующего суммарные затраты на реализацию газопровода в разрезе его фрагментов.

Приведённый алгоритм повторяется для значений и, соответственно, левой и правой границ оценки суммарных затрат для фрагментов маршрута газопровода. В результате определяются оптимальные по критерию минимизации суммарных затрат на строительство и эксплуатацию газопровода с учётом компенсационных выплат населению (1) варианты маршрутов.

На основе полученных оптимальных значений суммарных затрат следует построить риск-функцию. В соответствии с постановкой задачи, были скорректированы формулы построения риск-функции суммарных затрат и построен алгоритм, состоящий из следующих трёх шагов:

Шаг 1. Расчёт промежуточного значения а при варьировании О от до ^^

с априорно заданным интервалом Д (О = ^^ - Д,^Д,...^^):

а = <

Zmax G Z*,< G < Zl

ze - z* ' av

max av 7*

(7).

g - zm

7* _ 7* ' mm Zav Zmin

zl < G < Z*

Шаг 2. Расчёт промежуточного значения R при варьировании G от zЕ до zE с априорно заданным интервалом Д:

R =

Zmax - G Zmax Zmin

ZInn ^ G < Zmax

(8).

Шаг 3. Расчёт значений риск-функции Risk(G) при варьировании G от Zm;ln до Zm;ax с априорно заданным интервалом Д:

Risk (G) =

R[1 + i-^ln(1 -а)), Z1 < G < Z\

1 -а

1 -(1 - r) 1 + — 1П (1 -а)1, Zl < G < Z

а

ji ^min

(9).

1, G = Zl

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

Численный пример

На рис. 2 представлены маршруты газопровода от месторождения (вершина 1) до терминала отгрузки (вершина 6). На дугах представленного ориентированного графа указаны нечёткие оценки суммарных затрат.

(2;4;6)

(3;5;6)

(2;3;4)

Рис. 2. Ориентированный граф, представляющий совокупность маршрутов, составленных из отдельных фрагментов (у) е О снечёткимиоценкамисуммарных затрат

В результате скорректированного авторами статьи алгоритма Шимбела-Оттер-мана были получены результаты для границ нечётких оценок суммарных затрат (табл. 1). На рис. 3 представлено графическое отображение треугольного значения минимальных суммарных затрат для найденных маршрутов газопровода.

Таблица 1

Варианты оптимального маршрута газопровода для нечётких оценок суммарных затрат

Вариант расчёта Минимальные суммарные затраты Выбранный маршрут

Левая граница нечёткойоценки затрат при значении функции принадлежности ц = 0 9,4 L = {1 - 2 - 4 - 5 - 6}

Среднее гначение нечёткой оценки затрат при значении функции принадлежности ц = 1 16,0 L = {1 - 2 - 4 - 6}

Правая граница нечёткой оцетки затрат при значении функции принадлежности ц = 0 21,0 L = {1 - 2 - 4 - 6}

Рис. 3. Графическое представление нечётких оптимальных суммарных затрат при выборе оптимального маршрута газопровода

Для детализации полученного решения рассчитывается риск-функция затрат (рис. 4).

Риск-функция оценки затрат на строите л ьсгео

газопровода I,{l-2-4-S-6| ¿ = [1-2-4-6}

9А 16,0

Затраты на строительство газопровода, млрд руб.

Рис. 4. График риск-функции затрат на реализацию газопровода

На графике выделена точка, соответствующая среднему значению нечёткого треугольного числа. Слева от этой точки оптимальным является маршрут L = {1 -2-4-5-6}, а справа - L = {1 -2-4-6} . Полученный график риск-функции следует использовать для определения оптимального варианта маршрута прокладки при заданном уровне риска. Например, при уровне риска 0,2 (20%)оптимальный маршрут газопровода L = {1 - 2 - 4 - 5 - 6}.

Заключение

Предложенные алгоритмы были реализованы в виде программного комплекса в среде MS-Excel с использованием авторских макросов, написанных на языке VBA. Программный комплекс позволяет выявить все маршруты для найденного нечёткого значения минимальных суммарных затрат на строительство и эксплуатацию газопровода и провести визуализацию всех этапов поиска решения поставленной задачи. Апробация разработанных алгоритмов с использованием данного программного комплекса доказывает его высокую эффективность.

ЛИТЕРАТУРА

1. Кормен Т.Х., Лейзерсон Ч.И., Ривест Р.Л., Штайн К. Алгоритмы: построение и анализ. 2-е изд. М.: Вильямс, 2006. 1296 с.

2. Левит Б.Ю. Алгоритмы поиска кратчайших путей на графе. Труды института гидродинамики СО АН СССР // Моделирование процессов управления. 1971. Вып. 4. С. 117-148.

3. Левитин А.В. Алгоритмы: введение в разработку и анализ. М.: Вильямс, 2006. 576 с.

4. Новоселов А.Л., Медведева О.Е., Новоселова И.Ю. Экономика, организация и управление в области недропользования: учеб. и практикум. М.: Юрайт, 2014. 625 с.

5. Порфирьев Б.Н., Тулупов А.С. Оценка экологической опасности и прогноз экономического ущерба от аварийных ситуаций на промышленных предприятиях // Проблемы прогнозирования. 2017. № 6. С. 37-46.

6. Bellman R. On a Routing Problem // Quarterly of Applied Mathematics. 1958. Vol. 16. No. 1. P. 87-90.

7. Krupnik I., Mason R. Northern Ethnographic Landscapes // Perspectives from Circumpolar Nations / ed. T. Horton. Washington, 2004. P. 1-13.

8. Murashko O.A. What is the Etnologicheskaia Ekspertiza in Russia? // Special issue of Si-birica: the Interdisciplinary Journal of Siberian. Studies. 2006. Vol. 5. No. 2. P. 77-94.

9. Otterman J. Matrix Multiplication in Search for Alternate Routes // Electronica Com. IT. 1963. Vol. 38. No. 2. P. 156-164.

10. Shimbel A. Structural Parameters of Communication networks // Bulletin of Mathematical Biophysics. 1953. Vol. 15. No. 4. P. 501-507.

REFERENCES

1. Cormen T.H., Leiserson Ch.E., Rivest R.L., Stein C. Algoritmy:postroenie i analiz [Introduction

to Algorithms]. Moscow, Williams Publ., 2006. 1296 p.

2. Levit B.Yu. [Algorithms for Finding the Shortest Paths in a Graph. Proceedings of the Institute

of Hydrodynamics of the Siberian Branch of the Academy of Sciences of the USSR]. In: Modelirovanie protsessov upravleniya [Modeling of Management Processes], 1971, no. 4, pp. 117-148.

3. Levitin A.V. Algoritmy: vvedenie v razrabotku i analiz [Algorithms: Introduction to Design

and Analysis]. Moscow, Williams Publ., 2006. 576 p.

4. Novoselov A.L., Medvedeva O.E., Novoselova I.Yu. Ekonomika, organizatsiya i upravlenie

v oblasti nedropol'zovaniya [Economics, Organization and Management in the Field of Mineral Resources]. Moscow, Yurait Publ., 2014. 625 p.

5. Porfir'ev B.N., Tulupov A.S. [Environmental Risk Assessment and Forecast of Economic

Damage in Emergency Cases at Industrial Enterprises]. In: Problemyprognozirovaniya [The Problems of Forecasting], 2017, no. 6, pp. 37-46.

6. Bellman R. On a Routing Problem. In: Quarterly of Applied Mathematics, 1958, vol. 16, no. 1,

pp. 87-90.

7. Krupnik I., Mason R., Northern Ethnographic Landscapes. In: Horton T., ed. Perspectives from Circumpolar Nations. Washington, 2004, pp. 1-13.

8. Murashko O.A. What is the Ethnological Expertise in Russia? In: Special issue of Sibirica: the

Interdisciplinary Journal of Siberian Studies, 2006, vol. 5, no. 2, pp. 77-94.

9. Otterman J. Matrix Multiplication in Search for Alternate Routes. In: Electronica Com. IT,

1963, vol. 38, no. 2, pp. 156-164.

10. Shimbel A. Structural Parameters of Communication Networks. In: Bulletin of Mathematical Biophysics, 1953, vol. 15, no. 4, pp. 501-507.

БЛАГОДАРНОСТИ

Исследование выполнено при финансовой поддержке РФФИ, проект «Разработка экономического механизма согласования интересов государства и бизнеса для реализации региональных природоохранных проектов и программ» № 1702-00010- ОГН

ACKNOWLEDGEMENTS

The study was carried out with the financial support of the Russian Foundation for Basic Research, the project "Development of an economic mechanism for coordinating the interests of the state and business for the implementation of regional environmental projects and programs" No. 17-02-00010-0GN

ИНФОРМАЦИЯ ОБ АВТОРАХ

Новоселов Андрей Леонидович - доктор экономических наук, профессор кафедры математических методов в экономике Российского экономического университета им. Г.В. Плеханова; e-mail: alnov2004@yandex.ru

Новоселова Ирина Юрьевна - доктор экономических наук, профессор кафедры математических методов в экономике Российского экономического университета им. Г.В. Плеханова; e-mail: iunov2010@yandex.ru

Vioy

Желтенков Александр Владимирович - доктор экономических наук, профессор, заведующий кафедрой менеджмента и государственного управления Московского государственного областного университета; e-mail: kaf-menedg@mgou.ru

INFORMATION ABOUT AUTHORS

Andrey L. Novoselov - Doctor of Economics, professor at the Department of Mathematical Methods in Economics, Plekhanov Russian University of Economics; e-mail: alnov2004@yandex.ru

Irina Yu. Novoselova - Doctor of Economics, professor at the Department of Mathematical Methods in Economics, Plekhanov Russian University of Economics; e-mail: iunov20l0@yandex.ru

Alexander V. Zheltenkov - Doctor of Economics, professor, head of the Department of Management and Public Administration, Moscow Region State University; e-mail: kaf-menedg@mgou.ru

ПРАВИЛЬНАЯ ССЫЛКА НА СТАТЬЮ Новоселов А.Л., Новоселова И.Ю., Желтенков А.В. Разработка эффективного маршрута газопровода с учётом социально-культурных факторов // Вестник Московского государственного областного университета. Серия: Экономика. 2018. № 3. С. 94-103. DOI: 10.18384/2310-6646-2018-3-94-103

FOR CITATION

Novoselov A.L., Novoselova I.Y., Zheltenkov A.V. Optimal Design of a Gas Transmisión Network Regarding Socio-Cultural Factors In: Bulletin of Moscow Region State University. Series: Economics, 2018, no. 3, рр. 94-103. DOI: 10.18384/2310-6646-2018-3-94-103