WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!


 

На правах рукописи

Свеженцева Ольга Владимировна

МЕТОДЫ И АЛГОРИТМЫ ОБОСНОВАНИЯ РАЦИОНАЛЬНОЙ КОНФИГУРАЦИИ СИСТЕМ ЭЛЕКТРОСНАБЖЕНИЯ

Специальность 05.14.02 –Электрические станции и

  электроэнергетические системы

АВТОРЕФЕРАТ

диссертации на соискание ученой степени

кандидата технических наук

Иркутск 2012

Работа выполнена на кафедре электроснабжения и электротехники Национального исследовательского Иркутского государственного технического университета

Научные руководители:                чл. - корр. РАН, доктор технических наук,

                                               профессор Воропай Николай Иванович

                                               доцент, кандидат технических наук

                               

Официальные оппоненты:                доктор технических наук, профессор

Иркутского государственного

университета путей сообщения

                                               Крюков Андрей Васильевич

               кандидат технических наук, доцент

       Иркутского государственного

технического университета

Акишин Леонид Александрович        

Ведущая организация:                Национальный исследовательский Томский политехнический университет        

Защита состоится «27» ноября 2012 г. в 9-00 на заседании диссертационного совета Д003. 017.01 при Институте систем энергетики им. Л.А. Мелентьева СО РАН по адресу: 664033, Иркутск, ул. Лермонтова, 130, к. 355.

С диссертацией можно ознакомиться в библиотеке Института систем энергетики им. Л.А. Мелентьева СО РАН.

Отзывы на автореферат в двух экземплярах с подписью составителя, заверенные печатью организации, просим отправлять по адресу: 664033, Иркутск, ул. Лермонтова, 130, на имя ученого секретаря диссертационного совета.

Автореферат разослан «  »  2012 года.

Ученый секретарь диссертационного совета,                

доктор технических наук, профессор                                        А.М. Клер        

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы. Современные электроэнергетические системы (в том числе системы электроснабжения) являются сложными территориально протяженными системами, имеющими неоднородную структуру электрических сетей. Ввиду сложности и многомерности современных систем электроснабжения, многовариантности и многокритериальности, наличия различных предпочтений при выборе решений, проблема обоснования развития этих систем в виде общей задачи исследования операций является громоздкой и с практической точки зрения непреодолимой.

С учетом сложности проблемы ее решение рассматривается как система задач, поэтапно уточняющих и детализирующих решения по развитию систем электроснабжения.

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

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

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

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

Задачи исследования. В соответствии с целью были поставлены и решены следующие задачи:

  1. Разработка  и содержательное толкование иерархической модели по обоснованию развития систем электроснабжения городов;
  2. Формализация комплекса задач оптимального выбора решений при обосновании рациональной конфигурации систем электроснабжения в зависимости от территориальных уровней;
  3. Разработка математических моделей решения следующих задач:
    • размещения единичного ИП в радиальной электрической сети с учетом произвольных ограничений на местности;
    • оптимального закрепления потребителей за несколькими ИП;
    • оптимального размещения нескольких ИП и закрепления за ними потребителей в системе электроснабжения (два случая: одинаковые типоразмеры ИП,  разные типоразмеры ИП);
  4. Разработка методов решения следующих задач:
    • оптимальное закрепление потребителей за ИП с использованием, адаптированного эвристического алгоритма на основе методов комбинаторного анализа;
    • оптимальное размещение нескольких ИП и закрепления за ними потребителей в системе электроснабжения (два случая: одинаковые типоразмеры ИП разные типоразмеры ИП);
    • оптимальный выбор конфигурации распределительной электрической сети с помощью алгоритмов оптимизации на графах:
        1. выбор оптимальной конфигурации распре-делительной электрической сети петлевой структуры;
        2. построение оптимальной трассы линии электропередачи, связывающей потребителя с ИП с учетом ограничений на местности;
  1. Исследование процесса размещения источников питания и закрепления за ними потребителей на примере решения задачи синтеза рациональной конфигурации распределительной электрической сети микрорайона города.

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

Объект исследования. Объектом исследования являются системы электроснабжения различного назначения.

Предмет исследования. Методы и алгоритмы выбора оптимальных решений при обосновании рациональной конфигурации систем электроснабжения.

Научную новизну диссертации представляют следующие основные результаты, которые выносятся на защиту:

