Способ синхрофазорного измерения для использования в устройстве измерения фазоров (pmu) р-класса

Изобретение относится к способу синхрофазорного измерения для использования в устройстве измерения фазоров (PMU) Р-класса. Упомянутый способ измерения основывают на математической модели динамического фазора. Конструируют цифровой фильтр низких частот для коэффициентов фазора, объединенный с DFT. Этот фильтр устраняет утечку спектра, вызванную входными сигналами динамического фазора, причем после устранения утечки спектра могут провести измерения исходного фазора. Динамический фазор аппроксимируют с использованием ряда Тейлора второго порядка. Исследуют линейную зависимость между ошибками измерения, вызванными усредняющим эффектом DFT, и коэффициентами ряда Тейлора второго порядка. Затем используют упомянутую линейную зависимость для компенсации исходных ошибок измерения в динамических условиях. Техническим результатом при реализации заявленного способа измерения является возможность точного и быстрого измерения фазора как в статических, так и в динамических условиях. Точность упомянутого способа измерения не только удовлетворяет техническим требованиям соответствующих стандартов, но и на порядок превышает требования этих стандартов. 4 з.п. ф-лы, 3 табл., 7 ил.

 

ОБЛАСТЬ ТЕХНИКИ, К КОТОРОЙ ОТНОСИТСЯ ИЗОБРЕТЕНИЕ

Настоящее изобретение относится к области техники, связанной с измерением фазоров, в частности к разновидности способа синхрофазорного измерения, применяемого к устройству измерения фазоров (PMU) Р-класса.

УРОВЕНЬ ТЕХНИКИ

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

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

С выпуском и постепенным усовершенствованием стандартов IEEE и Китая для PMU точность измерений PMU в динамических условиях получает признание все большего количества исследовательских организаций, при этом в способе измерения применяют много новых технологий. Например, в некоторых литературных источниках предлагается моделировать трехфазный фазор способом смещаемого окна времени, чтобы устранить утечку спектра DFT в условиях отклонения частоты. Однако этот способ подходит только для статического входного сигнала и не может решать задачу утечки спектра, которая вызвана динамическими входными сигналами. Таким образом, в существующем техническом решении отсутствует единый алгоритм оценки, основанный на модели динамического фазора. Кроме того, существующие алгоритмы не могут решить проблемы утечки спектра и усредняющего эффекта, которые вызваны динамическими входными сигналами.

РАСКРЫТИЕ ИЗОБРЕТЕНИЯ

Задача данного изобретения заключается в создании способа синхрофазорного измерения для использования в устройстве измерения фазоров (PMU) Р-класса. Этот способ измерения обеспечивает возможность точного и быстрого измерения фазора как в статических, так и в динамических условиях. Точность способа не только удовлетворяет техническим требованиям соответствующих стандартов, но и на порядок выше требований стандартов.

Способ синхрофазорного измерения для использования в устройстве измерения фазоров (PMU) Р-класса выполняют на основе математической модели динамического фазора, при этом конструируют цифровой фильтр низких частот для коэффициентов фазора, объединенный с DFT и устраняющий утечку спектра, вызванную входными сигналами динамического фазора, причем после устранения утечки спектра могут проводить измерения исходного фазора.

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

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

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

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

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

Упомянутый способ измерения после выполнения точных измерений динамического фазора дополнительно включает в себя:

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

Частота среза цифрового фильтра низких частот, объединенного с дискретным преобразованием Фурье (DFT), составляет 2,5 Гц, при этом окно времени составляет два номинальных периода динамического входного сигнала и расчетная частота исходного фазора составляет 400 Гц.

Частота среза второго цифрового фильтра низких частот составляет 2,5 Гц, при этом окно времени составляет 27,5 миллисекунд.

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

КРАТКОЕ ОПИСАНИЕ ЧЕРТЕЖЕЙ

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

На фиг. 1 показана блок-схема способа синхрофазорного измерения для использования в устройстве измерения фазоров (PMU) Р-класса, который выполняют в соответствии с вариантом осуществления изобретения.

На фиг. 2 показано соответствие ошибок измерения амплитуды и частот, которое получено в соответствии с вариантом осуществления изобретения.

На фиг. 3 показана максимальная суммарная векторная ошибка (TVE) для различных основных частот при моделировании теста частотного сканирования, которая получена в соответствии с вариантом осуществления изобретения.

