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

СПОСОБ И СИСТЕМА ДЛЯ ВСЕСТОРОННЕГО АНАЛИЗА НИЗКОЧАСТОТНЫХ КОЛЕБАНИЙ Способ и система для всестороннего анализа низкочастотных колебаний. Способ включает в себя сначала активизацию широкопространственной системы (71) измерений выполнения контроля в реальном масштабе времени; периодическую передачу в реальном масштабе времени информации о низкочастотных колебаниях контролируемой системы в динамическую систему (72) раннего предупреждения; затем преобразование информации о низкочастотных колебаниях, полученной динамической системой (72) раннего предупреждения во входной файл задачи для расчета устойчивости с малыми возмущениями; и выполнение связанных с этим расчетов, передачу вычисленного результата назад в широкопространственную систему (71) измерений и на интерфейс"человек-машина" для выходного воспроизведения. Система содержит широкопространственную систему (71) измерений, динамическую систему (72) раннего предупреждения, систему (73) для расчета устойчивости с малыми возмущениями и выходную систему (74).(71)(73) Заявитель и патентовладелец: СТЕЙТ ГРИД КОРПОРЕЙШН ОФ ЧАЙНА (CN) По заявке на настоящий патент испрашивается приоритет по дате подачи заявки на патент Китая 200810168193.6, поданной 27 сентября 2008 г. в Патентное ведомство Китая, озаглавленной "Способ и система для всестороннего анализа низкочастотного колебания", которая во всей своей полноте включена здесь в качестве ссылки. Область техники, к которой относится изобретение Настоящее изобретение относится к области анализа устойчивости энергетической системы и, в частности, к способу и системе для всестороннего анализа низкочастотных колебаний в энергетической сети. Уровень техники В современной энергетической системе с увеличением масштаба единой энергетической сети возрастает объем передачи энергии на дальние расстояния, при этом в энергетической сети Китая часто возникают низкочастотные колебания. Критическим вопросом, который должен быть решен в крупномасштабной единой энергетической сети, является быстрая и точная идентификация низкочастотных колебаний в сети, а также выполнение онлайновых измерений этих колебаний в реальном масштабе времени посредством анализа с малыми возмущениями. Обычным и классическим способом анализа низкочастотных колебаний является способ анализа собственных значений с малыми возмущениями. Для ожидаемого режима работы системы информация, включающая в себя подгруппы низкочастотных колебаний, взаимосвязь колебаний, а также коэффициенты корреляции и т.д., извлекается из анализа собственных значений и собственных векторов уравнения состояния рабочих точек. Способ анализа устойчивости с малыми возмущениями обычно позволяет производить анализ в автономном режиме. Однако режим колебаний энергетической системы обычно с изменением рабочих условий изменяется, поэтому точные онлайновые вычисления не могут быть обеспечены на основании результата, полученного в автономном режиме. Средством контроля динамических процессов в энергетической системе является широко пространственная система измерений (WAMS). Она может производить быстрые измерения переменных, которые тесно связаны с переходными электромеханическими процессами генератора, такими как внутренний потенциал, угол мощности, угловая скорость, шинное напряжение и т.д., и передавать эту информацию на централизованную станцию для обеспечения возможности выполнения онлайнового анализа и управления низкочастотными колебаниями. Однако в настоящее время работа широко пространственной системы измерений, а также документация, касающаяся детекции низкочастотных колебаний на основе этой системы, ограничены лишь функцией динамического контроля. Т.е. динамические кривые относительного фазного угла напряжения энергетической сети, частоты и мощности постоянно контролируются для проведения в реальном масштабе времени вычислений и анализа спектра этих динамических кривых. При детекции наличия относительно сильной компоненты слабых демпфирующих колебаний в диапазоне от 0,2 до 2,5 Гц управляющему оператору передается предупредительная информация, при этом на региональной карте энергетической сети отмечается эта аномальная зона и активизируется таблица технических данных для высокоскоростной записи в реальном масштабе времени соответствующих параметров. Однако по финансовым соображениям невозможно на все узлы энергетической системы установить блоки измерения вектора. Обычно осцилляционные характеристики мощности на выходном конце линии или станции получаются на основании функции контроля широко пространственной системы измерений, поэтому невозможно определить, какой конкретный генератор вызывает колебания. Т.е. невозможно принятие никаких управляющих мер против колебаний ввиду отсутствия дополнительной аналитической информации относительно механизма низкочастотных колебаний. Сущность изобретения Задачей данного изобретения является обеспечение способа и системы для всестороннего анализа низкочастотных колебаний в комбинации с анализом с малыми возмущениями, а также проведения измерений в реальном масштабе времени и идентификации энергетической сети. Для достижения этой цели настоящее изобретение обеспечивает способ всестороннего анализа низкочастотных колебаний в комбинации с анализом с малыми возмущениями, а также проведения измерений в реальном масштабе времени и идентификации энергетической сети, который включает в себя: 1) активизацию широко пространственной системы измерений на предмет контроля динамических процессов в энергетической системе в реальном масштабе времени; 2) периодическую передачу широко пространственной системой измерений информации о низкочастотных колебаниях контролируемой системы в реальном масштабе времени в динамическую систему раннего предупреждения; 3) прием динамической системой раннего предупреждения информации о низкочастотных колебаниях, переданной широко пространственной системой измерений, и преобразование ее во входной файл задачи для расчета устойчивости с малыми возмущениями; 4) выполнение расчета устойчивости с малыми возмущениями для получения собственного значения и выдачи результирующего файла расчета устойчивости с малыми возмущениями; 5) вычисление коэффициента демпфирования и частоты низкочастотных колебаний в реальном масштабе времени по содержанию вычисленного результирующего файла для создания результирующе-1 018309 го файла онлайновых вычислений с малыми возмущениями; и 6) возврат результата онлайновых вычислений с малыми возмущениями в широко пространственную систему измерений и на интерфейс "человек-машина" для выходного воспроизведения. Кроме того, настоящее изобретение обеспечивает систему для всестороннего анализа низкочастотных колебаний, которая включает в себя: широко пространственную систему измерений, выполненную с возможностью контроля динамических процессов в энергетической системе в реальном масштабе времени и периодической передачи информации о низкочастотных колебаниях контролируемой системы в реальном масштабе времени в динамическую систему раннего предупреждения; динамическую систему раннего предупреждения, связанную с широко пространственной системой измерений и выполненную с возможностью приема контролируемой информации о низкочастотных колебаниях в реальном масштабе времени, переданной широко пространственной системой измерений, и преобразования ее во входной файл задач для вычисления устойчивости с малыми возмущениями; систему для вычисления устойчивости с малыми возмущениями, связанную динамической системой раннего предупреждения и выполненную с возможностью выполнения вычисления устойчивости с малыми возмущениями, выдачи результирующего файла вычисления устойчивости с малыми возмущениями и вычисления коэффициента демпфирования и частоты низкочастотных колебаний в реальном масштабе времени по вычислению результирующего файла для создания результирующего файла онлайновых вычислений с малыми возмущениями; и выходную систему, выполненную с возможностью выдачи результата онлайновых вычислений с малыми возмущениями. Данный патент объединяет в себе детальные преимущества динамического контроля параметров,производимого широко пространственной системой измерений, и модель системы для анализа устойчивости с малыми возмущениями. Сначала обычный расчет устойчивости с малыми возмущениями основан на имитации офлайновых параметров, а ввод динамических параметров, полученных от широко пространственной системы измерений, может эффективно исключить ошибки данных вследствие моделирования динамических элементов. Кроме того, обеспечиваемая широко пространственной системой измерений соответствующая информация по частоте может способствовать быстрой локализации диапазона поиска низкочастотных колебаний и исключению проблем, связанных со значительным объемом вычислительных работ, загрузкой памяти и медленной скорости вычислений вследствие вычисления всех собственных величин. Далее, по экономическим ограничениям в широко пространственной системе измерений невозможно произвести конфигурацию всех блоков измерения вектора на всех станциях сети, но обычно внимание обращается на линии передач и каналы передачи энергии. Даже если возникшее в линии низкочастотное колебание будет детектировано, все еще трудно принять соответствующие меры управления в отсутствии дополнительной информации, касающейся механизма колебаний. Вышеупомянутые проблемы могут быть устранены вводом расчета устойчивости с малыми возмущениями. Расчет устойчивости с малыми возмущениями может дать детальную аналитическую информацию относительно собственных величин, соответствующих частоте колебаний и группе (группам) генераторов, к которым будет привлечено внимание. Очевидно, расчет устойчивости с малыми возмущениями в комбинации с функцией контроля широко пространственной системы измерений дает возможность производить онлайновое исследование причин возникновения колебаний, влияющих на колебания первичных факторов и способа эффективного подавления этих колебаний, а также дает более полную информацию для предупредительных мер и анализа низкочастотных колебаний системы. Краткое описание чертежей Фиг. 1 представляет собой блок-схему, иллюстрирующую способ всестороннего анализа низкочастотных колебаний в соответствии с изобретением. Фиг. 2 представляет собой блок-схему, иллюстрирующую способ преобразования информации о низкочастотных колебаниях, переданной широко пространственной системой измерений, во входной файл задачи для расчета устойчивости с малыми возмущениями на этапе 3) способа всестороннего анализа низкочастотных колебаний в соответствии с изобретением. Фиг. 3 представляет собой блок-схему, иллюстрирующую расчет устойчивости на этапе 4) способа всестороннего анализа низкочастотных колебаний в соответствии с изобретением. Фиг. 4 представляет собой условную диаграмму, иллюстрирующую диапазон поиска набора собственных величин в соответствии с изобретением. Фиг. 5 представляет собой условную блок-схему, иллюстрирующую систему для всестороннего анализа низкочастотных колебаний в соответствии с изобретением. Подробное описание изобретения Далее будут подробно описаны вместе с чертежами предпочтительные варианты исполнения настоящего изобретения. Фиг. 1 иллюстрирует рабочую последовательность способа всестороннего анализа низкочастотных колебаний в комбинации с анализом с малыми возмущениями, а также измерениями в реальном масштабе времени и идентификацией энергетической сети в соответствии с изобретением. Эта рабочая последовательность включает в себя следующие этапы:S101 - активизация широко пространственной системы измерений на предмет контроля динамических процессов в энергетической системе;S102 - периодическая передача широко пространственной системой измерений в реальном масштабе времени информации о низкочастотных колебаниях контролируемой системы в динамическую систему раннего предупреждения;S103 - прием динамической системой раннего предупреждения информации о низкочастотных колебаниях, переданной широко пространственной системой измерений, и преобразование ее во входной файл задачи для расчета устойчивости с малыми возмущениями;S104 - выполнение расчета устойчивости с малыми возмущениями для получения собственной величины и выдачи результирующего файла расчета устойчивости с малыми возмущениями;S105 - вычисление отношения демпфирования и частоты низкочастотных колебаний в реальном масштабе времени по содержанию результирующего файла вычислений для создания результирующего файла онлайновых вычислений с малыми возмущениями;S106 - возврат результата онлайновых вычислений с малыми возмущениями в широко пространственную систему измерений и на интерфейс "человек-машина" для выходного воспроизведения. При этом информация о низкочастотных колебаниях контролируемой системы в реальном масштабе времени на этапе S102 включает в себя информацию о частоте колебаний, о линии с колебаниями и т.д. в соответствии с содержанием табл. 1. Таблица 1 Информация о низкочастотных колебаниях, передаваемая широко пространственной системой измерений Далее, как показано на фиг. 2, этап S103 преобразования информации о низкочастотных колебаниях, переданной широко пространственной системой измерений, во входной файл задачи для расчета устойчивости с малыми возмущениями, включает в себя следующие этапы: 30) установка действительных составляющих начальной и конечной точек собственного значения на ноль; 31) ввод диапазона поиска частоты колебаний; 32) ввод порядкового номера алгоритма, допустимой погрешности, количества начальных значений итераций, величины шага вычисления, общего времени вычисления, количества полученных за каждый раз собственных величин, а также верхнего предела количества итераций во входном файле с малыми возмущениями как значений по умолчанию. При этом диапазон поиска частоты колебаний обычно выше или равен 0,5f и меньше или равен l,5f. Входной файл SST.e0 вычисления с малыми возмущениями этапа 32) обычно включает в себя порядковый номер алгоритма, допустимую погрешность, действительную и мнимую составляющую начальной и конечной точек, количество итерационных начальных значений, величину шага вычисления,общее время вычисления, количество полученных за каждый раз собственных значений, а также верхний предел количества итераций. Порядковый номер алгоритма может быть выбран из предварительно введенного в систему перечня способов вычисления собственных значений (представленных в табл. 2). По умолчанию для получения собственных значений система может автоматически выбрать "метод параллельных итераций", т.е. порядковый номер алгоритма есть 2. Более того, по умолчанию, допустимая погрешность, количество итерационных начальных значений, величина шага вычисления, общее время вычисления, количество полученных за каждый раз собственных значений, а также верхний предел количества итераций во входном файле SST.e0 могут быть обеспечены значениями по умолчанию без ввода каким-либо иным образом, и только информация, касающаяся действительной и мнимой составляющей начальной и конечной точек во входном файле SST.e0, будет модифицироваться перед каждой активизацией программы расчета устойчивости с малыми возмущениями. В предпочтительном варианте исполнения настоящего изобретения содержания входного файла SST.e0 составляют 2; 8; -0,000000; 1,5f;-0,000000; 0,5f; 3; 0,010000; 10,000000; 2; 200. Таблица 2 Перечень методов вычисления собственных значений Как указано в этапе S104, диапазон поиска собственных значений определен из созданного на этапеS103 входного файла SST.e0, так что для поиска собственного значения выполняется расчет устойчивости с малыми возмущениями, как показано на фиг. 3. Значокна фиг. 3, введенный самой программой,определяет диапазон [0,01; 0,05]. В предпочтительном варианте исполнения настоящего изобретения авторы получают собственное значение методом параллельных итераций следующим образом. Взята матрица АCnn с собственными величинами и собственными векторами соответственно i и q1, где i=1, 2, , n, и где Матрица прямых собственных векторов относительно матрицы А имеет вид где q1 есть собственный вектор, соответствующий собственному значению i. Число m взаимно независимых начальных векторов полагается равным Над вышеупомянутыми начальными векторами производится следующая операция: Поскольку U может быть преобразовано в Относительная величина первого члена в правой стороне уравнения (2) больше, чем его относительная величина в уравнении (1). Другими словами, компоненты в соответствующих столбцах матрицыV, относящиеся к Qb, в некоторой степени подавлены. Собственное значение в а далее получается в результате последующего процесса итерации. Сначала дается определение: где верхний индекс H представляет собой операцию сопряженного транспонирования комплексной матрицы. Произведение UHQa полагаем невырожденным, тогда Другими словами, если матрица В - это решение выражения GB=H, тогда Предыдущее уравнение показывает, что матрицы а и Са включают в себя приблизительные собственные значения и оставшиеся собственные вектора матрицы В. Если матрица Р это матрица прямых собственных векторов матрицы В, т.е. РСа-1, тогда В предыдущем уравнении W дает более точно вычисленный набор прямых собственных векторов матрицы А. Полагаем, что набор начальных векторов U=W, и повторяем предыдущее вычисление до схождения результата. Метод параллельной итерации может позволить каждый раз получать одновременно несколько первичных собственных значений и собственных векторов. Более того, расчет устойчивости с малыми возмущениями на этапе 4 позволяет получить также параметры, включая отношение корреляции электромеханического контура, наименование шины генератора, коэффициент корреляции, прямой собственный вектор и т.д., относящиеся к настоящему изобретению. Отношение корреляции электромеханических контуров. Коэффициент корреляции 1 электромеханических контуров собственной величины i определяется как Коэффициент корреляции i электромеханических контуров отражает степень, с которой собственная величина i связана с переменнымии . На практике, если для собственной величины i нижеследующая система справедлива: тогда i находится в режиме низкочастотного колебания, т.е. в электромеханическом режиме. Коэффициент корреляции. Коэффициент корреляции pk1 это физическая величина, которая является мерой корреляции k-го состояния переменной xk с i-й собственной величиной i: Коэффициент корреляции pki есть всеобъемлющий показатель, который отражает управляемость ki и различимость uki переменной xk и собственной величиной i. На практике коэффициент корреляции pki действует в качестве определяющего принципа при выборе места для установки PSS и в значительной степени может отражать переменную состояния, генератор которой тесно связан с ее колебательным режимом, тем самым подавляя соответствующий колебательный режим при предпочтительной установкеPSS на этом генераторе. Мода колебания и форма моды. Для матрицы ACnn ее собственные величины (i) и собственные векторы (ui) удовлетворяют следующему уравнению: Устойчивость системы может быть определена по собственной величине:(1) действительная собственная величина соответствует режиму отсутствия колебаний; а(2) комплексная собственная величина всегда проявляется в форме сопряженной пары, т.е. =j. Действительная составляющая собственной величины указывает противодействие системы колеба-5 018309 ниям, а ее мнимая составляющая указывает частоту колебания. 0 означает колебание с увеличивающейся амплитудой, при котором система неустойчива; 0 означает колебание с уменьшающейся амплитудой, при котором система устойчива; и =0 означает колебание с постоянной амплитудой, при котором система находится в состоянии критической устойчивости. Собственные векторы u1 и u2 отражают соотношение между величиной относительной амплитуды и фазой, когда наблюдается соответствующее колебание вектора состояния X. Физически, пара сопряженных собственных величин называется модой колебания, и их соответствующий собственный вектор называется формой моды колебания. Физическое значение правого собственного вектора. Определенный в уравнении (5) собственный вектор (ui) называется также правым собственным векявляется линейно преобразованным с использованием матрицы,тором. Уравнение состояния составленной из собственных векторов ui, и может быть разделено. Собственный вектор ui, соответствующий собственной величине i (i=1, 2, , n), отражает относительную амплитуду и фазу моды i, наблюдаемой на соответствующем векторе состояния. Чем больше модуль uki, тем более значительно соотношение между xk и i, что тем самым отражает различимость xk относительно i. На основании этого свойства прямого собственного вектора авторы изобретения смогли непосредственно из формы моды колебания, соответствующей моде колебания i (прямой собственный вектор ui), сделать вывод, что мода колебания i соответствует неустойчивой моде среди групп генераторов. Определение и физическое значение левого собственного вектора. Вектор (i), удовлетворяющий приведенному ниже уравнению, называется левым собственным вектором: Т.е. i это правый собственный вектор той же самой собственной величины i матрицы АТ и может быть получен на основании этого свойства. Левый и правый собственные векторы удовлетворяют следующему соотношению: Левый собственный вектор отражает управляемость вектора состояния относительно собственной величины. На этапе S104 по результатам расчета устойчивости с малыми возмущениями могут быть созданы два выходных файла, т.е. соответственно sst.eg1 и sst.eg2, в которых: содержание выходного файла sst.eg1 включает в себя порядковый номер собственного значения,собственное значение и отношение корреляции электромеханических контуров; содержание выходного файла sst.eg2 включает в себя порядковый номер собственного значения,наименование шины генератора, коэффициент корреляции и правое собственное значение. Далее, как указано в этапе S105, по содержанию вычисленных результирующих файлов sst.eg1 иsst.eg2 производится вычисление отношения демпфирования и частоты колебаний для создания результирующего файла онлайновых вычислений с малыми возмущениями PMU.OUT; структура этого файла указана в табл. 3. Результирующий файл PMU.OUT содержит подробную информацию относительно отношения демпфирования, о группе генераторов колебаний подгрупп и т.д. в соответствии с воспринимаемой частотой, измеренной широко пространственной системой измерений. Таблица 3 Результаты онлайновых вычислений с малыми возмущениями Отношение демпфирования. Низкочастотные колебания обычно возникают в слабокоррелированной системе большой энергетической сети, и такая дестабилизация обычно является результатом слабого демпфирования и отрицательного демпфирования, причем степень демпфирования поддается определению, исходя из отношения демпфирования нескольких доминантных мод колебаний. Полагаем, что все собственные величины системы имеют вид Частота колебаний равна fi=i/2, и отношение демпфирования (коэффициент демпфирования) i определяется как где i 0,1 указывает на сильное демпфирование системы;i0,03 указывает на слабое демпфирование системы иi0 указывает на отрицательное демпфирование системы, что может привести к возникновению колебаний с увеличивающейся амплитудой. Далее, как указано в этапе S106, результирующий файл PMU.OUT, полученный из онлайновых вычислений с малыми возмущениями, возвращается и в широко пространственную систему измерений, и в В варианте исполнения по настоящему изобретению широко пространственная система измерений в каждый момент времени обеспечивает указанную в табл.4 информацию раннего оповещения о низкочастотных колебаниях в объединенной энергетической сети (извлекается только указанная в табл. 1 информация о колебаниях в соединительных линиях и о частоте колебаний, необходимая для проведения вычислений с малыми возмущениями). Таблица 4 Информация по низкочастотным колебаниям, передаваемая широко пространственной системой измерений Вышеприведенная информация о частоте записывается в файл SST.e0, а вышеуказанные частотные диапазоны берутся в качестве объектов детального онлайнового анализа методом онлайнового анализа с малыми возмущениями на предмет анализа собственных величин и получения результата вычислений с малыми возмущениями (см. табл. 5). Результат вычислений демонстрирует явление низкочастотных колебаний, возникающих в вышеуказанных областях (включая внутри- и "межпровинциальные" колебания), и, кроме того, указывает первичные группы "участвующих" генераторов, вызывающих низкочастотные колебания, моды колебаний и отношения демпфирования колебаний соответствующих групп генераторов. Таблица 5 Результаты анализа и вычислений с малыми возмущениями Как показано на фиг. 5, система для всестороннего анализа низкочастотных колебаний в комбинации с анализом с малыми возмущениями, а также с измерениями в реальном масштабе времени и идентификацией энергетической сети в соответствии с настоящим изобретением включает в себя: широко пространственную систему 71 измерений, адаптированную для контроля в реальном масштабе времени динамических процессов в энергетической системе и периодической передачи "снятой" таким образом информации о низкочастотных колебаниях в реальном масштабе времени; динамическую систему 72 раннего предупреждения, связанную с широко пространственной системой 71 измерений и адаптированную для приема информации о низкочастотных колебаниях в реальном масштабе времени, переданной широко пространственной системой 71 измерений, и преобразование ее во входной файл задачи для расчета устойчивости с малыми возмущениями; систему 73 для расчета устойчивости с малыми возмущениями, связанную с динамической системой 72 раннего предупреждения и адаптированную для проведения расчета устойчивости с малыми возмущениями для выдачи результирующего файла расчета устойчивости с малыми возмущениями и для вычисления отношения демпфирования и частоты колебаний низкочастотных колебаний в реальном масштабе времени по вычисленному результирующему файлу для создания результирующего файла онлайновых вычислений с малыми возмущениями; и выходную систему 74, адаптированную для вывода результата онлайновых вычислений системы 73 для расчета устойчивости с малыми возмущениями. ФОРМУЛА ИЗОБРЕТЕНИЯ 1. Способ всестороннего анализа низкочастотных колебаний, включающий этапы:S101 активизации широко пространственной системы измерений (WAMS) на предмет контроля динамических процессов в энергетической системе;S102 периодической передачи широко пространственной системой измерений в реальном масштабе времени информации о низкочастотных колебаниях контролируемой системы в динамическую систему раннего предупреждения;S103 приема динамической системой раннего предупреждения информации о низкочастотных колебаниях, переданной широко пространственной системой измерений, и преобразование ее во входной файл задачи для расчета устойчивости с малыми возмущениями;S104 выполнения расчета устойчивости с малыми возмущениями для получения собственной величины и выдачи результирующего файла расчета устойчивости с малыми возмущениями;S105 вычисления отношения демпфирования и частоты низкочастотных колебаний в реальном масштабе времени по содержанию результирующего файла вычислений для создания результирующего файла онлайновых вычислений с малыми возмущениями;S106 возврата результата онлайновых вычислений с малыми возмущениями в широко пространственную систему измерений и на человекомашинный интерфейс для выходного воспроизведения. 2. Способ всестороннего анализа низкочастотных колебаний по п.1, в котором информация в реальном масштабе времени о низкочастотных колебаниях контролируемой системы этапа S102 содержит информацию о частоте колебаний и о линии, в которых колебание происходит. 3. Способ всестороннего анализа низкочастотных колебаний по п.2, в котором преобразование во входной файл задачи для расчета устойчивости с малыми возмущениями на этапе S103 также содержит: 30) установку действительных составляющих начальной и конечной точек собственного значения на ноль; 31) ввод диапазона поиска частоты колебаний; 32) ввод порядкового номера алгоритма, допустимой погрешности, количества начальных значений итераций, величины шага вычисления, общего времени вычисления, количества полученных за каждый раз собственных величин, а также верхнего предела количества итераций во входном файле с малыми возмущениями как значений по умолчанию. 4. Способ всестороннего анализа низкочастотных колебаний по п.3, в котором диапазон поиска частоты больше или равен 0,5f и меньше или равен l,5f. 5. Способ всестороннего анализа низкочастотных колебаний по п.4, в котором этап 4) также содержит этапы: 40) определения диапазона поиска собственной величины; 41) поиска собственной величины; 42) вычисления отношения корреляции электромеханических контуров, коэффициента корреляции и прямого собственного вектора; 43) создания и вывода результирующего файла расчета устойчивости с малыми возмущениями. 6. Система для всестороннего анализа низкочастотных колебаний, содержащая широко пространственную систему измерений, адаптированную для контроля в реальном масштабе времени динамических процессов в энергетической системе и периодической передачи полученной таким образом информации о низкочастотных колебаниях в реальном масштабе времени; динамическую систему раннего предупреждения, связанную с широко пространственной системой измерений и адаптированную для приема информации о низкочастотных колебаниях в реальном масштабе времени, переданной широко пространственной системой измерений, и преобразование ее во входной файл задачи для расчета устойчивости с малыми возмущениями; систему для расчета устойчивости с малыми возмущениями, связанную с динамической системой раннего предупреждения и адаптированную для проведения расчета устойчивости с малыми возмущениями для выдачи результирующего файла расчета устойчивости с малыми возмущениями и для вычисления отношения демпфирования и частоты колебаний низкочастотных колебаний в реальном масштабе времени по вычисленному результирующему файлу для создания результирующего файла онлайновых вычислений с малыми возмущениями; выходную систему, адаптированную для вывода результата онлайновых вычислений для расчета устойчивости с малыми возмущениями.