1. Методический подход :

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

Все перечисленные результаты диссертационной работы, представляющие научную новизну, получены впервые.

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

  1. Выбора оптимальной мощности ИП в распределительных электрических сетях;
  2. Равномерного распределения потребителей электрической мощности между ИП;
  3. Сокращения расстояний от ИП до потребителей;
  4. Оптимизации размещения ИП с учетом ограничений на местности;
  5. Оптимизации прокладки трасс линий с учетом ограничений на местности;
  6. Оптимизации конфигурации распределительной электрической сети петлевой структуры.

Апробация работы. Основные результаты по различным разделам диссертационной работы докладывались и обсуждались на 9 Всесоюзной научной конференции «Моделирование электроэнергетических систем», г.Рига, 1987 г.; на Всероссийских научно- технических конференциях «Повышение эффективности производства и использования энергии в условиях Сибири», г. Иркутск, 1990, 2005, 2007, 2009, 2011 гг.

Личный вклад автора. Результаты, составляющие новизну и выносимые на защиту, получены лично автором.

Публикации. По материалам диссертации опубликовано 16 работ, в том числе пять статей в журналах из списка изданий, рекомендованных ВАК.

Структура и объем работы. Диссертация состоит из введения, четырех глав, заключения, списка литературы, двух приложений. Работа представлена на 167 страницах машинописного текста, включает 22 рисунка, 42 таблицы. Библиографический список включает 102 наименования.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

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

В работе систематизирован комплекс задач по обоснованию развития систем электроснабжения в виде нескольких этапов.

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

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

  1. Определение оптимального расположения ИП в системе электроснабжения;
  2. Оптимального закрепления потребителей за ИП;
  3. Оптимального выбора мощности ИП;
  4. Оптимизация прокладки линий электропередачи, связывающей потребителя с ИП с учетом ограничений на местности;
  5. Оптимального соединения элементов системы электроснабжения в схемы определенной структуры.

2. На втором этапе уточняются найденные решения и формируется оптимальная схема электрической сети с учетом требований по качеству электрической энергии и надежности электроснабжения потребителей.

3. На третьем этапе проводят анализ функционирования сформированной системы электроснабжения путем расчетов различных режимов работы системы.

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

Известные методики (Каждан А.И., Каялов Г.М., Балабанян Г.А., Дульзон Н.А., Закиров Р.И., Гладилин Л.В. и др.), применяемые при решении вышеперечисленных задач и описанные в русскоязычной специальной литературе, не исчерпывают того многообразия проблем, которые приходится решать при синтезе рациональной конфигурации систем электроснабжения. Все постановки задач по оптимальному размещению ИП (с учетом и без учета ограничений на местности) в русскоязычной специальной литературе носят непрерывный характер, т.е. решение ищется на всей плоскости , и могут быть применены к решению задач в очень ограниченном числе случаев. Как правило, задача ставится в контексте размещения единичного ИП, для размещения нескольких ИП не существует математических моделей и методик решения задачи, позволяющих осуществить такое размещения в самом общем случае неравномерно распределенных электрических нагрузок и произвольной формы территории, на которой проектируется система электроснабжения.

Достаточно интенсивно ведутся исследования по синтезу оптимальной конфигурации распределительных электрических сетей за рубежом  (Khator S. K., Diaz-Dorado E., Aoki K., Navarro A., Ramirez-Rosado I. J. и др.), о чем свидетельствует большое количество публикаций по данной тематике. Уже в течение более сорока лет различные постановки задач оптимизации распределительных электрических сетей решаются с помощью самых разнообразных моделей. Показано, что эти задачи могут решаться как в непрерывной, так и дискретной постановках, к решению этих задач могут быть применены алгоритмы оптимизации на графах, эволюционное моделирование, генетические и эвристические алгоритмы. Нигде не решается эта проблема в общем виде: она всегда разбивается на множество частных подзадач.

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

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

Первый уровень -  питающая сеть систем электроснабжения городов включает в себя сети 35, 110 кВ, в крупных городах 220, 330 кВ, схема электроснабжения в этом случае зависит от размера города.

Второй уровень для систем электроснабжения городов включает в себя питающие сети 6, 10 кВ. Эти сети представляют собой линии электропередачи, отходящие от центров питания и подающие электроэнергию к распределительным пунктам (РП) без присоединения к ним других потребителей.