На фиг. 4 показаны максимальные ошибки частоты и скорости изменения частоты (ROCOF) для различных основных частот при моделировании теста частотного сканирования, полученные в соответствии с вариантом осуществления изобретения.

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

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

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

ОСУЩЕСТВЛЕНИЕ ИЗОБРЕТЕНИЯ

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

Упомянутый вариант осуществления настоящего изобретения подробно описан ниже с привлечением чертежей. На фиг. 1 изображена блок-схема способа синхрофазорного измерения для использования в устройстве измерения фазоров (PMU) Р-класса, который выполняют в соответствии с вариантом осуществления изобретения. Упомянутый способ включает в себя:

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

На этом этапе формулируют математическую модель динамического фазора. Для выполнения вычислений предложена математическая модель динамического фазора, основанная на ряде Тейлора второго порядка. В этой модели все параметры фазора, включая амплитуду, угол сдвига фаз, частоту и ROCOF, изменяются с течением времени в расчетном окне времени.

Например, в этом случае прежде всего выражают динамический сигнал основной частоты в виде:

где x(t) - текущее значение сигнала, Xm(t) - амплитуда, f(t) - частота сигнала, ϕ0 - исходный угол сдвига фаз, f0 - номинальная частота, Δf(t) - отклонение частоты. Амплитуда и частота сигнала являются функцией времени.

Выражение (1) можно перезаписать следующим образом:

где Re - вещественная часть формулы. В этом случае в вышеприведенном выражении обычно можно пренебречь величиной , учитывая при этом, что базовая система координат синхронно вращается с угловой частотой 2πf0. Тогда формула (1) выражается величиной и формула математической модели динамического фазора представляет собой:

ROCOF имеет вид

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

где m2=d2Xm(t)/dt2|t=0, m1=dXm(t)/dt|t=0, m0=Xm(0); n2=d2ϕ(t)/dt2|t=0, n1=dϕ(t)/dt|t=0, n0=ϕ(0); p2=d2f(t)/dt2|t=0, p1=df(t)/dt|t=0, p0=f(0); q2=d2ROCOF(t)/dr2|t=0, q1=dROCOF(t)/dt|t=0, q0=ROCOF(0).

Затем, на основе математической модели динамического фазора, анализируют явление утечки спектра, вызванное DFT в случае динамического фазора, и его влияние на измерения фазора. Предложено DFT с цифровым фильтром низких частот для устранения утечки спектра и получения исходного фазора ( и ). Частота среза цифрового фильтра низких частот, объединенного с DFT, составляет 2,5 Гц, при этом окно времени составляет 2 номинальных периода входного сигнала динамического фазора и расчетная частота исходного фазора составляет 400 Гц.

Например, формула стандартного алгоритма DFT в отношении сигнала электрического тока выглядит как

где - измеренный фазор, N - объем выборки, x(k) - k-е значение выборки, f0 - номинальная частота и Δt- интервал выборки.

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

Вещественную и мнимую части γ следующим образом при этом выделяют:

где x={diag[x(k); n=0, 1, …, N-1]}, CT=[cos[2πf0(0)Δt], cos[2πf0(1)Δt], …, cos[2nf0(N-1)Δt]], ST=[sin[2πf0(0)Δt], sin[2πf0(1)Δt], …, sin[2πf0(N-1)Δt]]. Здесь Cx и Sx относятся к коэффициентам фазора, которые существенны для характеристики измерения DFT. При этом формула (1) может быть перезаписана с коэффициентами фазора следующим образом:

Если входной сигнал статичен и имеет номинальную частоту

то Cx и Sx

В соответствии с выражениями (15) и (16) при увеличении k коэффициенты фазора состоят из компоненты 2f0 и компоненты постоянного тока. Путем усреднения коэффициентов фазора в одном цикле сигнала компонента 2f0 может быть исключена и затем может быть получен действительный фазор, как показано в формуле (17). Это алгоритм DFT с одним циклом.

Алгоритм DFT широко используют из-за того, что он позволяет исключать гармоники и извлекать основную волну. Однако если частота смещается от номинального значения, то будет происходить утечка спектра. В частности, если частота изменяется в зависимости от времени, то трудно устранить утечку спектра. Кроме того, при поступлении динамического сигнала также будет происходить утечка спектра.

Прежде всего, используют формулу (1) в качестве входного сигнала, тогда Cx и Sx будут

