УДК 621.316.1.05
05.00.00 Технические науки
К ВОПРОСУ НАХОЖДЕНИЯ ОПТИМАЛЬНЫХ ПОТЕРЬ МОЩНОСТИ В ЭЛЕКТРОЭНЕРГЕТИЧЕСКИХ СЕТЯХ НА ОСНОВЕ РЕШЕНИЯ ТРАНСПОРТНОЙ ЗАДАЧИ МЕТОДОМ ПОТЕНЦИАЛОВ
Лаптев Владимир Николаевич к.т.н, доцент
Степанов Владимир Васильевич д.т.н., профессор
Атрощенко Валерий Александрович д.т.н., профессор
Кабанков Юрий Андреевич к.т.н, профессор
Степанова Марина Валерьевна
преподаватель
vvs04367@mail.ru
Кубанский государственный аграрный университет, г. Краснодар Кубанский государственный технологический университет, Россия г. Краснодар, ул. Красная 135, к. 123
Краснодарское высшее военное авиационное училище им. А.К. Серова, г. Краснодар
Работа носит научный и практический характер, так как рассматриваемый в ней процесс оптимизации потерь мощности в электроэнергетических сетях основан на использовании транспортной задачи, а именно правила минимального элемента транспортной матрицы, с помощью которой находится допустимое решение. Механизм улучшения полученного допустимого решения поставленной задачи ориентирован на использование метода потенциалов, присваивая каждой строке и столбцу матрицы соответствующий потенциал. Используя базисную переменную и манипулируя элементами транспортной матрицы по определённому алгоритму, осуществляем построение нового
UDC 621.316.1.05 Technical science
TO THE QUESTION OF FINDING THE OPTIMAL POWER LOSS ELECTRIC POWER NETWORKS BASED UPON SOLUTIONS OF TRANSPORT PROBLEMS BY THE POTENTIAL METHOD
Laptev Vladimir Nikolaevich Cand.Tech.Sci., Associate Professor
Stepanov Vladimir Vasilievich Doctor of engineering sciences, professor
Atroschenko Valery Aleksandrovich Doctor of engineering sciences, professor
Kabankov Yuri Andreevich Cand.Tech.Sci., Professor
Stepanova Marina Valerievna lecturer
vvs04367@mail.ru
Kuban State Agrarian University, Krasnodar Kuban State Technological University, Krasnodar, Russia, Krasnodar
Krasnodar Higher Military Aviation School named after A.K. Serov, Krasnodar, Russia
The work is scientific and practical in nature, same as the discussed process of optimization of power losses in power networks based on the use of the transportation problem, namely the rules of the minimum element of the transport matrix, which is a valid solution. The mechanism of the improvement of the obtained feasible solutions to the tasks focuses on using the potential method, assigning each row and column of the matrix their corresponding potential. Using a basic variable and manipulating the elements of the transport matrix according to a specific algorithm, we implement the construction of a new support plan. Based on the plan and its acyclic, we build two cycles and determine the values of the estimates between direct and indirect transportation
опорного плана. На основе представленного плана и его ацикличности, производится построение двух циклов и определяются значения оценок между прямыми и косвенными затратами на транспортировку единицы мощности. Полученный алгоритм позволяет легко запрограммировать все вычислительные операции и быстро получить результат оптимизации потерь мощности в проектируемой электроэнергетической сети, что обеспечивает снижение экономических и технологических затрат. В результате совместного использования транспортной задачи и метода потенциалов удается более эффективно найти оптимальное решение, на основе улучшенного решения, и получить схему оптимальной электрической сети. В качестве потребителей электроэнергии на практике могут выступать городские и сельскохозяйственные, промышленные предприятия, краевого, так и районного подчинения, включая промышленных, сельскохозяйственных и индивидуальных потребителей
costs per unit of output. The resulting algorithm makes it easy to program all computational operations and quickly obtain the result of the optimization of power losses in the design of the electricity network, thereby reducing the economic and technological costs. As a result of joint use of transportation problems and potential method, we could find the optimal solution based on improved solutions more efficiently, and a scheme of optimal electric network. As consumers, in practice, there may be urban and agricultural, industrial enterprises, regional and district subordination, including industrial, agricultural and individual customers
КЛЮЧЕВЫЕ СЛОВА: МЕТОД ПОТЕНЦИАЛОВ, МАТЕМАТИЧЕСКАЯ МОДЕЛЬ, ИСТОЧНИКИ ЭНЕРГИИ, ИСТОЧНИКИ ПОТРЕБИТЕЛЕЙ, МОЩНОСТЬ, ПРОЕКТИРУЕМЫЕ ЛИНИИ ЭЛЕКТРОПЕРЕДАЧ, ПРЯМЫЕ И КОСВЕННЫЕ ЗАТРАТЫ, УДЕЛЬНЫЕ ЗАТРАТЫ
Keywords: METHOD OF POTENTIALS, MATHEMATICAL MODEL, ENERGY SOURCES, SOURCES OF CONSUMERS, POWER OF THE DESIGNED POWER LINES, DIRECT AND INDIRECT COSTS, UNIT COSTS
В настоящее время все больше и больше ужесточаются требования к внедрению энергосберегающих технологий практически во всех областях производства и распределительных сетей электроэнергии. Возникающий дефицит мощностей, связанный с потерями электроэнергии, рост цен на услуги ее транспортировки вызывают необходимость поиска экономии как за счет снижения потерь мощностей, так и за счет эффективного ее распределения.
Одним из перспективных направлений в области оптимизации в электроэнергетических сетях являются пути усовершенствования систем
эксплуатации и оптимизации схем распределительных сетей, в основу которого положены задача линейного программирования и ее модификации - такие как транспортная задача.
Известно, что транспортная задача - это задача, связанная с нахождением таких путей доставки готового продукта от пунктов производства к пунктам потребления, при которых общая стоимость транспортировки будет оптимальной. Несомненно, что такая математическая модель этой задачи может быть использована и в электроэнергетике, где под продуктом следует понимать электрическую мощность, передаваемую от источников энергии к потребителям по линиям электропередач.
В качестве потребителей электроэнергии выступают как городские, так и сельскохозяйственные, промышленные предприятия как краевого, так и районного подчинения, включая промышленных, сельскохозяйственных и индивидуальных потребителей.
В общем случае оптимизации подлежать затраты на формирование (проектирование) схемы электрической сети [1, 2, 3], которая состоит из линий электропередачи, связывающей узлы источников электроэнергии и узлы ее потребителей.
В настоящее время на базе использования математической модели транспортной задачи решено достаточно много задач электроэнергетики как в части оптимизации схем электроэнергетической сети, так и в части нахождения оптимального пути ее транспортировки, а также в области
проектирования возможного варианта расположения узлов и варианта сооружения линий электропередач.
Основополагающим фактором в этих задачах являлось введение удельных затрат на транспортировку мощностей по линиям между узлами источников и потребителей, которые подлежат оптимизации и представление этой задачи в виде матрицы.
В работе [3] проведен полный расчет проектируемой электрической сети межрегиональной сетевой компании с использованием транспортной задачи и представлена матрица, рассчитанная по определенному алгоритму. Полученное решение является допустимым, а найденные в результате расчетов свободные переменные: Х14 = 0, = 0-, Х21 = 0, Х22 = 0, Х22 = 0, Хзз = 0, Х34 = 0 позволили значительно упростить электрическую сеть компании, найти оптимальные пути транспортировки электроэнергии, существенно уменьшить потребление дорогостоящего проводника (рисунок 1).
Ми 2
I
I
Мщ
Рисунок 1 - Схема электрической сети, отвечающая допустимому
решению
Последнее явилось исходным материалом для решения поставленной в этой работе задачи оптимизации потери мощностей в сети электропередачи на базе использования метода потенциалов (метод улучшения допустимого решения).
С целью более детального рассмотрения метода потенциалов сформулируем поставленную задачу следующим образом, применительно к полученному допустимому решению.
Из трех источников электроэнергии, с соответствующей мощностью Мш ^где 1=1,3), необходимо эффективно транспортировать электро-
энергию, снижая потери между узлами источников энергии и узлами потребителей Мщ ^где / '15).
Составим математическую модель, вводя, как и ранее, удельные стоимости транспортировки электроэнергии от источников к потребителям, обозначив их через в виде матрицы {2^}.
Запишем математическую модель задачи:
Составим опорный план по правилу минимального элемента, и запишем в виде матрицы (таблица 1).
МИ1 =55 Ми2=35 Миз=20
Мп1=15 Мп2=20 Мш=30 МП4=20 Мп5=25
Источники энергии Потребители энергии Мощность источника
Мп1 Мп2 МП3 Мп4 Мп5
Ми1 1,1 0 1,5 0 1,3 0 1,7 0 1,8 0 55
Ми2 1,8 0 2,2 0 2,3 0 2,0 0 2,1 0 35
Ми3 1,7 0 2,1 0 2,5 0 2,1 0 2,3 10 20
Мощность потребителей 15 20 30 20 25 110
*
Обозначим излишек нераспределенной мощности от поставщика М т ,
*
и недостаток в поставке мощности потребителю через М П.
Находим первую строку и первый столбец таблицы 1 с минимальной удельной стоимостью на транспортировку единицы мощности - в нашем
случае она равна 1,1. Помещаем в эту клетку вместо нуля меньшую из
* *
величин М И1=55, М П1=15. В результате в первом столбце установился баланс мощностей, но в первой строке данной таблицы баланс не установлен.
Поступаем аналогично, находим следующую клетку с минимальной удельной стоимостью транспортировки по этой же строке - в нашем
*
случае она равна 1,3 и помещаем в нее меньшее из чисел М И1=40 *
(излишек), М П1 =30.
Повторяем перечисленные выше операции для первой строки второго
* *
столбца в рабочей матрице, внося меньшую из величин М П2=20 и М И1=10 (остаток).
В результате проведенных манипуляций установится баланс мощностей по первой строке и первому столбцу транспортной матрицы таблицы 1, а в оставшиеся клетки первой строки и столбца вносим нулевые значения.
Для оставшихся клеток матрицы повторяем аналогичную процедуру,
начиная со второй строки и второго столбца, находим клетку - 214=2,0,
* *
вносим меньшее из значений М П4=20 и М И2=35.
Повторяя операцию, находим клетку с минимальным значением
* *
225=2,1 и вносим в нее меньшее из М И2=15 иМ П5=20.
В результате установлен баланс по всем строкам и столбцам транспортной матрицы.
Источники энергии Потребители энергии Мощность источника
МП1 МП2 МПз МП4 МП5
МИ1 1,1 15 1,5 10 1,3 30 1,7 0 1,8 0 55
МИ2 1,8 0 2,2 0 2,3 0 2,0 20 2,1 15 35
Миз 1,7 0 2,1 10 2,5 0 2,1 0 2,3 10 20
Мощность потребителей 15 20 30 20 25 110
Транспортные расходы при транспортировке электрической энергии от источников к потребителям, согласно таблице 2, составят:
7=711X11+Ъ12Х12+Ъ13Х13+Ъ24Х24+Ъ25Х25+Ъ32Х32+Ъ35Х35= =15*1,1+10*1,5+30*1,3+20*2,0+15*2,1+10*2,1+10*2,3=186
Как видно из приведенного выше алгоритма, метод нахождения допустимого решения достаточно трудоемок, поскольку нахождения его для соответствующей свободной переменной необходимо строить цикл, и каждый раз определять изменение целевой функции.
С целью улучшения полученного допустимого решения воспользуемся методом потенциалов, учитывая, что загруженных клеток семь, а свободных неизвестных восемь, приходим к выводу, что план является нецикличным (ацикличный).
Введем для этого в рассматриваемую матрицу таблицы 2 каждой строке - и и столбцам -У]-, соответствующие потенциалы и и у-, связанные с базисной переменной 2^ соотношением:
2ц= и + V (1.1)
где 1=1 ,т;
7=1 ,п
2у - удельные затраты на транспортировку единицы мощности по линиям, между узлами источников и потребителей.
Полагая, что и1 =0 определим остальные значения потенциалов по строкам и столбцам - и2, и3, У1, У2, У3, У4, У5, используя соотношение (1.1), учитывая, что нам известны и1 =0, а 2д в занятых клетках матрицы тоже известны.
Получим: Ul =0
У= 2ll- Ul = 1,1 - 0=1,1 У2= 212- и! = 1,5 - 0=1,5 У3= 213- и1 = 1,3 - 0=1,3 и2= 225- У5 = 2,1 -1,7=0,4 и3= 232- У2 = 2,1 -1,5=0,6 У5= 235- и3 = 2,3 - 0,6=1,7
У4= 224- и2 = 2,0 - 0,4=1,6
Определим значения оценок, то есть разницу между прямыми и косвенными удельными затратами на транспортировку единицы мощности идущего от источника к потребителям, обозначив их соответственно Zг^, для свободных клеток матрицы (таблица 2):
г14= 2±а - {и1+У4) = 1,7 - (0+1,6)=0,1;
^15=^15 " {и&У5) =1,8-(0+1,7)=0,1;
Г =221 - ■;и+У1) = 1,8 - (0,4+1,1)=0,3;
222 = 222 - ■ и2+У2) = 2,2 -(0,4+1,5)=0,3;
¿22 = 222 - и+Уз) = 2,3 - (0,4+1,3)=0,6
^31 = ^31 - ■Лз+У1) = 1,7 - (0,6+1,1)=0
^33 = ^33 - ■ из+Уз) = 2,5 - (0,6+1,3)=0,6
^34 = 2^4 ■ - и+Уд = = 2,1 - (0,6+1,6) =-0,1.
Как видно из полученных расчетов имеет место клетка таблицы с отрицательной оценкой, а именно третья строка и четвертый столбец -134= 2,1. Строим для нее цикл на основе минимального значения
Vi У2 У3 У4 У5
Источники энергии Потребители энергии Мощность источника
Мп1 МП2 Мп3 Мп4 Мп5
и Ми1 1,1 15 1,5 10 1,3 30 1,7 0 1,8 0 55
и Ми2 1,8 0 2,2 0 2,3 0 2,0 10 2,1 25 35
и Ми3 1,7 0 2,1 10 2,5 0 2,1 10 2,3 0 20
Мощность потребителей 15 20 30 20 25 110
Найдем функцию цели:
2=211Х11 +212X12+213X1 з+224X24+225X25+232X32+234X34= =15*1,1+10*1,5+30*1,3+10*2,0+25*2,1 + 10*2,1 + 10*2.1=186
Выбираем опять клетку таблицы с наименьшей оценкой и строим для нее цикл, перемещая мощность в 10 удельных единиц, прибавляя там, где стоит знак «+» и вычитаем там, где стоит знак «-».
В результате перемещений по построенному циклу получаем новый план:
Проверим полученный план на оптимальность. Для этого подсчитаем потенциалы для клеток таблицы, содержащих свободные переменные:
и1 =0
г., = и+ У4; У4= 1,7- и=1,7 = и+ У5; У5= 1,8- и=1,8 = и2+ У1; 71=1,8- и=1,5 Г:: = и2+ У2; У2=2,2- и2=0,4 !:, = и2+ У3; и2=2,3- У3=0,3 = и3+ У1; У=1,7- и3=1,2 277= и3+ У3; У3=2,5- и3=2,0 11== и3+ У5; и3=2,3- У5=0,5
Определяем значение оценок для всех свободных клеток матрицы таблицы 4:
г14Г214-(и1+У4)= 1,7-1,7=0;
=г15 - (и,+ У5) =1,8-1,8=0; Лп7=221 - (1/2+У0 = 1,7 - (0,3+2,1)=0,2; Т72 =122 - (и2+У2) = 2,2 -(0,3+0,4)=!,5; = ¿23 - и+У3) = 2,3 - (0,3+2,0)=0;
= ^31 - из+Уг) = 1,7- (0,3+1,5)=-0,1;
777 = - . = 2,5 - (0,5+2,0)=0; = {из+У5} = 2,3 - (0,5+1,5)=0,3.
Как видно из полученных расчетов имеет место клетка таблицы с отрицательной оценкой, а именно третья строка и первый столбец - 231 = 1,7, то есть план для нее не является оптимальным. Строим для нее цикл, перемещая по циклу мощность в 10 удельных единиц, прибавляя там, где стоит знак «+» и вычитаем там, где стоит знак «-».
В результате проведенных операций получаем новый опорный план (таблица 4).
На основе полученного нового плана найдем уже для него значение целевой функции:
^=^11Х11+^12Х12+^13Х13+^24Х24+^25Х25+^31Х31+^35Х35 =
=5*1,1+20*1,5+30*1,3+10*2,0+25*2,1+10*1,7+10*2,1=185
Таблица 4
Источники энергии Потребители энергии Мощность источника
Мп1 МП2 Мп3 Мп4 Мп5
Ми1 1,1 5 1,5 20 1,3 30 1,7 0 1,8 0 55
Ми2 1,8 0 2,2 0 2,3 0 2,0 10 2,1 25 35
Миз 1,7 10 2,1 0 2,5 0 2,1 10 2,3 0 20
Мощность потребителей 15 20 30 20 25 110
Мощность, передаваемая по проектируемым линиям электропередач, связывающая потребителей с источниками питания является искомой величиной и обозначается, как и ранее, Х^.
Таким образом, можно сделать заключение, что в результате использования метода потенциалов удается более эффективно определить оптимальное решение на основе улучшенного решения (таблица 4) и получить схему оптимальной электрической сети.
Учитывая условие (1.1) перевод любой другой переменной в базис приведет к увеличению Следовательно, полученное решение
оптимально, а схема электрической сети, соответствующая этому решению, представлена на рисунке 2.
Ми 2
МП2
Рисунок 2 - Схема оптимальной электрической сети
Обсуждая полученный результат, можно сделать вывод, что решение задачи оптимизации потерь мощности в электроэнергетических сетях базируется на использовании транспортной задачи и представления транспортной матрицы размерности п*т (где п - количество источников питания, а т - количество потребителей). На основе представленной матрицы находятся допустимые решения по правилу минимального элемента удельной стоимости.
Следующим шагом на основе полученного допустимого решения, используя метод потенциалов (1.1) находим улучшенное допустимое решение, а затем проверяем его на оптимальность.
Полученный таким образом алгоритм позволяет легко запрограммировать все вычисленные операции и быстро получить
результат оптимизации потерь мощности в проектируемой электроэнергетической сети,
что обеспечивает снижение экономических и технологических затрат.
Список литературы:
1. Аввакумов В.Г. Постановка и решение электроэнергетических задач исследования операций. - Киев: Вища школа, 1983.
2. Костин В.Н. Оптимизационные задачи электроэнергетики Учебное пособие -СПб.: СЗ ГЗТУ. Издательско-полиграфическая ассоциация вузов Санкт-Петербурга, 2003.
3. Степанов В.В. Степанов. Оптимизация поставки электроэнергии на примере использования транспортной задачи. - / В.В. Степанов, М.В. Степанова, В.В. Степанов. Сборник научных статей XVII всероссийская научно-практическая конференция (заочная) «Инновационные технологии в образовательном процессе, Краснодар, 2015.
4. Степанов В.В. Применение транспортной задачи для оптимизации в межрегиональных сетевых компаниях передачи электроэнергии. - / В.В. Степанов, М.В. Степанова, Ю.А. Кабанков. «YI Международная научно-практическая конференция молодых ученых, посвященная 55-й годовщине полета Ю.А.Гагарина в космос», Краснодар, 2016.
5. Курицкий Б .Я. Поиск оптимальных решений средствами EXCEL 7.0.- СПб.: ВНУ, 1997.
References
1. Avvakumov V.G. Postanovka i reshenie jelektrojenergeticheskih zadach issledovanija operacij. - Kiev: Vishha shkola, 1983.
2. Kostin V.N. Optimizacionnye zadachi jelektrojenergetiki Uchebnoe posobie - SPb.: SZ GZTU. Izdatel'sko-poligraficheskaja associacija vuzov Sankt-Peterburga, 2003.
3. Stepanov V.V. Stepanov. Optimizacija postavki jelektrojenergii na primere ispol'zovanija transportnoj zadachi. - / V.V. Stepanov, M.V. Stepanova, V.V. Stepanov. Sbornik nauchnyh statej XVII vserossijskaja nauchno-prakticheskaja konferencija (zaochnaja) «Innovacionnye tehnologii v obrazovatel'nom processe, Krasnodar, 2015.
4. Stepanov V.V. Primenenie transportnoj zadachi dlja optimizacii v mezhregional'nyh setevyh kompanijah peredachi jelektrojenergii. - / V.V. Stepanov, M.V. Stepanova, Ju.A. Kabankov. «YI Mezhdunarodnaja nauchno-prakticheskaja konferencija molodyh uchenyh, posvjashhennaja 55-j godovshhine poleta Ju.A.Gagarina v kosmos», Krasnodar, 2016.
5. Kurickij B.Ja. Poisk optimal'nyh reshenij sredstvami EXCEL 7.0.- SPb.: VNV, 1997.