Третий уровень. В системах электроснабжения городов к указанному уровню относятся распределительные сети 6, 10 кВ, выполняемые в основном кольцевыми с трансформаторными подстанциями (ТП) транзитного типа.

Четвертый уровень включает в себя распределительную сеть напряжением 0,4 кВ, которая выполняется в различных модификациях: радиальная с присоединением потребителей на нерезервируемых ответвлениях; резервируемая и нерезервируемая; петлевая и магистральная. В городах замкнутые сети на напряжении 0,4 кВ не применяются.

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

1. На каждом территориальном уровне можно выделить элементы, которые играют роль источников питания, а также элементы, которые играют роль потребителей;

2. Можно выделить два способа соединения источника питания с потребителем электрической мощности (пусть координаты ИП, координаты потребителя, расстояние от ИП до потребителя), первый способ - на территории проектируемого объекта нет ограничений в виде прямоугольной сетки проездов, второй способ – линия, связывающая потребитель с ИП проходит параллельно осям координат:

а)- расстояние между ИП и потребителем вычисляется с помощью метрики Евклида,

б) расстояние между ИП и потребителем вычисляется с помощью Манхэттенской метрики;

3. Можно выделить различные способы соединения элементов в схему определенной структуры

       а) радиальная схема,

       б) магистральная схема,

       в) схема петлевой структуры,

       г) схема кольцевой структуры и т.д.

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

На рис.1 представлена формализация решаемых задач в зависимости от территориальных уровней системы электроснабжения городов.

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

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

Математически данная задача формулируется следующим образом: в

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

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

       (1)

Здесь - количество сторон в том выпуклом многоугольнике.

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

За расстояние между точками в  сети принимается Евклидова метрика. Тогда задача, требующая решения, формулируется следующим образом: необходимо найти такую точку с координатами , которая доставляла бы минимум  целевой функции (ЦФ) вида

                               (2)

на множестве .

Область определения целевой функции  представляет собой многосвязную область.

Рис. 1 Формализация решаемых задач в зависимости от территориальных уровней системы электроснабжения городов

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

Воспользуемся следующей схемой решения данной задачи.

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

       .                                        (3)

Минимум функции (2) находится как решение системы уравнений вида

.                        (4)

Для заданной функции можно показать, что в ее критической точке выполняются условия теоремы о достаточных условиях минимума функции

,

где , , .

Функция строго выпуклая функция, когда не все точки лежат на одной прямой, и в этом случае минимум ЦФ достигается в единственной точке.

Функция вида (2) представляет собой негладкую функцию, так как ее первые производные не являются непрерывными функциями в области определения функции (2). При (т.е. при совпадении искомой точки с одной из данных) в левых частях системы уравнений (4) имеем для одного из слагаемых неопределенность вида .

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

Выполнение условия Ламе :

       (5)

является необходимым и достаточным для того, чтобы одна из данных точек доставляла бы абсолютный минимум функции (2). В случае, когда условие Ламе не выполняется ни для одной из точек , система уравнений (4) решается каким- либо численным итерационным методом.

  1. После того, как решение найдено на всем множестве , проверяется, попадает ли найденная точка внутрь какого-либо многоугольника . Если нет, то найдено окончательное решение задачи.
  2. Пусть найденная точка попадает в многоугольник . Нужно найти минимум функции на каждом из отрезков, из которых образован выпуклый многоугольник , и из этих значений выбрать наименьшее.

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

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

Данная задача является классической комбинаторной задачей - определения числа способов размещения некоторых объектов в каком-то количестве ящиков так, чтобы были выполнены заданные ограничения. Формально эту задачу можно сформулировать так: даны множества , причем , ,  здесь - мощность множества , т.е. число его элементов. Сколько существует функций , удовлетворяющих заданным ограничениям? Элементы множества соответствуют объектам, элементы множества соответствуют ящикам, а каждая функция определяет некоторое размещение, указывая для каждого объекта ящик , в котором наш объект находится.

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