Реальный фазор должен походить на формулу (3), которая составлена из C2 и S2 согласно формуле (18) и формуле (19). Однако C1 и S1 не могут быть устранены, так как их частоты больше не являются кратными f0. Поэтому происходит утечка спектра.

Как показано в формулах (18) и (19), все частоты частотных компонент, приводящих к утечке спектра, находятся вблизи 2f0, в то время как частоты частотных компонент, составляющих реальный фазор, относительно низки. Согласно требованиям для PMU Р-класса, содержащимся в IEEE стандарте С37.118.1, все значения Δf и f1 составляют менее 2 Гц, и для сохранения низкочастотных компонент может быть осуществлен фильтр низких частот, в то время как высокочастотные компоненты фильтруются в коэффициентах фазора. По этой причине в варианте осуществления настоящего изобретения предложен усовершенствованный алгоритм DFT с цифровым фильтром низких частот для устранения утечки спектра и получения исходного фазора ( и ). Формула выглядит следующим образом:

где - рассчитываемый фазор, - коэффициенты фильтра нижних частот и . Следует отметить, что k начинается с -(N-1)/2, потому что положение метки времени установлено в середине окна времени.

Переписав формулу (24) с коэффициентом фазора, получают:

Можно заметить, что, несмотря на то, что утечка спектра устранена, амплитуда зависит от Xm(t), Δf(t) и коэффициентов фильтра. Точный фазор, который выражается как , может быть получен, только если Xm(t) является константой и Δf(t)=0. Иначе существует ошибка измерения амплитуды. Тем не менее, коэффициенты фильтра считают величинами, фиксированными для выполняемого алгоритма. Ошибка вычисления амплитуды связана только с Δf(t). Зависимость между ними может быть принята и аппроксимирована в виде квадратичной функции путем подбора кривой для результатов моделирования, как показано на фиг.2, где

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

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

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

Во-вторых, уточняют параметры исходного динамического фазора для получения , и .

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

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

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

Заменяя Xm(t) в формуле (25) выражением (5) и предполагая для упрощения процесса получения производной Δf(t)=0, формулу (25) можно перезаписать в следующем виде:

где , , и

Поскольку коэффициенты фильтра являются симметричными, H может быть выражен как

.

Таким образом,

где - средневзвешенное значение K2. Заменив относительные переменные в формуле (26) выражениями (27), (28), (29), получим

Поскольку метка времени установлена в середине окна времени, действительный фазор будет . Таким образом, ошибка измерения амплитуды фазора

где em - ошибка измерения амплитуды, . Для алгоритма, выполняемого в устройстве, интервал выборки и коэффициенты фильтра являются фиксированными, то есть Δt и постоянны. Следовательно, ошибка измерения амплитуды фазора имеет линейную зависимость от m2. Кроме того, коэффициент наклона линейной зависимости зависит от фильтра низких частот.

Для компенсации ошибки измерения сначала необходимо вычислить m2. В предположении, что фазор “М” измерений амплитуды, вычисленный по формуле (20), имеет вид

где - r-е измерение амплитуды фазора, Δtc - интервал расчета фазора. Вышеупомянутая формула может быть записана в матричном представлении как

где Т - матрица коэффициентов. Неизвестную матрицу Cm вычисляют по способу взвешенных наименьших квадратов (WLS):

Опять же, поскольку метка времени установлена в середине, повторно вычисляют для получения аппроксимированной кривой:

Затем вычисленный m2 используют для калибровки . Компенсирующий коэффициент α может быть определен модельными экспериментами.

Аналогично, этим способом с использованием измеренного фазора () может быть вычислен также n2. Затем выполняют измерения первоначальной частоты и ROCOF:

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

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

Для рассматриваемого конкретного случая частота среза второго цифрового фильтра низких частот составляет 2,5 Гц. Для обеспечения высокой скорости динамической реакции окно времени устанавливают величиной 27,5 миллисекунд.

Для проведения испытания способом моделирования для вышеупомянутого способа измерения использованы конкретные варианты осуществления настоящего изобретения. В частности:

В стандарте IEEE С37.118.1 четко определены максимальные ошибки при измерении фазора, а также статические и динамические тесты для имитирования статических и динамических процессов в энергетической системе. Предложенный вариант выполнения способа синхрофазорного измерения в соответствии с настоящим изобретением был подвергнут имитационному тесту в условиях отклонения частоты, гармонических искажений, колебаний, рассинхронизации и замыканий. Показано сравнение с предельными ошибками согласно стандарту, чтобы продемонстрировать преимущество предложенного способа. Для выполнения алгоритма номинальная частота составляет 50 Гц, скорость выдачи отчетов составляет 50 Гц и частота выборки составляет 1200 Гц. N=48, M=7, N'=10. Матрица Cm в формуле (34) может быть вычислена предварительно, при этом при оценке фазора может использоваться способ таблиц подстановки. Полная вычислительная нагрузка при расчете фазора, частоты и ROCOF составляет приблизительно 296 операций умножения и 237 операций суммирования.

1) Тест частотного сканирования

Различные режимы работы энергетической системы могут вызывать дрейф частоты относительно ее номинального значения. Кроме того, замыкания могут приводить к большим отклонениям частоты. В этой секции исследована точность измерения предложенного алгоритма в отношении сигналов с отклонением частоты. Диапазон изменения частоты входного сигнала составляет 48-52 Гц и остается неизменным во время каждого состояния. В тесте частотного сканирования, в соответствии с IEEE С37.118.1а, максимальная суммарная векторная ошибка (TVE) ограничена в пределах 1%, максимальная ошибка частоты (FE) составляет 0,005 Гц и максимальная ошибка ROCOF (RFE) равна 0,1 Гц/с. Соответствующие измерения показаны на фиг. 3 и 4. На фиг. 3 изображена максимальная TVE для различных основных частот при моделировании теста частотного сканирования, выполненного в соответствии с этим вариантом осуществления изобретения. На фиг. 4 изображены максимальные ошибки частоты и ROCOF для различных основных частот при моделировании теста частотного сканирования, выполненного в соответствии с этим вариантом осуществления изобретения. Можно заметить, что предложенный алгоритм достигает намного лучшей точности, чем это определено требованиями IEEE С37.118.1а.

2) Тест влияния гармоник

В этой секции для исследования характеристики восприимчивости предложенного алгоритма к гармоникам использован сигнал энергетической системы с частотой 50 Гц, содержащий 2-ю, 3-ю, 8-ю и 13-ю гармоники в количестве 1%. В Таблице 1 также показаны допустимые ошибки согласно стандарту (станд.).

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

3) Тест модуляции

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

В этом тесте считается, что глубина модуляции амплитуды составляет 10%, глубина модуляции угла сдвига фаз составляет 0,1 рад, частота модуляции изменяется от 0,1 до 2 Гц. В тесте модуляции, согласно IEEE С37.118.1а, максимальная TVE ограничена величиной 3%, максимальная ошибка частоты (FE) составляет 3 Гц. Ошибки измерения показаны на фиг. 5, 6 и 7. На фиг. 5 показана максимальная TVE для различных частот модуляции при моделировании теста модуляции, выполненного в соответствии с этим вариантом осуществления изобретения. На фиг. 6 показана максимальная ошибка частоты для различных частот модуляции при моделировании теста модуляции, выполненного в соответствии с этим вариантом осуществления настоящего изобретения. На фиг. 7 показана максимальная ошибка ROCOF для различных частот модуляции при моделировании теста модуляции, выполненного в соответствии с этим вариантом осуществления настоящего изобретения. Можно заметить, что все показанные на фиг.5-7 ошибки измерений непрерывно увеличиваются при увеличении частоты модуляции, потому что в пределах одного окна времени сигнал изменяется все быстрее. Тем не менее, предложенный алгоритм показывает возможность точного динамического отслеживания по сравнению с предельными ошибками, указанными в IEEE С37.118.1а.

4) Тест линейного изменения частоты

Тест линейного изменения частоты имитирует рассинхронизацию энергетической системы, что отличается от теста частотного сканирования, в котором основная частота сигнала изменяется от 48 до 52 Гц со скоростью 1 Гц/с. Из этого теста можно заметить, что способ измерения, предложенный в настоящем изобретении, обеспечивает возможность точного измерения фазора, частоты и скорости изменения частоты при непрерывном изменении частоты, как показано в Таблице 2.

5) Тест ступенчатого изменения

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

В этом тесте входной сигнал изменяли путем скачкообразного изменения амплитуды на 10% и путем скачкообразного изменения угла сдвига фаз на 10°. Время реакции в каждом случае показано в Таблице 3. Можно видеть, что время реакции фазора удовлетворяет требованиям стандартов.

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

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