.                                                (6)

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

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

  1. Упорядочим  мощности (веса) потребителей в порядке убывания . Также упорядочим мощности используемых для данного варианта задачи ИП в порядке убывания  .
  2. Полагаем .
  3. Закрепляем -тый потребитель за ИП при  одновременном выполнении двух условий: а) кратчайшего расстояния от потребителя до ИП; б) суммарная мощность потребителей, закрепленных за ИП, не превосходит его мощности. Если ИП, для которых будут выполняться оба этих условия, несколько, то выбираем ИП с наименьшей мощностью. Здесь, говоря о мощности ИП, мы рассматриваем эту величину не постоянной, а изменяющейся динамически в процессе работы алгоритма; считаем, что мощность ИП на момент времени есть мощность ИП минус суммарная мощность потребителей, закрепленных за этим ИП на данном этапе выполнения алгоритма.
  4. Полагаем . Если , то останов в выполнении программы, иначе  переход на шаг .

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

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

Задача оптимального размещения нескольких ИП и закрепления за ними потребителей в системе электроснабжения. После того, как решена задача об оптимальном закреплении потребителей за источниками питания, можно перейти к более сложной общей задаче оптимального размещения нескольких источников питания в системе электроснабжения (указывается избыточное число мест возможного размещения ИП) и рационального закрепления за ними потребителей электрической мощности. В этой постановке возникает необходимость многократного решения задачи об оптимальном закреплении потребителей за фиксированными источниками питания. Это необходимо для вычисления целевой функции для каждого возможного варианта решения задачи. В содержательном плане данная постановка задачи также опосредовано учитывает ограничения генерального плана на выбор местоположения ИП: указывается избыточное число мест возможного размещения ИП (эти точки ищутся на генплане объекта проектирования с учетом площади, отводимой под установку ИП), число таких мест может быть достаточно большим, как правило, эти точки должны покрывать всю площадь проектируемого объекта.

Эта, более общая, задача состоит из следующего ряда подзадач:

  • выбор типоразмеров источников питания из предложенного типового ряда;
  • выбор местоположения отдельных источников питания;
  • рациональное закрепление потребителей электрической мощности за источниками питания.

Математическая постановка задачи. Пусть в декартовой системе координат (на плоскости) задана система точек c координатами, , которые будем именовать в дальнейшем точками потребления электроэнергии. Множество множество пунктов потребления электроэнергии.

Каждому потребителю сопоставляется некоторый вес , равный потребляемой мощности этим потребителем. Множество однозначно определяет объемы потребления. Таким образом, задана также суммарная потребляемая мощность всеми потребителями в рассматриваемой задаче.

                                               (7)

Пусть также в декартовой системе координат задана система точек с координатами возможных мест размещения источников питания, где количество этих возможных мест; априори предполагается, что это число заведомо больше практически необходимого числа мест размещения ИП. Так, например, в задаче следует определить не более мест размещения ИП, а возможных мест размещения можно указать несколько десятков. Множество множество пунктов возможного размещения источников питания.

Зададим типовой ряд используемых для решения задачи источников электрической мощности. Пусть этот типовой ряд представляет собой множество , где количество различных типоразмеров ИП. Очевидно, что суммарная отдаваемая мощность определяется как:

.                                                (8)

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

Пусть, например, суммарная потребляемая мощность в задаче равняется МВт, предполагаемый типовой ряд включает в себя источники питания следующих мощностей: МВт, МВт, МВт. При заданном типовом ряде  получаем, что МВт.

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

Здесь первый сомножитель в произведении указывает мощность в МВт источника питания, второй сомножитель указывает количество источников питания данного типоразмера. Как следует из описанных выше вариантов, в первом случае для размещения ИП понадобится возможных мест, во втором и т.д., наконец, в последнем варианте понадобится возможных места размещения ИП.

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

Обобщая задачу, исходя из найденной величины и заданных типоразмеров источников питания, найдем все возможные варианты электроснабжения при условии размещения в одной возможной точке одного ИП (это очень важное условие, в противном случае просто нужно добавить к используемому ряду другие ИП, мощность которых может быть удвоена, утроена и т.д. в зависимости от типоразмеров, используемых первоначально). Пронумеруем все эти варианты и получим в результате множество . Каждый элемент этого множества однозначно определяет количество мест размещения для каждого типоразмера ИП, участвующего в данном варианте электроснабжения.

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

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

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

Целевая функция, минимум которой будем находить, может быть представлена в виде:

.

(9)

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

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

Второй отличительной особенностью данной задачи является то, что возникает необходимость оптимизировать одновременно сразу несколько параметров, т.е. данная задача относится к классу задач многопараметрической оптимизации.

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

Как показывает опыт исследований, подобная постановка задачи позволяет формализовать поиск решения на основе теории «мягких вычислений», в частности на основе генетических алгоритмов.

Для любого генетического алгоритма (ГА) выделяется четыре основных этапа:

  1. формирование начальной популяции;
  2. синтез новых хромосом (операторы скрещивания и мутации);
  3. целенаправленное изменение вновь полученных хромосом (операторы инверсии);
  4. селекция текущей популяции.

Для выше сформулированной оптимизационной задачи был разработан  ГА, максимально адаптированный к специфике решаемой задачи. При этом были разработаны следующие его составные части и процедуры:

  1. способы кодировки возможных решений для двух случаев одинаковых типоразмеров используемых ИП и разных типоразмеров используемых ИП;
  2. алгоритм создания начальной популяции, учитывающий весь генетический материал задачи (все возможные места размещения ИП, все типоразмеры источников электрической мощности);
  3. алгоритм операции скрещивания для случая одинаковых типоразмеров источников питания;
  4. алгоритм операции скрещивания для случая разных типоразмеров источников питания;
  5. алгоритм вычисления целевой функции на основе эффективного эвристического алгоритма закрепления потребителей за источниками питания;
  6. разработанные алгоритмы реализованы в виде программ на встроенном языке программирования системы MATLAB.

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

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

Модифицируем известный алгоритм Краскала применительно к решению данной задачи. Исходные данные задачи:  трансформаторных подстанций, для каждой из них известно месторасположение ТП. Строим матрицу всевозможных расстояний от каждой ТП до каждой ТП. Так как эта матрица всевозможных расстояний является симметричной относительно главной диагонали, то для решения задачи следует взять элементы, стоящие в верхней треугольной матрице. При числе вершин в графе равном , число элементов в верхней треугольной матрице всевозможных расстояний вычисляется по формуле . Исходные данные зададим в виде трех одномерных массивов с числом элементов, равным . В первом массиве зададим номер вершины (начало ребра), во втором массиве номер вершины (конец ребра), в третьем – соответствующую длину ребра. Нумерацию ребер организуем по следующему правилу: номер вершины – начала ребра всегда меньше номера вершины – конца ребра. Этим  ребрам соответствует одномерный массив расстояний. Обозначим этот массив R. Далее упорядочим этот массив в порядке возрастания соответствующих длин. Преобразуем также первый и второй массивы, образующие соответствующие соединения. Соединяем две ТП с наиболее коротким соединяющим ребром. На каждом из следующих шагов добавляем самое короткое из оставшихся ребер , при присоединении которых к уже имеющимся ребрам не образовывалось бы никакого цикла, и при этом степень каждой вершины не превышала бы числа 2. Любое ребро, которое хотим присоединить к искомому минимальному остову, проверяем по двум условиям: первое – образуется ли с помощью этого ребра цикл в строящемся подграфе исходного графа, второе – приобретает ли какая-нибудь вершина в подграфе степень, равную трем. Если хотя бы одно из этих условий будет выполнено, то анализируемое ребро должно быть отброшено. Этот процесс продолжается до тех пор, пока впервые для трансформаторных подстанций не выберем ребро. Построенный по этому алгоритму остов графа будет оптимальным в смысле решаемой нами задачи.

Третья глава посвящена описанию комплекса программ по обоснованию рациональной конфигурации систем электроснабжения.  Данные программы реализованы в среде программирования пакета MATLAB.

Описание программной реализации генетического алгоритма. Данная программа, реализующая генетический алгоритм, состоит из основных функциональных модулей:

  • binIntObjects_MM.m- головная программа, осуществляет формирование исходных данных задачи, вызов основных функциональных модулей, управляет режимом расчета;
  • binFitness_MM.m- функциональный модуль, осуществляющий формирование текущей популяции, выполнение для каждой особи популяции процедуры оптимального закрепления потребителей за ИП, вычисление значения целевой функции для каждого возможного решения;
  • binCrossover_MM.m- функциональный модуль, осуществляющий процедуру скрещивания;
  • Binplot_MM.m- функциональный модуль, осуществляющий процедуру вывода результатов.

Программа Wprog_MM.m реализует  специальный алгоритм решения задачи нахождения минимума выпуклой функции в многосвязной области.

Программа Floid_MM.m, реализующая алгоритм Флойда нахождения кратчайшего пути между двумя вершинами связного неориентированного графа, предназначена для нахождения оптимальной трассы, связывающей потребителя с ИП.