1. Способ синхрофазорного измерения для использования в устройстве измерения фазоров (PMU) Р-класса со следующими характеристиками:

упомянутый способ измерения основывают на математической модели динамического фазора, конструируют цифровой фильтр низких частот для коэффициентов фазора, объединенный с дискретным преобразованием Фурье (DFT) и устраняющий утечку спектра, вызванную входными сигналами динамического фазора, при этом после устранения утечки спектра могут проводить измерения исходного фазора;

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

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

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

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

3. Способ по п. 1, который, после выполнения точных измерений динамического фазора, дополнительно включает в себя:

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

4. Способ по п. 1, в котором частота среза цифрового фильтра низких частот, объединенного с дискретным преобразованием Фурье (DFT), составляет 2,5 Гц, окно времени составляет два номинальных периода динамического входного сигнала и расчетная частота исходного фазора составляет 400 Гц.

5. Способ по п. 3, в котором частота среза второго цифрового фильтра низких частот составляет 2,5 Гц и окно времени составляет 27,5 мс.



 

Похожие патенты:

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

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

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

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

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

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

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

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

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

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

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

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

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

Измеритель содержит источник света и установленные последовательно многомодовое оптическое волокно, первый поляризатор, активный элемент ячейки Фарадея, второй поляризатор, плоскость пропускания которого составляет угол ±45° с плоскостью поляризации первого, собирающую линзу, второе многомодовое оптическое волокно и фотоприемник, а также линейный усилитель сигнала фотоприемника, блок преобразования сигналов и индикатор результатов измерения. Активный элемент ячейки Фарадея выполнен в виде четырехугольной призмы высотой h, одна пара боковых граней которой имеет ширину не менее диаметра D коллимированного пучка света, а противоположные боковые грани имеют ширину не менее 3D, первое основание призмы, на которое падает свет, полировано, и на его поверхности в центре нанесено зеркальное покрытие в виде прямоугольной полоски шириной D, другое основание призмы разделено на три равные прямоугольные зоны, по обе стороны от центральной прямоугольной зоны содержит две полированные поверхности с зеркальными покрытиями, составляющие с плоскостью центрально зоны углы γ = arctg(0,5D/h). Технический результат – уменьшение искажений состояния поляризации света, повышение чувствительности и точности измерений. 3 з.п. ф-лы, 5 ил.

Группа изобретений относится к области электрических измерений, в частности к высокоточным устройствам измерения постоянного и переменного напряжения на основе резистивных делителей. Способ измерения высокого напряжения с помощью резистивного делителя заключается в измерении значений напряжения на выходе низковольтного плеча делителя при периодической калиброванной вариации сопротивления его высоковольтного плеча. При этом выполняется формирование первого и второго наборов дискретных измеренных значений напряжения соответственно при первом значении сопротивления высоковольтного плеча на нечетных периодах вариации и при втором значении сопротивления высоковольтного плеча на четных периодах вариации. Далее производится формирование первого и второго наборов дискретных значений для всех нечетных и четных периодов вариации путем добавления к сформированным первому и второму наборам измеренных значений дискретных значений, вычисленных методом интерполяции соответственно при первом значении сопротивления высоковольтного плеча на четных периодах вариации и при втором значении сопротивления высоковольтного плеча на нечетных периодах вариации. На основании сформированных наборов дискретных значений напряжения рассчитывается набор дискретных значений высоковольтного сопротивления, после чего производится вычисление расчетного значения высоковольтного сопротивления путем интегрирования полученного набора дискретных значений сопротивлений с отбрасыванием аномальных значений. Определение значения высокого напряжения выполняют по расчетным формулам соответственно для нечетных и четных периодов вариации с учетом измеренных значений напряжений на выходе делителя, значений калиброванных сопротивлений и расчетного значения высоковольтного сопротивления. Для реализации способа предлагается устройство измерения высокого напряжения. Технический результат изобретения заключается в повышении точности измерения высокого напряжения резистивным делителем за счет исключения влияния изменения значения высокого напряжения в процессе измерений, возможности измерения мгновенных значений напряжения и уменьшении воздействия температурных и других медленно меняющихся процессов, в том числе старения, на точностные характеристики делителя. Кроме того, обеспечивается возможность использования в конструкции высоковольтного плеча делителя резисторов с большими температурными коэффициентами сопротивления. 2 н. и 2 з.п. ф-лы, 3 ил.

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