Программа Kraskal_MM.m, реализует модифицированный алгоритм Краскала соединения подстанций в схему петлевой структуры.

В четвертой главе диссертации проведено исследование предложенных методов решения задачи обоснования рациональной конфигурации системы электроснабжения на примере системы электроснабжения микрорайона города (на рис.2 представлен его генплан).

В задаче рассмотрен жилой микрорайон города, в котором в основном преобладает жилая застройка, кроме этого в данном микрорайоне имеются школы, детские сады, пансионат и другие общественные здания. Жилые дома имеют пристроенные торговые точки, которые уже учтены при расчете нагрузок жилых домов. Всего на территории данного микрорайона находится 51 потребитель электрической мощности. Суммарная нагрузка потребителей составляет величину . Проведя предварительный анализ проектируемой системы электроснабжения, предлагается к установке 8 двухтрансформаторных подстанций  мощностью . Предполагается, что трансформаторные подстанции загружены равномерно. Мощность нагрузки для одной подстанции принимает значение

.

В качестве возможных мест размещения ТП выберем 30  возможных мест, указанных на генплане микрорайона. Проведем тестирование по вышеописанной программе, реализующей ГА при следующих параметрах алгоритма: число особей в начальной популяции число особей = 500, число

Рис. 2. Расположение потребителей электроэнергии и ТП на

территории микрорайона города

итераций генетического алгоритма число итераций = 100, . При заданных параметрах генетического алгоритма получили следующие результаты, представленные на рис.3, 4.

Наилучшее значение целевой функции при заданных параметрах генетического алгоритма: 1.0000e+000  7.8845e+005. Как видно из приведенных графиков, разработанный алгоритм выходит на установившиеся оптимальные значения в период выполнения первых 20-30 итераций (рис. 3).

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

Рис. 3. График изменения целевой  функции (число особей = 500, число итераций = 100,

первые 30 итераций)

Рис. 4.  Экран вывода результата:

Оптимальное закрепление потребителей

(число особей = 500, число итераций = 100, )

Проведем тестовый расчет по разработанной программе, изменив один параметр в исходных данных, а именно приняв мощность ИП При этом получили следующие результаты. График изменения целевой функции при параметрах: число особей = 500, число итераций = 100 представлен на рис.5, на рис.6 представлено оптимальное закрепление потребителей за выбранными ИП.

Наилучшее значение целевой функции при заданных параметрах генетического алгоритма: 1.0000e+000  7.5693e+005.

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

Рис. 5.  График изменения целевой  функции (число особей = 500, число итераций = 100,

первые 30 итераций)

Рис. 6.  Экран вывода результата:

Оптимальное закрепление потребителей

(число особей = 500, число итераций = 100,

Проведем еще один тестовый расчет по разработанной программе, приняв мощность ИП При этом получили следующие результаты. График изменения целевой функции при параметрах: число особей = 500, число итераций = 100 представлен на рис. 7, закрепление потребителей на рис. 8.

Рис. 7.  График изменения целевой  функции (число особей = 500, число итераций = 100, первые 60 итераций)

Рис. 8  Экран вывода результата:

Оптимальное закрепление потребителей

(число особей = 500, число итераций = 100, )

Наилучшее значение целевой функции при заданных параметрах генетического алгоритма: 1.0000e+000  8.4225e+005.

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

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

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

Касаясь коэффициентов загрузки ИП, можно сделать следующие выводы:

  • Если суммарная мощность ИП превышает суммарную мощность потребителей на 10-15 %, то в этом случае практически всегда при заданных мощностях потребителей можно закрепить их за ближайшими ИП, но в этом случае разбиение потребителей за ИП не будет оптимальным: некоторые ИП будут загружены до «предела», а некоторые будут иметь достаточно большую свободную мощность.
  • Если взять суммарную мощность ИП практически равной суммарной мощности потребителей, то в этом случае существует вероятность того, что вообще не будет найдено оптимальное решение (это зависит от нагрузок потребителей: чем больше нагрузка отдельных потребителей, тем больше эта вероятность) и найденное по данному алгоритму оптимальное решение будет в числовом выражении больше оптимального решения, полученного в предыдущем случае.
  • Если взять суммарную мощность ИП больше на 5-7 % суммарной мощности потребителей, то в этом случае получается достаточно хорошее компромиссное решение: целевая функция получает значение близкое к оптимальному значению и разбивка потребителей за ИП происходит более равномерно.

После того, как на территории жилого микрорайона было размещено 8 ТП, возникает задача соединения этих подстанций в схему петлевой структуры. Исходными данными в этой задаче являются координаты размещения ТП, найденные на предыдущем этапе. Воспользуемся для решения этой задачи программой Kraskal_MM.m, реализующей модифицированный алгоритм Краскала. Результат работы программы  представлен в табл. 1.

Таблица 1

Результат работы программы Kraskal _MM.m

Номер начальной вершины ребра

Номер конечной вершины ребра

Вес ребра

Степень нача- льной вершины

Степень конеч- ной вершины

5

6

159

1

1

1

2

186

1

1

4

8

188

1

1

Продолжение табл. 1

4

6

208

2

2

1

7

248

2

1

2

3

258

2

1

3

8

334

2

2

Визуально соединение заданных 8 ТП в схему петлевой структуры представлено на рис. 9. Пунктирной линией на рисунке показаны связи в сторону РП. Каждому ребру на рисунке сопоставлена его длина.

Рис. 9. Соединение трансформаторных подстанций

в схему петлевой структуры

На следующем территориальном уровне ставится задача размещения РП на напряжении 10 кВ для электроснабжения определенных групп потребителей электрической энергии. На территории жилого микрорайона указаны места, где размещение данного РП недопустимо. Ограничения на плоскости заданы в виде выпуклых многоугольников, каждый из которых описывается своей системой  линейных неравенств (либо координатами вершин соответствующих многоугольников). Визуально эти ограничения представлены на рис.10.

Проведем расчеты по программе Wprog_MM.m. Результат расчета представлен в табл. 2.

Таблица 2

Результат работы программы Wprog_MM.m

№ п.п.

Отрезок

Минимум ЦФ на заданном отрезке

Координаты точки минимума ЦФ(м)

1

Min З=2854286

2

Min З=2831330

3

Min З=3257336

4

Min З=3401509

Минимум целевой функции в допустимой области достигается на отрезке в точке и равен Min З=2831330.

Рис. 10.  Ограничения на плоскости для  размещения РП

После того, как была решена задача оптимального размещения РП с учетом ограничений на местности, возникает задача об оптимальной прокладке линии, связывающей потребителя электрической энергии, в данном случае - это ТП, с ИП, на данном уровне системы электроснабжения - это РП на напряжении 10 кВ.

На рис. 11 представлен план местности, разбитый на зоны, на рис.12 показан кратчайший путь, найденный с помощью алгоритма Флойда.

Рис. 11. Граф, представляющий план местности, разбитый на зоны

Рис. 12. Кратчайший путь на заданном графе

ЗАКЛЮЧЕНИЕ

В диссертационной работе получены следующие основные результаты:

  1. Разработана  и получила содержательное толкование иерархия задач по обоснованию развития систем электроснабжения городов;
  2. Проведена формализация комплекса задач оптимального выбора решений при обосновании рациональной конфигурации систем электроснабжения в зависимости от территориальных уровней;
  3. Разработаны математические модели решения следующих задач:
    • размещения единичного ИП в радиальной электрической сети с учетом произвольных ограничений на местности;
    • оптимального закрепления потребителей за несколькими ИП;
    • оптимального размещения нескольких ИП и закрепления за ними потребителей в системе электроснабжения (два случая: одинаковые типоразмеры ИП,  разные типоразмеры ИП);
  1. Разработаны методы решения следующих задач:
    • оригинальный метод решения задачи, нахождения минимума выпуклой функции в многосвязной области;
    • адаптированный эвристический алгоритм решения задачи оптимального закрепления потребителей за ИП, основанный на методах комбинаторного анализа;
    • метод решения задачи оптимального размещения нескольких ИП и закрепления за ними потребителей в системе электроснабжения, основанный на принципах эволюционного моделирования и генетического программирования (два случая:  одинаковые типоразмеры ИП,  разные типоразмеры ИП);
    • методы решения задач об оптимальном выборе конфигурации распределительной электрической сети с помощью алгоритмов оптимизации на графах:
        1. выбор оптимальной конфигурации распреде-лительной электрической сети петлевой структуры;
        2. построение оптимальной трассы линии электро-передачи, связывающей потребителя с ИП с учетом ограничений на местности;
  1. Выполнено исследование процесса размещения источников питания и закрепления за ними потребителей на примере решения задачи синтеза рациональной конфигурации распределительной электрической сети микрорайона города.

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

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

Научные статьи, опубликованные в изданиях по списку ВАК

  1. Коновалов Ю.С., Свеженцева О.В. Задача оптимального размещения источников питания в многосвязной области Свеженцева // Изв. ВУЗов СССР. Энергетика, 1987,  №6.
  2. Свеженцева, О.В. Оптимальное размещение источников питания в системах электроснабжения методом генетического программирования//

Вестник Иркутского государственного технического университета, 2005, №4.

  1. Свеженцева О.В. Решение задачи об оптимальном размещении источника питания с учетом ограничений на местности // Вестник Иркутского государственного технического университета, 2011,  №7.
  2. Свеженцева О.В. Разработка и тестирование генетического алгоритма оптимального размещения источников питания в распределительной электрической сети // Вестник Иркутского государственного технического университета, 2012, №4.
  3. Воропай Н. И., Свеженцева О.В. Оптимизация размещения источников питания при формировании рациональной конфигурации систем электроснабжения // Электричество, 2012, №10.

Публикации в других изданиях

  1. Коновалов Ю.С., Свеженцева О.В. Задача оптимального размещения источников питания в многосвязной области  // Тезисы докладов на 9-ой Всесоюзной научной конференции «Моделирование электроэнергетических систем», Рига: ФЭИ, 1987.
  2. Коновалов, Ю.С., Свеженцева О.В. Размещение источников электрической мощности в системах электроснабжения// Тезисы докладов на региональной научно-технической конференции «Повышение эффективности производства и использования энергии в условиях Сибири»,  Иркутск: ИПИ, 1990.
  3. Коновалов Ю.С., Свеженцева О.В. Оптимизация размещения источников питания в системах электроснабжения // Тезисы докладов на научно-практической конференции «Региональные проблемы повышения качества электрической энергии», Астрахань, 1991.
  4. Свеженцева О.В. Решение задачи оптимального размещения источников питания в распределительной электрической сети методами генетического программирования // Всероссийская научно-практическая конференция «Повышение эффективности производства и использования энергии в условиях Сибири»,  Иркутск: ИрГТУ,  2005.
  5. Свеженцева, О.В. Решение задачи оптимального закрепления множества потребителей за источниками питания методом комбинаторного анализа // Материалы научно-практической конференции «Технико-экономические проблемы развития регионов»,  Иркутск, 2006,  Вып. 3.
  6. Свеженцева, О.В. Формализация задачи оптимального закрепления множества потребителей за источниками питания // Всероссийская научно-практическая конференция «Повышение эффективности производства и использования энергии в условиях Сибири», Иркутск: ИрГТУ, 2007.
  7. Свеженцева, О.В. Эвристические алгоритмы оптимального закрепления множества потребителей за источниками питания // Всероссийская научно-практическая конференция «Повышение эффективности производства и использования энергии в условиях Сибири», Иркутск: ИрГТУ, 2007.
  8. Свеженцева, О.В. Моделирование электрической сети петлевой структуры  // Всероссийская научно-практическая конференция «Повышение эффективности производства и использования энергии в условиях Сибири», Иркутск: ИрГТУ, 2009.
  9. Свеженцева, О.В. Эвристический алгоритм закрепления множества потребителей за источниками питания // Всероссийская научно-практическая конференция «Повышение эффективности производства и использования энергии в условиях Сибири», Иркутск: ИрГТУ, 2009.
  10. Свеженцева, О.В. Построение оптимальной трассы линии электро-передачи с помощью алгоритмов оптимизации на графах  // Всероссийская научно-практическая конференция «Повышение эффективности производства и использования энергии в условиях Сибири», Иркутск: ИрГТУ, 2011.
  11. Voropai N. I., Svezhentseva O.V. Optimization of supply source allocation in the problem of rational configuration of electricity supply system // IEEE ISGT Europe Int. lonf., Berlin, Germany, October 14 – 17, 2012.

___________________________________________________________________

Отпечатано в ИСЭМ СО РАН

664033, Иркутск, ул. Лермонтова, 130.

Заказ_____. Тираж 100 экз.







© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.