Группа изобретений относится к измерениям параметров электросетей, в частности к определению фазоров напряжения и тока в электрической сети среднего напряжения точным образом без необходимости в усложненных датчиках, и к определению и мониторингу мощности, развиваемой каждым из проводников, с использованием средств, обычно имеющихся в электрических сетях среднего напряжения. Раскрыты способ и соответствующее устройство для мониторинга параметров электрической сети среднего напряжения, включая определение силы тока, напряжения и мощности каждой фазы для электрической сети среднего напряжения. Текущие параметры электрической сети среднего напряжения, то есть фазоры тока и напряжения, а также мощности, определяются на основе измерений, выполненных датчиками (12, 14, 16), обычно установленными в электрической сети (5, 7) на уровне трансформатора (8). Конкретно определение фазора напряжения на каждом проводнике электрической сети (5) среднего напряжения выполняется с помощью амплитуды, выведенной из измеренной в электрической сети (7) низкого напряжения, и фазового угла, измеренного в электрической сети (5) среднего напряжения. Составления пар между фазорами тока среднего напряжения, углом, измеренным на среднем напряжении и выведенной амплитудой низкого напряжения выполняются с помощью сравнения с коэффициентом мощности cos ϕ электрической сети. Технический результат заключается в обеспечении приемлемой точности измерений мощности без применения усложненных датчиков за счет измерений трехфазных напряжений и мощностей в подстанциях MV/LV с особенностью обращения к информации о напряжениях, измеренных на стороне LV. 4 н. и 11 з.п. ф-лы, 2 ил.

Изобретение относится к контрольно-измерительной технике и может быть использовано при построении устройств для бесконтактного измерения мгновенных значений токов в симметричном трехжильном кабеле без металлических покровов с круглыми жилами. В предлагаемом способе определения мгновенных значений токов в трехжильном кабеле без металлических покровов, включающем измерение напряженности магнитного поля над жилой кабеля, вычисление мгновенного значения тока в жиле, измерения проводятся тремя однокомпонентными датчиками напряженности магнитного поля, расположенными на расстоянии d над жилами кабеля перпендикулярно радиусу и трассе кабеля, а вычисление мгновенных значений токов в каждой жиле осуществляется по формулам: ; ; ;где Н1, Н2, Н3 - напряженности магнитного поля, измеряемые с помощью однокомпонентных датчиков; ; ;r1 - расстояние от центра жилы до поверхности кабеля; r2 - расстояние между центрами жил трехжильного кабеля; d - расстояние от поверхности кабеля до оси датчика. Техническим результатом является обеспечение возможности определения мгновенных значений токов в трехжильном кабеле без металлических покровов при расположении датчиков на расстоянии от кабеля. 1 ил.

Использование – в области электротехники. Технический результат – снижение потерь и повышение надежности контактора. Настоящее изобретение относится к способу 60 в блоке 12 управления для замыкания контактора 1. Блок 12 управления выполнен с возможностью обеспечения перемещения между замкнутым положением и разомкнутым положением за счет возбуждения катушки 6 электромагнитной схемы. Способ 60 содержит подачу 61 напряжения на катушку 6; определение 62, в течение первого периода времени, тока, протекающего через катушку 6, и напряжения на катушке 6 и оценку на этом основании параметров модели для модели, предсказывающей поведение тока, протекающего через катушку 6, как если бы контактор 1 оставался в разомкнутом положении, и измерение 63, после окончания первого периода времени, тока, протекающего через катушку 6, и определение 64 разности между, с одной стороны, измеренным током и, с другой стороны, предсказанным током модели, и повторение измерения 63 и определения 64 до тех пор, пока не будет обнаружено изменение состояния с разомкнутого положения на замкнутое положение с помощью разности токов. 2 н. и 16 з.п. ф-лы, 10 ил.

Изобретение относится к крепежному элементу для сенсора тока и направлено на сокращение ручного труда при монтаже. Крепежный элемент имеет стопорное устройство, а также фланцевую область для крепления сенсора тока в вертикальном положении на крепежной поверхности. Фланцевая область имеет опорную поверхность и/или опорную планку, которые лежат в первой плоскости, крепежный элемент имеет также зажимную поверхность, на которой отформовано стопорное устройство. Зажимная поверхность лежит во второй плоскости, и при этом эти две плоскости имеют общий угол (α) пересечения, который составляет больше 90,25°. 2 н. и 13 з.п. ф-лы, 16 ил.
Наверх